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:）