I vectorized the function as there are many applications where is the useful.

Cheers, 

Brian Larsen
MSU - Physics
larsen@ssel.montana.edu