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