%ex0804a.m 用ode23 得到微分方程解并計算出該算法運行時間
fun =inline('-3*y^2+2*x.^2+3*x','x','y'); %用inline構(gòu)造函數(shù)f(x,y)
tic %計時開始
[x,y]=ode23(fun,[0,1],1); %可得到x,y輸出向量值
t=toc %得到ode23的運行時間
ode23(fun,[0,1],1) %可得到輸出得函數(shù)圖