level 1
fhn2021
楼主
tmr是输出函数，afa是变量，输出曲线为什么不光滑？
代码如下：
Clear["`*"]
vh = 2*(1.88 - 0.59*Cos[Pi/(1 + d1/(0.65*10^(-9)))]);
u = vh*0.89;
mu = 4.6;
gama = -4*10^(-29)*(1 - Cos[6*afa])*(1 - 1*10^(-8)*I);
d1 = 1.95*10^(-9);
d2 = 1.95*10^(-9);
HBAR = 1.05457266*10^(-34);
ME = 9.1093897*10^(-31);
ELEC = 1.60217733*10^(-19);
Kh = 2.29*10^(10);
kc = Sqrt[2*ME*ELEC/HBAR^2];
k := kc*Sqrt[mu]
k0 := Sqrt[k^2 - k^2 Sin[x]^2]
kh := Sqrt[k^2 - (Kh - k Sin[x] Cos[y])^2 - k^2 Sin[x]^2 Sin[y]^2]
khg := Sqrt[
k^2 - (2*Kh*Sin[afa/2]*Sin[afa/2] -
k Sin[x] Cos[y])^2 - (2*Kh*Sin[afa/2]*Cos[afa/2] +
k Sin[x] Sin[y])^2]
kg1 := Sqrt[
k^2 - (Kh*Cos[Pi/3 - afa] -
k Sin[x] Cos[y])^2 - (Kh*Sin[Pi/3 - afa] + k Sin[x] Sin[y])^2]
kg2 := Sqrt[
k^2 - (Kh*Cos[afa] - k Sin[x] Cos[y])^2 - (k Sin[x] Sin[y] -
Kh*Sin[afa])^2]
k0pl := Sqrt[k^2 - k^2 Sin[x]^2 + kc^2 vh - kc^2*u]
k0mi := Sqrt[k^2 - k^2 Sin[x]^2 - kc^2 vh - kc^2*u]
khpl := Sqrt[
k^2 - (Kh - k Sin[x] Cos[y])^2 - k^2 Sin[x]^2 Sin[y]^2 + kc^2 vh -
kc^2*u]
khmi := Sqrt[
k^2 - (Kh - k Sin[x] Cos[y])^2 - k^2 Sin[x]^2 Sin[y]^2 - kc^2 vh -
kc^2*u]
k0plpl := Sqrt[k^2 - k^2 Sin[x]^2 + 2*kc^2 vh]
k0mimi := Sqrt[k^2 - k^2 Sin[x]^2 - 2*kc^2 vh]
khplpl :=
Sqrt[k^2 - (Kh - k Sin[x] Cos[y])^2 - k^2 Sin[x]^2 Sin[y]^2 +
2*kc^2 vh]
khmimi :=
Sqrt[k^2 - (Kh - k Sin[x] Cos[y])^2 - k^2 Sin[x]^2 Sin[y]^2 -
2*kc^2 vh]
khgplpl :=
Sqrt[k^2 - (2*Kh*Sin[afa/2]*Sin[afa/2] -
k Sin[x] Cos[y])^2 - (2*Kh*Sin[afa/2]*Cos[afa/2] +
k Sin[x] Sin[y])^2 + 2*kc^2 vh]
khgmimi :=
Sqrt[k^2 - (2*Kh*Sin[afa/2]*Sin[afa/2] -
k Sin[x] Cos[y])^2 - (2*Kh*Sin[afa/2]*Cos[afa/2] +
k Sin[x] Sin[y])^2 - 2*kc^2 vh]
kg1plpl :=
Sqrt[k^2 - (Kh*Cos[Pi/3 - afa] -
k Sin[x] Cos[y])^2 - (Kh*Sin[Pi/3 - afa] + k Sin[x] Sin[y])^2 +
2*kc^2 vh]
kg1mimi :=
Sqrt[k^2 - (Kh*Cos[Pi/3 - afa] -
k Sin[x] Cos[y])^2 - (Kh*Sin[Pi/3 - afa] + k Sin[x] Sin[y])^2 -
2*kc^2 vh]
kg2plpl :=
Sqrt[k^2 - (Kh*Cos[afa] - k Sin[x] Cos[y])^2 - (k Sin[x] Sin[y] -
Kh*Sin[afa])^2 + 2*kc^2 vh]
kg2mimi :=
Sqrt[k^2 - (Kh*Cos[afa] - k Sin[x] Cos[y])^2 - (k Sin[x] Sin[y] -
Kh*Sin[afa])^2 - 2*kc^2 vh]
A1 := 1/(1 + I*ME*gama/(HBAR^2*k0pl))*Exp[I*k0pl*d1]/2
B1 := 1/(1 + I*ME*gama/(HBAR^2*k0pl))*Exp[I*k0pl*d1]/2
A2 := 1/(1 + I*ME*gama/(HBAR^2*khpl))*Exp[I*khpl*d1]/2
B2 := 1/(1 + I*ME*gama/(HBAR^2*khpl))*Exp[I*khpl*d1]/2
A3 := 1/(1 + I*ME*gama/(HBAR^2*k0mi))*Exp[I*k0mi*d1]/2
B3 := 1/(1 + I*ME*gama/(HBAR^2*k0mi))*Exp[I*k0mi*d1]/2
A4 := 1/(1 + I*ME*gama/(HBAR^2*khmi))*Exp[I*khmi*d1]/2
B4 := 1/(1 + I*ME*gama/(HBAR^2*khmi))*Exp[I*khmi*d1]/2
T1 := 1/4 Re[
Abs[A1]^2*k0plpl*Exp[I*(k0plpl - Conjugate[k0plpl])*d2] +
Abs[A1]^2*kg1plpl*Exp[I*(kg1plpl - Conjugate[kg1plpl])*d2] +
Abs[A2]^2*khplpl*Exp[I*(khplpl - Conjugate[khplpl])*d2] +
Abs[A2]^2*khgplpl*Exp[I*(khgplpl - Conjugate[khgplpl])*d2] +
Abs[B3]^2*k0mimi*Exp[I*(k0mimi - Conjugate[k0mimi])*d2] +
Abs[B3]^2*kg1mimi*Exp[I*(kg1mimi - Conjugate[kg1mimi])*d2] +
Abs[B4]^2*khmimi*Exp[I*(khmimi - Conjugate[khmimi])*d2] +
Abs[B4]^2*khgmimi*Exp[I*(khgmimi - Conjugate[khgmimi])*d2] +
Abs[B1 + A3]^2*k0*Exp[I*(k0 - Conjugate[k0])*d2] +
Abs[B1 - A3]^2*kg1*Exp[I*(kg1 - Conjugate[kg1])*d2] +
Abs[A4 - B2]^2*kh*Exp[I*(kh - Conjugate[kh])*d2] +
Abs[A4 + B2]^2*khg*Exp[I*(khg - Conjugate[khg])*d2] +
Conjugate[A1]*B3*k0mimi*Exp[I*(k0mimi - Conjugate[k0plpl])*d2] +
A1*Conjugate[B3]*k0plpl*Exp[I*(k0plpl - Conjugate[k0mimi])*d2] +
Conjugate[A1]*(B1 + A3)*k0*Exp[I*(k0 - Conjugate[k0plpl])*d2] +
A1*Conjugate[B1 + A3]*k0plpl*Exp[I*(k0plpl - Conjugate[k0])*d2] +
Conjugate[B3]*(B1 + A3)*k0*Exp[I*(k0 - Conjugate[k0mimi])*d2] +
B3*Conjugate[B1 + A3]*k0mimi*Exp[I*(k0mimi - Conjugate[k0])*d2] -
Conjugate[A1]*B3*kg1mimi*
Exp[I*(kg1mimi - Conjugate[kg1plpl])*d2] -
A1*Conjugate[B3]*kg1plpl*
Exp[I*(kg1plpl - Conjugate[kg1mimi])*d2] +
Conjugate[A1]*(A3 - B1)*kg1*
Exp[I*(kg1 - Conjugate[kg1plpl])*d2] +
A1*Conjugate[A3 - B1]*kg1plpl*
Exp[I*(kg1plpl - Conjugate[kg1])*d2] +
Conjugate[B3]*(B1 - A3)*kg1*
Exp[I*(kg1 - Conjugate[kg1mimi])*d2] +
B3*Conjugate[B1 - A3]*kg1mimi*
Exp[I*(kg1mimi - Conjugate[kg1])*d2] -
Conjugate[A2]*B4*khmimi*Exp[I*(khmimi - Conjugate[khplpl])*d2] -
A2*Conjugate[B4]*khplpl*Exp[I*(khplpl - Conjugate[khmimi])*d2] +
Conjugate[A2]*(B2 - A4)*kh*Exp[I*(kh - Conjugate[khplpl])*d2] +
A2*Conjugate[B2 - A4]*khplpl*Exp[I*(khplpl - Conjugate[kh])*d2] +
Conjugate[B4]*(A4 - B2)*kh*Exp[I*(kh - Conjugate[khmimi])*d2] +
B4*Conjugate[A4 - B2]*khmimi*Exp[I*(khmimi - Conjugate[kh])*d2] +
Conjugate[A2]*B4*khgmimi*
Exp[I*(khgmimi - Conjugate[khgplpl])*d2] +
A2*Conjugate[B4]*khgplpl*
Exp[I*(khgplpl - Conjugate[khgmimi])*d2] -
Conjugate[A2]*(A4 + B2)*khg*
Exp[I*(khg - Conjugate[khgplpl])*d2] -
A2*Conjugate[A4 + B2]*khgplpl*
Exp[I*(khgplpl - Conjugate[khg])*d2] -
Conjugate[B4]*(A4 + B2)*khg*
Exp[I*(khg - Conjugate[khgmimi])*d2] -
B4*Conjugate[A4 + B2]*khgmimi*
Exp[I*(khgmimi - Conjugate[khg])*d2]];
T2 := 1/4 Re[
Abs[A1]^2*k0plpl*Exp[I*(k0plpl - Conjugate[k0plpl])*d2] +
Abs[A1]^2*kg2plpl*Exp[I*(kg2plpl - Conjugate[kg2plpl])*d2] +
Abs[A2]^2*khplpl*Exp[I*(khplpl - Conjugate[khplpl])*d2] +
Abs[A2]^2*khgplpl*Exp[I*(khgplpl - Conjugate[khgplpl])*d2] +
Abs[B3]^2*k0mimi*Exp[I*(k0mimi - Conjugate[k0mimi])*d2] +
Abs[B3]^2*kg2mimi*Exp[I*(kg2mimi - Conjugate[kg2mimi])*d2] +
Abs[B4]^2*khmimi*Exp[I*(khmimi - Conjugate[khmimi])*d2] +
Abs[B4]^2*khgmimi*Exp[I*(khgmimi - Conjugate[khgmimi])*d2] +
Abs[B1 + A3]^2*k0*Exp[I*(k0 - Conjugate[k0])*d2] +
Abs[B1 - A3]^2*kg2*Exp[I*(kg2 - Conjugate[kg2])*d2] +
Abs[A4 - B2]^2*kh*Exp[I*(kh - Conjugate[kh])*d2] +
Abs[A4 + B2]^2*khg*Exp[I*(khg - Conjugate[khg])*d2] +
Conjugate[A1]*B3*k0mimi*Exp[I*(k0mimi - Conjugate[k0plpl])*d2] +
A1*Conjugate[B3]*k0plpl*Exp[I*(k0plpl - Conjugate[k0mimi])*d2] +
Conjugate[A1]*(B1 + A3)*k0*Exp[I*(k0 - Conjugate[k0plpl])*d2] +
A1*Conjugate[B1 + A3]*k0plpl*Exp[I*(k0plpl - Conjugate[k0])*d2] +
Conjugate[B3]*(B1 + A3)*k0*Exp[I*(k0 - Conjugate[k0mimi])*d2] +
B3*Conjugate[B1 + A3]*k0mimi*Exp[I*(k0mimi - Conjugate[k0])*d2] -
Conjugate[A1]*B3*kg2mimi*
Exp[I*(kg2mimi - Conjugate[kg2plpl])*d2] -
A1*Conjugate[B3]*kg2plpl*
Exp[I*(kg2plpl - Conjugate[kg2mimi])*d2] +
Conjugate[A1]*(A3 - B1)*kg2*
Exp[I*(kg2 - Conjugate[kg2plpl])*d2] +
A1*Conjugate[A3 - B1]*kg2plpl*
Exp[I*(kg2plpl - Conjugate[kg2])*d2] +
Conjugate[B3]*(B1 - A3)*kg2*
Exp[I*(kg2 - Conjugate[kg2mimi])*d2] +
B3*Conjugate[B1 - A3]*kg2mimi*
Exp[I*(kg2mimi - Conjugate[kg2])*d2] -
Conjugate[A2]*B4*khmimi*Exp[I*(khmimi - Conjugate[khplpl])*d2] -
A2*Conjugate[B4]*khplpl*Exp[I*(khplpl - Conjugate[khmimi])*d2] +
Conjugate[A2]*(B2 - A4)*kh*Exp[I*(kh - Conjugate[khplpl])*d2] +
A2*Conjugate[B2 - A4]*khplpl*Exp[I*(khplpl - Conjugate[kh])*d2] +
Conjugate[B4]*(A4 - B2)*kh*Exp[I*(kh - Conjugate[khmimi])*d2] +
B4*Conjugate[A4 - B2]*khmimi*Exp[I*(khmimi - Conjugate[kh])*d2] +
Conjugate[A2]*B4*khgmimi*
Exp[I*(khgmimi - Conjugate[khgplpl])*d2] +
A2*Conjugate[B4]*khgplpl*
Exp[I*(khgplpl - Conjugate[khgmimi])*d2] -
Conjugate[A2]*(A4 + B2)*khg*
Exp[I*(khg - Conjugate[khgplpl])*d2] -
A2*Conjugate[A4 + B2]*khgplpl*
Exp[I*(khgplpl - Conjugate[khg])*d2] -
Conjugate[B4]*(A4 + B2)*khg*
Exp[I*(khg - Conjugate[khgmimi])*d2] -
B4*Conjugate[A4 + B2]*khgmimi*
Exp[I*(khgmimi - Conjugate[khg])*d2]];
G := Re[NIntegrate[
k Sin[x] T1, {x, 0, Pi/2}, {y, -Pi/6, -Pi/6 + afa}] +
NIntegrate[k Sin[x] T2, {x, 0, Pi/2}, {y, -Pi/6 + afa, Pi/6}]];
tmr := {afa, G};
Export["D://ceshi.txt", Table[tmr, {afa, 0, Pi/6, 0.001}], "Table"];
2022年07月26日 03点07分
1
代码如下：
Clear["`*"]
vh = 2*(1.88 - 0.59*Cos[Pi/(1 + d1/(0.65*10^(-9)))]);
u = vh*0.89;
mu = 4.6;
gama = -4*10^(-29)*(1 - Cos[6*afa])*(1 - 1*10^(-8)*I);
d1 = 1.95*10^(-9);
d2 = 1.95*10^(-9);
HBAR = 1.05457266*10^(-34);
ME = 9.1093897*10^(-31);
ELEC = 1.60217733*10^(-19);
Kh = 2.29*10^(10);
kc = Sqrt[2*ME*ELEC/HBAR^2];
k := kc*Sqrt[mu]
k0 := Sqrt[k^2 - k^2 Sin[x]^2]
kh := Sqrt[k^2 - (Kh - k Sin[x] Cos[y])^2 - k^2 Sin[x]^2 Sin[y]^2]
khg := Sqrt[
k^2 - (2*Kh*Sin[afa/2]*Sin[afa/2] -
k Sin[x] Cos[y])^2 - (2*Kh*Sin[afa/2]*Cos[afa/2] +
k Sin[x] Sin[y])^2]
kg1 := Sqrt[
k^2 - (Kh*Cos[Pi/3 - afa] -
k Sin[x] Cos[y])^2 - (Kh*Sin[Pi/3 - afa] + k Sin[x] Sin[y])^2]
kg2 := Sqrt[
k^2 - (Kh*Cos[afa] - k Sin[x] Cos[y])^2 - (k Sin[x] Sin[y] -
Kh*Sin[afa])^2]
k0pl := Sqrt[k^2 - k^2 Sin[x]^2 + kc^2 vh - kc^2*u]
k0mi := Sqrt[k^2 - k^2 Sin[x]^2 - kc^2 vh - kc^2*u]
khpl := Sqrt[
k^2 - (Kh - k Sin[x] Cos[y])^2 - k^2 Sin[x]^2 Sin[y]^2 + kc^2 vh -
kc^2*u]
khmi := Sqrt[
k^2 - (Kh - k Sin[x] Cos[y])^2 - k^2 Sin[x]^2 Sin[y]^2 - kc^2 vh -
kc^2*u]
k0plpl := Sqrt[k^2 - k^2 Sin[x]^2 + 2*kc^2 vh]
k0mimi := Sqrt[k^2 - k^2 Sin[x]^2 - 2*kc^2 vh]
khplpl :=
Sqrt[k^2 - (Kh - k Sin[x] Cos[y])^2 - k^2 Sin[x]^2 Sin[y]^2 +
2*kc^2 vh]
khmimi :=
Sqrt[k^2 - (Kh - k Sin[x] Cos[y])^2 - k^2 Sin[x]^2 Sin[y]^2 -
2*kc^2 vh]
khgplpl :=
Sqrt[k^2 - (2*Kh*Sin[afa/2]*Sin[afa/2] -
k Sin[x] Cos[y])^2 - (2*Kh*Sin[afa/2]*Cos[afa/2] +
k Sin[x] Sin[y])^2 + 2*kc^2 vh]
khgmimi :=
Sqrt[k^2 - (2*Kh*Sin[afa/2]*Sin[afa/2] -
k Sin[x] Cos[y])^2 - (2*Kh*Sin[afa/2]*Cos[afa/2] +
k Sin[x] Sin[y])^2 - 2*kc^2 vh]
kg1plpl :=
Sqrt[k^2 - (Kh*Cos[Pi/3 - afa] -
k Sin[x] Cos[y])^2 - (Kh*Sin[Pi/3 - afa] + k Sin[x] Sin[y])^2 +
2*kc^2 vh]
kg1mimi :=
Sqrt[k^2 - (Kh*Cos[Pi/3 - afa] -
k Sin[x] Cos[y])^2 - (Kh*Sin[Pi/3 - afa] + k Sin[x] Sin[y])^2 -
2*kc^2 vh]
kg2plpl :=
Sqrt[k^2 - (Kh*Cos[afa] - k Sin[x] Cos[y])^2 - (k Sin[x] Sin[y] -
Kh*Sin[afa])^2 + 2*kc^2 vh]
kg2mimi :=
Sqrt[k^2 - (Kh*Cos[afa] - k Sin[x] Cos[y])^2 - (k Sin[x] Sin[y] -
Kh*Sin[afa])^2 - 2*kc^2 vh]
A1 := 1/(1 + I*ME*gama/(HBAR^2*k0pl))*Exp[I*k0pl*d1]/2
B1 := 1/(1 + I*ME*gama/(HBAR^2*k0pl))*Exp[I*k0pl*d1]/2
A2 := 1/(1 + I*ME*gama/(HBAR^2*khpl))*Exp[I*khpl*d1]/2
B2 := 1/(1 + I*ME*gama/(HBAR^2*khpl))*Exp[I*khpl*d1]/2
A3 := 1/(1 + I*ME*gama/(HBAR^2*k0mi))*Exp[I*k0mi*d1]/2
B3 := 1/(1 + I*ME*gama/(HBAR^2*k0mi))*Exp[I*k0mi*d1]/2
A4 := 1/(1 + I*ME*gama/(HBAR^2*khmi))*Exp[I*khmi*d1]/2
B4 := 1/(1 + I*ME*gama/(HBAR^2*khmi))*Exp[I*khmi*d1]/2
T1 := 1/4 Re[
Abs[A1]^2*k0plpl*Exp[I*(k0plpl - Conjugate[k0plpl])*d2] +
Abs[A1]^2*kg1plpl*Exp[I*(kg1plpl - Conjugate[kg1plpl])*d2] +
Abs[A2]^2*khplpl*Exp[I*(khplpl - Conjugate[khplpl])*d2] +
Abs[A2]^2*khgplpl*Exp[I*(khgplpl - Conjugate[khgplpl])*d2] +
Abs[B3]^2*k0mimi*Exp[I*(k0mimi - Conjugate[k0mimi])*d2] +
Abs[B3]^2*kg1mimi*Exp[I*(kg1mimi - Conjugate[kg1mimi])*d2] +
Abs[B4]^2*khmimi*Exp[I*(khmimi - Conjugate[khmimi])*d2] +
Abs[B4]^2*khgmimi*Exp[I*(khgmimi - Conjugate[khgmimi])*d2] +
Abs[B1 + A3]^2*k0*Exp[I*(k0 - Conjugate[k0])*d2] +
Abs[B1 - A3]^2*kg1*Exp[I*(kg1 - Conjugate[kg1])*d2] +
Abs[A4 - B2]^2*kh*Exp[I*(kh - Conjugate[kh])*d2] +
Abs[A4 + B2]^2*khg*Exp[I*(khg - Conjugate[khg])*d2] +
Conjugate[A1]*B3*k0mimi*Exp[I*(k0mimi - Conjugate[k0plpl])*d2] +
A1*Conjugate[B3]*k0plpl*Exp[I*(k0plpl - Conjugate[k0mimi])*d2] +
Conjugate[A1]*(B1 + A3)*k0*Exp[I*(k0 - Conjugate[k0plpl])*d2] +
A1*Conjugate[B1 + A3]*k0plpl*Exp[I*(k0plpl - Conjugate[k0])*d2] +
Conjugate[B3]*(B1 + A3)*k0*Exp[I*(k0 - Conjugate[k0mimi])*d2] +
B3*Conjugate[B1 + A3]*k0mimi*Exp[I*(k0mimi - Conjugate[k0])*d2] -
Conjugate[A1]*B3*kg1mimi*
Exp[I*(kg1mimi - Conjugate[kg1plpl])*d2] -
A1*Conjugate[B3]*kg1plpl*
Exp[I*(kg1plpl - Conjugate[kg1mimi])*d2] +
Conjugate[A1]*(A3 - B1)*kg1*
Exp[I*(kg1 - Conjugate[kg1plpl])*d2] +
A1*Conjugate[A3 - B1]*kg1plpl*
Exp[I*(kg1plpl - Conjugate[kg1])*d2] +
Conjugate[B3]*(B1 - A3)*kg1*
Exp[I*(kg1 - Conjugate[kg1mimi])*d2] +
B3*Conjugate[B1 - A3]*kg1mimi*
Exp[I*(kg1mimi - Conjugate[kg1])*d2] -
Conjugate[A2]*B4*khmimi*Exp[I*(khmimi - Conjugate[khplpl])*d2] -
A2*Conjugate[B4]*khplpl*Exp[I*(khplpl - Conjugate[khmimi])*d2] +
Conjugate[A2]*(B2 - A4)*kh*Exp[I*(kh - Conjugate[khplpl])*d2] +
A2*Conjugate[B2 - A4]*khplpl*Exp[I*(khplpl - Conjugate[kh])*d2] +
Conjugate[B4]*(A4 - B2)*kh*Exp[I*(kh - Conjugate[khmimi])*d2] +
B4*Conjugate[A4 - B2]*khmimi*Exp[I*(khmimi - Conjugate[kh])*d2] +
Conjugate[A2]*B4*khgmimi*
Exp[I*(khgmimi - Conjugate[khgplpl])*d2] +
A2*Conjugate[B4]*khgplpl*
Exp[I*(khgplpl - Conjugate[khgmimi])*d2] -
Conjugate[A2]*(A4 + B2)*khg*
Exp[I*(khg - Conjugate[khgplpl])*d2] -
A2*Conjugate[A4 + B2]*khgplpl*
Exp[I*(khgplpl - Conjugate[khg])*d2] -
Conjugate[B4]*(A4 + B2)*khg*
Exp[I*(khg - Conjugate[khgmimi])*d2] -
B4*Conjugate[A4 + B2]*khgmimi*
Exp[I*(khgmimi - Conjugate[khg])*d2]];
T2 := 1/4 Re[
Abs[A1]^2*k0plpl*Exp[I*(k0plpl - Conjugate[k0plpl])*d2] +
Abs[A1]^2*kg2plpl*Exp[I*(kg2plpl - Conjugate[kg2plpl])*d2] +
Abs[A2]^2*khplpl*Exp[I*(khplpl - Conjugate[khplpl])*d2] +
Abs[A2]^2*khgplpl*Exp[I*(khgplpl - Conjugate[khgplpl])*d2] +
Abs[B3]^2*k0mimi*Exp[I*(k0mimi - Conjugate[k0mimi])*d2] +
Abs[B3]^2*kg2mimi*Exp[I*(kg2mimi - Conjugate[kg2mimi])*d2] +
Abs[B4]^2*khmimi*Exp[I*(khmimi - Conjugate[khmimi])*d2] +
Abs[B4]^2*khgmimi*Exp[I*(khgmimi - Conjugate[khgmimi])*d2] +
Abs[B1 + A3]^2*k0*Exp[I*(k0 - Conjugate[k0])*d2] +
Abs[B1 - A3]^2*kg2*Exp[I*(kg2 - Conjugate[kg2])*d2] +
Abs[A4 - B2]^2*kh*Exp[I*(kh - Conjugate[kh])*d2] +
Abs[A4 + B2]^2*khg*Exp[I*(khg - Conjugate[khg])*d2] +
Conjugate[A1]*B3*k0mimi*Exp[I*(k0mimi - Conjugate[k0plpl])*d2] +
A1*Conjugate[B3]*k0plpl*Exp[I*(k0plpl - Conjugate[k0mimi])*d2] +
Conjugate[A1]*(B1 + A3)*k0*Exp[I*(k0 - Conjugate[k0plpl])*d2] +
A1*Conjugate[B1 + A3]*k0plpl*Exp[I*(k0plpl - Conjugate[k0])*d2] +
Conjugate[B3]*(B1 + A3)*k0*Exp[I*(k0 - Conjugate[k0mimi])*d2] +
B3*Conjugate[B1 + A3]*k0mimi*Exp[I*(k0mimi - Conjugate[k0])*d2] -
Conjugate[A1]*B3*kg2mimi*
Exp[I*(kg2mimi - Conjugate[kg2plpl])*d2] -
A1*Conjugate[B3]*kg2plpl*
Exp[I*(kg2plpl - Conjugate[kg2mimi])*d2] +
Conjugate[A1]*(A3 - B1)*kg2*
Exp[I*(kg2 - Conjugate[kg2plpl])*d2] +
A1*Conjugate[A3 - B1]*kg2plpl*
Exp[I*(kg2plpl - Conjugate[kg2])*d2] +
Conjugate[B3]*(B1 - A3)*kg2*
Exp[I*(kg2 - Conjugate[kg2mimi])*d2] +
B3*Conjugate[B1 - A3]*kg2mimi*
Exp[I*(kg2mimi - Conjugate[kg2])*d2] -
Conjugate[A2]*B4*khmimi*Exp[I*(khmimi - Conjugate[khplpl])*d2] -
A2*Conjugate[B4]*khplpl*Exp[I*(khplpl - Conjugate[khmimi])*d2] +
Conjugate[A2]*(B2 - A4)*kh*Exp[I*(kh - Conjugate[khplpl])*d2] +
A2*Conjugate[B2 - A4]*khplpl*Exp[I*(khplpl - Conjugate[kh])*d2] +
Conjugate[B4]*(A4 - B2)*kh*Exp[I*(kh - Conjugate[khmimi])*d2] +
B4*Conjugate[A4 - B2]*khmimi*Exp[I*(khmimi - Conjugate[kh])*d2] +
Conjugate[A2]*B4*khgmimi*
Exp[I*(khgmimi - Conjugate[khgplpl])*d2] +
A2*Conjugate[B4]*khgplpl*
Exp[I*(khgplpl - Conjugate[khgmimi])*d2] -
Conjugate[A2]*(A4 + B2)*khg*
Exp[I*(khg - Conjugate[khgplpl])*d2] -
A2*Conjugate[A4 + B2]*khgplpl*
Exp[I*(khgplpl - Conjugate[khg])*d2] -
Conjugate[B4]*(A4 + B2)*khg*
Exp[I*(khg - Conjugate[khgmimi])*d2] -
B4*Conjugate[A4 + B2]*khgmimi*
Exp[I*(khgmimi - Conjugate[khg])*d2]];
G := Re[NIntegrate[
k Sin[x] T1, {x, 0, Pi/2}, {y, -Pi/6, -Pi/6 + afa}] +
NIntegrate[k Sin[x] T2, {x, 0, Pi/2}, {y, -Pi/6 + afa, Pi/6}]];
tmr := {afa, G};
Export["D://ceshi.txt", Table[tmr, {afa, 0, Pi/6, 0.001}], "Table"];