Taylor de terceira ordem

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

function y=f(t,u)
 y=u*t
endfunction
function y=ft(t,u)
 y=(u*t)*t + u
endfunction
function y=ft2(t,u)
 y=(u*t^3)+3*u*t
endfunction
function [ufinal]=taylor(N)
// T tempo final
// N numero intervalos
t(1)=1;
u(1)=0.1;
T=2;
h=( T-t(1) )/N
for n=1:N
 t(n+1)=t(n)+h
 u(n+1)=u(n)+h*f(t(n),u(n))+(h^2/2)*ft( t(n),u(n))+(h^3/6)*ft2(t(n),u(n))
 
 //u(n+a)=u(n)+h*f(t(n),u(n))+(h^2/2!)*ft(t(n),u(n))+ ...+ (h^n/n!)*ftn(t(n),u(n))
 //n! = fatorial de n
 //ftn(t(n),u(n)) = derivada n de y
end
ufinal=u(n+1)
//disp(u)
//plot(t,u,'b.-') ; xgrid 
endfunction