function F = gravitational_attraction(m1, m2, r)
%Function to the gravitational attraction between a masses m1 and m2 a
%distance r apart

%Gravitational constant
G = 6.674e-11;

%Force as per Newton's law of universal gravitation
F = (G*m1*m2)/(r^2);