clc;
L=[-10?1?2;3?-12?0;1?2?-9];
M=[-1?0?3;1/2?1/2?1;1/2?3/2?0];
N=1/72*[-1?5?2;4?0?3;-2?4?1];?????????%生成三個已知矩陣
B=110;????????????????????????????????%計算步數
h=0.1;????????????????????????????????%計算步長
tau=1;????????????????????????????????%時間延遲
m=tau/h;??????????????????????????????%延時時間計算步數
for?i=1:m?????????????????????????????%初始化參數以及u
????k1（:i）=[0?0?0];
????k2（:i）=[0?0?0];
????k3（:i）=[0?0?0];
????k4（:i）=[0?0?0];
????u（:i）?=[1?1?1];
end
u（:11）=[1?1?1];
for?n=11:B????????????????????????????%計算k1k2k3k4和u在第n步的取值
????k1（:n）=h*L*u（:n）+h*M*u（:n