runge_kutta.m

用matlab实现Romberg积分

function [t,w]=Runge_Kutta(a,b,h,w0)
n=(b-a)/h;
t=a;
w=w0;
for i=1:n
    k1=h*feval('f',t,w);
    k2=h*feval('f',t+h/2,w+k1/2);
    k3=h*feval('f',t+h/2,w+k2/2);
    k4=h*feval('f',t+h,w+k3);
    w=w+(k1+2*k2+2*k3+k4)/6;
    t=a+i*h;
end