function [poinwert,isterminal,richt] = poinevent(t,y,epsilon)
% Eventfunktion
% Ereignis ist das Durchstossen der Poincareflaeche
% Die Poincareflaeche ist definiert durch y(2)=0

% W. Schweizer
% Simulation physikalischer Systeme - Computational Physics mit MATLAB -
% Beispiele zu Kapitel 1
% vers. 3.2 2013
% Eventfunktionen



%persistent np
%if isempty(np)
%    np=0
%end
%np=np+1
poinwert=y(2);
isterminal=0;
richt=0;
% if (abs(poinwert)<5.e-14)
%     xa=sqrt(-2*epsilon)*y(1);
%     %zwei=0.5*(y(3)*y(3)+y(4)*y(4)) -epsilon*(y(1)*y(1)+y(2)*y(2))+...
%     %1/8*(y(1)*y(1)*y(2)*y(2))*(y(1)*y(1)+y(2)*y(2))
%     plot(xa,y(3),'.',...
%     'MarkerSize',1,'MarkerFaceColor',farb,'MarkerEdgeColor',farb)
%     plot(-xa,y(3),'.',...
%     'MarkerSize',1,'MarkerFaceColor',farb,'MarkerEdgeColor',farb)
%     plot(xa,-y(3),'.',...
%     'MarkerSize',1,'MarkerFaceColor',farb,'MarkerEdgeColor',farb)
%     plot(-xa,-y(3),'.',...
%     'MarkerSize',1,'MarkerFaceColor',farb,'MarkerEdgeColor',farb)
% end

