New answers tagged computational-physics
1
Here is a working code and the result shown below, just analyze it. Note that the normalization for the probability density function $p(x)$ used here is $\int p(x) dx$=1.
import numpy as np
import matplotlib.pyplot as plt
def show_histogram(vals):
n, bins, patches = plt.hist(x=vals, bins='auto', density=True, color='#0504aa', alpha=0.7, rwidth=0.85)
...
0
What exactly do you think the formula
to_add=(rj-ri)*big_g*masses(i)*masses(j)/((abs(rj-ri))**3.0d0)
does, especially the denominator? For the correct physics it should be the third power of the Euclidean distance.
Top 50 recent answers are included
