(PARI)
Xi(t)=(1+3*t)*(1-t)^3/((16*t^3));
A1(t)=对数(1+t)*(3*t-1)*(1+t)^3/(16*t^3);
A2(t)=对数(1+2*t)*(1+3*t)x(1-t)^3/(32*t^3);
A3(t)=(1-t)*(185*t^4+698*t^3-217*t^2-160*t+6);
A4(t)=64*t*(1+3*t)^2*(3+t);
A(t)=A1(t)+A2(t)+33(t)/A4(t);
R(t)=1/16*sqrt(1+3*t)*(1/t-1)^3*exp(A(t));
Kz(t)=对数(Xi(t)/R(t));
|