aiuto MatLab

[noncrederci]

ho scritto il metodo di heun, definendo una funzione in questo modo:

function [x,y]=heun2(f,x0,y0,xf,N)

x(1) = x0 ;

y(1, = y0 ;

h = (xf-x0)/N ;

for k=1:N

x(k+1) = x(k) + h;

fk = feval(f,x(k),y(k,);

ytilde = y(k, + h*fk;

fk1 = feval(f,x(k+1),ytilde);

y(k+1, = y(k, + h*(fk+fk1)/2 ;

end;

end

devo usare questa function per risolvere questo sistema nell'intervallo 1,10 con

condizioni iniziali 1,2 con 1000 sottointervalli. Come devo fare?

function dx=kum(t,x)

dx=zeros(2,1);

dx(1)=x(2);

dx(2)=(5/4)*((x(2)^2)/x(1))-4*x(1)^2+4…

end

potreste scrivermi lo script per usarle che non ci riesco?