Introduction to computational physics

import numpy as np
import scipy as sp
# import astropy.units as u