def levypdf(x, gamma, alpha):
N = len(x)<EOL>if N%<NUM_LIT:2> != <NUM_LIT:0>:<EOL><INDENT>raise ValueError('<STR_LIT>' +'<STR_LIT>')<EOL><DEDENT>dx = x[<NUM_LIT:1>]-x[<NUM_LIT:0>]<EOL>f = <NUM_LIT:1>/(N*dx)*np.arange(-N/<NUM_LIT:2>, N/<NUM_LIT:2>, np.float_)<EOL>ind = np.concatenate([np.arange(N/<NUM_LIT:2>, N, int),<EOL>np.arange(<NUM_LIT:0>, N/<NUM_LIT:2>, int)])<EOL>df = f[<NUM_LIT:1>]-f[<NUM_LIT:0>]<EOL>cfl = exp(-gamma*np.absolute(<NUM_LIT:2>*pi*f)**alpha)<EOL>px = np.fft.fft(np.take(cfl,ind)*df).astype(np.float_)<EOL>return np.take(px, ind)<EOL>
Returm the levy pdf evaluated at *x* for params *gamma*, *alpha*