odefun3.m

《matlab使用详解》的部分源码

function f = odefun3(t,y,yp)
% 全隐式常微分方程 f(t,y,y') = 0

f = t*y^2*yp^3 - y^3*yp^2 + t*(t^2+1)*yp - t^2*y;