16.12.2012 Views

Computer Algebra Recipes

Computer Algebra Recipes

Computer Algebra Recipes

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

200 CHAPTER 4. NONLINEAR ODE MODELS<br />

OK, X, there's your curve shown in Figure 4.17. You can arrange for a window<br />

to be left open at the point where the curve intersects the museum roof.<br />

Now that we know the distance through which you will be dropping, we<br />

can easily calculate your velocity from the expression v = p 2 g y1. Takingthe<br />

gravitational acceleration to be g =9:81 m/s 2 ,<br />

> g:=9.81: vel:=sqrt(2*g*y1);<br />

vel := 49:98095625<br />

you will be traveling at almost 50 m/s as you pass through the window on<br />

the museum roof, if you haven't been braking. This is an upper bound on<br />

your speed, of course, because we have completely neglected friction and air<br />

resistance. I can see by the look on your face that you don't quite appreciate<br />

how fast this is. So, let me convert the velocity to kilometers per hour.<br />

> vel2:=convert(vel,units,m/s,km/h)*km/h;<br />

179:9314425 km<br />

vel2 :=<br />

h<br />

Your theoretical speed would be almost 180 km/hr but, as I said, in reality it<br />

would be somewhat less. The time of descent can now be calculated from<br />

Z £<br />

T =<br />

p<br />

(dx=dμ) 2 +(dy=dμ) 2<br />

dμ;<br />

2 gy(μ)<br />

viz:;<br />

μ=0<br />

> T:=int(sqrt(diff(x,theta)^2+diff(rhs(Y),theta)^2)<br />

/sqrt(2*g*rhs(Y)),theta=0..evalf(Theta));<br />

T := 5:659009626<br />

The minimum time of descent is about 5:7 seconds. Again, friction, air resistance,<br />

and any braking on your part would lengthen this time. Given the<br />

estimates that we have come up with, do you still want to go through with your<br />

bizarreescaperoute?"<br />

\Sure, Mike," Mr. X replies, \and I am counting on you to help attach the<br />

wire to the appropriate window in the museum roof and leave it open for me."<br />

\Oh no," Mike groans. \If I get caught, Vectoria's parents are sure to call<br />

o® our recently announced engagement."<br />

PROBLEMS:<br />

Problem 4-41: A di®erent escape route<br />

Determine Mr. X's escape route if the museum roof slanted at 30 ± to the horizontal,<br />

the street is 50 m wide, and the vertical section of the museum wall<br />

adjoining the slanted roof is 25 m tall. Make a plot similar to that in Figure<br />

4.17. How close to the edge of the roof would Mr. X land? What would be<br />

his speed at this point? Neglect friction and air resistance.<br />

Problem 4-42: In search of reality<br />

Suppose that Mr. X decides that sliding down a steel wire at nearly 180 km/hr<br />

is a little too insane even for his bizarre tastes. His friend Mike is unavailable for<br />

technical advice because he is busy with preparations for his upcoming wedding

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!