t = linspace(-4, 4, 1000);
x = 2*t.^3 - 3*t.^2 - 12*t;
y = 2*t.^3 + 6*t.^2 - 18*t;
figure(1);
plot(x, y, 'k')
hold on;
plot(-20, 4, 'ro', [-20 -20], ylim, 'r', 7, 22, 'bo', [7 7], ylim, 'b');
hold off;
legend('c', 'P_1', 'e_1', 'P_2', 'e_2');
xlabel('x');
ylabel('y');