EL/FL Calculations

SetDirectory["/Users/zachscofield/OneDrive - Wake Forest University/Research/Simpson's"];

<<"chifunctions5.m";

Variables

ClearAll;

xb = 1/10;

x0=(101*xb)/100

m=1;

t=0;

xc = ((-xb)/2)+Sqrt[1-(3*xb^2)/4];

f = -(1/x)*(x-xc)*(x-xb)*(x+xc+xb);

fp = D[f,x];

κb=(1/(2*xb))*(xc-xb)*(xc+(2*xb))

κc=(1/(2*xc))*(xc-xb)*((2*xc)+xb)

κN=(1/(2(xc+xb)))*((2*xc)+xb)*(xc+(2*xb))

xstar = (1/(2*κb))*Log[(x-xb)/(xc-xb)]-(1/(2*κc))*Log[-(x-xc)/(xc-xb)]+(1/(2*κN))*Log[(x+xc+xb)/(xc+(2*xb))]-(xc/(4*xb*κb))*Log[((2*xc)+xb)/(xc+(2*xb))]-(xb*xc)/(2*(xc-xb))*Log[(xb/xc)];

x1a=(((10^9)-1)*xc)/(10^9)

x1b=(9999*xc)/10000

Band 1

Length[chiL1]

ax = (11/10)*x0;

bx = (108/100)*x0;

cx = (106/100)*x0;

dx = (104/100)*x0;

ex = (102/100)*x0;

fx = (101/100)*x0;

gx = (1005/1000)*x0;

hx = (1003/1000)*x0;

ix = (1001/1000)*x0;

α1 = Table[{chi[ax]/.chiL1[[i]]},{i,1,991,1}];

