#include<iostream>
using namespace std;

/* Calcul du module de la force gravitationnelle exercee par une sphere */

double force(double M, double R, double m, double d){
  const double G = 6.67e-11;
  double res;

  if (d >= R)
    res = G*M*m/(d*d);
  else
    res = G*M*m*d/(R*R*R);
  return res;
}

int main(){
  double M, R, m, d;

  cout << "Valeur de M ? ";
  cin >> M;
  cout << "Valeur de R ? ";
  cin >> R;
  cout << "Valeur de m ? ";
  cin >> m;
  cout << "Valeur de d ? ";
  cin >> d;
  cout << "Force gravitationnelle : " << force(M,R,m,d) << endl;
  return 0;
}
