f=zero(x)

% Author: Housam Binous

% Velocity Profile of Laminar Flow of Carreau 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 lamda n nu0 nuI rend

X0=[0 x nu0];

opts = odeset('Mass','M','MassSingular','yes');

[r,X]=ode15s('goveq',[1e-5 rend],X0,opts);

f=X(end,2);

end