Singular Jacobian Error using bvp4c.

Hi,

I have been having problems with the following code:

function [] = IBVE()
    solinit = bvpinit(linspace(0,0.1,30),@schaetzungbewegungsgleichung);
    options = bvpset('stats','on');
    sol = bvp4c(@bewegungsgl,@residuum,solinit,options);

    function dydt = bewegungsgl(t,y)

    %   M1= (J1+l1^2*(0.25*m1+m2+m3);
    %   M2= (J2+l2^2*(0.25*m2+m3);  
    %   M3= l1*l2*(0.5*m2+m3); 
    %   cd == Federkonstante

        M1=0.70336;
        M2=0.6594;
        M3=0.03768;
        deltaphi0 = 0;
        cd=1;  

        C1 = (y(2)*sin(y(1)-y(3))*y(2));
        C2 = ((y(2))^2*sin(y(1)-y(3)));
        C3 = (cd*(y(1)-y(3)-deltaphi0));
        C4 = cos (y(1)-y(3));
        D1 = 1/(M1*M3*M3/M2*(cos(y(1)-y(3)))^2);
        D2 = (M3)^2/M2*C1+M3/M2*C3*C4+C2*M3+C3;

        dydt = [    y(2)
                D1*D2-D1*(t^3+t^+t+2)
                y(4)
                M3/M2*(C1-C4*D1*(D2+C3/C2))+M3/M2*C4*D1*(t^3+t^+t+2)  ];

    function res = residuum(ya,yb)
        res = [ ya(1)
                ya(2)
                yb(1)-2.3562
                yb(2)     ];

    function v = schaetzungbewegungsgleichung(t)
        v = [4*t^3+3*t^2+2*t+5
            12*t^2+6*t+2
            9*t^3+7*t^2+5*t+7
            27*t^2+14*t+5   ];

When I try to run it, this error appears:

Error using bvp4c (line 252)
Unable to solve the collocation
equations -- a singular Jacobian
encountered.

Error in IBVE (line 4)
    sol =
    bvp4c(@bewegungsgl,@residuum,solinit,options);

Is my guess simply too poor, or is there something that I am missing?

Thanks!

Answer 1

[Diablo]

Hi

Were you able to solve your problem. I have a similar issue.