function y = yinitfcn(x,region)

global n eta lambda kappa

   switch region
   case 1            % x in [0, 1]
      y = [1;1];     % initial guess for y(x) 

   case 2 % x in [1, lambda]
      y = [1;1];     % initial guess for y(x), 1<=x<=lambda

end