rugekutta2

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

N = 20
h = 2/N
t = [0:h:2]'
y = zeros(N+1,1)
//rhs
deff('y = f(t,y)','y = y - t^2 + 1')
//c.i.
y(1) = 0.5
for k = 1:N
 y(k+1) = y(k) + h*f(t(k) + h/2), y(k) + h*f(t(k),y(k))/2
end
disp([t,y])