    %   Function for falkner.m
%   Y=[f f' f''],  dY/deta=[f' f'' f''']=[Y(2) Y(3) -f*f''/2]
    function dY=falknerf(eta,Y,beta)
    dY(1)=Y(2);
    dY(2)=Y(3);
    dY(3)=-Y(1)*Y(3)-beta*(1-Y(2)^2);
    dY=dY';  
    
    