UNISAM - UNIformly SAMple a function y(x) specified as x,y pairs unisam xin= yin= nout= [optional parameters] >binaryfile ... or ... unisam xfile= yfile= npairs= nout= [optional parameters] >binaryfile ... or ... unisam xyfile= npairs= nout= [optional parameters] >binaryfile Required Parameters: xin=,,, array of x values (number of xin = number of yin) yin=,,, array of y values (number of yin = number of xin) ... or xfile= binary file of x values yfile= binary file of y values ... or xyfile= binary file of x,y pairs npairs= number of pairs input (active only if xfile= and yfile= or xyfile= are set) nout= number of y values output to binary file Optional Parameters: dxout=1.0 output x sampling interval fxout=0.0 output first x method=linear =linear for linear interpolation (continuous y) =mono for monotonic cubic interpolation (continuous y') =akima for Akima's cubic interpolation (continuous y') =spline for cubic spline interpolation (continuous y'') isint=,,, where these sine interpolations to apply amp=,,, amplitude of sine interpolations phase0=,,, starting phase (defaults: 0,0,0,...,0) totalphase=,,, total phase (default pi,pi,pi,...,pi.) nwidth=0 apply window smoothing if nwidth>0 AUTHOR: Dave Hale, Colorado School of Mines, 07/07/89 Zhaobo Meng, Colorado School of Mines, added sine interpolation and window smoothing, 09/16/96 CWP: John Stockwell, added file input options, 24 Nov 1997 Remarks: In interpolation, suppose you need 2 pieces of sine interpolation before index 3 to 4, and index 20 to 21 then set: isint=3,20. The sine interpolations use a sine function with starting phase being phase0, total phase being totalphase (i.e. ending phase being phase0+totalphase for each interpolation).