%   function for falkner.m
    function f=fpeq1(A,beta)
    global range
    Vinit=[0 0 A];
    options=odeset('Events',@events);
    [eta,Y]=ode45(@falknerf,range,Vinit,options,beta);
    f=Y(end,2)-1;
