

dy = @(t, x)[-3*x(1)+3*x(2); -2*x(1)+4*x(2)-2*x(2).^3];

[T, Y] = ode45(dy, [-10, 10], [1, -1]);
plot(Y(:, 1), Y(:, 2))
hold on

[T, Y] = ode45(dy, [-10, 10], [-1, 1]);
plot(Y(:, 1), Y(:, 2))
hold on

[T, Y] = ode45(dy, [-10, 10], [-2, -2]);
plot(Y(:, 1), Y(:, 2))

[T, Y] = ode45(dy, [-10, 10], [2, 2]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [-0.1 -0.1]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [0.1 -0.1]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [-0.1 0.1]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [0.1 0.1]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [0.5, 0.5]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [-0.5, -0.5]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [4, -4]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [-4, -4]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [5, -5]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [-5, -5]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [-5, 5]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [5, 5]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [1, -4]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [-4, 1]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [2, 4]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [1, 0]);
plot(Y(:, 1), Y(:, 2))
hold on;

[T, Y] = ode45(dy, [-10, 10], [2, 0.5]);
plot(Y(:, 1), Y(:, 2), 'r')
hold on;




quiver([0 0], [0 0], [3 1],[1 2],0,'Red')

plot(1, 1, 'or')
plot(-1,-1, 'or')
plot(0,0, 'or')
text(0.5,0, 'GGP1')
text(1,1.5, 'GGP2')
text(-1,-2, 'GGP3')

title('Phasenportrait')
xlabel('u_C')
ylabel('i_L')


