function [x,f, fp, fs]=fBlasiusBVP

a=0;  %limita dreapta 
b=7;  %limita stanga 


solinit.x=linspace(0,7,701);
solinit.y=[solinit.x; exp(-solinit.x);exp(-solinit.x)];%initializarea solutiei

options=bvpset('AbsTol', 1e-15);
sol = bvp4c(@BlasiusEqs,@BlasiusBC,solinit);

x=sol.x;f=sol.y(1,:);fp=sol.y(2,:);fs=sol.y(3,:);
plot(x,fp)
fs(1)

function dydx = BlasiusEqs(x,y)
dydx = [  y(2)
          y(3)
          -y(3)*y(1)];

function res = BlasiusBC(ya,yb)
res = [ ya(1)
        ya(2)
        yb(2)-1];
