欧拉法解常微分方程组数值解的MATLAB程序

function y=Eular（funcsabny0）
%用Eular法解常微分方程组的数值解，其中funcs为方程组表达式的右侧项
%如dy1/dx=y1+2*y2+3*y3，则取等式后面项
%其形式是一个细胞数组
%ab为区间左右端点
%n为节点总数
%y0为方程组在自变量为a时的取值形式为一行向量或一列向量
m=max（size（funcs））;
if m~=max（size（y0））;
    error（‘方程组个数与初值个数不匹配‘）;
end
for i=1:m
    funcs{i}=inline（funcs{i}）;
end
y=zeros（nm）;
y（1:）