主题:同一个程序运行得到不同的结果,怎么办?
一个程序,运行了几次,结果发现并不是每一次的结果都相同,请问高手,这怎么办?下面是程序,和摘录的两个结果,请帮帮忙!先行谢过啦!
clear all;
syms C x beta3 beta9 tau delta delta2 beta4 beta8
C=(4*beta3)^(tau)*sqrt(beta9/pi)*gamma(tau+1/2)/gamma(tau);
eq1=(beta9*x^2+4*beta3)^(-tau-1/2)/(1+x^2);
i1=2*C*int(eq1,x,0,inf);
i11=simple(i1)
eq2=(beta9*x^2+4*beta3)^(-tau-1/2)/(1+x^2)^2;
i2=2*C*int(eq2,x,0,inf);
i22=simple(i2)
下面是两个得到的结果,除了这两个结果还有其他的形式
结果一:
i11 =
1/2*(4*beta3)^tau*(beta9/pi)^(1/2)/gamma(tau)*2^(-2*tau)*beta3^(-tau-1/2)*(-pi*gamma(tau+1/2)/(1/4*beta9/beta3-1)*(1-1/4*beta9/beta3)^(-tau+1/2)-(pi*beta9/beta3)^(1/2)*gamma(1+tau)*hypergeom([1, 1+tau],[3/2],1/4*beta9/beta3))
i22 =
2*(4*beta3)^tau*(beta9/pi)^(1/2)/gamma(tau)*beta3^(-tau-1/2)*2^(-2*tau-1)/(beta9/beta3)^(1/2)*(pi^(1/2)*gamma(tau)*hypergeom([1/2, 2],[1-tau],4*beta3/beta9)+pi*csc(pi*tau)*(-1-tau)*gamma(tau+1/2)*2^(2*tau)*(beta9/beta3)^(-tau)*hypergeom([tau+1/2, 2+tau],[1+tau],4*beta3/beta9))
结果二:
i11 =
2*(4*beta3)^tau*(beta9/pi)^(1/2)/gamma(tau)*beta3^(-tau-1/2)*2^(-1-2*tau)/(beta9/beta3)^(1/2)*(pi^(1/2)*gamma(tau)*hypergeom([1/2, 1,[1-tau],4*beta3/beta9)+pi*csc(pi*tau)*gamma(tau+1/2)*2^(2*tau)*(beta9/beta3)^(-tau)/(4*beta3/beta9-1)*(1-4*beta3/beta9)^(-tau+1/2))
i22 =
1/3*(4*beta3)^tau*(beta9/pi)^(1/2)*2^(-2-2*tau)*beta3^(-tau-1/2)*(3*pi*gamma(tau+1/2)*hypergeom([1/2, tau+1/2],[-1/2],1/4*beta9/beta3)+pi^(1/2)*gamma(2+tau)*(beta9/beta3)^(3/2)*hypergeom([2, 2+tau],[5/2],1/4*beta9/beta3))/gamma(tau)
clear all;
syms C x beta3 beta9 tau delta delta2 beta4 beta8
C=(4*beta3)^(tau)*sqrt(beta9/pi)*gamma(tau+1/2)/gamma(tau);
eq1=(beta9*x^2+4*beta3)^(-tau-1/2)/(1+x^2);
i1=2*C*int(eq1,x,0,inf);
i11=simple(i1)
eq2=(beta9*x^2+4*beta3)^(-tau-1/2)/(1+x^2)^2;
i2=2*C*int(eq2,x,0,inf);
i22=simple(i2)
下面是两个得到的结果,除了这两个结果还有其他的形式
结果一:
i11 =
1/2*(4*beta3)^tau*(beta9/pi)^(1/2)/gamma(tau)*2^(-2*tau)*beta3^(-tau-1/2)*(-pi*gamma(tau+1/2)/(1/4*beta9/beta3-1)*(1-1/4*beta9/beta3)^(-tau+1/2)-(pi*beta9/beta3)^(1/2)*gamma(1+tau)*hypergeom([1, 1+tau],[3/2],1/4*beta9/beta3))
i22 =
2*(4*beta3)^tau*(beta9/pi)^(1/2)/gamma(tau)*beta3^(-tau-1/2)*2^(-2*tau-1)/(beta9/beta3)^(1/2)*(pi^(1/2)*gamma(tau)*hypergeom([1/2, 2],[1-tau],4*beta3/beta9)+pi*csc(pi*tau)*(-1-tau)*gamma(tau+1/2)*2^(2*tau)*(beta9/beta3)^(-tau)*hypergeom([tau+1/2, 2+tau],[1+tau],4*beta3/beta9))
结果二:
i11 =
2*(4*beta3)^tau*(beta9/pi)^(1/2)/gamma(tau)*beta3^(-tau-1/2)*2^(-1-2*tau)/(beta9/beta3)^(1/2)*(pi^(1/2)*gamma(tau)*hypergeom([1/2, 1,[1-tau],4*beta3/beta9)+pi*csc(pi*tau)*gamma(tau+1/2)*2^(2*tau)*(beta9/beta3)^(-tau)/(4*beta3/beta9-1)*(1-4*beta3/beta9)^(-tau+1/2))
i22 =
1/3*(4*beta3)^tau*(beta9/pi)^(1/2)*2^(-2-2*tau)*beta3^(-tau-1/2)*(3*pi*gamma(tau+1/2)*hypergeom([1/2, tau+1/2],[-1/2],1/4*beta9/beta3)+pi^(1/2)*gamma(2+tau)*(beta9/beta3)^(3/2)*hypergeom([2, 2+tau],[5/2],1/4*beta9/beta3))/gamma(tau)