import jovian
%pylab inline
import numpy as np
#from sincfit import *
import matplotlib.pyplot as plt
import matplotlib as mpl
Populating the interactive namespace from numpy and matplotlib
#plt.style.use('e21')
from matplotlib.colors import LogNorm
#Determining Kappa, Janoschek PRB 2013
#Eq. 3:

κ\kappa = κGI\kappa_{GI} \cdot τ+(1τ3+12τ3)1/321/3[1τ3+(12τ3)1/6]\frac{\tau + (1-\tau^3+ \sqrt{1-2\tau^3} )^{1/3}}{\sqrt{ 2^{1/3}\cdot\left[ 1- \tau^3 + (1-2\tau^3)^{1/6} \right]}}