double calculateSectionalDensity(double weight, double diameter) => weight / pow(diameter, 2) / 7000;