Matlab equadiff

[invite642d375a]

Bonjour, je suis embêté par matlab qui me renvoit une erreur bizarre:
je cherche à résoudre:
y1''=-2*y1'/x+2*y1/x²+sin(ln(x))/x² pour x entre 1 et 2
avec les conditions: y1(1)=1 ; y1'(1)=0

voilà les codes:
clear all;
xspan=[1 2];
y1ini=[1 0]';
y2ini=[0 1]';
[x,y1]=ode45(@func1,xspan,y1ini);

function ydot1=func1(x,y1)
ydot(1)=y1(2);
ydot(2)=-2*y1(2)/x+2*y1(1)/x^2+sin(log(x))/x^2;
ydot=ydot(: );


et ça me renvoie:
Error using ==> C:\matlab\toolbox\matlab\funfu n\private\odearguments
Solving FUNC1 requires an initial condition vector of length 0.



-----

[invite5e1117d5]

Je vois pas l'erreur. Peut-être que le nom de ta fonction est ydot et non ydot1. J'ai des doutes sur le fait que xspan soit en ligne et non en colonne. Bref, c'est de la syntaxe.

Je te propose le code suivant, qui fonctionne apparemment.




f=inline('[x(2) ; -2.*x(2)./(t+2)+2.*x(1)./t^2+sin(log(t))/t^2 ]','t','x');

[tsol,ysol]=ode45(f,[1 2]',[1 0]');

plot(tsol,ysol)
figure
a=ysol(:,1);
b=ysol(:,2);
plot(a,b)

[invite642d375a]

ok merci pour l'info; problème, jamais vu inline en cours, ca risque de ne pas passer dans un contrôle...
enfin je vais faire avec ça parce qu'au moins là ca marche