xclear; clc; close all;%%f = @(x, y) [ y(2); (1 - y(1)^2)*y(2) - y(1) ];[x, y] = ode45(f, [0 10], [2 0]);plot(x, y, '-*')
clear; clc; close all;
%%
f = @(x, y) [
y(2);
(1 - y(1)^2)*y(2) - y(1)
];
[x, y] = ode45(f, [0 10], [2 0]);
plot(x, y, '-*')
效果