function yp=Blasius(x,y)
yp=zeros(3,1);
yp(1)=-0.5*y(3)*y(1);
yp(2)=y(1);
yp(3)=y(2);

