curr_eng = 50   ;set 50 keV for electron flux computation

;set the three energy fitting parameters 
lowcut = 37.9659 & breakeng = 139.870 & highcut = 5000.  

p = 3.67 & q = 4.5     ;set low and high electron indices
flux = 0.87d35                ;total electron flux, i.e. a(2) * 1.0d35
energy_arr = [ curr_eng, lowcut, breakeng, highcut]
energy_arr = energy_arr / 511.0d0 + 1.0                  ; convert to the ? energies
brm_distrn, energy_arr[0], energy_arr[1], energy_arr[2],  energy_arr[3], p, q, fcn
fcn = flux * fcn / 511.0d0                    ;get the final result, i.e. electron flux for 10 keV electrons
print, fcn
end