04-56.txt

来自「MATLABR2006a基础教程」· 文本 代码 · 共 13 行

>> t0 = 1;
>> y0 = sqrt(3/2);
>> yp0 = 0;
>> [y0,yp0] = decic(@weissinger,t0,y0,1,yp0,0);
>> yp0


>> [t,y] = ode15i(@weissinger,[1 10],y0,yp0);


>> ytrue = sqrt(t.^2 + 0.5);
>> plot(t,y,t,ytrue,'o');