Algoritmo do método de Heun para Matlab

%este programa calcula o metodo de heun para um PVI


%para a equaçao f(x,y)=x^2+Y^2

%entrada de dados



x0=3 ;y0=1 ;h=0.5;



%inicializaçao



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



for i=1:4

k1=(x(i))*(y(i)^2-1);

k2=((x(i)+h)*((y(i)+h*k1)^2)-1);

y(i+1)=y(i)+(h/2)*(k1+k2);

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

end



%saida de dados

for i=1:4

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

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

end



%Saida de dados

y

yexato

erro



plot (x,y)