Lösa en ODE i matlab - Mekanik

Hej, jag skall lösa ett svängningsproblem med matlab. Här är frågan jag sitter med:

Detta är härledningen som jag utförde för ekvationerna, en för skivan och en för hjulen:

Massan, längden, radien och tröghetsmomentet är konstanter som vi fått angivna, detta gäller även: $\overset{. .}{x} = r * \overset{. .}{φ}, T = T_{0} \cos Ω t$

När jag försökte lösa denna differential ekvation i matlab testade jag följande kod:

[t,x]=ode113(@deriv, [0 100], [1;0;0.2;0]);

function thetder=deriv(t,x)
G = (-m*L*(x(2))^2*r*sin(x(1))+m*L*r*(x(4))*cos(x(1))+T)/(I_hjul+m*r^2+m_hjul*r^2);
F = (-L^2*m*(x(1))^2*(cos(x(1)))^2*sin(x(1))+L*m*g*sin(x(1))*cos(x(1))+L*m*r*G+L^2*m*(x(2))^2*sin(x(1)))/((cos(x(1))*(I_skiva+L^2*m+L^2*m*(sin(x(1)))^2)));

thetder = [x(2);F;x(4);G];
end

*Jag har också programmerat in värdena för konstanterna men tog inte med i inlägget här så det inte blir för långt. m = massan för skivan

$θ = x (1), \dot{θ} = x (2) φ = x (3), \dot{φ} = x (4) F = \overset{. .}{θ}, G = \overset{. .}{φ}$

När jag försöker köra detta program får jag dock error, någon som kan se vart jag gjort fel?

Svara