#NEWLINE# def cat( x , y , alpha , eta = 1 , purity = 1 ) : #NEWLINE# #INDENT# ' Multiphoton shrodinger cat. eta is the fidelity, alpha the number\n        of photons' #NEWLINE# cos = numpy . cos #NEWLINE# exp = numpy . exp #NEWLINE# return ( ( 1 + ( ( eta * ( ( exp ( ( ( - ( x ** 2 ) ) - ( ( y - alpha ) ** 2 ) ) ) + exp ( ( ( - ( x ** 2 ) ) - ( ( y + alpha ) ** 2 ) ) ) ) + ( ( ( 2 * purity ) * exp ( ( ( - ( x ** 2 ) ) - ( exp ** 2 ) ) ) ) * cos ( ( ( 2 * alpha ) * x ) ) ) ) ) / ( 2 * ( 1 + exp ( ( - ( alpha ** 2 ) ) ) ) ) ) ) / 2 )