diff --git a/benchmarks/benchmark_MathematicsSuite.py b/benchmarks/benchmark_MathematicsSuite.py
index 737ca89..bbaf9b9 100644
--- a/benchmarks/benchmark_MathematicsSuite.py
+++ b/benchmarks/benchmark_MathematicsSuite.py
@@ -1,7 +1,7 @@
 """Functions from mathematics to benchmark with aerospeed velocity """
 
 import astropy.units as u
-from plasmapy.mathematics import (
+from plasmapy.formulary import (
     plasma_dispersion_func,
     plasma_dispersion_func_deriv,
     Fermi_integral
diff --git a/benchmarks/benchmark_PhysicsSuite.py b/benchmarks/benchmark_PhysicsSuite.py
index 3ec51a0..6fa47ae 100644
--- a/benchmarks/benchmark_PhysicsSuite.py
+++ b/benchmarks/benchmark_PhysicsSuite.py
@@ -3,21 +3,21 @@
 import astropy.constants as const
 import astropy.units as u
 from numpy import pi
-from plasmapy.physics import parameters
+from plasmapy.formulary import parameters
 
 #For dielectric class
-from plasmapy.physics import (
+from plasmapy.formulary import (
     cold_plasma_permittivity_SDP,
     cold_plasma_permittivity_LRP,
     permittivity_1D_Maxwellian
 )
 
 #For dimensionless class
-from plasmapy.physics.dimensionless import quantum_theta, beta
+from plasmapy.formulary.dimensionless import quantum_theta, beta
 
 
 #For distribution class
-from plasmapy.physics import (
+from plasmapy.formulary import (
     Maxwellian_1D,
     Maxwellian_velocity_2D,
     Maxwellian_velocity_3D,
@@ -29,7 +29,7 @@
 )
 
 #For parameters class
-from plasmapy.physics.parameters import (
+from plasmapy.formulary.parameters import (
     mass_density,
     Alfven_speed,
     ion_sound_speed,
@@ -50,7 +50,7 @@
 )
 
 #For quantum class
-from plasmapy.physics import(
+from plasmapy.formulary import(
     deBroglie_wavelength,
     thermal_deBroglie_wavelength,
     Fermi_energy,
@@ -60,7 +60,7 @@
 )
 
 #For relativity class
-from plasmapy.physics import Lorentz_factor
+from plasmapy.formulary import Lorentz_factor
 
 
 class dielectric: