f=zero(x)

% 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: [e-mail]

% 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