New answers tagged

7

We can use the original code after small correction as follows Clear[gab, \[Rho]0, \[HBar], \[Omega]0, mb, \[Tau], k, r0, t, s, drc] drc = 0; \[HBar] = 1; \[Omega]1 = 10; \[Omega]0 = 10; gab = 1; mb = 1; ma = 1; \[Rho]0 = 1; drc = 0; r0 = Sqrt[(\[HBar])/(ma*\[Omega]0)]; f[\[Tau]_] := (Sqrt[ 2] E^(I (\[Tau]) \[Omega]0) gab^2 \[Pi]^(3/ 2) r0^2 \[...


2

To have a coarse view of $c_1,c_2$ we can proceed as follows Clear["Global`*"] n = 100; tmax = 5; dT = tmax/n; rc = 1; \[HBar] = 1; \[Omega] = 10; gab = 1; mb = 1; ma = 1; \[Rho]0 = 1; r0 = 1; \[CapitalDelta] = 1; f[k_] := (E^(I k dT \[Omega]) gab^2 \[Pi]^(3/2) r0^2\[Rho]0)/(Sqrt[2] \[HBar]^2 (r0^2 + (I k dT \[HBar])/mb)^(3/2)); g[k_] := (Sqrt[2] E^...


Top 50 recent answers are included