diff options
Diffstat (limited to 'xmds2/realistic_Rb_and_fields/RbPropEquations.cdata')
| -rw-r--r-- | xmds2/realistic_Rb_and_fields/RbPropEquations.cdata | 12 |
1 files changed, 12 insertions, 0 deletions
diff --git a/xmds2/realistic_Rb_and_fields/RbPropEquations.cdata b/xmds2/realistic_Rb_and_fields/RbPropEquations.cdata new file mode 100644 index 0000000..f098f3d --- /dev/null +++ b/xmds2/realistic_Rb_and_fields/RbPropEquations.cdata @@ -0,0 +1,12 @@ +<![CDATA[ +//---------------- RbPropEquations.cdata starts ------------------ +dEp1_dz = -((i*Ndens*(rt2*conj(r0309) + rt6*(conj(r0410) + rt2*conj(r0511)))*eta1)/(rt2*rt6)) - Lt[Ep1]; +dEm1_dz = -((i*Ndens*(rt2*rt6*conj(r0109) + rt6*conj(r0210) + rt2*conj(r0311))*eta1)/(rt2*rt6)) - Lt[Em1]; +dEp2_dz = (i*Ndens*(conj(r0709) + conj(r0810))*eta1)/rt6 - Lt[Ep2]; +dEm2_dz = -((i*Ndens*(conj(r0610) + conj(r0711))*eta1)/rt6) - Lt[Em2]; +dEp3_dz = -((i*Ndens*(rt2*rt3*conj(r0612) + rt3*conj(r0713) + conj(r0814))*eta2)/rt3) - Lt[Ep3]; +dEm3_dz = -((i*Ndens*(conj(r0614) + rt3*(conj(r0715) + rt2*conj(r0816)))*eta2)/rt3) - Lt[Em3]; +dEp4_dz = -((i*Ndens*(rt2*conj(r0212) + rt3*conj(r0313) + rt3*conj(r0414) + rt2*conj(r0515))*eta2)/rt3) - Lt[Ep4]; +dEm4_dz = (i*Ndens*(rt2*conj(r0113) + rt3*conj(r0214) + rt3*conj(r0315) + rt2*conj(r0416))*eta2)/rt3 - Lt[Em4]; +//---------------- RbPropEquations.cdata ends ------------------ +]]> |