www.gusucode.com > symbolic 源码程序 matlab案例代码 > symbolic/OdeVectorFieldSolveSecondOrderDifferentialEquationExample.m

    %% Solve Second-Order Differential Equation
%% 1

% Copyright 2015 The MathWorks, Inc.

syms y(t)
V = odeToVectorField(diff(y, 2) == (1 - y^2)*diff(y) - y)
M = matlabFunction(V,'vars', {'t','Y'})
sol = ode45(M,[0 20],[2 0]);
%% 2
x = linspace(0,20,100);
y = deval(sol,x,1);
plot(x,y)