Algortimo do método de RK4 para Matlab.

%este programa calcula o metodo de rk4


%para um PVI de primeira ordem

%para a equaçao f(x,y)=-y

%entrada de dados



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



%inicializaçao



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



for i=1:10

k1=-y(i);

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

k3=-(y(i)+h*(k2/2));

k4=-(y(i)+h*(k3));

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

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

end



%Saida de dados

for i=1:11

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

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

end



%Saida de dados

%

y

x

yexato

erro



plot (x,y)