Clear[lbk]"`*"[rbk]a[lbk]x_[rbk] := Cos[lbk]Pi*x[rbk]; A = E^Sqrt[lbk]1/2[rbk]; B = 0; j = 1; M = 7; (*u=\[lbk]ExponentialE[rbk]^Sqrt[lbk]1/2-x[rbk]*)uuu[lbk]x[rbk] = E^Sqrt[lbk]1/2[rbk] + b0;f[lbk]x_[rbk] := -(Pi*Sin[lbk]Pi*x[rbk] + Cos[lbk]Pi*x[rbk]*1/(2 Sqrt[lbk](1/2 - x)[rbk]))*uuu[lbk]x[rbk];step = 1;(*ai*)fs = Normal[lbk]Series[lbk](f[lbk]x[rbk]) /. x -> (1/2 - x), {x, 0, 2}[rbk][rbk];gamma = Exponent[lbk]fs, x, List[rbk];b = gamma[lbk][lbk]j[rbk][rbk];c = Coefficient[lbk]fs, x, b[rbk];imax = IntegerPart[lbk]M - 2 - b[rbk];alpj = Table[lbk] 2 + b + 2*i, {i, 0, imax}[rbk];(*\[lbk]Alpha[rbk]i=2+i+\[lbk]CapitalUpsilon[rbk]j(b),i=0,1,...,M-(1 - \\[lbk]Alpha[rbk] + \[lbk]CapitalUpsilon[rbk]j(b)*)a = Table[lbk]((-1)^i EulerE[lbk]2 i[rbk]* Pi^(2 i))/(2 i)! (c/((1 + b) alpj[lbk][lbk]i + 1[rbk][rbk])), {i, 0, imax}[rbk];uuu[lbk]x[rbk] = uuu[lbk]x[rbk] + a.(1/2 - x)^alpj;While[lbk] b < M - 2, j++;nx = Ceiling[lbk]b + step[rbk];(*ceiling为取整函数,级数展开次数小于等于nx*) fs = Normal[lbk]Series[lbk](f[lbk]x[rbk]) /. x -> (1/2 - x), {x, 0, nx}[rbk][rbk];While[lbk] Length[lbk]fs[rbk] < j, nx++; fs = Normal[lbk]Series[lbk](f[lbk]x[rbk]) /. x -> (1/2 - x), {x, 0, nx}[rbk][rbk] // Simplify [rbk];(*右端项的展开项数足够多,保证取到第j项的值*) gamma = Exponent[lbk]fs, x, List[rbk];(*取展开项x的指数所构成数组*) (*step=Max[lbk]beta[lbk][lbk]j[rbk][rbk]-beta[lbk][lbk]j-1[rbk][rbk],beta[lbk][lbk]j+1[rbk][rbk]-beta[lbk][lbk]j[rbk][rbk][rbk];*)\
b = gamma[lbk][lbk]j[rbk][rbk];(*第j项指数*) c = Coefficient[lbk]fs, x, b[rbk];(*第j项的系数*)imax = IntegerPart[lbk]M - 2 - b[rbk]; alpj = Table[lbk]2 + b + 2*i, {i, 0, imax}[rbk]; a = Table[lbk]((-1)^i EulerE[lbk]2 i[rbk]* Pi^(2 i))/(2 i)! (c/((1 + b) alpj[lbk][lbk]i + 1[rbk][rbk])), {i, 0, imax}[rbk]; uuu[lbk]x[rbk] = uuu[lbk]x[rbk] + a.(1/2 - x)^alpj;[rbk](uuu[lbk]x[rbk] /. b0 -> (1 - E^(1/Sqrt[lbk]2[rbk]))) /. x -> (1/2 - x)
得到的结果与预设的结果不一样尤其是次数小于等于1的没有,想问问怎么修改,已经调试很久了代码
b = gamma[lbk][lbk]j[rbk][rbk];(*第j项指数*) c = Coefficient[lbk]fs, x, b[rbk];(*第j项的系数*)imax = IntegerPart[lbk]M - 2 - b[rbk]; alpj = Table[lbk]2 + b + 2*i, {i, 0, imax}[rbk]; a = Table[lbk]((-1)^i EulerE[lbk]2 i[rbk]* Pi^(2 i))/(2 i)! (c/((1 + b) alpj[lbk][lbk]i + 1[rbk][rbk])), {i, 0, imax}[rbk]; uuu[lbk]x[rbk] = uuu[lbk]x[rbk] + a.(1/2 - x)^alpj;[rbk](uuu[lbk]x[rbk] /. b0 -> (1 - E^(1/Sqrt[lbk]2[rbk]))) /. x -> (1/2 - x)
得到的结果与预设的结果不一样尤其是次数小于等于1的没有,想问问怎么修改,已经调试很久了代码