
    \h                         S SK JrJrJr  S SKJr  S SKJr  S SKJ	r	  S SK
Jr  S SKJr  S SKJr  S SKJrJrJr  S S	KJrJrJr  S S
KJr  \\-  \-  rS rSS jrSS jrSS jrSS jrg)    )Rationaloopi)S)Symbol)exp)sqrt)	integrate)simplify)omegamx)psi_nE_ncoherent_state)hbarc            	         [         [        -  [        SS5      -  [        [         * [        S-  -  S-  5      -  [         [        -  [        SS5      -  [        S[         -  5      -  [        -  [        [         * [        S-  -  S-  5      -  [         [        -  [        SS5      -  S[         -  [        S-  -  S-
  -  [        S5      -  [        [         * [        S-  -  S-  5      -  [         [        -  [        SS5      -  [        [         S-  5      -  S[         -  [        S-  -  S[        -  -
  -  [        [         * [        S-  -  S-  5      -  S.n U  H1  n[        [        U[        [        [        5      X   -
  5      S:X  a  M1   e   g )N            )r   r   r   r   r   )
nur   r   r   r   r	   r   r   r   r   )Psins     W/var/www/auris/envauris/lib/python3.13/site-packages/sympy/physics/tests/test_qho_1d.pytest_wavefunctionr      sD   rEHQN"S"q!tQ%77rEHQN"T!B$Z/!3c2#1*a-6HHrEHQN"a"fq!tma&78a@3sQPQTzST}CUUrEHQN"T"Q$Z/1r6AqD=1q53HICQSPSVWYZVZPZ\]P]L^^	C aAu-671<<<     c           	          [        U S-   5       H8  n[        [        U[        SS5      S-  [        [        * [        45      S:X  a  M8   e   g )Nr   r   ranger
   r   r   r   r   is     r   	test_normr#      sA    1q5\q!Q*A-B3|<AAA r   c                     [        U S-   5       Ha  n[        US-   U S-   5       HH  n[        [        U[        SS5      [        U[        SS5      -  [        [        * [        45      S:X  a  MH   e   Mc     g )Nr   r   r   )r   r"   js      r   test_orthogonalityr&       ss    1q5\q1ua!e$AaAq!%1a"33a"b\CFGH H H % r   c                     [        U S-   5       H8  n[        U[        5      [        [        -  U[        R
                  -   -  :X  a  M8   e   g )Nr   )r    r   r   r   r   Halfr!   s     r   test_energiesr)   (   s9    1q5\1e}uAFF
 ;;;; r   c           	          [        S5      n[        U S-   5       HE  n[        [        U S-   5      [	        U S-   U5      -  5      [        U[	        X5      -  5      :X  a  ME   e   g )Nalphar   )r   r    r   r	   r   )r   r+   r"   s      r   test_coherent_stater,   -   s\     7OE1q5\QUnQUE&BBCxPUXfghXpPpGqqqq r   N)r   )
   ) sympy.core.numbersr   r   r   sympy.core.singletonr   sympy.core.symbolr   &sympy.functions.elementary.exponentialr   (sympy.functions.elementary.miscellaneousr	   sympy.integrals.integralsr
   sympy.simplify.simplifyr   	sympy.abcr   r   r   sympy.physics.qho_1dr   r   r   sympy.physics.quantum.constantsr   r   r   r#   r&   r)   r,    r   r   <module>r9      sV    1 1 " $ 6 9 / , ! ! ; ; 0Y=BH<
rr   