R=0.3/100L=0.6/100M=0.5/100l=0.3/100m=0.3/100k=0.3/100m1=4.9m2=4.9m0=4*Pi*(10^-7)g=9.8h=0.3/100sigma=beta'[t]omega=((1/Abs[Cos[alpha[t]]])-Abs[Cos[alpha[t]]])*sigmaa=R*r*Abs[Sin[ceta]]*Abs[Sin[alpha[t]]]b=R*pi-R*L-(R^2)*Abs[Sin[alpha[t]]]+R*L*Cos[alpha[t]]+R*pi*Abs[Tan[alpha[t]/2]]*Abs[Tan[alpha[t]]]-L*R*Abs[Tan[alpha[t]/2]]*Abs[Tan[alpha[t]]]-(R^2)Abs[Tan[alpha[t]/2]]*Abs[Tan[alpha[t]]]*Abs[Sin[alpha[t]]]-R*r*Abs[Sin[alpha[t]/2]]*Abs[Cos[ceta]]c=-r*R*Abs[Sin[ceta]]+r*R*Abs[Sin[ceta]]*Abs[Cos[alpha[t]]]-R*r*Abs[Tan[alpha[t]/2]]*Abs[Tan[alpha[t]]]+Abs[Sin[ceta]]d=(a^2+b^2+c^2)*(M/(R^2)*Pi*L)*rH=(R*Abs[Tan[alpha[t]/2]]+L/2)*Abs[Cos[alpha[t]]]+R*Abs[Tan[alpha[t]/2]]+l/2+hx=(l/2)+(l/2*Abs[Cos[alpha[t]]])+k*Abs[Tan[alpha[t]]]+hI1=(1/2)*M*(R^2)I2=Integrate[d,{pi,0,L},{ceta,0,2*Pi},{r,0,R}]I3=(5/4)*M*(R^2)+(1/3)*M*(L^2)T1=(1/2)*I1*(omega^2)T2=(1/2)*I2*(sigma^2)T3=(1/2)*I3*(alpha'[t])^2T=T1+T2+T3V1=m*g*HV2=-((m0*m1*m2)/(2*Pi))*(1/x^3)*Abs[Cos[alpha[t]]]V=V1+V2p=T-Vsol=NDSolve[{D[D[p,alpha'[t]],t]-D[p,alpha[t]]==0,D[D[p,beta'[t]],t]-D[p,beta[t]]==0,alpha[0]==0,beta[0]==0,alpha'[0]==0,beta'[0]==0},{alpha[t],beta[t]}{t,0,0.000001}]Plot[alpha[t]/.sol,{t,0,0.000001}]