Advanced Search

Question: error: fsolve unrecognized form

Question:error: fsolve unrecognized form

susanhl 5 Maple 16
fsolve rootof
October 30 2019
1

Hello,

I don't understand why Maple gives me the following error:

Error, (in fsolve/Iterators) unrecognized form

The full maple file:


 

> 

restart
> 

mbal := 0.28e-2:
> 

rbal := 0.2e-1:
> 

Cw := .47:
> 

g := 9.81:
> 

A := Pi*rbal^2:
> 

beta := 3.33*Pi*(1/180):
> 

s0x := .2:
> 

s0y := .25:
> 

rho := 1.293:
> 

l := 2:
> 

vx := diff(sx(t), t);

diff(sx(t), t)

 (1)
> 

vy := diff(sy(t), t);

diff(sy(t), t)

 (2)
> 

ax := diff(sx(t), `$`(t, 2));

diff(diff(sx(t), t), t)

 (3)
> 

ay := diff(sy(t), `$`(t, 2));

diff(diff(sy(t), t), t)

 (4)
> 

v0x := cos(beta)*v0;

cos(0.1850000000e-1*Pi)*v0

 (5)
> 

v0y := sin(beta)*v0;

sin(0.1850000000e-1*Pi)*v0

 (6)
> 

`ΣFx` := -Fdx = mbal*ax;

-Fdx = 0.28e-2*(diff(diff(sx(t), t), t))

 (7)
> 

`ΣFy` := -Fz-Fdy = mbal*ay;

-Fz-Fdy = 0.28e-2*(diff(diff(sy(t), t), t))

 (8)
> 

Fz := mbal*g;

0.27468e-1

 (9)
> 

Fdx := .5*rho*vx^2*A*Cw;

0.1215420e-3*(diff(sx(t), t))^2*Pi

 (10)
> 

Fdy := .5*rho*vy^2*A*Cw;

0.1215420e-3*(diff(sy(t), t))^2*Pi

 (11)
> 

ics1 := sx(0) = s0x, (D(sx))(0) = v0x;

sx(0) = .2, (D(sx))(0) = cos(0.1850000000e-1*Pi)*v0

 (12)
> 

ics2 := sy(0) = s0y, (D(sy))(0) = v0y;

sy(0) = .25, (D(sy))(0) = sin(0.1850000000e-1*Pi)*v0

 (13)
> 

verplaatsingx := dsolve({ics1, `ΣFx`}, sx(t));

sx(t) = (1400000/60771)*ln((60771/1400000)*cos((37/2000)*Pi)*v0*exp((60771/7000000)*Pi)*t*Pi+exp((60771/7000000)*Pi))/Pi

 (14)
> 

verplaatsingy := dsolve({ics2, `ΣFy`});

sy(t) = (1/243084)*(60771*Pi+2800000*ln((1/92736546000)*(20257*Pi^(1/2)*sin((37/2000)*Pi)*v0*sin((3/70000)*231841365^(1/2)*Pi^(1/2)*t)+20*231841365^(1/2)*cos((3/70000)*231841365^(1/2)*Pi^(1/2)*t))^2))/Pi

 (15)
> 

eindy := rhs(verplaatsingy) = rbal;

(1/243084)*(60771*Pi+2800000*ln((1/92736546000)*(20257*Pi^(1/2)*sin((37/2000)*Pi)*v0*sin((3/70000)*231841365^(1/2)*Pi^(1/2)*t)+20*231841365^(1/2)*cos((3/70000)*231841365^(1/2)*Pi^(1/2)*t))^2))/Pi = 0.2e-1

 (16)
> 

eindx := rhs(verplaatsingx) = l;

(1400000/60771)*ln((60771/1400000)*cos((37/2000)*Pi)*v0*exp((60771/7000000)*Pi)*t*Pi+exp((60771/7000000)*Pi))/Pi = 2

 (17)
> 

t := solve(eindy, t);

RootOf(1397733*Pi+70000000*ln((1/92736546000)*(20257*Pi^(1/2)*sin((37/2000)*Pi)*v0*sin((3/70000)*231841365^(1/2)*Pi^(1/2)*_Z)+20*231841365^(1/2)*cos((3/70000)*231841365^(1/2)*Pi^(1/2)*_Z))^2))

 (18)
> 

t := evalf(t[1]);

RootOf(1397733*Pi+70000000*ln((1/92736546000)*(20257*Pi^(1/2)*sin((37/2000)*Pi)*v0*sin((3/70000)*231841365^(1/2)*Pi^(1/2)*_Z)+20*231841365^(1/2)*cos((3/70000)*231841365^(1/2)*Pi^(1/2)*_Z))^2))[1]

 (19)
> 

v0 := fsolve(eindx, v0)

Error, (in fsolve/Iterators) unrecognized form

  
> 

``
> 

``


 

Download ttbuddymaple.mw

 

Thank you!

 
Please Wait...