import math u = 100000000000 # distance to detectors, which is chosen to be about twice b, b = distance to neutral surface s = .00011 # distance between apertures in meters d = "{:.2e}".format(u) # distance to gamma ray detectors D = 16000 # size of mirror, or distance between two gamma ray detectors in meters print('\n') print('This is for the Moon \n') print(s, ' = distance between apertures in meters \n') print(d, ' = distance to gamma ray detectors in meters \n') AU = u * .0000000000067 print(AU, 'AU = distance to gamma ray detectors in astronomical units \n') print(D, ' = size of mirror, or smallest distance between two gamma ray detectors in meters \n') d = u R = math.atan(s/d) * (3600*180) * 1/3.14159 print("R =", R, 'arcseconds, the angular distance between the two apertures in seconds of arc from distance = d \n') w = (R * D)/0.21 * .000001 # wavelength of required ray in meters print('wavelength =', w, 'meters \n') f = 3e+8/w # frequency in Hertz print('frequency =', f, 'Hz \n') keV = f * 4.14e-18 print(keV, 'keV') print(keV/1000, ' MeV')