## Singular Jacobian Error using bvp4c.

on 20 Jun 2014

### Diablo (view profile)

on 5 Jan 2015
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!