def bpMagnitudeErrorEoM(G, vmini, nobs=70):
  """
  Calculate the end-of-mission photometric standard error in the BP band as a function of G and (V-I).
  Note: this refers to the integrated flux from the BP spectrophotometer. A margin of 20% is included.

  Parameters
  ----------

  G     - Value(s) of G-band magnitude.
  vmini - Value(s) of (V-I) colour.

  Keywords
  --------

  nobs  - Number of observations collected (default 70).

  Returns
  -------

  The BP band photometric standard error in units of magnitude.
  """
  return sqrt( (power(bpMagnitudeError(G, vmini)/_scienceMargin,2) +
      _eomCalibrationFloorBP*_eomCalibrationFloorBP)/nobs ) * _scienceMargin
Calculate the end-of-mission photometric standard error in the BP band as a function of G and (V-I).
  Note: this refers to the integrated flux from the BP spectrophotometer. A margin of 20% is included.

  Parameters
  ----------

  G     - Value(s) of G-band magnitude.
  vmini - Value(s) of (V-I) colour.

  Keywords
  --------

  nobs  - Number of observations collected (default 70).

  Returns
  -------

  The BP band photometric standard error in units of magnitude.