(* Plots for E(t)/E0 for t1 = - t2 *)

EtE0 = (1/2)*(Tanh[b*(t+t2)]-Tanh[b*(t-t2)])

(* Here is a plot for t2 = 20, 40, 60 and b = 0.5. Blue is 20, orange is 40, green is 60. *)

Plot[{EtE0/.{b->1,t2->20},EtE0/.{b->1,t2->40},EtE0/.{b->1,t2->60}},{t,-70,70}]

(* Now fix t2 = 60 and try various values of b.  b = 0.1 is blue, 0.5 is orange, 1 is green *)

Plot[{EtE0/.{b->1/10,t2->60},EtE0/.{b->1/4,t2->60},EtE0/.{b->1,t2->60}},{t,-100,100}]

(* OK now I need to generate tables for the plots. *)

T1 = Table[{N[i/100],N[EtE0/.{b->1,t2->20,t->i/100},17]},{i,-5000,5000}];

T2 = Table[{N[i/100],N[EtE0/.{b->1,t2->40,t->i/100},17]},{i,-5000,5000}];

T3 = Table[{N[i/100],N[EtE0/.{b->1,t2->60,t->i/100},17]},{i,-10000,10000}];

ListPlot[{T1,T2,T3}]

T4 = Table[{N[i/100],N[EtE0/.{b->1/10,t2->60,t->i/100},17]},{i,-10000,10000}];

T5 = Table[{N[i/100],N[EtE0/.{b->1/4,t2->60,t->i/100},17]},{i,-10000,10000}];

T6 = Table[{N[i/100],N[EtE0/.{b->1,t2->60,t->i/100},17]},{i,-10000,10000}];

ListPlot[{T4,T5,T6}]

Export["C:\\Users\\anderson\\Research\\de Sitter States\\Dillon\\E_field_2_parameter_model\\Efield_b_1_t2_20.dat",T1]

Export["C:\\Users\\anderson\\Research\\de Sitter States\\Dillon\\E_field_2_parameter_model\\Efield_b_1_t2_40.dat",T2]

Export["C:\\Users\\anderson\\Research\\de Sitter States\\Dillon\\E_field_2_parameter_model\\Efield_b_1_t2_60.dat",T3]

Export["C:\\Users\\anderson\\Research\\de Sitter States\\Dillon\\E_field_2_parameter_model\\Efield_b_0p1_t2_60.dat",T4]

Export["C:\\Users\\anderson\\Research\\de Sitter States\\Dillon\\E_field_2_parameter_model\\Efield_b_0p25_t2_60.dat",T5]

(* Note that T3 and T6 are the same.  The plot below confirms this. *)

ListPlot[{T3,T6}]