kutta-2

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

N=4
h=0.25
t=[1:h:2]'
y=zeros(N+1,1)
y(1) = 2
deff('y=f(y,t)','y=t^(-2) * (sin(2*t) - 2*t*y)')
//kutta-4
for k=1:N
 k1 = h*(f(y(k),t(k)))
 k2 = h*(f(y(k) + k1/2, t(k) + h/2))
 k3 = h*(f(y(k) + k2/2, t(k) + h/2))
 k4 = h*(f(y(k) + k3, t(k) + h))
 
 y(k+1) = y(k) + (k1 + 2*k2 + 2*k3 + k4)/6
end
disp([t,y])