euler (2)

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

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