β1 = Table[{chi'[ax]/.chiL1[[i]]},{i,1,991,1}];

α2 = Table[{chi[bx]/.chiL1[[i]]},{i,1,991,1}];

β2 = Table[{chi'[bx]/.chiL1[[i]]},{i,1,991,1}];

α3 = Table[{chi[cx]/.chiL1[[i]]},{i,1,991,1}];

β3 = Table[{chi'[cx]/.chiL1[[i]]},{i,1,991,1}];

α4 = Table[{chi[dx]/.chiL1[[i]]},{i,1,991,1}];

β4 = Table[{chi'[dx]/.chiL1[[i]]},{i,1,991,1}];

α5 = Table[{chi[ex]/.chiL1[[i]]},{i,1,991,1}];

β5 = Table[{chi'[ex]/.chiL1[[i]]},{i,1,991,1}];

α6 = Table[{chi[fx]/.chiL1[[i]]},{i,1,991,1}];

β6 = Table[{chi'[fx]/.chiL1[[i]]},{i,1,991,1}];

α7 = Table[{chi[gx]/.chiL1[[i]]},{i,1,991,1}];

β7 = Table[{chi'[gx]/.chiL1[[i]]},{i,1,991,1}];

α8 = Table[{chi[hx]/.chiL1[[i]]},{i,1,991,1}];

β8 = Table[{chi'[hx]/.chiL1[[i]]},{i,1,991,1}];

α9 = Table[{chi[ix]/.chiL1[[i]]},{i,1,991,1}];

β9 = Table[{chi'[ix]/.chiL1[[i]]},{i,1,991,1}];

α10 = Table[{chi[x0]/.chiL1[[i]]},{i,1,991,1}];

β10 = Table[{chi'[x0]/.chiL1[[i]]},{i,1,991,1}];

omegatable1 = Table[i,{i,1/10000,1/100,1/100000}];

EL1 = Table[{((α1[[i]]*I*omegatable1[[i]]+f*(β1[[i]]))/(2*I*omegatable1[[i]]*Exp[I*omegatable1[[i]]*xstar]))/.x-> ax},{i,1,991,1}];

FL1 = Table[{((α1[[i]]*I*omegatable1[[i]]-f*(β1[[i]]))/(2*I*omegatable1[[i]]*Exp[-I*omegatable1[[i]]*xstar]))/.x-> ax},{i,1,991,1}];

EL2 = Table[{((α2[[i]]*I*omegatable1[[i]]+f*(β2[[i]]))/(2*I*omegatable1[[i]]*Exp[I*omegatable1[[i]]*xstar]))/.x-> bx},{i,1,991,1}];

FL2 = Table[{((α2[[i]]*I*omegatable1[[i]]-f*(β2[[i]]))/(2*I*omegatable1[[i]]*Exp[-I*omegatable1[[i]]*xstar]))/.x-> bx},{i,1,991,1}];

EL3 = Table[{((α3[[i]]*I*omegatable1[[i]]+f*(β3[[i]]))/(2*I*omegatable1[[i]]*Exp[I*omegatable1[[i]]*xstar]))/.x-> cx},{i,1,991,1}];

FL3 = Table[{((α3[[i]]*I*omegatable1[[i]]-f*(β3[[i]]))/(2*I*omegatable1[[i]]*Exp[-I*omegatable1[[i]]*xstar]))/.x-> cx},{i,1,991,1}];

EL4 = Table[{((α4[[i]]*I*omegatable1[[i]]+f*(β4[[i]]))/(2*I*omegatable1[[i]]*Exp[I*omegatable1[[i]]*xstar]))/.x-> dx},{i,1,991,1}];

FL4 = Table[{((α4[[i]]*I*omegatable1[[i]]-f*(β4[[i]]))/(2*I*omegatable1[[i]]*Exp[-I*omegatable1[[i]]*xstar]))/.x-> dx},{i,1,991,1}];

EL5 = Table[{((α5[[i]]*I*omegatable1[[i]]+f*(β5[[i]]))/(2*I*omegatable1[[i]]*Exp[I*omegatable1[[i]]*xstar]))/.x-> ex},{i,1,991,1}];

FL5 = Table[{((α5[[i]]*I*omegatable1[[i]]-f*(β5[[i]]))/(2*I*omegatable1[[i]]*Exp[-I*omegatable1[[i]]*xstar]))/.x-> ex},{i,1,991,1}];

EL6 = Table[{((α6[[i]]*I*omegatable1[[i]]+f*(β6[[i]]))/(2*I*omegatable1[[i]]*Exp[I*omegatable1[[i]]*xstar]))/.x-> fx},{i,1,991,1}];

FL6 = Table[{((α6[[i]]*I*omegatable1[[i]]-f*(β6[[i]]))/(2*I*omegatable1[[i]]*Exp[-I*omegatable1[[i]]*xstar]))/.x-> fx},{i,1,991,1}];

EL7 = Table[{((α7[[i]]*I*omegatable1[[i]]+f*(β7[[i]]))/(2*I*omegatable1[[i]]*Exp[I*omegatable1[[i]]*xstar]))/.x-> gx},{i,1,991,1}];

FL7 = Table[{((α7[[i]]*I*omegatable1[[i]]-f*(β7[[i]]))/(2*I*omegatable1[[i]]*Exp[-I*omegatable1[[i]]*xstar]))/.x-> gx},{i,1,991,1}];

EL8 = Table[{((α8[[i]]*I*omegatable1[[i]]+f*(β8[[i]]))/(2*I*omegatable1[[i]]*Exp[I*omegatable1[[i]]*xstar]))/.x-> hx},{i,1,991,1}];

FL8 = Table[{((α8[[i]]*I*omegatable1[[i]]-f*(β8[[i]]))/(2*I*omegatable1[[i]]*Exp[-I*omegatable1[[i]]*xstar]))/.x-> hx},{i,1,991,1}];

EL9 = Table[{((α9[[i]]*I*omegatable1[[i]]+f*(β9[[i]]))/(2*I*omegatable1[[i]]*Exp[I*omegatable1[[i]]*xstar]))/.x-> ix},{i,1,991,1}];

FL9 = Table[{((α9[[i]]*I*omegatable1[[i]]-f*(β9[[i]]))/(2*I*omegatable1[[i]]*Exp[-I*omegatable1[[i]]*xstar]))/.x-> ix},{i,1,991,1}];

EL10 = Table[{((α10[[i]]*I*omegatable1[[i]]+f*(β10[[i]]))/(2*I*omegatable1[[i]]*Exp[I*omegatable1[[i]]*xstar]))/.x-> x0},{i,1,991,1}];

FL10 = Table[{((α10[[i]]*I*omegatable1[[i]]-f*(β10[[i]]))/(2*I*omegatable1[[i]]*Exp[-I*omegatable1[[i]]*xstar]))/.x-> x0},{i,1,991,1}];

Tce = Table[{{ax,Part[EL1[[i]],1]},{bx,Part[EL2[[i]],1]},{cx,Part[EL3[[i]],1]},{dx,Part[EL4[[i]],1]},{ex,Part[EL5[[i]],1]}, {fx, Part[EL6[[i]],1]}, {gx, Part[EL7[[i]],1]},{hx,Part[EL8[[i]],1]},{ix,Part[EL9[[i]],1]},{x0,Part[EL10[[i]],1]}},{i,1,991,1}];

TceCorrected= Tce/.{x_,{{{y_}}}}->{x,y};

Tre=Table[{Part[TceCorrected[[i]],j,1],Re[Part[TceCorrected[[i]],j,2]]},{i,1,991,1},{j,10}];

ELFitr = Table[Fit[Tre[[i]],{1,y,y^2,y^3},y],{i,1,991,1}];

Tie=Table[{Part[TceCorrected[[i]],j,1],Im[Part[TceCorrected[[i]],j,2]]},{i,1,991,1},{j,10}];

ELFiti = Table[Fit[Tie[[i]],{1,y,y^2,y^3},y],{i,1,991,1}];

Tcf = Table[{{ax,Part[FL1[[i]],1]},{bx,Part[FL2[[i]],1]},{cx,Part[FL3[[i]],1]},{dx,Part[FL4[[i]],1]},{ex,Part[FL5[[i]],1]}, {fx, Part[FL6[[i]],1]}, {gx, Part[FL7[[i]],1]},{hx,Part[FL8[[i]],1]},{ix,Part[FL9[[i]],1]},{x0,Part[FL10[[i]],1]}},{i,1,991,1}];

TcfCorrected= Tcf/.{x_,{{{y_}}}}->{x,y};

Trf=Table[{Part[TcfCorrected[[i]],j,1],Re[Part[TcfCorrected[[i]],j,2]]},{i,1,991,1},{j,10}];

FLFitr = Table[Fit[Trf[[i]],{1,y,y^2,y^3},y],{i,1,991,1}];

Tif=Table[{Part[TcfCorrected[[i]],j,1],Im[Part[TcfCorrected[[i]],j,2]]},{i,1,991,1},{j,10}];

FLFiti = Table[Fit[Tif[[i]],{1,y,y^2,y^3},y],{i,1,991,1}];

ReELfin = Table[{omegatable1[[i]],N[ELFitr[[i]]/.{y-> 1/10}]},{i,1,991,1}];

ImELfin = Table[{omegatable1[[i]],N[ELFiti[[i]]/.{y-> 1/10}]},{i,1,991,1}];

ReFLfin = Table[{omegatable1[[i]],N[FLFitr[[i]]/.{y-> 1/10}]},{i,1,991,1}];

ImFLfin = Table[{omegatable1[[i]],N[FLFiti[[i]]/.{y-> 1/10}]},{i,1,991,1}];

EL1rp = ListPlot[ReELfin]

EL1ip = ListPlot[ImELfin, PlotRange -> All]

FL1rp = ListPlot[ReFLfin]

FL1ip = ListPlot[ImFLfin, PlotRange-> All]

Band 2

Length[chiL2]

α1b = Table[{chi[ax]/.chiL2[[i]]},{i,1,991,1}];

β1b = Table[{chi'[ax]/.chiL2[[i]]},{i,1,991,1}];

α2b = Table[{chi[bx]/.chiL2[[i]]},{i,1,991,1}];

β2b = Table[{chi'[bx]/.chiL2[[i]]},{i,1,991,1}];

α3b = Table[{chi[cx]/.chiL2[[i]]},{i,1,991,1}];

β3b = Table[{chi'[cx]/.chiL2[[i]]},{i,1,991,1}];

α4b = Table[{chi[dx]/.chiL2[[i]]},{i,1,991,1}];

β4b = Table[{chi'[dx]/.chiL2[[i]]},{i,1,991,1}];

α5b = Table[{chi[ex]/.chiL2[[i]]},{i,1,991,1}];

β5b = Table[{chi'[ex]/.chiL2[[i]]},{i,1,991,1}];

α6b = Table[{chi[fx]/.chiL2[[i]]},{i,1,991,1}];

β6b = Table[{chi'[fx]/.chiL2[[i]]},{i,1,991,1}];

α7b = Table[{chi[gx]/.chiL2[[i]]},{i,1,991,1}];

β7b = Table[{chi'[gx]/.chiL2[[i]]},{i,1,991,1}];

α8b = Table[{chi[hx]/.chiL2[[i]]},{i,1,991,1}];

β8b = Table[{chi'[hx]/.chiL2[[i]]},{i,1,991,1}];

α9b = Table[{chi[ix]/.chiL2[[i]]},{i,1,991,1}];

β9b = Table[{chi'[ix]/.chiL2[[i]]},{i,1,991,1}];

α10b = Table[{chi[x0]/.chiL2[[i]]},{i,1,991,1}];

β10b = Table[{chi'[x0]/.chiL2[[i]]},{i,1,991,1}];

omegatable2 = Table[i,{i,1/100,1,1/1000}];

EL1b = Table[{((α1b[[i]]*I*omegatable2[[i]]+f*(β1b[[i]]))/(2*I*omegatable2[[i]]*Exp[I*omegatable2[[i]]*xstar]))/.x-> ax},{i,1,991,1}];

FL1b = Table[{((α1b[[i]]*I*omegatable2[[i]]-f*(β1b[[i]]))/(2*I*omegatable2[[i]]*Exp[-I*omegatable2[[i]]*xstar]))/.x-> ax},{i,1,991,1}];

EL2b = Table[{((α2b[[i]]*I*omegatable2[[i]]+f*(β2b[[i]]))/(2*I*omegatable2[[i]]*Exp[I*omegatable2[[i]]*xstar]))/.x-> bx},{i,1,991,1}];

FL2b= Table[{((α2b[[i]]*I*omegatable2[[i]]-f*(β2b[[i]]))/(2*I*omegatable2[[i]]*Exp[-I*omegatable2[[i]]*xstar]))/.x-> bx},{i,1,991,1}];

EL3b= Table[{((α3b[[i]]*I*omegatable2[[i]]+f*(β3b[[i]]))/(2*I*omegatable2[[i]]*Exp[I*omegatable2[[i]]*xstar]))/.x-> cx},{i,1,991,1}];

FL3b = Table[{((α3b[[i]]*I*omegatable2[[i]]-f*(β3b[[i]]))/(2*I*omegatable2[[i]]*Exp[-I*omegatable2[[i]]*xstar]))/.x-> cx},{i,1,991,1}];

EL4b= Table[{((α4b[[i]]*I*omegatable2[[i]]+f*(β4b[[i]]))/(2*I*omegatable2[[i]]*Exp[I*omegatable2[[i]]*xstar]))/.x-> dx},{i,1,991,1}];

FL4b = Table[{((α4b[[i]]*I*omegatable2[[i]]-f*(β4b[[i]]))/(2*I*omegatable2[[i]]*Exp[-I*omegatable2[[i]]*xstar]))/.x-> dx},{i,1,991,1}];

EL5b = Table[{((α5b[[i]]*I*omegatable2[[i]]+f*(β5b[[i]]))/(2*I*omegatable2[[i]]*Exp[I*omegatable2[[i]]*xstar]))/.x-> ex},{i,1,991,1}];

FL5b= Table[{((α5b[[i]]*I*omegatable2[[i]]-f*(β5b[[i]]))/(2*I*omegatable2[[i]]*Exp[-I*omegatable2[[i]]*xstar]))/.x-> ex},{i,1,991,1}];

EL6b = Table[{((α6b[[i]]*I*omegatable2[[i]]+f*(β6b[[i]]))/(2*I*omegatable2[[i]]*Exp[I*omegatable2[[i]]*xstar]))/.x-> fx},{i,1,991,1}];

FL6b = Table[{((α6b[[i]]*I*omegatable2[[i]]-f*(β6b[[i]]))/(2*I*omegatable2[[i]]*Exp[-I*omegatable2[[i]]*xstar]))/.x-> fx},{i,1,991,1}];

EL7b = Table[{((α7b[[i]]*I*omegatable2[[i]]+f*(β7b[[i]]))/(2*I*omegatable2[[i]]*Exp[I*omegatable2[[i]]*xstar]))/.x-> gx},{i,1,991,1}];

FL7b = Table[{((α7b[[i]]*I*omegatable2[[i]]-f*(β7b[[i]]))/(2*I*omegatable2[[i]]*Exp[-I*omegatable2[[i]]*xstar]))/.x-> gx},{i,1,991,1}];

EL8b = Table[{((α8b[[i]]*I*omegatable2[[i]]+f*(β8b[[i]]))/(2*I*omegatable2[[i]]*Exp[I*omegatable2[[i]]*xstar]))/.x-> hx},{i,1,991,1}];

FL8b = Table[{((α8b[[i]]*I*omegatable2[[i]]-f*(β8b[[i]]))/(2*I*omegatable2[[i]]*Exp[-I*omegatable2[[i]]*xstar]))/.x-> hx},{i,1,991,1}];

EL9b = Table[{((α9b[[i]]*I*omegatable2[[i]]+f*(β9b[[i]]))/(2*I*omegatable2[[i]]*Exp[I*omegatable2[[i]]*xstar]))/.x-> ix},{i,1,991,1}];

FL9b = Table[{((α9b[[i]]*I*omegatable2[[i]]-f*(β9b[[i]]))/(2*I*omegatable2[[i]]*Exp[-I*omegatable2[[i]]*xstar]))/.x-> ix},{i,1,991,1}];

EL10b = Table[{((α10b[[i]]*I*omegatable2[[i]]+f*(β10b[[i]]))/(2*I*omegatable2[[i]]*Exp[I*omegatable2[[i]]*xstar]))/.x-> x0},{i,1,991,1}];

FL10b = Table[{((α10b[[i]]*I*omegatable2[[i]]-f*(β10b[[i]]))/(2*I*omegatable2[[i]]*Exp[-I*omegatable2[[i]]*xstar]))/.x-> x0},{i,1,991,1}];

Tceb = Table[{{ax,Part[EL1b[[i]],1]},{bx,Part[EL2b[[i]],1]},{cx,Part[EL3b[[i]],1]},{dx,Part[EL4b[[i]],1]},{ex,Part[EL5b[[i]],1]}, {fx, Part[EL6b[[i]],1]}, {gx, Part[EL7b[[i]],1]},{hx,Part[EL8b[[i]],1]},{ix,Part[EL9b[[i]],1]},{x0,Part[EL10b[[i]],1]}},{i,1,991,1}];

TceCorrectedb= Tceb/.{x_,{{{y_}}}}->{x,y};

Treb=Table[{Part[TceCorrectedb[[i]],j,1],Re[Part[TceCorrectedb[[i]],j,2]]},{i,1,991,1},{j,10}];

ELFitrb = Table[Fit[Treb[[i]],{1,y,y^2,y^3},y],{i,1,991,1}];

Tieb=Table[{Part[TceCorrectedb[[i]],j,1],Im[Part[TceCorrectedb[[i]],j,2]]},{i,1,991,1},{j,10}];

ELFitib = Table[Fit[Tieb[[i]],{1,y,y^2,y^3},y],{i,1,991,1}];

Tcfb = Table[{{ax,Part[FL1b[[i]],1]},{bx,Part[FL2b[[i]],1]},{cx,Part[FL3b[[i]],1]},{dx,Part[FL4b[[i]],1]},{ex,Part[FL5b[[i]],1]}, {fx, Part[FL6b[[i]],1]}, {gx, Part[FL7b[[i]],1]},{hx,Part[FL8b[[i]],1]},{ix,Part[FL9b[[i]],1]},{x0,Part[FL10b[[i]],1]}},{i,1,991,1}];

TcfCorrectedb= Tcfb/.{x_,{{{y_}}}}->{x,y};

Trfb=Table[{Part[TcfCorrectedb[[i]],j,1],Re[Part[TcfCorrectedb[[i]],j,2]]},{i,1,991,1},{j,10}];

FLFitrb = Table[Fit[Trfb[[i]],{1,y,y^2,y^3},y],{i,1,991,1}];

Tifb=Table[{Part[TcfCorrectedb[[i]],j,1],Im[Part[TcfCorrectedb[[i]],j,2]]},{i,1,991,1},{j,10}];

FLFitib = Table[Fit[Tifb[[i]],{1,y,y^2,y^3},y],{i,1,991,1}];

ReELfinb = Table[{omegatable2[[i]],N[ELFitrb[[i]]/.{y-> 1/10}]},{i,1,991,1}];

ImELfinb = Table[{omegatable2[[i]],N[ELFitib[[i]]/.{y-> 1/10}]},{i,1,991,1}];

ReFLfinb = Table[{omegatable2[[i]],N[FLFitrb[[i]]/.{y-> 1/10}]},{i,1,991,1}];

ImFLfinb = Table[{omegatable2[[i]],N[FLFitib[[i]]/.{y-> 1/10}]},{i,1,991,1}];

EL2rp = ListPlot[ReELfinb]

EL2ip = ListPlot[ImELfinb, PlotRange -> All]

FL2rp = ListPlot[ReFLfinb]

FL2ip = ListPlot[ImFLfinb, PlotRange-> All]

Band 3

Length[chiL3]

α1c = Table[{chi[ax]/.chiL3[[i]]},{i,1,1001,1}];

β1c = Table[{chi'[ax]/.chiL3[[i]]},{i,1,1001,1}];

α2c = Table[{chi[bx]/.chiL3[[i]]},{i,1,1001,1}];

β2c = Table[{chi'[bx]/.chiL3[[i]]},{i,1,1001,1}];

α3c = Table[{chi[cx]/.chiL3[[i]]},{i,1,1001,1}];

β3c = Table[{chi'[cx]/.chiL3[[i]]},{i,1,1001,1}];

α4c = Table[{chi[dx]/.chiL3[[i]]},{i,1,1001,1}];

β4c = Table[{chi'[dx]/.chiL3[[i]]},{i,1,1001,1}];

α5c = Table[{chi[ex]/.chiL3[[i]]},{i,1,1001,1}];

β5c = Table[{chi'[ex]/.chiL3[[i]]},{i,1,1001,1}];

α6c = Table[{chi[fx]/.chiL3[[i]]},{i,1,1001,1}];

β6c = Table[{chi'[fx]/.chiL3[[i]]},{i,1,1001,1}];

α7c = Table[{chi[gx]/.chiL3[[i]]},{i,1,1001,1}];

β7c = Table[{chi'[gx]/.chiL3[[i]]},{i,1,1001,1}];

α8c = Table[{chi[hx]/.chiL3[[i]]},{i,1,1001,1}];

β8c = Table[{chi'[hx]/.chiL3[[i]]},{i,1,1001,1}];

α9c = Table[{chi[ix]/.chiL3[[i]]},{i,1,1001,1}];

β9c = Table[{chi'[ix]/.chiL3[[i]]},{i,1,1001,1}];

α10c = Table[{chi[x0]/.chiL3[[i]]},{i,1,1001,1}];

β10c = Table[{chi'[x0]/.chiL3[[i]]},{i,1,1001,1}];

omegatable3 = Table[i,{i,1,6,5/1000}];

EL1c = Table[{((α1c[[i]]*I*omegatable3[[i]]+f*(β1c[[i]]))/(2*I*omegatable3[[i]]*Exp[I*omegatable3[[i]]*xstar]))/.x-> ax},{i,1,1001,1}];

FL1c = Table[{((α1c[[i]]*I*omegatable3[[i]]-f*(β1c[[i]]))/(2*I*omegatable3[[i]]*Exp[-I*omegatable3[[i]]*xstar]))/.x-> ax},{i,1,1001,1}];

EL2c = Table[{((α2c[[i]]*I*omegatable3[[i]]+f*(β2c[[i]]))/(2*I*omegatable3[[i]]*Exp[I*omegatable3[[i]]*xstar]))/.x-> bx},{i,1,1001,1}];

FL2c = Table[{((α2c[[i]]*I*omegatable3[[i]]-f*(β2c[[i]]))/(2*I*omegatable3[[i]]*Exp[-I*omegatable3[[i]]*xstar]))/.x-> bx},{i,1,1001,1}];

EL3c = Table[{((α3c[[i]]*I*omegatable3[[i]]+f*(β3c[[i]]))/(2*I*omegatable3[[i]]*Exp[I*omegatable3[[i]]*xstar]))/.x-> cx},{i,1,1001,1}];

FL3c = Table[{((α3c[[i]]*I*omegatable3[[i]]-f*(β3c[[i]]))/(2*I*omegatable3[[i]]*Exp[-I*omegatable3[[i]]*xstar]))/.x-> cx},{i,1,1001,1}];

EL4c = Table[{((α4c[[i]]*I*omegatable3[[i]]+f*(β4c[[i]]))/(2*I*omegatable3[[i]]*Exp[I*omegatable3[[i]]*xstar]))/.x-> dx},{i,1,1001,1}];

FL4c = Table[{((α4c[[i]]*I*omegatable3[[i]]-f*(β4c[[i]]))/(2*I*omegatable3[[i]]*Exp[-I*omegatable3[[i]]*xstar]))/.x-> dx},{i,1,1001,1}];

EL5c = Table[{((α5c[[i]]*I*omegatable3[[i]]+f*(β5c[[i]]))/(2*I*omegatable3[[i]]*Exp[I*omegatable3[[i]]*xstar]))/.x-> ex},{i,1,1001,1}];

FL5c = Table[{((α5c[[i]]*I*omegatable3[[i]]-f*(β5c[[i]]))/(2*I*omegatable3[[i]]*Exp[-I*omegatable3[[i]]*xstar]))/.x-> ex},{i,1,1001,1}];

EL6c = Table[{((α6c[[i]]*I*omegatable3[[i]]+f*(β6c[[i]]))/(2*I*omegatable3[[i]]*Exp[I*omegatable3[[i]]*xstar]))/.x-> fx},{i,1,1001,1}];

FL6c = Table[{((α6c[[i]]*I*omegatable3[[i]]-f*(β6c[[i]]))/(2*I*omegatable3[[i]]*Exp[-I*omegatable3[[i]]*xstar]))/.x-> fx},{i,1,1001,1}];

EL7c = Table[{((α7c[[i]]*I*omegatable3[[i]]+f*(β7c[[i]]))/(2*I*omegatable3[[i]]*Exp[I*omegatable3[[i]]*xstar]))/.x-> gx},{i,1,1001,1}];

FL7c = Table[{((α7c[[i]]*I*omegatable3[[i]]-f*(β7c[[i]]))/(2*I*omegatable3[[i]]*Exp[-I*omegatable3[[i]]*xstar]))/.x-> gx},{i,1,1001,1}];

EL8c = Table[{((α8c[[i]]*I*omegatable3[[i]]+f*(β8c[[i]]))/(2*I*omegatable3[[i]]*Exp[I*omegatable3[[i]]*xstar]))/.x-> hx},{i,1,1001,1}];

FL8c = Table[{((α8c[[i]]*I*omegatable3[[i]]-f*(β8c[[i]]))/(2*I*omegatable3[[i]]*Exp[-I*omegatable3[[i]]*xstar]))/.x-> hx},{i,1,1001,1}];

EL9c = Table[{((α9c[[i]]*I*omegatable3[[i]]+f*(β9c[[i]]))/(2*I*omegatable3[[i]]*Exp[I*omegatable3[[i]]*xstar]))/.x-> ix},{i,1,1001,1}];

FL9c = Table[{((α9c[[i]]*I*omegatable3[[i]]-f*(β9c[[i]]))/(2*I*omegatable3[[i]]*Exp[-I*omegatable3[[i]]*xstar]))/.x-> ix},{i,1,1001,1}];

EL10c = Table[{((α10c[[i]]*I*omegatable3[[i]]+f*(β10c[[i]]))/(2*I*omegatable3[[i]]*Exp[I*omegatable3[[i]]*xstar]))/.x-> x0},{i,1,1001,1}];

FL10c = Table[{((α10c[[i]]*I*omegatable3[[i]]-f*(β10c[[i]]))/(2*I*omegatable3[[i]]*Exp[-I*omegatable3[[i]]*xstar]))/.x-> x0},{i,1,1001,1}];

Tcec = Table[{{ax,Part[EL1c[[i]],1]},{bx,Part[EL2c[[i]],1]},{cx,Part[EL3c[[i]],1]},{dx,Part[EL4c[[i]],1]},{ex,Part[EL5c[[i]],1]}, {fx, Part[EL6c[[i]],1]}, {gx, Part[EL7c[[i]],1]},{hx,Part[EL8c[[i]],1]},{ix,Part[EL9c[[i]],1]},{x0,Part[EL10c[[i]],1]}},{i,1,1001,1}];

TceCorrectedc = Tcec/.{x_,{{{y_}}}}->{x,y};

Trec =Table[{Part[TceCorrectedc[[i]],j,1],Re[Part[TceCorrectedc[[i]],j,2]]},{i,1,1001,1},{j,10}];

ELFitrc = Table[Fit[Trec[[i]],{1,y,y^2,y^3},y],{i,1,1001,1}];

Tiec =Table[{Part[TceCorrectedc[[i]],j,1],Im[Part[TceCorrectedc[[i]],j,2]]},{i,1,1001,1},{j,10}];

ELFitic = Table[Fit[Tiec[[i]],{1,y,y^2,y^3},y],{i,1,1001,1}];

Tcfc = Table[{{ax,Part[FL1c[[i]],1]},{bx,Part[FL2c[[i]],1]},{cx,Part[FL3c[[i]],1]},{dx,Part[FL4c[[i]],1]},{ex,Part[FL5c[[i]],1]}, {fx, Part[FL6c[[i]],1]}, {gx, Part[FL7c[[i]],1]},{hx,Part[FL8c[[i]],1]},{ix,Part[FL9c[[i]],1]},{x0,Part[FL10c[[i]],1]}},{i,1,1001,1}];

TcfCorrectedc = Tcfc/.{x_,{{{y_}}}}->{x,y};

Trfc =Table[{Part[TcfCorrectedc[[i]],j,1],Re[Part[TcfCorrectedc[[i]],j,2]]},{i,1,1001,1},{j,10}];

FLFitrc = Table[Fit[Trfc[[i]],{1,y,y^2,y^3},y],{i,1,1001,1}];

Tifc =Table[{Part[TcfCorrectedc[[i]],j,1],Im[Part[TcfCorrectedc[[i]],j,2]]},{i,1,1001,1},{j,10}];

FLFitic = Table[Fit[Tifc[[i]],{1,y,y^2,y^3},y],{i,1,1001,1}];

ReELfinc = Table[{omegatable3[[i]],N[ELFitrc[[i]]/.{y-> 1/10}]},{i,1,1001,1}];

ImELfinc = Table[{omegatable3[[i]],N[ELFitic[[i]]/.{y-> 1/10}]},{i,1,1001,1}];

ReFLfinc = Table[{omegatable3[[i]],N[FLFitrc[[i]]/.{y-> 1/10}]},{i,1,1001,1}];

ImFLfinc = Table[{omegatable3[[i]],N[FLFitic[[i]]/.{y-> 1/10}]},{i,1,1001,1}];

EL3rp = ListPlot[ReELfinc]

EL3ip = ListPlot[ImELfinc, PlotRange -> All]

FL3rp = ListPlot[ReFLfinc]

FL3ip = ListPlot[ImFLfinc, PlotRange-> All]

Band 4

Length[chiL4]

α1d = Table[{chi[ax]/.chiL4[[i]]},{i,1,1101,1}];

β1d = Table[{chi'[ax]/.chiL4[[i]]},{i,1,1101,1}];

α2d = Table[{chi[bx]/.chiL4[[i]]},{i,1,1101,1}];

β2d = Table[{chi'[bx]/.chiL4[[i]]},{i,1,1101,1}];

α3d = Table[{chi[cx]/.chiL4[[i]]},{i,1,1101,1}];

β3d = Table[{chi'[cx]/.chiL4[[i]]},{i,1,1101,1}];

α4d = Table[{chi[dx]/.chiL4[[i]]},{i,1,1101,1}];

β4d = Table[{chi'[dx]/.chiL4[[i]]},{i,1,1101,1}];

α5d = Table[{chi[ex]/.chiL4[[i]]},{i,1,1101,1}];

β5d = Table[{chi'[ex]/.chiL4[[i]]},{i,1,1101,1}];

α6d = Table[{chi[fx]/.chiL4[[i]]},{i,1,1101,1}];

β6d = Table[{chi'[fx]/.chiL4[[i]]},{i,1,1101,1}];

α7d = Table[{chi[gx]/.chiL4[[i]]},{i,1,1101,1}];

β7d = Table[{chi'[gx]/.chiL4[[i]]},{i,1,1101,1}];

α8d = Table[{chi[hx]/.chiL4[[i]]},{i,1,1101,1}];

β8d = Table[{chi'[hx]/.chiL4[[i]]},{i,1,1101,1}];

α9d = Table[{chi[ix]/.chiL4[[i]]},{i,1,1101,1}];

β9d = Table[{chi'[ix]/.chiL4[[i]]},{i,1,1101,1}];

α10d = Table[{chi[x0]/.chiL4[[i]]},{i,1,1101,1}];

β10d = Table[{chi'[x0]/.chiL4[[i]]},{i,1,1101,1}];

omegatable4 = Table[i,{i,6,5000/100,4/100}];

EL1d = Table[{((α1d[[i]]*I*omegatable4[[i]]+f*(β1d[[i]]))/(2*I*omegatable4[[i]]*Exp[I*omegatable4[[i]]*xstar]))/.x-> ax},{i,1,1101,1}];

FL1d = Table[{((α1d[[i]]*I*omegatable4[[i]]-f*(β1d[[i]]))/(2*I*omegatable4[[i]]*Exp[-I*omegatable4[[i]]*xstar]))/.x-> ax},{i,1,1101,1}];

EL2d = Table[{((α2d[[i]]*I*omegatable4[[i]]+f*(β2d[[i]]))/(2*I*omegatable4[[i]]*Exp[I*omegatable4[[i]]*xstar]))/.x-> bx},{i,1,1101,1}];

FL2d = Table[{((α2d[[i]]*I*omegatable4[[i]]-f*(β2d[[i]]))/(2*I*omegatable4[[i]]*Exp[-I*omegatable4[[i]]*xstar]))/.x-> bx},{i,1,1101,1}];

EL3d = Table[{((α3d[[i]]*I*omegatable4[[i]]+f*(β3d[[i]]))/(2*I*omegatable4[[i]]*Exp[I*omegatable4[[i]]*xstar]))/.x-> cx},{i,1,1101,1}];

FL3d = Table[{((α3d[[i]]*I*omegatable4[[i]]-f*(β3d[[i]]))/(2*I*omegatable4[[i]]*Exp[-I*omegatable4[[i]]*xstar]))/.x-> cx},{i,1,1101,1}];

EL4d = Table[{((α4d[[i]]*I*omegatable4[[i]]+f*(β4d[[i]]))/(2*I*omegatable4[[i]]*Exp[I*omegatable4[[i]]*xstar]))/.x-> dx},{i,1,1101,1}];

FL4d = Table[{((α4d[[i]]*I*omegatable4[[i]]-f*(β4d[[i]]))/(2*I*omegatable4[[i]]*Exp[-I*omegatable4[[i]]*xstar]))/.x-> dx},{i,1,1101,1}];

EL5d = Table[{((α5d[[i]]*I*omegatable4[[i]]+f*(β5d[[i]]))/(2*I*omegatable4[[i]]*Exp[I*omegatable4[[i]]*xstar]))/.x-> ex},{i,1,1101,1}];

FL5d = Table[{((α5d[[i]]*I*omegatable4[[i]]-f*(β5d[[i]]))/(2*I*omegatable4[[i]]*Exp[-I*omegatable4[[i]]*xstar]))/.x-> ex},{i,1,1101,1}];

EL6d = Table[{((α6d[[i]]*I*omegatable4[[i]]+f*(β6d[[i]]))/(2*I*omegatable4[[i]]*Exp[I*omegatable4[[i]]*xstar]))/.x-> fx},{i,1,1101,1}];

FL6d = Table[{((α6d[[i]]*I*omegatable4[[i]]-f*(β6d[[i]]))/(2*I*omegatable4[[i]]*Exp[-I*omegatable4[[i]]*xstar]))/.x-> fx},{i,1,1101,1}];

EL7d = Table[{((α7d[[i]]*I*omegatable4[[i]]+f*(β7d[[i]]))/(2*I*omegatable4[[i]]*Exp[I*omegatable4[[i]]*xstar]))/.x-> gx},{i,1,1101,1}];

FL7d = Table[{((α7d[[i]]*I*omegatable4[[i]]-f*(β7d[[i]]))/(2*I*omegatable4[[i]]*Exp[-I*omegatable4[[i]]*xstar]))/.x-> gx},{i,1,1101,1}];

EL8d = Table[{((α8d[[i]]*I*omegatable4[[i]]+f*(β8d[[i]]))/(2*I*omegatable4[[i]]*Exp[I*omegatable4[[i]]*xstar]))/.x-> hx},{i,1,1101,1}];

FL8d = Table[{((α8d[[i]]*I*omegatable4[[i]]-f*(β8d[[i]]))/(2*I*omegatable4[[i]]*Exp[-I*omegatable4[[i]]*xstar]))/.x-> hx},{i,1,1101,1}];

EL9d = Table[{((α9d[[i]]*I*omegatable4[[i]]+f*(β9d[[i]]))/(2*I*omegatable4[[i]]*Exp[I*omegatable4[[i]]*xstar]))/.x-> ix},{i,1,1101,1}];

FL9d = Table[{((α9d[[i]]*I*omegatable4[[i]]-f*(β9d[[i]]))/(2*I*omegatable4[[i]]*Exp[-I*omegatable4[[i]]*xstar]))/.x-> ix},{i,1,1101,1}];

EL10d = Table[{((α10d[[i]]*I*omegatable4[[i]]+f*(β10d[[i]]))/(2*I*omegatable4[[i]]*Exp[I*omegatable4[[i]]*xstar]))/.x-> x0},{i,1,1101,1}];

FL10d = Table[{((α10d[[i]]*I*omegatable4[[i]]-f*(β10d[[i]]))/(2*I*omegatable4[[i]]*Exp[-I*omegatable4[[i]]*xstar]))/.x-> x0},{i,1,1101,1}];

Tced = Table[{{ax,Part[EL1d[[i]],1]},{bx,Part[EL2d[[i]],1]},{cx,Part[EL3d[[i]],1]},{dx,Part[EL4d[[i]],1]},{ex,Part[EL5d[[i]],1]}, {fx, Part[EL6d[[i]],1]}, {gx, Part[EL7d[[i]],1]},{hx,Part[EL8d[[i]],1]},{ix,Part[EL9d[[i]],1]},{x0,Part[EL10d[[i]],1]}},{i,1,1101,1}];

TceCorrectedd = Tced/.{x_,{{{y_}}}}->{x,y};

Tred =Table[{Part[TceCorrectedd[[i]],j,1],Re[Part[TceCorrectedd[[i]],j,2]]},{i,1,1101,1},{j,10}];

ELFitrd = Table[Fit[Tred[[i]],{1,y,y^2,y^3},y],{i,1,1101,1}];

Tied =Table[{Part[TceCorrectedd[[i]],j,1],Im[Part[TceCorrectedd[[i]],j,2]]},{i,1,1101,1},{j,10}];

ELFitid = Table[Fit[Tied[[i]],{1,y,y^2,y^3},y],{i,1,1101,1}];

Tcfd = Table[{{ax,Part[FL1d[[i]],1]},{bx,Part[FL2d[[i]],1]},{cx,Part[FL3d[[i]],1]},{dx,Part[FL4d[[i]],1]},{ex,Part[FL5d[[i]],1]}, {fx, Part[FL6d[[i]],1]}, {gx, Part[FL7d[[i]],1]},{hx,Part[FL8d[[i]],1]},{ix,Part[FL9d[[i]],1]},{x0,Part[FL10d[[i]],1]}},{i,1,1101,1}];

TcfCorrectedd = Tcfd/.{x_,{{{y_}}}}->{x,y};

Trfd =Table[{Part[TcfCorrectedd[[i]],j,1],Re[Part[TcfCorrectedd[[i]],j,2]]},{i,1,1101,1},{j,10}];

FLFitrd = Table[Fit[Trfd[[i]],{1,y,y^2,y^3},y],{i,1,1101,1}];

Tifd =Table[{Part[TcfCorrectedd[[i]],j,1],Im[Part[TcfCorrectedd[[i]],j,2]]},{i,1,1101,1},{j,10}];

FLFitid = Table[Fit[Tifd[[i]],{1,y,y^2,y^3},y],{i,1,1101,1}];

ReELfind = Table[{omegatable4[[i]],N[ELFitrd[[i]]/.{y-> 1/10}]},{i,1,1101,1}];

ImELfind = Table[{omegatable4[[i]],N[ELFitid[[i]]/.{y-> 1/10}]},{i,1,1101,1}];

ReFLfind = Table[{omegatable4[[i]],N[FLFitrd[[i]]/.{y-> 1/10}]},{i,1,1101,1}];

ImFLfind = Table[{omegatable4[[i]],N[FLFitid[[i]]/.{y-> 1/10}]},{i,1,1101,1}];

EL4rp = ListPlot[ReELfind]

EL4ip = ListPlot[ImELfind, PlotRange -> All]

FL4rp = ListPlot[ReFLfind]

FL4ip = ListPlot[ImFLfind, PlotRange-> All]

Band 5

Length[chiL5]

α1e = Table[{chi[ax]/.chiL5[[i]]},{i,1,1001,1}];

β1e = Table[{chi'[ax]/.chiL5[[i]]},{i,1,1001,1}];

α2e = Table[{chi[bx]/.chiL5[[i]]},{i,1,1001,1}];

β2e = Table[{chi'[bx]/.chiL5[[i]]},{i,1,1001,1}];

α3e = Table[{chi[cx]/.chiL5[[i]]},{i,1,1001,1}];

β3e = Table[{chi'[cx]/.chiL5[[i]]},{i,1,1001,1}];

α4e = Table[{chi[dx]/.chiL5[[i]]},{i,1,1001,1}];

β4e = Table[{chi'[dx]/.chiL5[[i]]},{i,1,1001,1}];

α5e = Table[{chi[ex]/.chiL5[[i]]},{i,1,1001,1}];

β5e = Table[{chi'[ex]/.chiL5[[i]]},{i,1,1001,1}];

α6e = Table[{chi[fx]/.chiL5[[i]]},{i,1,1001,1}];

β6e = Table[{chi'[fx]/.chiL5[[i]]},{i,1,1001,1}];

α7e = Table[{chi[gx]/.chiL5[[i]]},{i,1,1001,1}];

β7e = Table[{chi'[gx]/.chiL5[[i]]},{i,1,1001,1}];

α8e = Table[{chi[hx]/.chiL5[[i]]},{i,1,1001,1}];

β8e = Table[{chi'[hx]/.chiL5[[i]]},{i,1,1001,1}];

α9e = Table[{chi[ix]/.chiL5[[i]]},{i,1,1001,1}];

β9e = Table[{chi'[ix]/.chiL5[[i]]},{i,1,1001,1}];

α10e = Table[{chi[x0]/.chiL5[[i]]},{i,1,1001,1}];

β10e = Table[{chi'[x0]/.chiL5[[i]]},{i,1,1001,1}];

omegatable5 = Table[i,{i,5000/100,10000/100,5/100}];

EL1e = Table[{((α1e[[i]]*I*omegatable5[[i]]+f*(β1e[[i]]))/(2*I*omegatable5[[i]]*Exp[I*omegatable5[[i]]*xstar]))/.x-> ax},{i,1,1001,1}];

FL1e = Table[{((α1e[[i]]*I*omegatable5[[i]]-f*(β1e[[i]]))/(2*I*omegatable5[[i]]*Exp[-I*omegatable5[[i]]*xstar]))/.x-> ax},{i,1,1001,1}];

EL2e = Table[{((α2e[[i]]*I*omegatable5[[i]]+f*(β2e[[i]]))/(2*I*omegatable5[[i]]*Exp[I*omegatable5[[i]]*xstar]))/.x-> bx},{i,1,1001,1}];

FL2e = Table[{((α2e[[i]]*I*omegatable5[[i]]-f*(β2e[[i]]))/(2*I*omegatable5[[i]]*Exp[-I*omegatable5[[i]]*xstar]))/.x-> bx},{i,1,1001,1}];

EL3e = Table[{((α3e[[i]]*I*omegatable5[[i]]+f*(β3e[[i]]))/(2*I*omegatable5[[i]]*Exp[I*omegatable5[[i]]*xstar]))/.x-> cx},{i,1,1001,1}];

FL3e = Table[{((α3e[[i]]*I*omegatable5[[i]]-f*(β3e[[i]]))/(2*I*omegatable5[[i]]*Exp[-I*omegatable5[[i]]*xstar]))/.x-> cx},{i,1,1001,1}];

EL4e = Table[{((α4e[[i]]*I*omegatable5[[i]]+f*(β4e[[i]]))/(2*I*omegatable5[[i]]*Exp[I*omegatable5[[i]]*xstar]))/.x-> dx},{i,1,1001,1}];

FL4e = Table[{((α4e[[i]]*I*omegatable5[[i]]-f*(β4e[[i]]))/(2*I*omegatable5[[i]]*Exp[-I*omegatable5[[i]]*xstar]))/.x-> dx},{i,1,1001,1}];

EL5e = Table[{((α5e[[i]]*I*omegatable5[[i]]+f*(β5e[[i]]))/(2*I*omegatable5[[i]]*Exp[I*omegatable5[[i]]*xstar]))/.x-> ex},{i,1,1001,1}];

FL5e = Table[{((α5e[[i]]*I*omegatable5[[i]]-f*(β5e[[i]]))/(2*I*omegatable5[[i]]*Exp[-I*omegatable5[[i]]*xstar]))/.x-> ex},{i,1,1001,1}];

EL6e = Table[{((α6e[[i]]*I*omegatable5[[i]]+f*(β6e[[i]]))/(2*I*omegatable5[[i]]*Exp[I*omegatable5[[i]]*xstar]))/.x-> fx},{i,1,1001,1}];

FL6e = Table[{((α6e[[i]]*I*omegatable5[[i]]-f*(β6e[[i]]))/(2*I*omegatable5[[i]]*Exp[-I*omegatable5[[i]]*xstar]))/.x-> fx},{i,1,1001,1}];

EL7e = Table[{((α7e[[i]]*I*omegatable5[[i]]+f*(β7e[[i]]))/(2*I*omegatable5[[i]]*Exp[I*omegatable5[[i]]*xstar]))/.x-> gx},{i,1,1001,1}];

FL7e = Table[{((α7e[[i]]*I*omegatable5[[i]]-f*(β7e[[i]]))/(2*I*omegatable5[[i]]*Exp[-I*omegatable5[[i]]*xstar]))/.x-> gx},{i,1,1001,1}];

EL8e = Table[{((α8e[[i]]*I*omegatable5[[i]]+f*(β8e[[i]]))/(2*I*omegatable5[[i]]*Exp[I*omegatable5[[i]]*xstar]))/.x-> hx},{i,1,1001,1}];

FL8e = Table[{((α8e[[i]]*I*omegatable5[[i]]-f*(β8e[[i]]))/(2*I*omegatable5[[i]]*Exp[-I*omegatable5[[i]]*xstar]))/.x-> hx},{i,1,1001,1}];

EL9e = Table[{((α9e[[i]]*I*omegatable5[[i]]+f*(β9e[[i]]))/(2*I*omegatable5[[i]]*Exp[I*omegatable5[[i]]*xstar]))/.x-> ix},{i,1,1001,1}];

FL9e = Table[{((α9e[[i]]*I*omegatable5[[i]]-f*(β9e[[i]]))/(2*I*omegatable5[[i]]*Exp[-I*omegatable5[[i]]*xstar]))/.x-> ix},{i,1,1001,1}];

EL10e = Table[{((α10e[[i]]*I*omegatable5[[i]]+f*(β10e[[i]]))/(2*I*omegatable5[[i]]*Exp[I*omegatable5[[i]]*xstar]))/.x-> x0},{i,1,1001,1}];

FL10e = Table[{((α10e[[i]]*I*omegatable5[[i]]-f*(β10e[[i]]))/(2*I*omegatable5[[i]]*Exp[-I*omegatable5[[i]]*xstar]))/.x-> x0},{i,1,1001,1}];

Tcee = Table[{{ax,Part[EL1e[[i]],1]},{bx,Part[EL2e[[i]],1]},{cx,Part[EL3e[[i]],1]},{dx,Part[EL4e[[i]],1]},{ex,Part[EL5e[[i]],1]}, {fx, Part[EL6e[[i]],1]}, {gx, Part[EL7e[[i]],1]},{hx,Part[EL8e[[i]],1]},{ix,Part[EL9e[[i]],1]},{x0,Part[EL10e[[i]],1]}},{i,1,1001,1}];

TceCorrectede = Tcee/.{x_,{{{y_}}}}->{x,y};

Tree =Table[{Part[TceCorrectede[[i]],j,1],Re[Part[TceCorrectede[[i]],j,2]]},{i,1,1001,1},{j,10}];

ELFitre = Table[Fit[Tree[[i]],{1,y,y^2,y^3},y],{i,1,1001,1}];

Tiee =Table[{Part[TceCorrectede[[i]],j,1],Im[Part[TceCorrectede[[i]],j,2]]},{i,1,1001,1},{j,10}];

ELFitie = Table[Fit[Tiee[[i]],{1,y,y^2,y^3},y],{i,1,1001,1}];

Tcfe = Table[{{ax,Part[FL1e[[i]],1]},{bx,Part[FL2e[[i]],1]},{cx,Part[FL3e[[i]],1]},{dx,Part[FL4e[[i]],1]},{ex,Part[FL5e[[i]],1]}, {fx, Part[FL6e[[i]],1]}, {gx, Part[FL7e[[i]],1]},{hx,Part[FL8e[[i]],1]},{ix,Part[FL9e[[i]],1]},{x0,Part[FL10e[[i]],1]}},{i,1,1001,1}];

TcfCorrectede = Tcfe/.{x_,{{{y_}}}}->{x,y};

Trfe =Table[{Part[TcfCorrectede[[i]],j,1],Re[Part[TcfCorrectede[[i]],j,2]]},{i,1,1001,1},{j,10}];

FLFitre = Table[Fit[Trfe[[i]],{1,y,y^2,y^3},y],{i,1,1001,1}];

Tife =Table[{Part[TcfCorrectede[[i]],j,1],Im[Part[TcfCorrectede[[i]],j,2]]},{i,1,1001,1},{j,10}];

FLFitie = Table[Fit[Tife[[i]],{1,y,y^2,y^3},y],{i,1,1001,1}];

ReELfine = Table[{omegatable5[[i]],N[ELFitre[[i]]/.{y-> 1/10}]},{i,1,1001,1}];

ImELfine = Table[{omegatable5[[i]],N[ELFitie[[i]]/.{y-> 1/10}]},{i,1,1001,1}];

ReFLfine = Table[{omegatable5[[i]],N[FLFitre[[i]]/.{y-> 1/10}]},{i,1,1001,1}];

ImFLfine = Table[{omegatable5[[i]],N[FLFitie[[i]]/.{y-> 1/10}]},{i,1,1001,1}];

EL5rp = ListPlot[ReELfine]

EL5ip = ListPlot[ImELfine, PlotRange -> All]

FL5rp = ListPlot[ReFLfine]

FL5ip = ListPlot[ImFLfine, PlotRange-> All]

Logarithmic Plots of All Bands

EL Real:

p1b = ListLogLinearPlot[{ReELfin,ReELfinb,ReELfinc, ReELfind,ReELfine}, PlotRange-> {0,1}, PlotStyle-> Blue, Frame ->  True,FrameLabel-> {"ω'"}, PlotLabel-> Style[Framed["Scattering Parameters Re("<>ToString[Subscript["E",L],StandardForm]<>"), Re("<>ToString[Subscript["F",L],StandardForm]<>") versus frequency"],Black], LabelStyle->Directive[Black],
PlotLabels-> Placed[Style[Subscript["E",L],Black, FontSize-> 14],Above]]

EL Imaginary:

pi = ListLogLinearPlot[{ImELfin,ImELfinb,ImELfinc, ImELfind,ImELfine}, PlotRange-> All, PlotStyle-> Blue, Frame-> True,FrameLabel-> {"ω'"},PlotLabel-> Style[Framed["Scattering Parameters Im("<>ToString[Subscript["E",L],StandardForm]<>"), Im("<>ToString[Subscript["F",L],StandardForm]<>") versus frequency"],Black], PlotLabels-> Placed[Style[Subscript["E",L],Black, FontSize-> 14],Right]]

FL Real:

p2b = ListLogLinearPlot[{ReFLfin,ReFLfinb,ReFLfinc, ReFLfind,ReFLfine}, PlotRange-> {0,1}, PlotStyle-> Red, Frame ->  True,FrameLabel-> {"ω'"},LabelStyle->Directive[Black, Medium],PlotLabels-> Placed[Style[Subscript["F",L],Black, FontSize-> 14],Above]]

Show[{p1b,p2b}]

Show[%839,ImageSize->Large]

FL Imaginary:

p3b = ListLogLinearPlot[{ImFLfin,ImFLfinb,ImFLfinc, ImFLfind,ImFLfine},PlotRange-> All,PlotStyle-> Red, Frame-> True,FrameLabel-> {"ω'"}, PlotLabels-> Placed[Style[Subscript["F",L],Black, FontSize-> 14],Right]]

Show[pi,p3b,PlotRange-> All]

Show[%829,ImageSize->Large]

(* Save["ELFLvalues.m",{ReELfin,ReELfinb,ReELfinc, ReELfind, ReELfine, ImELfin, ImELfinb,ImELfinc,ImELfind,ImELfine,ReFLfin, ReFLfinb,ReFLfinc,ReFLfind,ReFLfine,ImFLfin,ImFLfinb,ImFLfinc,ImFLfind,ImFLfine}] *)