# -*- coding: utf-8 -*-
"""
Created on Wed Oct 16 19:55:23 2019

@author: pc
"""

from scipy.constants import *

def E(n):
    return -13.6*1/n**2

def dE(n1, n2):
    return E(n1) - E(n2)

def l(E):
    return h*c/E

for n in range(2, 9):
    print("Lyman in eV: ", dE(n, 1), " ", n)
    print("Lyman in J: ", dE(n, 1)*1.602e-19, " ", n)
    print("Lyman corresponding wavelength in nm: ", l(dE(n, 1)*1.602e-19)*1e9)
print("limit Lyman in eV: ", E(1))
print("limit Lyman in J: ", E(1)*1.602e-19)

for n in range(3, 10):
    print("Balmer in eV: ", dE(n, 2), " ", n)
    print("Balmer in J: ", dE(n, 2)*1.602e-19, " ", n)
    print("Balmer corresponding wavelength in nm: ", l(dE(n, 2)*1.602e-19)*1e9)
print("limit Balmer in eV: ", E(2))
print("limit Balmer in J: ", E(2)*1.602e-19)

for n in range(4, 11):
    print("Paschen in eV: ", dE(n, 3), " ", n)
    print("Paschen in J: ", dE(n, 3)*1.602e-19, " ", n)
    print("Paschen corresponding wavelength in nm: ", l(dE(n, 3)*1.602e-19)*1e9)
print("limit Paschen in eV: ", E(3))
print("limit Paschen in J: ", E(3)*1.602e-19)