TRABALHO 3 MÉTODO DE EULER

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

clear
clc
//definir função
function saida=f(t,y)
 saida=2*(t**2)*(y**2)
endfunction
//definir condição inicial
t(1)=0
y(1)=1
//definir deltaT
dt=0.2
//definir o n° de passos
tf=1 //tempo final
N=int((tf-t(1))/dt)
//avaliar ao longo dos passos
for i=1:N
 k1=f(t(i),y(i))
 y(i+1)=y(i)+dt*k1
 t(i+1)=t(i)+dt
end
//letra C
//definindo valor inicial
yc(1)=y(1)
tc(1)=t(1)
//avaliando ao longo dos passos
for i=1:N
 k1=f(tc(i),yc(i))
 k2=f(tc(i)+dt,yc(i)+k1*dt)
 yc(i+1)=yc(i)+dt*((k1+k2)/2)
 tc(i+1)=tc(i)+dt
end
erro=(3-yc(N+1))/3*100 //calculo do erro
disp(y(N+1))
disp(yc(N+1))
disp(erro)
plot(t(1:N+1),y(1:N+1)) //plotando as soluções para letra B e C
plot(tc(1:N+1),yc(1:N+1))