使用计算机绘制出动力系统的轨线

典型常微分方程，使用ode45求解，代码：

df = @(t,x) [x(1)+x(2); x(1)-x(2)];
[T,X] = ode45(df,[0 1],[1 1]);
figure
plot(T,X)
legend('x(t)','y(t)')
figure
plot(X(:,1),X(:,2))
xlabel('x')
ylabel('y')

结果：