matlab开发-数字集成功能RK1RK2RK3RK4

matlab开发-数字集成功能RK1RK2RK3RK4。数值积分函数（RK1、RK2、RK3、RK4）

function [ sol t ] = RK1（ ABX0dtt_initialt_finalorder ）
%% Coded by
% Mohamed Mohamed El-Sayed Atyya
%% Inputs
% A              : state matix of the diffrential equations  （nxn）
% B              : B matix of the diffrential equations （nx1）
% X0            : initial condition matix （nx1）
% dt             : time step
% t_initial    : initial time
% t_final      : final time
% order       : n
%% Outputs
% sol              : the solution of diffrential equation [x dx/dt d2x/dt2 ...  d（order）x/dt（order） ; t]
% t                  : solution time
%% Function body
n=（t_final-t_initial）/dt;
X（1:order1）=X0;
t（1）=t_initial;
for m=2:ceil（n）+1
    k1=（A*X（1:orderm-1）+B）*dt;
    X（1:orderm）=X（1:orderm-1）+k1;
    t（m）=t（m-1）+dt;
end
sol=[X];
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         779  2016-03-21 04:55  numerical integration functions\RK1.m
     文件         826  2016-03-21 04:54  numerical integration functions\RK2.m
     文件         878  2016-03-21 04:54  numerical integration functions\RK3.m
     文件         916  2016-03-21 04:55  numerical integration functions\RK4.m
     文件        1164  2016-03-21 04:52  numerical integration functions\RK4low.m
     文件        1174  2016-03-21 04:54  numerical integration functions\RK5.m
     文件        1315  2016-07-25 16:50  license.txt