Taylor de segunda ordem

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

function y=f(t,u)
 y=-sin(u)*sin(cos(u))*sin(cos(cos(u)))
endfunction
function y=f_t(t,u)
 y=-cos(cos(cos(u)))*(sin(u)^2)*(sin(cos(u)))^2 + cos(cos(u))*sin(cos(cos(u))*(sin(u))^2-cos(u)*sin(cos(u))*sin(cos(cos(u)))
 
endfunction
function [ufinal]=taylor(N)
// T tempo final
// Nint 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)*f_t( t(n),u(n))
end
ufinal=u(n+1)
//disp(u)
//plot(t,u,'b.-') ; xgrid 
endfunction