getting a singular jaobian error

hi all, i am getting singular jacobian error in code please help me as early as possible its important function Bvvp4c

xlow=10^(-6);xhigh=0.6;

solinit=bvpinit(linspace(xlow,xhigh,10),[2;4;9;12;3;16;5;88;7;8;11;25;34],10^(-2));

sol=bvp4c(@bvp4ode,@bvp4bc,solinit);

xint=linspace(xlow,xhigh);

Sxint=deval(sol,xint);

plot(xint,Sxint(1,:))

hold on %---------------------------------------------------

function dydx = bvp4ode(x,y,muo)

dydx =[(-y(2)*(y(3)+y(6))/9-y(3)*y(5));

    (y(3)*y(4)+y(1)*(y(3)+y(6)/9));

    (y(1)*y(5)-y(2)*y(4));

    (-(y(5)+y(2))*(y(3)+y(6))-(y(6)+y(3))*y(5));

    ((y(6)+y(3))*y(4)+(y(4)+y(1))*(y(3)+y(6))/9);

    ((y(4)+y(1))*y(5)-(y(5)+y(2))*y(4));

    (((y(10))^2+(y(11))^2-1/2)*y(4)+
(y(11)*y(12)-y(10)*y(13)*y(5)-(y(11)*y(13)+y(10)*y(12))*(y(3)+y(6))/9));

    ((y(11)*y(12)+y(10)*y(13))*y(4)+((y(10))^2+
(y(12))^2-1/2)*y(5)-(y(12)*y(13)-y(10)*y(13))*((y(3)+y(6))/9));

    ((y(11)*y(13)-y(10)*y(12))*y(4)+(y(12)*y(13)+y(10)*y(11))*y(5)-((y(10))^2+
(y(13))^2-1/2)*((y(3)+y(6))/9));

    (-y(12)/2*(1-y(3)/9-y(6)/9)+muo*y(10));

    (y(12)/2*(1-y(3)/9-y(6)/9)+muo*y(11));

    (-y(11)/2*(1-y(3)/9-y(6)/9)+muo*y(12));

    (y(10)/2*(1-y(3)/9-y(6)/9)+muo*y(13))];
  %----------------------------------------------------
function res = bvp4bc(ya,yb,muo)
    res=[ya(2);
        ya(5);
        ya(7);
        ya(9);
        ya(11);
        ya(12);
         yb(7)
         yb(8)
         yb(9)-0.6
         yb(10)-1
         yb(11)
         yb(12)
         yb(13)
         ya(10)^2+ya(11)^2+ya(12)^2+ya(13)^2-1];