euler-melhorado;runge-kutta

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

N = 10
h = 0.1
t = [0:h:1]'
y = zeros(N+1,1)
//Inicialização
y(1) = 1/3
deff('z = f(y,t)','z = -5*y + 5*t^2 +2*t')
//Euler Explícito
 y(2) = y(1) + h*(f(y(1),t(1)))
//Euler Melhorado
//for k = 1:N
 // y(k+1) = y(k) + h*f(y(k),t(k))
 //y(k+1) = y(k) + (h/2)*(f(y(k),t(k)) + f(y(k+1),t(k+1)))
//end
//Ruge-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 +1))
 
// y(k+1) = y(k) + (k1 + 2*k2 + 2*k3 + k4)/6 
//end
//Adams-Bashforth-2
for k=2:N
 y(k+1) = y(k) + (h/2)*(3*f(y(k),t(k)) - f(y(k-1),t(k-1)))
end
disp([t,y])