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: $\stackrel{..}{x} = r*\stackrel{..}{\phi }, T\hspace{0.17em}=\hspace{0.17em}{T}_0\cos \Omega 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 

$$\begin{gathered} \theta =\hspace{0.17em}x\left(1\right), \stackrel{.}{\theta } = x\left(2\right) \\ \phi = x\left(3\right), \stackrel{.}{\phi } =\hspace{0.17em}x\left(4\right) \\ F\hspace{0.17em}=\hspace{0.17em}\stackrel{..}{\theta }, G\hspace{0.17em}=\hspace{0.17em}\stackrel{..}{\phi } \end{gathered}$$

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