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]
2022-05-24