problem with bvp4c

Can anybody help me for the following code? I'm trying to solve a second degree differential equation by bvp4c. I encounter singular jacobian error.

function out= test100(t,y) clc

h=.01; for ii=1:5

    if ii==1
        q= @(t) -t;    
    end;

    sol=bvpinit(linspace(0,20,20),[1  1]);
    sol=bvp4c(@bvpexample,@bc1,sol);
    xx=sol.x;yy=sol.y;
    ht= @(t) spline(xx,yy(1,:),t);
    q= @(t) -numder(ht,t,h);

end

function yprime = bvpexample(t,y) a=10;b=1;c=5;d=5;p=.05;

yprime=[y(2);-1./(2*a.*exp(y(2)-q(t)).*(1+y(2))).*(p*c.*exp(-y(1))+p*d)+y(2).*(p+numder(q,t,h))./(1+y(2))] end

function res=bc1(ya,yb) %BC: Evaluates the residue of the boundary condition res=[ya(1)-10;yb(1)]; end

function der = numder(fun, x, h, varargin) der = (feval(fun,x+h,varargin{:})-feval(fun,x-h,varargin{:}))/(2*h); end

end

[Walter Roberson]

http://www.mathworks.com/matlabcentral/answers/13205-tutorial-how-to-format-your-question-with-markup

Answer 1

[Walter Roberson]

Please see http://www.mathworks.com/matlabcentral/answers/6369-bvp4c