Heun Arquivo

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

function y=f(t,u)
 y=u^2
endfunction
function [u]=heun(N,cor) //N é o tamanho do intervalo
u(1)=0.1;
h=10^-3;
T=4; // t o qual se quer saber
t(1)=T-(h*N);
for n=1:N
 t(n+1)=t(n)+h;
 util=u(n)+h*f(t(n),u(n)) //+h*f1
 F1=f(t(n),u(n))
 F2=f(t(n+1),util) //heun
 u(n+1)=u(n)+(h/2)*(F1+F2)
end
ultimo=u(n+1)
plot(t,u,cor) ; xgrid 
endfunction