thermal[T_, x_, p_] := If[T === 0, coherent[0, x, p], 1/Pi Tanh[1/(2 T)]*Exp[-(x^2/Coth[1/(2 T)]) - p^2/Coth[1/(2 T)]] ]; rq = 4; rp = 4; t = 0.003; FontS...
thermal[T_, x_, p_] := If[T === 0, coherent[0, x, p], 1/Pi Tanh[1/(2 T)]*Exp[-(x^2/Coth[1/(2 T)]) - p^2/Coth[1/(2 T)]] ]; rq = 4; rp = 4; t = 0.003; FontS...
thermal[T_, x_, p_] := If[T === 0, coherent[0, x, p], 1/Pi Tanh[1/(2 T)]*Exp[-(x^2/Coth[1/(2 T)]) - p^2/Coth[1/(2 T)]] ]; rq = 4; rp = 4; t = 0.003; FontS...
thermal[T_, x_, p_] := If[T === 0, coherent[0, x, p], 1/Pi Tanh[1/(2 T)]*Exp[-(x^2/Coth[1/(2 T)]) - p^2/Coth[1/(2 T)]] ]; rq = 4; rp = 4; t = 0.003; FontS...