# Integral of quantum wavefunction

**URL:** <https://discourse.julialang.org/t/integral-of-quantum-wavefunction/61104>\
**Category:** General Usage\
**Created:** [May 13, 2021, 8:07pm UTC](https://discourse.julialang.org/t/integral-of-quantum-wavefunction/61104 "2021-05-13T20:07:09Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![Rasoul\_Jafarzadeh](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/rasoul_jafarzadeh/32/9967_2.png) [@Rasoul\_Jafarzadeh](https://discourse.julialang.org/u/Rasoul_Jafarzadeh)\
**Post date:** [May 13, 2021, 8:07pm UTC](https://discourse.julialang.org/t/integral-of-quantum-wavefunction/61104/1 "2021-05-13T20:07:09Z")

</div>

Hi. Using QuantumOptics package my code is in this form  
using QuantumOptics  
b\_c = FockBasis(10); #Fockbasis (cavity field)  
b\_a = SpinBasis(1//2); #2-level basis (atom)  
a = destroy(bc) ⊗ one(ba);  
σ⁻ = one(bc) ⊗ sigmam(ba);  
H = Δ\*dagger(a)_a + g_(dagger(a)_σ⁻ + a_dagger(σ⁻));  
ψ₀ = coherentstate(bc, α) ⊗ spindown(ba);  
T = [0:0.1:max\_time;];  
tout, psi = timeevolution.schroedinger(T, ψ₀, H);  
becasue result of psi is numeracal,max\_time should be a fixed number. but I need to integrate on psi such as Integrate( \<psi(t’)| \sigma\_z |psi(t’)\> , t’=0:10). How can I do this? any help would be appreciated.

---

<div class="post-metadata">

**Author:** ![rafael.guerra](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/rafael.guerra/32/216610_2.png) [@rafael.guerra](https://discourse.julialang.org/u/rafael.guerra)\
**Post date:** [May 13, 2021, 9:54pm UTC](https://discourse.julialang.org/t/integral-of-quantum-wavefunction/61104/2 "2021-05-13T21:54:01Z")

</div>

Repeating same advice to you as [here](https://discourse.julialang.org/t/save-figures-as-pdf/57503/5) and [here](https://discourse.julialang.org/t/parallel-run-in-for-loop/58858/2).