section i of routines in random.i

functions in random.i - i

ipq_compute


    ipq_compute  


  
     Interpreted function, defined at i/random.i   line 362

ipq_function


    ipq_function  


  
     Interpreted function, defined at i/random.i   line 450

ipq_setup


    model= ipq_setup(x, u)  
 or model= ipq_setup(x, u, power=[pleft,prght])  
 or model= ipq_setup(x, u, power=[pleft,prght], slope=[sleft,srght])  


compute a model for the ipq_compute function, which computes the  
inverse of a piecewise quadratic function.  This function occurs  
when computing random numbers distributed according to a piecewise  
linear function.  The piecewise linear function is u(x), determined  
by the discrete points X and U input to ipq_setup.  None of the  
values of U may be negative, and X must be strictly increasing,  
X(i)0 while SRGHT<0.  If either power is greater than or equal to  
100, an exponential tail will be used.  As a convenience, you may  
also specify PLEFT or PRGHT of 0 to get an exponential tail.  
Note: ipq_function(model, xp) returns the function values u(xp) at  
the points xp, including the tails (if any).  ipq_compute(model, yp)  
returns the xp for which (integral from -infinity to xp) of u(x)  
equals yp; i.e.- the inverse of the piecewise quadratic.  
Interpreted function, defined at i/random.i   line 277