Some Basic Equations (See: X-Ray Data Booklet)

nForm[num_]:=Print[NumberForm[N[num,8],{10,7},NumberFormat->(SequenceForm[#1,"D",#3]&),ExponentFunction->(#&)]]

wx[yIn_]:=(3/(5Pi))NIntegrate[BesselK[5/3,yi],{yi,yIn,Infinity}]

dwx[yIn_]:=-(3/(5 Pi))BesselK[5/3,yIn]

LogLogPlot[wx[z],{z,10^-6,1}]

swx[yIn_]:=(3/(5 Pi))yIn^(2/3) NIntegrate[BesselK[5/3,yi],{yi,yIn,Infinity}]

dswx[yIn_]:=(3/(5 Pi)) yIn^(2/3) (-BesselK[5/3,yIn]+(2/(3yIn)) NIntegrate[BesselK[5/3,yi],{yi,yIn,Infinity}])

dswxLog[yIn_]:=dswx[yIn]yIn Log[10]

nForm[Map[swx,PowerRange[10^-7,10^-2,10]]]

nForm[Map[dswxLog,PowerRange[10^-7,10^-2,10]]]

nForm[Map[swx,Range[0.01,0.1,0.01]]]

nForm[Map[dswx,Range[0.01,0.1,0.01]]]

nForm[Map[swx,Range[0.1,1,0.02]]]

nForm[Map[dswx,Range[0.1,1,0.02]]]

ewx[yIn_]:=(3/(5 Pi))Exp[yIn] Sqrt[yIn]NIntegrate[BesselK[5/3,yi],{yi,yIn,Infinity}]

dewx[yIn_]:=(3/(5 Pi))Exp[yIn]Sqrt[yIn](-BesselK[5/3,yIn]+(1+1/(2yIn))NIntegrate[BesselK[5/3,yi],{yi,yIn,Infinity}])

nForm[Map[ewx,Range[1,3,0.05]]]

nForm[Map[dewx,Range[1,3,0.05]]]

nForm[Map[ewx,Range[3,10,0.2]]]

nForm[Map[dewx,Range[3,10,0.2]]]

nForm[Map[ewx,Range[10,40,1]]]

nForm[Map[dewx,Range[10,40,1]]]