首先D2y（0）能等于0二阶微程能两初始值想想Dy（0）=0,y（0）=1带进D2y（0）-Dy（0）+y（0）=0 D2y（0）自先定义微程组函数建立rigid.m文件function dy = rigid(t,y)dy = zeros(2,1); % a column vectordy(1) = y(2) ;dy(2) = y(2) - y(1);保存另脚本文件！
