We fitted the model for acceleration driven by radiation pressure described in the main text using a Markov chain Monte Carlo, solving for the differential equations $$a(r)=\frac{{{\rm{d}}}^{2}r}{{\rm{d}}{t}^{2}}=\left\{\begin{array}{ll}0 & {\rm{if}}\,0\le r < {r}...
importnumpyasnpimportpolarsaspldf=pl.DataFrame({"value":np.random.standard_normal(7)})defzscore(series,n=3):mean=series.mean()std=series.std()high=mean+n*stdlow=mean-n*stdseries=series.clip(low,high)minval=series.min()return(series[-1]-minval)/(series.max()-minval)print(df.with_colu...