adams-2

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

N=10
h=.1
t=[0:h:1]'
y = zeros(N+1,1)
y(1)=1/3
deff('y=f(y,t)','y= -5*y + 5*t^2 + 2*t')
//euler
y(2) = y(1) + h*(f(y(1),t(1)))
//adams
for k = 2:N
 y(1+k) = y(k) + (h/2)*(3*f(y(k),t(k)) - f(y(k-1),t(k-1)))
end
disp([t,y])