Introduction to computational physics# 1import numpy as np 2import scipy as sp 3# import astropy.units as u