-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathaux_analysisOscDelta.py
32 lines (27 loc) · 46 KB
/
aux_analysisOscDelta.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
import utils
# 1. Turn off plotting in config --> plot: false
# 2. Run multiple simulations with different dts and all algos
# ./target/tp4-simu-1.0/damped-osc.sh -Dalgo=Beeman -Ddt=0.0001 --> BEEMAN-0.0001.txt
# python analysisOsc.py BEEMAN-0.0001.txt BEEMAN-0.001.txt BEEMAN-0.01.txt
# 3. Replace delta_t values and ecms below with obtained values
# 4. Run python aux_analysisOscDelta.py
# Different delta t step for each fragment
# ./multipleDtOsc.sh 0.00001 0.00001 0.00010; ./multipleDtOsc.sh 0.00011 0.00002 0.00505;
# ./multipleDtOsc.sh 0.00506 0.00005 0.0200; ./multipleDtOsc.sh 0.0201 0.0005 0.0601
# Replace delta_t values here (delta_t_sim = delta_t_print)
dts = [1.0000000000E-05, 2.0000000000E-05, 3.0000000000E-05, 4.0000000000E-05, 5.0000000000E-05, 6.0000000000E-05, 7.0000000000E-05, 8.0000000000E-05, 9.0000000000E-05, 1.0000000000E-04, 1.1000000000E-04, 1.3000000000E-04, 1.5000000000E-04, 1.7000000000E-04, 1.9000000000E-04, 2.1000000000E-04, 2.3000000000E-04, 2.5000000000E-04, 2.7000000000E-04, 2.9000000000E-04, 3.1000000000E-04, 3.3000000000E-04, 3.5000000000E-04, 3.7000000000E-04, 3.9000000000E-04, 4.1000000000E-04, 4.3000000000E-04, 4.5000000000E-04, 4.7000000000E-04, 4.9000000000E-04, 5.1000000000E-04, 5.3000000000E-04, 5.5000000000E-04, 5.7000000000E-04, 5.9000000000E-04, 6.1000000000E-04, 6.3000000000E-04, 6.5000000000E-04, 6.7000000000E-04, 6.9000000000E-04, 7.1000000000E-04, 7.3000000000E-04, 7.5000000000E-04, 7.7000000000E-04, 7.9000000000E-04, 8.1000000000E-04, 8.3000000000E-04, 8.5000000000E-04, 8.7000000000E-04, 8.9000000000E-04, 9.1000000000E-04, 9.3000000000E-04, 9.5000000000E-04, 9.7000000000E-04, 9.9000000000E-04, 1.0100000000E-03, 1.0300000000E-03, 1.0500000000E-03, 1.0700000000E-03, 1.0900000000E-03, 1.1100000000E-03, 1.1300000000E-03, 1.1500000000E-03, 1.1700000000E-03, 1.1900000000E-03, 1.2100000000E-03, 1.2300000000E-03, 1.2500000000E-03, 1.2700000000E-03, 1.2900000000E-03, 1.3100000000E-03, 1.3300000000E-03, 1.3500000000E-03, 1.3700000000E-03, 1.3900000000E-03, 1.4100000000E-03, 1.4300000000E-03, 1.4500000000E-03, 1.4700000000E-03, 1.4900000000E-03, 1.5100000000E-03, 1.5300000000E-03, 1.5500000000E-03, 1.5700000000E-03, 1.5900000000E-03, 1.6100000000E-03, 1.6300000000E-03, 1.6500000000E-03, 1.6700000000E-03, 1.6900000000E-03, 1.7100000000E-03, 1.7300000000E-03, 1.7500000000E-03, 1.7700000000E-03, 1.7900000000E-03, 1.8100000000E-03, 1.8300000000E-03, 1.8500000000E-03, 1.8700000000E-03, 1.8900000000E-03, 1.9100000000E-03, 1.9300000000E-03, 1.9500000000E-03, 1.9700000000E-03, 1.9900000000E-03, 2.0100000000E-03, 2.0300000000E-03, 2.0500000000E-03, 2.0700000000E-03, 2.0900000000E-03, 2.1100000000E-03, 2.1300000000E-03, 2.1500000000E-03, 2.1700000000E-03, 2.1900000000E-03, 2.2100000000E-03, 2.2300000000E-03, 2.2500000000E-03, 2.2700000000E-03, 2.2900000000E-03, 2.3100000000E-03, 2.3300000000E-03, 2.3500000000E-03, 2.3700000000E-03, 2.3900000000E-03, 2.4100000000E-03, 2.4300000000E-03, 2.4500000000E-03, 2.4700000000E-03, 2.4900000000E-03, 2.5100000000E-03, 2.5300000000E-03, 2.5500000000E-03, 2.5700000000E-03, 2.5900000000E-03, 2.6100000000E-03, 2.6300000000E-03, 2.6500000000E-03, 2.6700000000E-03, 2.6900000000E-03, 2.7100000000E-03, 2.7300000000E-03, 2.7500000000E-03, 2.7700000000E-03, 2.7900000000E-03, 2.8100000000E-03, 2.8300000000E-03, 2.8500000000E-03, 2.8700000000E-03, 2.8900000000E-03, 2.9100000000E-03, 2.9300000000E-03, 2.9500000000E-03, 2.9700000000E-03, 2.9900000000E-03, 3.0100000000E-03, 3.0300000000E-03, 3.0500000000E-03, 3.0700000000E-03, 3.0900000000E-03, 3.1100000000E-03, 3.1300000000E-03, 3.1500000000E-03, 3.1700000000E-03, 3.1900000000E-03, 3.2100000000E-03, 3.2300000000E-03, 3.2500000000E-03, 3.2700000000E-03, 3.2900000000E-03, 3.3100000000E-03, 3.3300000000E-03, 3.3500000000E-03, 3.3700000000E-03, 3.3900000000E-03, 3.4100000000E-03, 3.4300000000E-03, 3.4500000000E-03, 3.4700000000E-03, 3.4900000000E-03, 3.5100000000E-03, 3.5300000000E-03, 3.5500000000E-03, 3.5700000000E-03, 3.5900000000E-03, 3.6100000000E-03, 3.6300000000E-03, 3.6500000000E-03, 3.6700000000E-03, 3.6900000000E-03, 3.7100000000E-03, 3.7300000000E-03, 3.7500000000E-03, 3.7700000000E-03, 3.7900000000E-03, 3.8100000000E-03, 3.8300000000E-03, 3.8500000000E-03, 3.8700000000E-03, 3.8900000000E-03, 3.9100000000E-03, 3.9300000000E-03, 3.9500000000E-03, 3.9700000000E-03, 3.9900000000E-03, 4.0100000000E-03, 4.0300000000E-03, 4.0500000000E-03, 4.0700000000E-03, 4.0900000000E-03, 4.1100000000E-03, 4.1300000000E-03, 4.1500000000E-03, 4.1700000000E-03, 4.1900000000E-03, 4.2100000000E-03, 4.2300000000E-03, 4.2500000000E-03, 4.2700000000E-03, 4.2900000000E-03, 4.3100000000E-03, 4.3300000000E-03, 4.3500000000E-03, 4.3700000000E-03, 4.3900000000E-03, 4.4100000000E-03, 4.4300000000E-03, 4.4500000000E-03, 4.4700000000E-03, 4.4900000000E-03, 4.5100000000E-03, 4.5300000000E-03, 4.5500000000E-03, 4.5700000000E-03, 4.5900000000E-03, 4.6100000000E-03, 4.6300000000E-03, 4.6500000000E-03, 4.6700000000E-03, 4.6900000000E-03, 4.7100000000E-03, 4.7300000000E-03, 4.7500000000E-03, 4.7700000000E-03, 4.7900000000E-03, 4.8100000000E-03, 4.8300000000E-03, 4.8500000000E-03, 4.8700000000E-03, 4.8900000000E-03, 4.9100000000E-03, 4.9300000000E-03, 4.9500000000E-03, 4.9700000000E-03, 4.9900000000E-03, 5.0100000000E-03, 5.0300000000E-03, 5.0500000000E-03, 5.0600000000E-03, 5.1100000000E-03, 5.1600000000E-03, 5.2100000000E-03, 5.2600000000E-03, 5.3100000000E-03, 5.3600000000E-03, 5.4100000000E-03, 5.4600000000E-03, 5.5100000000E-03, 5.5600000000E-03, 5.6100000000E-03, 5.6600000000E-03, 5.7100000000E-03, 5.7600000000E-03, 5.8100000000E-03, 5.8600000000E-03, 5.9100000000E-03, 5.9600000000E-03, 6.0100000000E-03, 6.0600000000E-03, 6.1100000000E-03, 6.1600000000E-03, 6.2100000000E-03, 6.2600000000E-03, 6.3100000000E-03, 6.3600000000E-03, 6.4100000000E-03, 6.4600000000E-03, 6.5100000000E-03, 6.5600000000E-03, 6.6100000000E-03, 6.6600000000E-03, 6.7100000000E-03, 6.7600000000E-03, 6.8100000000E-03, 6.8600000000E-03, 6.9100000000E-03, 6.9600000000E-03, 7.0100000000E-03, 7.0600000000E-03, 7.1100000000E-03, 7.1600000000E-03, 7.2100000000E-03, 7.2600000000E-03, 7.3100000000E-03, 7.3600000000E-03, 7.4100000000E-03, 7.4600000000E-03, 7.5100000000E-03, 7.5600000000E-03, 7.6100000000E-03, 7.6600000000E-03, 7.7100000000E-03, 7.7600000000E-03, 7.8100000000E-03, 7.8600000000E-03, 7.9100000000E-03, 7.9600000000E-03, 8.0100000000E-03, 8.0600000000E-03, 8.1100000000E-03, 8.1600000000E-03, 8.2100000000E-03, 8.2600000000E-03, 8.3100000000E-03, 8.3600000000E-03, 8.4100000000E-03, 8.4600000000E-03, 8.5100000000E-03, 8.5600000000E-03, 8.6100000000E-03, 8.6600000000E-03, 8.7100000000E-03, 8.7600000000E-03, 8.8100000000E-03, 8.8600000000E-03, 8.9100000000E-03, 8.9600000000E-03, 9.0100000000E-03, 9.0600000000E-03, 9.1100000000E-03, 9.1600000000E-03, 9.2100000000E-03, 9.2600000000E-03, 9.3100000000E-03, 9.3600000000E-03, 9.4100000000E-03, 9.4600000000E-03, 9.5100000000E-03, 9.5600000000E-03, 9.6100000000E-03, 9.6600000000E-03, 9.7100000000E-03, 9.7600000000E-03, 9.8100000000E-03, 9.8600000000E-03, 9.9100000000E-03, 9.9600000000E-03, 1.0010000000E-02, 1.0060000000E-02, 1.0110000000E-02, 1.0160000000E-02, 1.0210000000E-02, 1.0260000000E-02, 1.0310000000E-02, 1.0360000000E-02, 1.0410000000E-02, 1.0460000000E-02, 1.0510000000E-02, 1.0560000000E-02, 1.0610000000E-02, 1.0660000000E-02, 1.0710000000E-02, 1.0760000000E-02, 1.0810000000E-02, 1.0860000000E-02, 1.0910000000E-02, 1.0960000000E-02, 1.1010000000E-02, 1.1060000000E-02, 1.1110000000E-02, 1.1160000000E-02, 1.1210000000E-02, 1.1260000000E-02, 1.1310000000E-02, 1.1360000000E-02, 1.1410000000E-02, 1.1460000000E-02, 1.1510000000E-02, 1.1560000000E-02, 1.1610000000E-02, 1.1660000000E-02, 1.1710000000E-02, 1.1760000000E-02, 1.1810000000E-02, 1.1860000000E-02, 1.1910000000E-02, 1.1960000000E-02, 1.2010000000E-02, 1.2060000000E-02, 1.2110000000E-02, 1.2160000000E-02, 1.2210000000E-02, 1.2260000000E-02, 1.2310000000E-02, 1.2360000000E-02, 1.2410000000E-02, 1.2460000000E-02, 1.2510000000E-02, 1.2560000000E-02, 1.2610000000E-02, 1.2660000000E-02, 1.2710000000E-02, 1.2760000000E-02, 1.2810000000E-02, 1.2860000000E-02, 1.2910000000E-02, 1.2960000000E-02, 1.3010000000E-02, 1.3060000000E-02, 1.3110000000E-02, 1.3160000000E-02, 1.3210000000E-02, 1.3260000000E-02, 1.3310000000E-02, 1.3360000000E-02, 1.3410000000E-02, 1.3460000000E-02, 1.3510000000E-02, 1.3560000000E-02, 1.3610000000E-02, 1.3660000000E-02, 1.3710000000E-02, 1.3760000000E-02, 1.3810000000E-02, 1.3860000000E-02, 1.3910000000E-02, 1.3960000000E-02, 1.4010000000E-02, 1.4060000000E-02, 1.4110000000E-02, 1.4160000000E-02, 1.4210000000E-02, 1.4260000000E-02, 1.4310000000E-02, 1.4360000000E-02, 1.4410000000E-02, 1.4460000000E-02, 1.4510000000E-02, 1.4560000000E-02, 1.4610000000E-02, 1.4660000000E-02, 1.4710000000E-02, 1.4760000000E-02, 1.4810000000E-02, 1.4860000000E-02, 1.4910000000E-02, 1.4960000000E-02, 1.5010000000E-02, 1.5060000000E-02, 1.5110000000E-02, 1.5160000000E-02, 1.5210000000E-02, 1.5260000000E-02, 1.5310000000E-02, 1.5360000000E-02, 1.5410000000E-02, 1.5460000000E-02, 1.5510000000E-02, 1.5560000000E-02, 1.5610000000E-02, 1.5660000000E-02, 1.5710000000E-02, 1.5760000000E-02, 1.5810000000E-02, 1.5860000000E-02, 1.5910000000E-02, 1.5960000000E-02, 1.6010000000E-02, 1.6060000000E-02, 1.6110000000E-02, 1.6160000000E-02, 1.6210000000E-02, 1.6260000000E-02, 1.6310000000E-02, 1.6360000000E-02, 1.6410000000E-02, 1.6460000000E-02, 1.6510000000E-02, 1.6560000000E-02, 1.6610000000E-02, 1.6660000000E-02, 1.6710000000E-02, 1.6760000000E-02, 1.6810000000E-02, 1.6860000000E-02, 1.6910000000E-02, 1.6960000000E-02, 1.7010000000E-02, 1.7060000000E-02, 1.7110000000E-02, 1.7160000000E-02, 1.7210000000E-02, 1.7260000000E-02, 1.7310000000E-02, 1.7360000000E-02, 1.7410000000E-02, 1.7460000000E-02, 1.7510000000E-02, 1.7560000000E-02, 1.7610000000E-02, 1.7660000000E-02, 1.7710000000E-02, 1.7760000000E-02, 1.7810000000E-02, 1.7860000000E-02, 1.7910000000E-02, 1.7960000000E-02, 1.8010000000E-02, 1.8060000000E-02, 1.8110000000E-02, 1.8160000000E-02, 1.8210000000E-02, 1.8260000000E-02, 1.8310000000E-02, 1.8360000000E-02, 1.8410000000E-02, 1.8460000000E-02, 1.8510000000E-02, 1.8560000000E-02, 1.8610000000E-02, 1.8660000000E-02, 1.8710000000E-02, 1.8760000000E-02, 1.8810000000E-02, 1.8860000000E-02, 1.8910000000E-02, 1.8960000000E-02, 1.9010000000E-02, 1.9060000000E-02, 1.9110000000E-02, 1.9160000000E-02, 1.9210000000E-02, 1.9260000000E-02, 1.9310000000E-02, 1.9360000000E-02, 1.9410000000E-02, 1.9460000000E-02, 1.9510000000E-02, 1.9560000000E-02, 1.9610000000E-02, 1.9660000000E-02, 1.9710000000E-02, 1.9760000000E-02, 1.9810000000E-02, 1.9860000000E-02, 1.9910000000E-02, 1.9960000000E-02, 2.0100000000E-02, 2.0600000000E-02, 2.1100000000E-02,
2.1600000000E-02, 2.2100000000E-02, 2.2600000000E-02, 2.3100000000E-02, 2.3600000000E-02, 2.4100000000E-02, 2.4600000000E-02, 2.5100000000E-02, 2.5600000000E-02, 2.6100000000E-02, 2.6600000000E-02, 2.7100000000E-02, 2.7600000000E-02, 2.8100000000E-02, 2.8600000000E-02, 2.9100000000E-02, 2.9600000000E-02, 3.0100000000E-02, 3.0600000000E-02, 3.1100000000E-02, 3.1600000000E-02, 3.2100000000E-02, 3.2600000000E-02, 3.3100000000E-02, 3.3600000000E-02, 3.4100000000E-02, 3.4600000000E-02, 3.5100000000E-02, 3.5600000000E-02, 3.6100000000E-02, 3.6600000000E-02, 3.7100000000E-02, 3.7600000000E-02, 3.8100000000E-02, 3.8600000000E-02, 3.9100000000E-02, 3.9600000000E-02, 4.0100000000E-02, 4.0600000000E-02, 4.1100000000E-02, 4.1600000000E-02, 4.2100000000E-02, 4.2600000000E-02, 4.3100000000E-02, 4.3600000000E-02, 4.4100000000E-02, 4.4600000000E-02, 4.5100000000E-02, 4.5600000000E-02, 4.6100000000E-02, 4.6600000000E-02, 4.7100000000E-02, 4.7600000000E-02, 4.8100000000E-02, 4.8600000000E-02, 4.9100000000E-02, 4.9600000000E-02, 5.0100000000E-02, 5.0600000000E-02, 5.1100000000E-02, 5.1600000000E-02, 5.2100000000E-02, 5.2600000000E-02, 5.3100000000E-02, 5.3600000000E-02, 5.4100000000E-02, 5.4600000000E-02, 5.5100000000E-02, 5.5600000000E-02, 5.6100000000E-02, 5.6600000000E-02, 5.7100000000E-02, 5.7600000000E-02, 5.8100000000E-02, 5.8600000000E-02, 5.9100000000E-02, 5.9600000000E-02, 6.0100000000E-02]
# Replace ecm values for verlet here
ecms_verlet = [4.9396572876E-10, 1.9757797432E-09, 4.4453069136E-09, 7.9024484683E-09, 1.2347052264E-08, 1.7779073374E-08, 2.4198030246E-08, 3.1604445600E-08, 3.9997371743E-08, 4.9377787151E-08, 5.9744015718E-08, 8.3436956085E-08, 1.1107785148E-07, 1.4265887083E-07, 1.7818536328E-07, 2.1765177138E-07, 2.6106860194E-07, 3.0841861905E-07, 3.5970058288E-07, 4.1494907859E-07, 4.7410883138E-07, 5.3719557710E-07, 6.0424150519E-07, 6.7520524728E-07, 7.5011117888E-07, 8.2899126919E-07, 9.1175152207E-07, 9.9847604712E-07, 1.0891344155E-06, 1.1836811068E-06, 1.2821579683E-06, 1.3845872115E-06, 1.4909295538E-06, 1.6012096278E-06, 1.7153403180E-06, 1.8334917081E-06, 1.9554847916E-06, 2.0816634264E-06, 2.2113846671E-06, 2.3454170194E-06, 2.4831231837E-06, 2.6248127591E-06, 2.7701314994E-06, 2.9195436720E-06, 3.0732325828E-06, 3.2304227935E-06, 3.3918054480E-06, 3.5571239138E-06, 3.7260668417E-06, 3.8989570022E-06, 4.0755070661E-06, 4.2569563183E-06, 4.4415407907E-06, 4.6297186579E-06, 4.8221141246E-06, 5.0195379652E-06, 5.2197937631E-06, 5.4235585042E-06, 5.6317183181E-06, 5.8441212242E-06, 6.0592457726E-06, 6.2794887070E-06, 6.5033380945E-06, 6.7305184877E-06, 6.9623711046E-06, 7.1988042378E-06, 7.4378819989E-06, 7.6811140915E-06, 7.9283362498E-06, 8.1793610839E-06, 8.4339771583E-06, 8.6942524788E-06, 8.9554205875E-06, 9.2219131242E-06, 9.4937475695E-06, 9.7681881231E-06, 1.0044879308E-05, 1.0329403338E-05, 1.0615874670E-05, 1.0903864005E-05, 1.1199650205E-05, 1.1496539915E-05, 1.1797669893E-05, 1.2102924684E-05, 1.2412175012E-05, 1.2725277233E-05, 1.3046309622E-05, 1.3366781176E-05, 1.3690586036E-05, 1.4017518905E-05, 1.4352244813E-05, 1.4689981838E-05, 1.5030483245E-05, 1.5373485458E-05, 1.5724314026E-05, 1.6077439413E-05, 1.6432552467E-05, 1.6789327600E-05, 1.7153807389E-05, 1.7519647660E-05, 1.7893277741E-05, 1.8267926337E-05, 1.8650430120E-05, 1.9033567259E-05, 1.9416902123E-05, 1.9807905555E-05, 2.0206810472E-05, 2.0605442941E-05, 2.1011964286E-05, 2.1417687167E-05, 2.1822078681E-05, 2.2243449393E-05, 2.2653706539E-05, 2.3081376832E-05, 2.3506998071E-05, 2.3940487543E-05, 2.4371231786E-05, 2.4809696070E-05, 2.5244672331E-05, 2.5698899668E-05, 2.6137422540E-05, 2.6595596849E-05, 2.7049258715E-05, 2.7510623827E-05, 2.7979891277E-05, 2.8443612420E-05, 2.8914958692E-05, 2.9394118830E-05, 2.9881292784E-05, 3.0361617855E-05, 3.0849623959E-05, 3.1345497302E-05, 3.1833257135E-05, 3.2328503832E-05, 3.2831406785E-05, 3.3342144833E-05, 3.3860906942E-05, 3.4369741248E-05, 3.4886147400E-05, 3.5410293499E-05, 3.5942356792E-05, 3.6462699611E-05, 3.7010692114E-05, 3.7546435624E-05, 3.8089916654E-05, 3.8641305662E-05, 3.9178688080E-05, 3.9745931381E-05, 4.0298562716E-05, 4.0858826960E-05, 4.1426882658E-05, 4.2002896108E-05, 4.2562016284E-05, 4.3128417438E-05, 4.3728236163E-05, 4.4310124593E-05, 4.4899736966E-05, 4.5497229258E-05, 4.6074590456E-05, 4.6687747303E-05, 4.7279985120E-05, 4.7909438505E-05, 4.8517135087E-05, 4.9132506222E-05, 4.9755698309E-05, 5.0354674032E-05, 5.0993332594E-05, 5.1640270342E-05, 5.2261590578E-05, 5.2890244416E-05, 5.3526356154E-05, 5.4170054958E-05, 5.4821475660E-05, 5.5480759264E-05, 5.6148053143E-05, 5.6823510898E-05, 5.7468031896E-05, 5.8159538621E-05, 5.8819011144E-05, 5.9485822882E-05, 6.0160090850E-05, 6.0884673589E-05, 6.1574954562E-05, 6.2228872491E-05, 6.2934232816E-05, 6.3647710870E-05, 6.4369455573E-05, 6.5052315892E-05, 6.5790189328E-05, 6.6487885062E-05, 6.7192810944E-05, 6.7955472850E-05, 6.8675960309E-05, 6.9404027921E-05, 7.0139794321E-05, 7.0883382603E-05, 7.1634920510E-05, 7.2339087895E-05, 7.3105997670E-05, 7.3881244119E-05, 7.4606880085E-05, 7.5398236811E-05, 7.6138458543E-05, 7.6946446930E-05, 7.7701692450E-05, 7.8464192191E-05, 7.9234044594E-05, 8.0011350918E-05, 8.0796215819E-05, 8.1655130923E-05, 8.2389059604E-05, 8.3197268280E-05, 8.4013495368E-05, 8.4837866734E-05, 8.5670512296E-05, 8.6439139225E-05, 8.7287581220E-05, 8.8144684912E-05, 8.8934934147E-05, 8.9731991507E-05, 9.0613421403E-05, 9.1425264216E-05, 9.2323953635E-05, 9.3151002367E-05, 9.3985304074E-05, 9.4826948332E-05, 9.5676027025E-05, 9.6532634806E-05, 9.7396869471E-05, 9.8268832256E-05, 9.9148628028E-05, 1.0003636539E-04, 1.0093215667E-04, 1.0183611783E-04, 1.0274836819E-04, 1.0357404844E-04, 1.0450186141E-04, 1.0543831315E-04, 1.0628477753E-04, 1.0723742225E-04, 1.0809782569E-04, 1.0906709510E-04, 1.0994176039E-04, 1.1092811572E-04, 1.1181738828E-04, 1.1271353499E-04, 1.1372556991E-04, 1.1463685485E-04, 1.1555524785E-04, 1.1659405097E-04, 1.1752817144E-04, 1.1846965819E-04, 1.1941859424E-04, 1.2037506474E-04, 1.2133915721E-04, 1.2231096185E-04, 1.2329057157E-04, 1.2427808221E-04, 1.2477483111E-04, 1.2728902091E-04, 1.2972184171E-04, 1.3220029084E-04, 1.3472558343E-04, 1.3729902053E-04, 1.3992200637E-04, 1.4259605336E-04, 1.4516510017E-04, 1.4794119228E-04, 1.5060625783E-04, 1.5331918608E-04, 1.5608127384E-04, 1.5871374633E-04, 1.6157316808E-04, 1.6429626437E-04, 1.6725827532E-04, 1.7007661153E-04, 1.7294203151E-04, 1.7585569002E-04, 1.7881882033E-04, 1.8183273871E-04, 1.8467274254E-04, 1.8778622206E-04, 1.9071734510E-04, 1.9393552047E-04, 1.9696211920E-04, 2.0003544717E-04, 2.0315654556E-04, 2.0632651904E-04, 2.0954654287E-04, 2.1281786403E-04, 2.1585608154E-04, 2.1922664058E-04, 2.2235387752E-04, 2.2582828449E-04, 2.2904824965E-04, 2.3231308320E-04, 2.3594860771E-04, 2.3931235593E-04, 2.4272389001E-04, 2.4618421083E-04, 2.4969437614E-04, 2.5325550298E-04, 2.5649898741E-04, 2.6015699471E-04, 2.6386935966E-04, 2.6724408208E-04, 2.7105959240E-04, 2.7452386945E-04, 2.7844706035E-04, 2.8200429356E-04, 2.8560580511E-04, 2.8925230792E-04, 2.9339423301E-04, 2.9714101690E-04, 3.0093540591E-04, 3.0477827168E-04, 3.0867053125E-04, 3.1261315095E-04, 3.1660714774E-04, 3.2013932799E-04, 3.2422852372E-04, 3.2837207770E-04, 3.3257121508E-04, 3.3627300314E-04, 3.4057469172E-04, 3.4436214031E-04, 3.4877019093E-04, 3.5264579922E-04, 3.5716447444E-04, 3.6113101822E-04, 3.6576501430E-04, 3.6982564531E-04, 3.7393014583E-04, 3.7807910783E-04, 3.8294184232E-04, 3.8719132647E-04, 3.9148740412E-04, 3.9583076275E-04, 4.0022211606E-04, 4.0466220995E-04, 4.0915182731E-04, 4.1369179142E-04, 4.1828296812E-04, 4.2292626631E-04, 4.2762263696E-04, 4.3237307011E-04, 4.3636245243E-04, 4.4120865205E-04, 4.4611168631E-04, 4.5107270858E-04, 4.5522314936E-04, 4.6028611735E-04, 4.6541037334E-04, 4.6968511483E-04, 4.7491653468E-04, 4.7927453239E-04, 4.8461655711E-04, 4.8905958852E-04, 4.9451604692E-04, 4.9904603116E-04, 5.0462113383E-04, 5.0924022508E-04, 5.1493851939E-04, 5.1964917931E-04, 5.2440176149E-04, 5.3028050066E-04, 5.3512823681E-04, 5.4001917484E-04, 5.4608847469E-04, 5.5107883016E-04, 5.5611382859E-04, 5.6119396370E-04, 5.6631971778E-04, 5.7149156911E-04, 5.7794266649E-04, 5.8322279536E-04, 5.8855093467E-04, 5.9392766028E-04, 5.9935356560E-04, 6.0482926535E-04, 6.1035539869E-04, 6.1593263173E-04, 6.2156165950E-04, 6.2724320727E-04, 6.3297803120E-04, 6.3876691827E-04, 6.4461068538E-04, 6.5051017751E-04, 6.5497541400E-04, 6.6096560970E-04, 6.6701363277E-04, 6.7312046985E-04, 6.7928712662E-04, 6.8551461942E-04, 6.9019156650E-04, 6.9651614209E-04, 7.0290431340E-04, 7.0935719353E-04, 7.1417972320E-04, 7.2073450118E-04, 7.2735717276E-04, 7.3229160747E-04, 7.3901985273E-04, 7.4581948314E-04, 7.5086723752E-04, 7.5777680691E-04, 7.6476148739E-04, 7.6992428565E-04, 7.7702398771E-04, 7.8226538431E-04, 7.8948289882E-04, 7.9678229516E-04, 8.0214216802E-04, 8.0956559540E-04, 8.1500599713E-04, 8.2255703827E-04, 8.2807891873E-04, 8.3576146602E-04, 8.4136585206E-04, 8.4918407952E-04, 8.5487213843E-04, 8.6060158539E-04, 8.6860355358E-04, 8.7441781819E-04, 8.8256629239E-04, 8.8846692840E-04, 8.9676691854E-04, 9.0275581125E-04, 9.0878652283E-04, 9.1729168635E-04, 9.2341251312E-04, 9.2957581198E-04, 9.3829558162E-04, 9.4455127488E-04, 9.5084991972E-04, 9.5979484134E-04, 9.6618871177E-04, 9.7262593260E-04, 9.7910717318E-04, 9.8834333562E-04, 9.9492302202E-04, 1.0015470649E-03, 1.0082159919E-03, 1.0177623189E-03, 1.0245339536E-03, 1.0313508789E-03, 1.0382135464E-03, 1.0451223538E-03, 1.0550622796E-03, 1.0620794442E-03, 1.0691433504E-03, 1.0762543852E-03, 1.0834129041E-03, 1.0906192360E-03, 1.1010550667E-03, 1.1083778834E-03, 1.1157494926E-03, 1.1231702677E-03, 1.1306405697E-03, 1.1381607517E-03, 1.1457311613E-03, 1.1533521441E-03, 1.1610240465E-03, 1.1687472181E-03, 1.1765220144E-03, 1.1879005026E-03, 1.1958090457E-03, 1.2037711171E-03, 1.2117871732E-03, 1.2198576813E-03, 1.2279831214E-03, 1.2361639873E-03, 1.2444007878E-03, 1.2526940481E-03, 1.2610443104E-03, 1.2694521346E-03, 1.2779180990E-03, 1.2864428006E-03, 1.2950268557E-03, 1.3036708992E-03, 1.3123755850E-03, 1.3211415854E-03, 1.3257817931E-03, 1.3346246919E-03, 1.3435298488E-03, 1.3524980336E-03, 1.3615300384E-03, 1.3706266771E-03, 1.3797887826E-03, 1.3890172055E-03, 1.3983128110E-03, 1.4076764756E-03, 1.4171090834E-03, 1.4266115215E-03, 1.4314972735E-03, 1.4410784149E-03, 1.4507312686E-03, 1.4604568636E-03, 1.4702562253E-03, 1.4801303678E-03, 1.4900802846E-03, 1.4951212146E-03, 1.5051510924E-03, 1.5152592871E-03, 1.5254469396E-03, 1.5357151495E-03, 1.5460649609E-03, 1.5512360501E-03, 1.5616676713E-03, 1.5721840644E-03, 1.5827864010E-03, 1.5934757553E-03, 1.5987624052E-03, 1.6095350748E-03, 1.6203984421E-03, 1.6313536727E-03, 1.6424017707E-03, 1.6477905708E-03, 1.6589250870E-03, 1.6701567071E-03, 1.6814864416E-03, 1.6869753968E-03, 1.6983928705E-03, 1.7099130571E-03, 1.7215368987E-03, 1.7271179267E-03, 1.7388318891E-03, 1.7506544736E-03, 1.7625864333E-03, 1.7682508508E-03, 1.7802765250E-03, 1.7924168872E-03, 1.7981777579E-03, 1.8104117183E-03, 1.8227658694E-03, 1.8352406666E-03, 1.8410690117E-03, 1.8536431890E-03, 1.8663437584E-03, 1.8722579702E-03, 1.8850592542E-03, 1.8979928322E-03, 1.9039890176E-03, 1.9170253223E-03, 1.9301999558E-03, 1.9362736404E-03, 1.9495536996E-03, 1.9629782130E-03, 1.9691245539E-03, 1.9826578994E-03, 1.9963418426E-03, 2.0025558358E-03, 2.0163527581E-03, 2.0303063251E-03, 2.0365830301E-03, 2.0506545073E-03, 2.0570358355E-03, 2.0907160203E-03, 2.2062340592E-03,
2.3284426238E-03, 2.4579662552E-03, 2.5843715986E-03, 2.7176082899E-03, 2.8581684864E-03, 2.9931412122E-03, 3.1489122417E-03, 3.2982626690E-03, 3.4551077424E-03, 3.6199080454E-03, 3.7750180934E-03, 3.9558195961E-03, 4.1260167830E-03, 4.3249134022E-03, 4.5116255158E-03, 4.7068772163E-03, 4.9110575616E-03, 5.1246505639E-03, 5.3482583956E-03, 5.5826235079E-03, 5.7959875052E-03, 6.0511476326E-03, 6.2835040304E-03, 6.5620846899E-03, 6.8148742943E-03, 7.0783240356E-03, 7.3527100617E-03, 7.6383585743E-03, 7.9821186069E-03, 8.2937510971E-03, 8.5672192216E-03, 8.9026384182E-03, 9.2520830721E-03, 9.6163969471E-03, 9.9965838347E-03, 1.0329403807E-02, 1.0738444626E-02, 1.1166165749E-02, 1.1539241800E-02, 1.2000356675E-02, 1.2403082074E-02, 1.2901137736E-02, 1.3335212593E-02, 1.3874501641E-02, 1.4341994861E-02, 1.4827538015E-02, 1.5431322534E-02, 1.5954123466E-02, 1.6496541231E-02, 1.7058677783E-02, 1.7760786779E-02, 1.8366884273E-02, 1.8995295237E-02, 1.9646585139E-02, 2.0321423403E-02, 2.1020586539E-02, 2.1744951417E-02, 2.2495483521E-02, 2.3273223760E-02, 2.4079276215E-02, 2.4914798367E-02, 2.5780994860E-02, 2.6679115799E-02, 2.7610460827E-02, 2.8576390645E-02, 2.9578348098E-02, 3.0617891035E-02, 3.1443675754E-02, 3.2553677593E-02, 3.3705164204E-02, 3.4899709290E-02, 3.6139144329E-02, 3.7104901640E-02, 3.8433458305E-02, 3.9810923595E-02, 4.1239231387E-02, 4.2332989270E-02, 4.3868148078E-02, 4.5458358631E-02]
# Replace ecm values for beeman here
ecms_beeman = [3.2795070765E-18, 5.2735174172E-17, 2.6694602117E-16, 8.4374045508E-16, 2.0598316186E-15, 4.2716021947E-15, 7.9131125050E-15, 1.3499375635E-14, 2.1623478218E-14, 3.2957631050E-14, 4.8252784502E-14, 9.4129093887E-14, 1.6684749502E-13, 2.7525927547E-13, 4.2949472007E-13, 6.4093182900E-13, 9.2226325907E-13, 1.2873551452E-12, 1.7513767325E-12, 2.3309579239E-12, 3.0435327337E-12, 3.9081364088E-12, 4.9452527447E-12, 6.1760914752E-12, 7.6236309956E-12, 9.3122932838E-12, 1.1266387258E-11, 1.3513442453E-11, 1.6080976802E-11, 1.8997393025E-11, 2.2293637334E-11, 2.6001761166E-11, 3.0153988051E-11, 3.4784970880E-11, 3.9928188497E-11, 4.5624182615E-11, 5.1906685468E-11, 5.8824139915E-11, 6.6399311595E-11, 7.4696170763E-11, 8.3738631109E-11, 9.3580541489E-11, 1.0425462865E-10, 1.1582422443E-10, 1.2834604681E-10, 1.4183803250E-10, 1.5637968187E-10, 1.7201245085E-10, 1.8877403132E-10, 2.0673378975E-10, 2.2593188145E-10, 2.4649441635E-10, 2.6838144676E-10, 2.9167520608E-10, 3.1647449839E-10, 3.4289854661E-10, 3.7086521874E-10, 4.0047979923E-10, 4.3187199005E-10, 4.6510099642E-10, 5.0011383973E-10, 5.3717252084E-10, 5.7622622938E-10, 6.1731940979E-10, 6.6064601695E-10, 7.0627913182E-10, 7.5410703253E-10, 8.0435114200E-10, 8.5707483070E-10, 9.1233804744E-10, 9.7019671674E-10, 1.0309725797E-09, 1.0941914394E-09, 1.1604577177E-09, 1.2298751368E-09, 1.3021861752E-09, 1.3774167897E-09, 1.4564168178E-09, 1.5384788994E-09, 1.6236094145E-09, 1.7128240945E-09, 1.8052245849E-09, 1.9013734847E-09, 2.0013534776E-09, 2.1052444128E-09, 2.2131228147E-09, 2.3258089562E-09, 2.4419227911E-09, 2.5622311340E-09, 2.6867926431E-09, 2.8166087539E-09, 2.9508900547E-09, 3.0896830667E-09, 3.2330274122E-09, 3.3821485139E-09, 3.5360117080E-09, 3.6946390687E-09, 3.8580443262E-09, 4.0277146079E-09, 4.2023188302E-09, 4.3834955798E-09, 4.5697492868E-09, 4.7628845767E-09, 4.9612236088E-09, 5.1647249898E-09, 5.3754608382E-09, 5.5936819575E-09, 5.8172986097E-09, 6.0486973638E-09, 6.2855507183E-09, 6.5277553773E-09, 6.7808799935E-09, 7.0366456140E-09, 7.3038630007E-09, 7.5766963379E-09, 7.8584019794E-09, 8.1456772626E-09, 8.4420706006E-09, 8.7439502936E-09, 9.0592575393E-09, 9.3760095230E-09, 9.7067580804E-09, 1.0043068190E-08, 1.0389493234E-08, 1.0746330894E-08, 1.1108612620E-08, 1.1481482932E-08, 1.1865241848E-08, 1.2260204459E-08, 1.2660473723E-08, 1.3072087262E-08, 1.3495358787E-08, 1.3923610183E-08, 1.4363618810E-08, 1.4815692411E-08, 1.5280152886E-08, 1.5757337324E-08, 1.6239106112E-08, 1.6733639504E-08, 1.7241263196E-08, 1.7762317589E-08, 1.8287319716E-08, 1.8835918067E-08, 1.9388430152E-08, 1.9955033194E-08, 2.0536079984E-08, 2.1120144340E-08, 2.1730737826E-08, 2.2344226257E-08, 2.2972747555E-08, 2.3616659034E-08, 2.4276332752E-08, 2.4937631110E-08, 2.5614472377E-08, 2.6322673172E-08, 2.7032145361E-08, 2.7758245700E-08, 2.8501356079E-08, 2.9244170429E-08, 3.0021875738E-08, 3.0798906715E-08, 3.1612795846E-08, 3.2425570963E-08, 3.3256505445E-08, 3.4105994545E-08, 3.4952330850E-08, 3.5839428853E-08, 3.6746323197E-08, 3.7649146813E-08, 3.8571133834E-08, 3.9512665140E-08, 4.0474133950E-08, 4.1455947274E-08, 4.2458526904E-08, 4.3482310368E-08, 4.4527751252E-08, 4.5564497973E-08, 4.6653659342E-08, 4.7733191990E-08, 4.8834310369E-08, 4.9957418551E-08, 5.1138441881E-08, 5.2307814186E-08, 5.3462916213E-08, 5.4678290948E-08, 5.5917884937E-08, 5.7182191303E-08, 5.8429630082E-08, 5.9743585340E-08, 6.1039287668E-08, 6.2359236256E-08, 6.3750629218E-08, 6.5121553150E-08, 6.6518068113E-08, 6.7940634365E-08, 6.9389726273E-08, 7.0865833248E-08, 7.2314576315E-08, 7.3844686005E-08, 7.5403311714E-08, 7.6931681902E-08, 7.8547264379E-08, 8.0130565039E-08, 8.1805186009E-08, 8.3445308663E-08, 8.5113821495E-08, 8.6811180439E-08, 8.8537852447E-08, 9.0294316956E-08, 9.2155187899E-08, 9.3898610443E-08, 9.5747470372E-08, 9.7628186073E-08, 9.9541313200E-08, 1.0148742353E-07, 1.0338135835E-07, 1.0539289087E-07, 1.0743914176E-07, 1.0942848718E-07, 1.1144944208E-07, 1.1359840311E-07, 1.1568583666E-07, 1.1790706836E-07, 1.2006311415E-07, 1.2225320445E-07, 1.2447783751E-07, 1.2673752266E-07, 1.2903278130E-07, 1.3136414817E-07, 1.3373217229E-07, 1.3613741792E-07, 1.3858046496E-07, 1.4106190953E-07, 1.4358236423E-07, 1.4614245776E-07, 1.4860796020E-07, 1.5124594825E-07, 1.5392545666E-07, 1.5650322026E-07, 1.5926419494E-07, 1.6191862578E-07, 1.6476352813E-07, 1.6749672979E-07, 1.7042815620E-07, 1.7324233761E-07, 1.7609643420E-07, 1.7916051177E-07, 1.8209900277E-07, 1.8507902104E-07, 1.8828197854E-07, 1.9135008080E-07, 1.9446144289E-07, 1.9761663710E-07, 2.0081624763E-07, 2.0406087163E-07, 2.0735112009E-07, 2.1068761850E-07, 2.1407100766E-07, 2.1578049077E-07, 2.2450916240E-07, 2.3331022516E-07, 2.4240850907E-07, 2.5181337561E-07, 2.6153468795E-07, 2.7158289118E-07, 2.8196906154E-07, 2.9239039020E-07, 3.0347132750E-07, 3.1457998817E-07, 3.2604125221E-07, 3.3786601076E-07, 3.4967281204E-07, 3.6224032802E-07, 3.7477822664E-07, 3.8813421125E-07, 4.0144629360E-07, 4.1515334089E-07, 4.2926655737E-07, 4.4379771584E-07, 4.5875920280E-07, 4.7358993274E-07, 4.8942484209E-07, 5.0510707945E-07, 5.2186612681E-07, 5.3844667910E-07, 5.5548254378E-07, 5.7298556081E-07, 5.9096810507E-07, 6.0944314144E-07, 6.2842425661E-07, 6.4707730239E-07, 6.6707704917E-07, 6.8671290534E-07, 7.0778502751E-07, 7.2845201006E-07, 7.4963797037E-07, 7.7240642419E-07, 7.9470231704E-07, 8.1755603831E-07, 8.4098101612E-07, 8.6499124611E-07, 8.8960132724E-07, 9.1352222715E-07, 9.3932730132E-07, 9.6577724425E-07, 9.9144349887E-07, 1.0191758301E-06, 1.0460577976E-06, 1.0751351599E-06, 1.1032869263E-06, 1.1320717346E-06, 1.1615023138E-06, 1.1933991812E-06, 1.2242154292E-06, 1.2557206424E-06, 1.2879294771E-06, 1.3208571351E-06, 1.3545194077E-06, 1.3889327050E-06, 1.4218543566E-06, 1.4577423510E-06, 1.4944307514E-06, 1.5319386077E-06, 1.5677266934E-06, 1.6068378934E-06, 1.6441171144E-06, 1.6848991157E-06, 1.7237246252E-06, 1.7662504274E-06, 1.8066810965E-06, 1.8510294869E-06, 1.8931288872E-06, 1.9360495006E-06, 1.9798054596E-06, 2.0279105610E-06, 2.0734673973E-06, 2.1199070359E-06, 2.1672452274E-06, 2.2154981589E-06, 2.2646825139E-06, 2.3148155253E-06, 2.3659150195E-06, 2.4179994515E-06, 2.4710879296E-06, 2.5252002259E-06, 2.5803567702E-06, 2.6317071381E-06, 2.6888640603E-06, 2.7471247463E-06, 2.8065123877E-06, 2.8616371672E-06, 2.9231720502E-06, 2.9859019359E-06, 3.0439960034E-06, 3.1089935057E-06, 3.1691188988E-06, 3.2364654757E-06, 3.2986825265E-06, 3.3684656365E-06, 3.4328375314E-06, 3.5051509088E-06, 3.5717447450E-06, 3.6466883058E-06, 3.7155760298E-06, 3.7855776040E-06, 3.8645144221E-06, 3.9369186852E-06, 4.0104845548E-06, 4.0936462180E-06, 4.1697349630E-06, 4.2470365605E-06, 4.3255680999E-06, 4.4053467397E-06, 4.4863897847E-06, 4.5783739235E-06, 4.6621939134E-06, 4.7473387005E-06, 4.8338275608E-06, 4.9216801921E-06, 5.0109167656E-06, 5.1015579707E-06, 5.1936250563E-06, 5.2871398660E-06, 5.3821248669E-06, 5.4786031716E-06, 5.5765985526E-06, 5.6761354463E-06, 5.7772389460E-06, 5.8667101033E-06, 5.9706891811E-06, 6.0763031597E-06, 6.1835800863E-06, 6.2925486331E-06, 6.4032380022E-06, 6.5006292535E-06, 6.6144395620E-06, 6.7300524620E-06, 6.8474999309E-06, 6.9504541776E-06, 7.0712010458E-06, 7.1938738237E-06, 7.3011528154E-06, 7.4272633691E-06, 7.5553985559E-06, 7.6671362170E-06, 7.7988615363E-06, 7.9327170882E-06, 8.0490551559E-06, 8.1866684502E-06, 8.3061424934E-06, 8.4476170843E-06, 8.5914122669E-06, 8.7157381150E-06, 8.8635882839E-06, 8.9912137768E-06, 9.1432402127E-06, 9.2742347683E-06, 9.4305659710E-06, 9.5650013952E-06, 9.7257729437E-06, 9.8637244691E-06, 1.0003400584E-05, 1.0170625904E-05, 1.0313921983E-05, 1.0485938904E-05, 1.0632939453E-05, 1.0809907588E-05, 1.0960704047E-05, 1.1113329280E-05, 1.1297477688E-05, 1.1454020976E-05, 1.1612446373E-05, 1.1804101368E-05, 1.1966577069E-05, 1.2130986680E-05, 1.2330503109E-05, 1.2499108072E-05, 1.2669698673E-05, 1.2842300571E-05, 1.3052403253E-05, 1.3229388236E-05, 1.3408436965E-05, 1.3589573401E-05, 1.3810919524E-05, 1.3996651585E-05, 1.4184526841E-05, 1.4374568587E-05, 1.4566799375E-05, 1.4802765147E-05, 1.4999871585E-05, 1.5199228508E-05, 1.5400858907E-05, 1.5604785390E-05, 1.5811030265E-05, 1.6065683077E-05, 1.6277179167E-05, 1.6491065459E-05, 1.6707365757E-05, 1.6926103833E-05, 1.7147303483E-05, 1.7370988587E-05, 1.7597183163E-05, 1.7825911421E-05, 1.8057197812E-05, 1.8291067074E-05, 1.8582712963E-05, 1.8822610110E-05, 1.9065187392E-05, 1.9310472474E-05, 1.9558493454E-05, 1.9809278904E-05, 2.0062857897E-05, 2.0319260041E-05, 2.0578515507E-05, 2.0840655053E-05, 2.1105710046E-05, 2.1373712487E-05, 2.1644695022E-05, 2.1918690957E-05, 2.2195734268E-05, 2.2475859601E-05, 2.2759102277E-05, 2.2973349533E-05, 2.3261655158E-05, 2.3553154403E-05, 2.3847885571E-05, 2.4145887779E-05, 2.4447200939E-05, 2.4751865733E-05, 2.5059923568E-05, 2.5371416524E-05, 2.5686387291E-05, 2.6004879084E-05, 2.6326935545E-05, 2.6566010969E-05, 2.6893598497E-05, 2.7224855260E-05, 2.7559829390E-05, 2.7898569344E-05, 2.8241123716E-05, 2.8587541023E-05, 2.8842099139E-05, 2.9194339758E-05, 2.9550568547E-05, 2.9910838955E-05, 3.0275203892E-05, 3.0643715355E-05, 3.0911964213E-05, 3.1286577483E-05, 3.1665482900E-05, 3.2048737150E-05, 3.2436394995E-05, 3.2716525043E-05, 3.3110522682E-05, 3.3509087000E-05, 3.3912276704E-05, 3.4320146915E-05, 3.4612209611E-05, 3.5026705723E-05, 3.5446064936E-05, 3.5870344238E-05, 3.6172538064E-05, 3.6603646740E-05, 3.7039872103E-05, 3.7481271105E-05, 3.7793642110E-05, 3.8242106755E-05, 3.8695956285E-05, 3.9155244410E-05, 3.9477822255E-05, 3.9944450774E-05, 4.0416743560E-05, 4.0747803458E-05, 4.1227567300E-05, 4.1713230529E-05, 4.2204841623E-05, 4.2546140440E-05, 4.3045569108E-05, 4.3551192863E-05, 4.3901035324E-05, 4.4414654495E-05, 4.4934724839E-05, 4.5293155286E-05, 4.5821420831E-05, 4.6356401955E-05, 4.6723448635E-05, 4.7266847816E-05, 4.7817234144E-05, 4.8192914710E-05, 4.8751966443E-05, 4.9318281955E-05, 4.9702609286E-05, 5.0277863788E-05, 5.0860660343E-05, 5.1253649238E-05, 5.1845686615E-05, 5.2245901342E-05, 5.3783254621E-05, 5.9348659568E-05,
6.5411378148E-05, 7.2017745312E-05, 7.8877757000E-05, 8.6300228506E-05, 9.4329870042E-05, 1.0254987089E-04, 1.1188810620E-04, 1.2141290909E-04, 1.3163633990E-04, 1.4260711614E-04, 1.5362297288E-04, 1.6616427813E-04, 1.7872006599E-04, 1.9305216136E-04, 2.0733788310E-04, 2.2253255148E-04, 2.3868798329E-04, 2.5586221930E-04, 2.7412099395E-04, 2.9353891698E-04, 3.1236134643E-04, 3.3413105885E-04, 3.5516293300E-04, 3.7957795195E-04, 4.0304341893E-04, 4.2776203696E-04, 4.5378683124E-04, 4.8117387089E-04, 5.1323044248E-04, 5.4376498025E-04, 5.7214472594E-04, 6.0564851731E-04, 6.4088338143E-04, 6.7794830860E-04, 7.1695401112E-04, 7.5277230547E-04, 7.9547015300E-04, 8.4044190526E-04, 8.8138670491E-04, 9.3062549331E-04, 9.7530093396E-04, 1.0292383024E-03, 1.0779070810E-03, 1.1370562740E-03, 1.1900192604E-03, 1.2451782829E-03, 1.3126080738E-03, 1.3725622151E-03, 1.4349277812E-03, 1.4997674263E-03, 1.5798101652E-03, 1.6502592499E-03, 1.7234612479E-03, 1.7995016558E-03, 1.8784701087E-03, 1.9604615511E-03, 2.0455768915E-03, 2.1339233208E-03, 2.2256144675E-03, 2.3207705387E-03, 2.4195185737E-03, 2.5219929058E-03, 2.6283359061E-03, 2.7386990447E-03, 2.8532442581E-03, 2.9721455254E-03, 3.0955904340E-03, 3.1931854547E-03, 3.3242524457E-03, 3.4602718769E-03, 3.6014737690E-03, 3.7481143103E-03, 3.8617971815E-03, 4.0173540617E-03, 4.1788330282E-03, 4.3465771575E-03, 4.4743792720E-03, 4.6522504577E-03, 4.8369646140E-03]
# Replace ecm values for gpc5 here
ecms_gear_5 = [5.5722832263E-23, 1.4449721781E-23, 1.2396949503E-23, 5.6103335292E-24, 2.8635243417E-24, 2.8922602250E-24, 2.4833890144E-25, 1.6255239213E-24, 7.8826497460E-25, 1.4362645934E-24, 7.8862578437E-25, 2.9332919713E-25, 6.4858677940E-25, 7.1207465417E-25, 3.5827140689E-25, 1.6401053548E-25, 2.0672015965E-25, 7.5457823056E-26, 1.9128713284E-26, 6.7305569034E-26, 1.3791156588E-25, 5.6149710402E-26, 5.1186697313E-26, 1.0664058654E-25, 2.5170390387E-25, 2.4450812743E-26, 7.0898755755E-26, 1.8789051202E-25, 2.2575474584E-25, 1.5108899373E-25, 4.9010393762E-25, 7.1163583505E-25, 5.5398119939E-25, 1.0812647924E-24, 1.5086811273E-24, 2.0864409887E-24, 2.9063356566E-24, 3.9825419779E-24, 5.8584341709E-24, 6.7322533839E-24, 9.2651097468E-24, 1.3488036150E-23, 1.6180671801E-23, 2.1149539752E-23, 2.8038816451E-23, 3.6334856641E-23, 4.4212570812E-23, 5.8651955271E-23, 7.4324557750E-23, 9.1560071042E-23, 1.1459426239E-22, 1.4270036340E-22, 1.7687994489E-22, 2.1811314928E-22, 2.6789639242E-22, 3.2989637918E-22, 3.9662178660E-22, 4.8145292850E-22, 5.8796190170E-22, 7.0200764667E-22, 8.4288970960E-22, 1.0116803600E-21, 1.2076133076E-21, 1.4357460221E-21, 1.7010671148E-21, 2.0128167376E-21, 2.3671825109E-21, 2.7858475800E-21, 3.2663461820E-21, 3.8212984342E-21, 4.4627064340E-21, 5.1958392289E-21, 6.0439185470E-21, 6.9915394830E-21, 8.1066240839E-21, 9.3580471319E-21, 1.0761410122E-20, 1.2386904148E-20, 1.4213864189E-20, 1.6284131055E-20, 1.8624338886E-20, 2.1258216499E-20, 2.4223657741E-20, 2.7540250454E-20, 3.1284436327E-20, 3.5506606544E-20, 4.0176816431E-20, 4.5424762109E-20, 5.1289308523E-20, 5.7809350875E-20, 6.5044502851E-20, 7.3165393120E-20, 8.2144501989E-20, 9.2042710919E-20, 1.0310293044E-19, 1.1531263290E-19, 1.2879491183E-19, 1.4366598284E-19, 1.6009864873E-19, 1.7821532969E-19, 1.9809796845E-19, 2.2001218388E-19, 2.4418560735E-19, 2.7053491327E-19, 2.9947300043E-19, 3.3124131060E-19, 3.6606235832E-19, 4.0397763698E-19, 4.4565268252E-19, 4.9098426916E-19, 5.4025191223E-19, 5.9430781915E-19, 6.5280358250E-19, 7.1682559369E-19, 7.8626825207E-19, 8.6184850360E-19, 9.4374295182E-19, 1.0325520476E-18, 1.1286682466E-18, 1.2335604600E-18, 1.3459991413E-18, 1.4685122228E-18, 1.6006069001E-18, 1.7434953597E-18, 1.8980014033E-18, 2.0641198045E-18, 2.2435494069E-18, 2.4370950909E-18, 2.6460201042E-18, 2.8700023681E-18, 3.1113086600E-18, 3.3711777641E-18, 3.6491838086E-18, 3.9482350791E-18, 4.2693433575E-18, 4.6147357872E-18, 4.9854252110E-18, 5.3805962450E-18, 5.8047477485E-18, 6.2592601659E-18, 6.7463731584E-18, 7.2646198329E-18, 7.8229763670E-18, 8.4165080380E-18, 9.0510027338E-18, 9.7294512208E-18, 1.0449109208E-17, 1.1223062498E-17, 1.2043412275E-17, 1.2918629755E-17, 1.3851840589E-17, 1.4846374769E-17, 1.5898660663E-17, 1.7018726153E-17, 1.8219976705E-17, 1.9489221905E-17, 2.0839096465E-17, 2.2274488832E-17, 2.3787427587E-17, 2.5408011023E-17, 2.7115074871E-17, 2.8942180718E-17, 3.0865572722E-17, 3.2906226413E-17, 3.5069247080E-17, 3.7341854788E-17, 3.9771641678E-17, 4.2346224130E-17, 4.5046724789E-17, 4.7905830418E-17, 5.0930256086E-17, 5.4128532305E-17, 5.7511606813E-17, 6.1087838232E-17, 6.4867783473E-17, 6.8862025817E-17, 7.3037274819E-17, 7.7491439099E-17, 8.2143183155E-17, 8.7050986697E-17, 9.2227645530E-17, 9.7745613568E-17, 1.0350262729E-16, 1.0950099180E-16, 1.1588990348E-16, 1.2262038765E-16, 1.2970928561E-16, 1.3708529246E-16, 1.4493841953E-16, 1.5310451126E-16, 1.6169138494E-16, 1.7083190450E-16, 1.8032720238E-16, 1.9030531521E-16, 2.0078969819E-16, 2.1180169964E-16, 2.2336934571E-16, 2.3535346679E-16, 2.4809387776E-16, 2.6146921113E-16, 2.7531372149E-16, 2.9003105472E-16, 3.0525613511E-16, 3.2143926716E-16, 3.3817017671E-16, 3.5569713930E-16, 3.7405662094E-16, 3.9328319104E-16, 4.1341498233E-16, 4.3480765585E-16, 4.5655144310E-16, 4.7963829884E-16, 5.0379447236E-16, 5.2907005505E-16, 5.5550849214E-16, 5.8272042458E-16, 6.1161231657E-16, 6.4182114631E-16, 6.7288658482E-16, 7.0532798381E-16, 7.3976566753E-16, 7.7515364782E-16, 8.1271861745E-16, 8.5129624874E-16, 8.9155219838E-16, 9.3354737717E-16, 9.7735707548E-16, 1.0230506969E-15, 1.0707016829E-15, 1.1203847496E-15, 1.1721863013E-15, 1.2261833988E-15, 1.2824607842E-15, 1.3411139314E-15, 1.4022305923E-15, 1.4646969550E-15, 1.5309671031E-15, 1.5999967959E-15, 1.6704924314E-15, 1.7452834096E-15, 1.8216252044E-15, 1.9026262839E-15, 1.9852541070E-15, 2.0729320898E-15, 2.1623240158E-15, 2.2552451684E-15, 2.3538485457E-15, 2.4542971859E-15, 2.5586731912E-15, 2.6694540596E-15, 2.7822083685E-15, 2.8993281331E-15, 3.0209676703E-15, 3.1472872337E-15, 3.2784505180E-15, 3.4146307854E-15, 3.5559973322E-15, 3.7027401731E-15, 3.7781872915E-15, 4.1771235551E-15, 4.6102882873E-15, 5.0844307172E-15, 5.6030658769E-15, 6.1700089507E-15, 6.7893798153E-15, 7.4656407531E-15, 8.1955001897E-15, 8.9993943873E-15, 9.8656782912E-15, 1.0808167247E-14, 1.1832980532E-14, 1.2933361740E-14, 1.4141684951E-14, 1.5437296377E-14, 1.6859018426E-14, 1.8381278844E-14, 2.0029174067E-14, 2.1812164582E-14, 2.3740414359E-14, 2.5824808136E-14, 2.8045582492E-14, 3.0474754205E-14, 3.3059553868E-14, 3.5885591362E-14, 3.8888677373E-14, 4.2121593588E-14, 4.5600427872E-14, 4.9342313193E-14, 5.3365535202E-14, 5.7689568405E-14, 6.2259329622E-14, 6.7241278475E-14, 7.2500335610E-14, 7.8231580740E-14, 8.4274283154E-14, 9.0743880324E-14, 9.7791347775E-14, 1.0520860825E-13, 1.1314110818E-13, 1.2162166986E-13, 1.3068504674E-13, 1.4036817740E-13, 1.5050952094E-13, 1.6153411576E-13, 1.7330139118E-13, 1.8560471714E-13, 1.9897655667E-13, 2.1294254472E-13, 2.2811972082E-13, 2.4395282047E-13, 2.6079188739E-13, 2.7869505568E-13, 2.9814506116E-13, 3.1839345634E-13, 3.3990235654E-13, 3.6274361963E-13, 3.8699306267E-13, 4.1273058823E-13, 4.4004046723E-13, 4.6831307769E-13, 4.9898357117E-13, 5.3150264688E-13, 5.6597418612E-13, 6.0158269659E-13, 6.4021655819E-13, 6.8008870154E-13, 7.2334547579E-13, 7.6794414219E-13, 8.1633299153E-13, 8.6616811180E-13, 9.2025239218E-13, 9.7588607057E-13, 1.0346020194E-12, 1.0965555480E-12, 1.1638006891E-12, 1.2328375185E-12, 1.3056304415E-12, 1.3823663808E-12, 1.4632408030E-12, 1.5484584510E-12, 1.6382335839E-12, 1.7327905509E-12, 1.8323641832E-12, 1.9372004499E-12, 2.0475569013E-12, 2.1637033297E-12, 2.2819134941E-12, 2.4102212394E-12, 2.5451855956E-12, 2.6871298587E-12, 2.8312973114E-12, 2.9878634409E-12, 3.1524464521E-12, 3.3193442221E-12, 3.5006882326E-12, 3.6844338666E-12, 3.8841083072E-12, 4.0862428692E-12, 4.3059576590E-12, 4.5281529219E-12, 4.7697719748E-12, 5.0138423363E-12, 5.2793972346E-12, 5.5473130883E-12, 5.8277084823E-12, 6.1329205899E-12, 6.4403999722E-12, 6.7620274809E-12, 7.1123873562E-12, 7.4647549324E-12, 7.8331472343E-12, 8.2182284945E-12, 8.6206879782E-12, 9.0412409059E-12, 9.4998287375E-12, 9.9597700749E-12, 1.0440172515E-11, 1.0941870430E-11, 1.1465729814E-11, 1.2012649899E-11, 1.2583564010E-11, 1.3179441511E-11, 1.3801289014E-11, 1.4450151836E-11, 1.5127115222E-11, 1.5833306647E-11, 1.6569896855E-11, 1.7338101954E-11, 1.8099830368E-11, 1.8932893564E-11, 1.9801372327E-11, 2.0706674862E-11, 2.1650263698E-11, 2.2633657362E-11, 2.3605565138E-11, 2.4670263543E-11, 2.5779501966E-11, 2.6935038645E-11, 2.8074663903E-11, 2.9324491424E-11, 3.0626091901E-11, 3.1908008185E-11, 3.3314810480E-11, 3.4779481767E-11, 3.6219668612E-11, 3.7801661243E-11, 3.9448309661E-11, 4.1064396262E-11, 4.2841833821E-11, 4.4585137617E-11, 4.6502934740E-11, 4.8498331580E-11, 5.0450902780E-11, 5.2602687091E-11, 5.4706327812E-11, 5.7025876729E-11, 5.9291188567E-11, 6.1790680484E-11, 6.4228952884E-11, 6.6921421010E-11, 6.9544680932E-11, 7.2262984152E-11, 7.5265137273E-11, 7.8187638076E-11, 8.1418905758E-11, 8.4559659681E-11, 8.8036509279E-11, 9.1410513394E-11, 9.4903755019E-11, 9.8773795813E-11, 1.0252414529E-10, 1.0640573657E-10, 1.1071040217E-10, 1.1487525823E-10, 1.1918443189E-10, 1.2396940764E-10, 1.2859054361E-10, 1.3337024653E-10, 1.3831353535E-10, 1.4380904058E-10, 1.4910632156E-10, 1.5458315752E-10, 1.6024513900E-10, 1.6654975924E-10, 1.7261318163E-10, 1.7887966128E-10, 1.8535543157E-10, 1.9204688854E-10, 1.9951189620E-10, 2.0667212377E-10, 2.1406871040E-10, 2.2170879695E-10, 2.2959971331E-10, 2.3774898367E-10, 2.4686349462E-10, 2.5557617980E-10, 2.6457177784E-10, 2.7385870185E-10, 2.8344558967E-10, 2.9334131099E-10, 3.0355497220E-10, 3.1409592553E-10, 3.2497377404E-10, 3.3619838057E-10, 3.4777987201E-10, 3.6079062637E-10, 3.7315527371E-10, 3.8591039163E-10, 3.9906733018E-10, 4.1263774904E-10, 4.2663362445E-10, 4.4106725904E-10, 4.5595128883E-10, 4.7129869682E-10, 4.8712281867E-10, 5.0343735275E-10, 5.2025637265E-10, 5.3759433504E-10, 5.5546609150E-10, 5.7388689640E-10, 5.9287242083E-10, 6.1243876064E-10, 6.3062699296E-10, 6.5132854661E-10, 6.7265873104E-10, 6.9463537102E-10, 7.1727677876E-10, 7.4060176322E-10, 7.6462964545E-10, 7.8938026674E-10, 8.1487400139E-10, 8.4113176383E-10, 8.6817502130E-10, 8.9602580266E-10, 9.2169775446E-10, 9.5111078332E-10, 9.8139669619E-10, 1.0125799210E-09, 1.0446855119E-09, 1.0777391570E-09, 1.1117671780E-09, 1.1429868194E-09, 1.1788882198E-09, 1.2158422783E-09, 1.2538781760E-09, 1.2930257821E-09, 1.3333156409E-09, 1.3701223627E-09, 1.4125968150E-09, 1.4563048398E-09, 1.5012802440E-09, 1.5475575269E-09, 1.5896906068E-09, 1.6384461562E-09, 1.6886072210E-09, 1.7402118862E-09, 1.7932988817E-09, 1.8414384186E-09, 1.8973361984E-09, 1.9548339133E-09, 2.0139741014E-09, 2.0674644670E-09, 2.1297061815E-09, 2.1937201348E-09, 2.2595528316E-09, 2.3189210546E-09, 2.3881757383E-09, 2.4593922597E-09, 2.5326211042E-09, 2.5984405253E-09, 2.6754449914E-09, 2.7546193400E-09, 2.8256884937E-09, 2.9089143304E-09, 2.9944791858E-09, 3.0824411723E-09, 3.1610853763E-09, 3.2535179953E-09, 3.3485332215E-09, 3.4333308153E-09, 3.5331430286E-09, 3.6357365491E-09, 3.7271171399E-09, 3.8348578010E-09, 3.9455923961E-09, 4.0440113065E-09, 4.1602688425E-09, 4.2797474904E-09, 4.3856878235E-09, 4.5110928666E-09, 4.6399612091E-09, 4.7539362913E-09, 4.8891643004E-09, 5.0281130286E-09, 5.1506693245E-09, 5.2964431073E-09, 5.4250389616E-09, 5.8452776275E-09, 7.5720654817E-09,
9.7616485630E-09, 1.2527616358E-08, 1.5938760265E-08, 2.0193084634E-08, 2.5480733325E-08, 3.1881441455E-08, 3.9929277214E-08, 4.9595330384E-08, 6.1390403926E-08, 7.5742897239E-08, 9.2679808218E-08, 1.1364053592E-07, 1.3820817178E-07, 1.6851556881E-07, 2.0380391586E-07, 2.4582995887E-07, 2.9576984931E-07, 3.5499035200E-07, 4.2507807339E-07, 5.0786905449E-07, 6.0174477833E-07, 7.1586409252E-07, 8.4453034489E-07, 1.0007245376E-06, 1.1758837760E-06, 1.3790251851E-06, 1.6142527695E-06, 1.8862149844E-06, 2.2155544442E-06, 2.5802637162E-06, 2.9787710174E-06, 3.4578253704E-06, 4.0079393029E-06, 4.6388890879E-06, 5.3616439237E-06, 6.1415043420E-06, 7.0791279616E-06, 8.1495663796E-06, 9.2966244438E-06, 1.0676313223E-05, 1.2147795365E-05, 1.3918755683E-05, 1.5799528292E-05, 1.8063587341E-05, 2.0459712136E-05, 2.3146504701E-05, 2.6385153926E-05, 2.9792374160E-05, 3.3604108995E-05, 3.7865967167E-05, 4.3011662577E-05, 4.8386779485E-05, 5.4384475752E-05, 6.1073118058E-05, 6.8528214709E-05, 7.6833146248E-05, 8.6079984751E-05, 9.6370413475E-05, 1.0781670790E-04, 1.2054292726E-04, 1.3468575413E-04, 1.5039692322E-04, 1.6784053659E-04, 1.8720820996E-04, 2.0867868922E-04, 2.3256682908E-04, 2.5895491563E-04, 2.8603913485E-04, 3.2146165149E-04, 3.6907631914E-04, 4.6691156526E-04, 7.2276321831E-04, 1.9136051785E-03, 6.0957763698E-03, 2.1096535677E-02, 7.4352636394E-02, 3.3574731916E-01, 1.1290934621E+00, 3.6697690114E+00]
# Initialize plotting
utils.init_plotter()
utils.plot_multiple_values([dts, dts, dts], 'dt (s)', [ecms_beeman, ecms_verlet, ecms_gear_5], 'error cuadrático medio (m^2)', ['beeman', 'verlet', 'gear'], sci_x=True, log_x=True, log_y=True, legend_loc='lower right')
# Hold execution
utils.hold_execution()