docs: def extract_parameters(Pressure, PressureErr, A, AErr, Gamma0, Gamma0Err, method="chang"):
      """
      Calculates the radius, mass and conversion factor and thier uncertainties.
      For values to be correct data must have been taken with feedback off and
      at pressures of around 1mbar (this is because the equations assume
      harmonic motion and at lower pressures the uncooled particle experiences
      anharmonic motion (due to exploring furthur outside the middle of the trap).
      When cooled the value of Gamma (the damping) is a combination of the
      enviromental damping and feedback damping and so is not the correct value
      for use in this equation (as it requires the enviromental damping).
      Environmental damping can be predicted though as A=const*Gamma0. By
      fitting to 1mbar data one can find the value of the const and therefore
      Gamma0 = A/const
  
      Parameters
      ----------
      Pressure : float
          Pressure in mbar when the data was taken
      PressureErr : float
          Error in the Pressure as a decimal (e.g. 15% error is 0.15) 
      A : float
          Fitting constant A
          A = γ**2*2*Γ_0*(K_b*T_0)/(π*m)
          where:
          γ = conversionFactor
          Γ_0 = Damping factor due to environment
          π = pi
      AErr : float
          Error in Fitting constant A
      Gamma0 : float
          The enviromental damping factor Gamma_0 = Γ_0
      Gamma0Err : float
          The error in the enviromental damping factor Gamma_0 = Γ_0
  
      Returns:
      Params : list
          [radius, mass, conversionFactor]
          The extracted parameters
      ParamsError : list
          [radiusError, massError, conversionFactorError]
          The error in the extracted parameters
      """
      Pressure = 100 * Pressure  # conversion to Pascals
  
      rho = 1800 # as quoted by Microspheres and Nanospheres  # kgm^3
      dm = 0.372e-9  # m O'Hanlon, 2003
      T0 = 300  # kelvin
      kB = Boltzmann  # m^2 kg s^-2 K-1
      eta = 18.27e-6  # Pa s, viscosity of air
  
      method = method.lower()
      if method == "rashid":
          radius = (0.619 * 9 * pi * eta * dm**2) / \
                   (_np.sqrt(2) * rho * kB * T0) * (Pressure/Gamma0)
          
      m_air = 4.81e-26 # molecular mass of air is 28.97 g/mol and Avogadro's Number 6.0221409^23
      if method == "chang":
          vbar = (8*kB*T0/(pi*m_air))**0.5
          radius = 16/(rho*pi*vbar)*(Pressure/Gamma0)/4 # CORRECTION FACTOR OF 4 APPLIED!!!!
      # see section 4.1.1 of Muddassar Rashid's 2016 Thesis for
      # derivation of this
      # see also page 132 of Jan Giesler's Thesis
      err_radius = radius * \
          _np.sqrt(((PressureErr * Pressure) / Pressure)
                   ** 2 + (Gamma0Err / Gamma0)**2)
      mass = rho * ((4 * pi * radius**3) / 3)
      err_mass = mass * 3 * err_radius / radius
      conversionFactor = _np.sqrt(A * mass / (4 * kB * T0 * Gamma0))
      err_conversionFactor = conversionFactor * \
          _np.sqrt((AErr / A)**2 + (err_mass / mass)
                   ** 2 + (Gamma0Err / Gamma0)**2)
      return [radius, mass, conversionFactor], [err_radius, err_mass, err_conversionFactor]