#ifndef MUON_H #define MUON_H double get_range(double T, double rho); double get_E(double T, double x, double rho); double get_dEdx(double T, double rho); #endif