function output=f_6b(theta1,theta2,omega1,omega2,t)
m1=1;
m2=2;
L1=1;
L2=2;
g=9.8;
output=omega2;