阶跃响应曲线 num=[23.15 39.65]; den=[2.835 7.726 20.92]; [y,x,t]=step(num, den);plot(t,y) 时域分析零输入响应曲线 A=[1 -1 0.5;2 -2 0.3;1 -4 -0.1] ; B=[0 0 1]'; C=[0 0 1]; D=0; x0=[1 0 0]'; t=0:0.1:20; initial(A,B,C,D,x0,t); title('The Initial Condition Reponse')　

计算流体力学程序，2阶精度龙哥库塔。包括CFL，网格读入，求解器编写，输出等多个源代码。-code of computational fluid dynamics, Runge-Kutta Scheme