Generation of random numbers using shape preserving quadratic splines

Abstract

Let F be an arbitrary continuous cumulative distribution function of a single variable specified by a finite set of points. A (smooth) increasing quadratic spline is constructed which interpolates the data points and preserves the convexity of the data [2]. The spline is compared with Akima's piecewise cubic approximation [1] for several common distributions F.