function myplotph2(x, y)
L=length(x);
plot(x,y)
hold on
plot(x(1),y(1),'r.')
plot(x(L),y(L),'k.')
xlabel('Решение y')
ylabel('Производная y\prime')
title('График фазовой кривой y=y(t), y\prime=y\prime(t)')