function [] = main()
format long;
z=0:0.2:10'; f0=[0 0 0.332057339]; 
[~,u0] = ode45(@blasius, z ,f0); 
u0 = u0(:,2)
end

function dy = blasius(t, y);
dy = zeros(3,1); 
dy(1) = y(2);
dy(2) = y(3);
dy(3) = -.5*y(1)*y(3);
end