Advanced Search

Question: Error, (in DEtools/convertsys) unable to convert to an explicit first-order system

Question:Error, (in DEtools/convertsys) unable to convert to an explicit first-order system

Cyrus Mehdipour 5 Maple 13
ode differential-equation dsolve numeric
February 13 2015
1

Hello i want to solve the differentiel equation but i have these problem i don't understund  why !?

drive.mw

restart

ode1 := (K[Q]*T*R[a]/K[kol]+R[a]*B[m]*sqrt(T/K[kol]))/K[i]+K[b]*sqrt(T/K[kol]) = 0

(K[Q]*T*R[a]/K[kol]+R[a]*B[m]*(T/K[kol])^(1/2))/K[i]+K[b]*(T/K[kol])^(1/2) = 0

 (1)

ode2 := (1/2)*(-(4*(diff(theta(t), t)+theta(t)))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/l[1]

(1/2)*(-4*(diff(theta(t), t)+theta(t))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/l[1]

 (2)

 

ode3 := subs(T = (1/2)*(-(4*(diff(theta(t), t)+theta(t)))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/l[1], ode1)

((1/2)*K[Q]*(-4*(diff(theta(t), t)+theta(t))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])*R[a]/(l[1]*K[kol])+(1/2)*R[a]*B[m]*2^(1/2)*((-4*(diff(theta(t), t)+theta(t))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))^(1/2))/K[i]+(1/2)*K[b]*2^(1/2)*((-4*(diff(theta(t), t)+theta(t))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))^(1/2) = 0

 (3)

simplify(ode3)

(1/2)*(-4*R[a]*K[Q]*M*(diff(theta(t), t))-4*R[a]*K[Q]*M*theta(t)+2*R[a]*K[Q]*B1*(diff(theta(t), t))-2*R[a]*K[Q]*w[2]*sin(theta(t))+R[a]*K[Q]*m1*g*sin(theta(t))*l[kol]+2*R[a]*K[Q]*w[1]*sin(theta(t))-2*R[a]*K[Q]*m1*g*sin(theta(t))*l[2]+R[a]*B[m]*2^(1/2)*((-4*M*(diff(theta(t), t))-4*M*theta(t)+2*B1*(diff(theta(t), t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))^(1/2)*l[1]*K[kol]+K[b]*2^(1/2)*((-4*M*(diff(theta(t), t))-4*M*theta(t)+2*B1*(diff(theta(t), t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))^(1/2)*l[1]*K[kol]*K[i])/(l[1]*K[kol]*K[i]) = 0

 (4)

eol := (1/2)*(-4*R[a]*K[Q]*M*(diff(theta(t), t))-4*R[a]*K[Q]*M*theta(t)+2*R[a]*K[Q]*B1*(diff(theta(t), t))-2*R[a]*K[Q]*w[2]*sin(theta(t))+R[a]*K[Q]*m1*g*sin(theta(t))*l[kol]+2*R[a]*K[Q]*w[1]*sin(theta(t))-2*R[a]*K[Q]*m1*g*sin(theta(t))*l[2]+R[a]*B[m]*sqrt(2)*sqrt((-4*M*(diff(theta(t), t))-4*M*theta(t)+2*B1*(diff(theta(t), t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))*l[1]*K[kol]+K[b]*sqrt(2)*sqrt((-4*M*(diff(theta(t), t))-4*M*theta(t)+2*B1*(diff(theta(t), t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))*l[1]*K[kol]*K[i])/(l[1]*K[kol]*K[i])

 

 (5)

R[a] := 1.42;

1.42

 (6)

K[b] := 0.649e-2;

0.649e-2

 (7)

K[i] := 0.649e-2

0.649e-2

 (8)

K[Q] := 0.1051618298e-6

0.1051618298e-6

 (9)

l[kol] := 1

1

 (10)

l[1] := .7

.7

 (11)

l[2] := .3

.3

 (12)

K[kol] := 0.1168464776e-5

0.1168464776e-5

 (13)

B1 := 0.955e-3

0.955e-3

 (14)

B[m] := 0.955e-3

0.955e-3

 (15)

J := 0.475e-6

0.475e-6

 (16)

M := 0.91e-2

0.91e-2

 (17)

m1 := 0.726e-1

0.726e-1

 (18)

w[1] := 0.72e-1

0.72e-1

 (19)

w[2] := .45

.45

 (20)

g := 9.81

9.81

 (21)

a1 := 0

0

 (22)

eol

-.4851223862*(diff(theta(t), t))-.5119876735*theta(t)-6.626549550*sin(theta(t))+.1077211171*2^(1/2)*(-42167.66273*(diff(theta(t), t))-44502.83918*theta(t)-575990.9557*sin(theta(t)))^(1/2)

 (23)

simplify(eol)

-.4851223862*(diff(theta(t), t))-.5119876735*theta(t)-6.626549550*sin(theta(t))+0.1523406647e-3*(-0.4216766273e11*(diff(theta(t), t))-0.4450283918e11*theta(t)-0.5759909557e12*sin(theta(t)))^(1/2)

 (24)

with(plots)

ic1 := theta(0) = a1

theta(0) = 0

 (25)

``

dsol1 := dsolve({eol, ic1}, numeric, output = listprocedure, range = 0 .. 10)

Error, (in DEtools/convertsys) unable to convert to an explicit first-order system

  

plots[odeplot](dsol1, [t, theta(t)], 0 .. 10)

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

  

``


thanks for your help

Download drive.mw
Please Wait...