% Author: Housam Binous
% Velocity Profile of Laminar Flow of Power Law fluid in a Pipe
% National Institute of Applied Sciences and Technology, Tunis, TUNISIA
% Email: binoushousam@yahoo.com
% We use the shooting method to get the no-slip boundary
% condition at the pipe wall
function f=zero(x)
global kappa n rend DPoverL
X0=[0,x];
[r X]=ode45(@goveq,[1e-5 rend],X0);
f=X(end,2);
end