Jovian
Sign In
3 years ago
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+1−2τ3)1/321/3⋅[1−τ3+(1−2τ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]}}