www.gusucode.com > symbolic 源码程序 matlab案例代码 > symbolic/FunctionHandleSuitableForODESolversExample.m
%% Function Handle Suitable for ODE Solvers %% 1 syms y(t) eqn = diff(y(t), t, t) == (1 - y(t)^2)*diff(y(t),t) - y(t); [eqs,vars] = reduceDifferentialOrder(eqn,y(t)); initConditions = [2,0]; [M,F] = massMatrixForm(eqs,vars); f = M\F %% 2 odefun = odeFunction(f,vars); ode15s(odefun, [0 10], initConditions)