projects / Analog_Engine.git / blame
| Commit | Line | Data |
|---|---|---|
| 6d26a6bd P |
1 | #PROGRAM-ID ChargeInEMfield |
| 2 | #VERSION 20240116 | |
| 3 | # A mass m with charge q>0 is located in a static electro-magnetic field [0,0,E], [0,-B,0]. | |
| 4 | # Find the trajetory | |
| 5 | # Differential equation (Lorentz force): | |
| 6 | # x'' = omega*z' | |
| 7 | # z'' = a*omega^2 - omega*x' | |
| 8 | # good start settings are omega=0,92, a*omega^2=0,05 | |
| 9 | ||
| 10 | # coefficient definition | |
| 11 | alias coefficient.1 -> omega_x # -omega for x' | |
| 12 | coefficient.2 (+1) -> a*omega^2 | |
| 13 | alias coefficient(3) -> omega_z # omega for z' | |
| 14 | ||
| 15 | # solving the DE | |
| 16 | iintegrate (x'') -> -x' | |
| 17 | iintegrate (-x') -> x | |
| 18 | ||
| 19 | coefficient.omega_x (-x') -> -omega*x' | |
| 20 | #cmultiply (-x', omega_x) -> omega*x' | |
| 21 | ||
| 22 | iintegrate (a*omega^2, -omega*x') -> -z' # input is z'' | |
| 23 | iintegrate (-z') -> z | |
| 24 | invert (-z') -> z' | |
| 25 | ||
| 26 | coefficient.omega_z (z') -> omega*z' | |
| 27 | #cmultiply (z', omega_z) -> omega*z' | |
| 28 | assign omega*z' -> x'' | |
| 29 | ||
| 30 | output(x) -> out.x | |
| 31 | output(z) -> out.z |