Advanced Search

Question: conversion to an explicit first-order system; Help pleeeeeeeease

Question:conversion to an explicit first-order system; Help pleeeeeeeease

MOLA 10 Maple 13
ode differential-equation dsolve numeric
May 31 2015
2

I'm currently having some difficulties in solving a system of differential equations numerically.

This is my code.

 

theta[0] := 3.14159*(1/4); R := .1; l := .2; mu := 0.1e-3; rho := 1000; w[0] := .3; M := .15; g := 9.81
 
L[cm] := (2/3*(2*sin(theta[0])^3/(2*theta[0]-sin(2*theta[0]))))*R; p := 2*mu*R^2*l*theta[0]; m := rho*(R^2*theta[0]-(1/2)*R^2*sin(2*theta[0]))
 
eq1 := m*(-L[cm]*(diff(alpha(t), `$`(t, 1)))*cos(alpha(t))+R*(diff(beta(t), `$`(t, 1))))+M*R*(diff(beta(t), `$`(t, 1))) = M*R*w[0]
 
eq2 := diff(alpha(t), `$`(t, 2))+p*(diff(alpha(t), `$`(t, 1))-(diff(beta(t), `$`(t, 1))))/(m*L[cm]^2)+(g*L[cm]*sin(alpha(t))-R*(diff(beta(t), `$`(t, 2)))*cos(alpha(t)))/L[cm]^2 = 0
 
CI := alpha(0) = 0, beta(0) = 0, (D(beta))(0) = w[0], (D(alpha))(0) = w[0]
 
solution := dsolve([eq1, eq2, CI], numeric, range = 0 .. 2)
 
 
Then, "Error, (in DEtools/convertsys) unable to convert to an explicit first-order system" comes out. What's wrong with my code? Hope for someone's help.
 
Please Wait...