Question: Plotting the Parameters that cause a D.E. Event to Execute

October 08 2012 jschulzb 65
Maple

1

Hi,

     So I've got a numerical solution to a D.E. and want to find which parameters (v, theta) cause the second event to occur. Here's my dsolve

> soln:= dsolve({eqx, eqy, inty, intx, x(0)=0, y(0)=0}, {x(t), y(t)}, numeric, method = rkf45, output=listprocedure, parameters=[v, theta], events=[ [ [diff(y(t),t)=0, y(t)

 The first event is "Does the trajectory reach its max at a height less than 3.05m? If so, stop computation" and second event is "when the height is 3.05m and derivative is negitive, stop the computation". 

I would ultimatly like to plot "v vs theta" to see blotches/areas showing parameters that worked. By worked I mean, did event #2 happen? If yes, plot the value (v, theta)  else continue. 

I've tried for loops, and seq's but I can't figure out how to 'store' the parameter to plot later, and it takes a long time. I'm sure there is an efficient way. Any ideas? 

Please Wait...