7716 Rate this article:
No rating

Dynamic Plots Using an Equation Function

Mark Bowersox

This blog post continues to explore the Dynamic Visualizations available using  the Equation argument on the PLOT  function. In this example, we explore the use of a function name rather than an equation string. See Dynamic  Plots Using an Equation String for the first article in this series.

Using an equation string has some limitations:

  • You can only have a single statement
  • You cannot easily change the equation unless youset a new string
  • You cannot pass parameters into your equation

As a different approach, create an IDL function containing your equation and then pass the function name to the Equation argument. This method allows you a bit more flexibility in your input equations.

For example, in the first line of the code above, we could have simply written:

p1= PLOT('LambertW', '2', DIMENSIONS=[400,400],$

  NAME='Upperbranch', $

  TITLE='LambertW Function', XRANGE=[-0.4, 2])

 

Notice that we no longer have an X variable in our equation, we just have the name of the function. We can also create our own routine which accepts our X vector and some optional user data.

First, create a new IDL routine called ex_plot_function andsave it in a file ex_plot_function.pro on IDL's path:

FUNCTIONex_plot_function, x, k

COMPILE_OPTIDL2

RETURN,LAMBERTW(DCOMPLEX(x), k)

END

 

Next, we create our plot visualization, passing in the nameof our equation along with our user data containing the "branch"parameter k:


p1= PLOT('ex_plot_function', '2', DIMENSIONS=[400,400],$

NAME='k= 0', EQN_USERDATA=0, $

TITLE='$\Re${LambertW}',XRANGE=[-1, 2])

p2r= PLOT('ex_plot_function', '2r', /OVERPLOT, $

NAME='k= -1', EQN_USERDATA=-1)

p3r= PLOT('ex_plot_function', '2g', /OVERPLOT, $

NAME='k= 1', EQN_USERDATA=1)

p4r= PLOT('ex_plot_function', '2b', /OVERPLOT, $

NAME='k= 2', EQN_USERDATA=2)

lg= LEGEND(/DATA, POSITION=[1.9, -4], LINESTYLE=6, SHADOW=0)

Our plot should now look like the following:

Again, we can pan and zoom around the plot, and IDL willautomatically update the equations.

Please login or register to post comments.