rk_wjl.m

(有源代码)数值分析作业,本文主要包括两个部分,第一部分是常微分方程(ODE)的三个实验题,第二部分是有关的拓展讨论,包括高阶常微分的求解和边值问题的求解(BVP).文中的算法和算例都是基于Matla

function [x,y]=rk_wjl(odefile,xi,xf,yi,h,varargin)

x = [xi:h:xf];		% Vector of x values
if x(end) ~= xf
   x(end+1) = xf;
end
y(1) = yi;		% Initial condition


d = diff(x);
% if no this sentence ,you  change d(i) to  h ,is OK
% this is a fine method 
for i = 1:length(x)-1
      k1 = d(i) * feval(odefile,x(i),y(i),varargin{:});
      k2 = d(i) * feval(odefile,x(i)+d(i)/2,y(i)+k1/2, ...
         varargin{:});
      k3 = d(i) * feval(odefile,x(i)+d(i)/2,y(i)+k2/2, ...
         varargin{:});
      k4 = d(i) * feval(odefile,x(i+1),y(i)+k3,varargin{:});
      y(i+1) = y(i) + (k1+2*k2+2*k3+k4)/6;
end
%  xx=x; 
% yy=y;