https://github.com/bluescarni/e3bp
Code for the solution of the three-dimensional problem of two fixed centres
https://github.com/bluescarni/e3bp
Last synced: 2 months ago
JSON representation
Code for the solution of the three-dimensional problem of two fixed centres
- Host: GitHub
- URL: https://github.com/bluescarni/e3bp
- Owner: bluescarni
- License: gpl-3.0
- Created: 2015-10-23T14:45:55.000Z (over 9 years ago)
- Default Branch: master
- Last Pushed: 2015-10-23T14:57:45.000Z (over 9 years ago)
- Last Synced: 2025-02-13T00:16:28.171Z (4 months ago)
- Language: Python
- Size: 160 KB
- Stars: 2
- Watchers: 3
- Forks: 0
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
- License: LICENSE
Awesome Lists containing this project
README
# e3bp
Python code for the solution of the three-dimensional problem of two fixed centres
The code has the following prerequisites:
* the w_elliptic library, for the computation of Weierstrass elliptic and related functions,
https://github.com/bluescarni/w_elliptic,
* the mpmath library, http://mpmath.org/,
* the matplotlib/numpy stack.
Usage example:
```python
import e3bp
e = e3bp.e3bp(1.1,1.2,1.3,[.1,.3,.4,.1,-.1,.2])
```
This will initialise an object ``e`` corresponding to a two fixed centres problem with the following parameters:
* ``a = 1.1``,
* ``mu_1 = 1.2``,
* ``mu_2 = 1.3``,
* initial cartesian position vector ``[.1,.3,.4]``,
* initial cartesian velocity vector ``[.1,-.1,.2]``.
The ``e`` object can then be used as follows:
```python
e.xi_tau(0.1) # Compute the xi coordinate at tau = 0.1
e.eta_tau(0.2) # Compute the eta coordinate at tau = 0.2
e.ell_ham_state_tau(0.3) # Compute the Hamiltonian state at tau = 0.3
```