euler

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

N = 10 //numero de subintervalos
dt = 1/N //passo no tempo
t = [0:dt:1]' //tempos discretos
w = zeros(N+1,1)//vetor solução
y = zeros(N+1,1)//vetor solução
//cond. inicial
w(1) = 0
y(1) = 1
//iteração de Euler
for k = 1:N
 w(k+1) = w(k) + dt*(-w(k) - y(k) + cos(t(k)))
 y(k+1) = y(k) + dt*w(k)
end
disp([t,y])
tt = linspace(0,1);
deff('y = ye(t)','y = sin(t) + exp(-t/2).*(cos(sqrt(3)*t/2) - sin(sqrt(3)*t/2)/sqrt(3))')
plot(t,w,'b.-',tt,ye(tt),'r-')
//plot(t,abs(y-se(t)))