3

We can solve this problem using FDM and implicit method of integration differential equations as follows: Clear[gab, \[Rho]0, \[HBar], mb, r01, r02, drc, \[Omega]1, \[Omega]2, \ \[Omega]] drc = 2*r01; \[HBar] = 1; \[Lambda] = 10; k = I \[Lambda]/\[HBar]; \[Omega]1 = 50; \[Omega]2 = 55; dom = \[Omega]2 - \[Omega]1; \[Omega] = (\[Omega]1*\[Omega]2)/(\[Omega]1 +...


2

You can not solve mixed time dependent and stationary equations. You'd have to make the first equation time dependent. Something like this starts to time integrate, though I did not wait for it to finish: pde = {D[σ[t, x, y, z], t] == ϵ Laplacian[σ[t, x, y, z], {x, y, z}] + s - δ ϕ[t, x, y, z] - γ σ[t, x, y, z], D[ϕ[t, x, y, z], t] == λ Laplacian[ϕ[...


Only top voted, non community-wiki answers of a minimum length are eligible