segunda-feira, 30 de agosto de 2010

Algoritmo do método de Euler para o Matlab.

Este problema soluciona um PVI atraves do metodo de euler


%entrada de dados

x0=0;y0=1;h=0.1;

%inicializaçao

x(1) = x0 ;y(1) = y0 ;

for i=1:10

x(i+1)=x(i)+h;

y(i+1)=y(i)+h*(-y(i));

end



%saida de dados

for i=1:11

yexato(i)=exp(-x(i));

erro(i)=abs(y(i)-yexato(i));

end

%Saida de dados

%

x

y

yexato

erro

Nenhum comentário:

Postar um comentário