function [u0]=inital(n,t)
%
% Function inital sets the initial condition vector
% for the 3 x 3 ODE problem
%
  Ca0=1.0;
  u0(1)=Ca0;
  u0(2)=0.0;
  u0(3)=0.0;
%
% End of inital
