Baixe o app para aproveitar ainda mais
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 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)))
Compartilhar