Python Sympy计算变限积分积分函数含有对数（log）时，遇到的问题。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
import sympy



Z, x, t ,H1 ,H2,H3,H4= sympy.symbols('Z,x,t,H1,H2,H3,H4')



z11 = Z ** 3 + 0.01



H1 = Z ** 2 * (- (
        -16 * 41.3399927 * Z ** 3 - 48 * 41.3399927 * sympy.log(Z) - Z * (16 * -1634.89996 - 144 * 41.3399927) - (
                -240 * 41.3399927 * Z ** 9 - 120 * 41.3399927 * Z ** 6 - 20 * 41.3399927 * Z ** 3 - 15 * 41.3399927 +
                Z ** 11 * (-60 * 190592.559 + 360 * -1634.89996 - 1620 * 41.3399927) + Z ** 10 * (
                            -120 * -1634.89996 + 1080 * 41.3399927) + Z ** 8 * (
                            -15 * 190592.559 + 90 * -1634.89996 - 405 * 41.3399927) + Z ** 7 * (
                            -24 * -1634.89996 + 216 * 41.3399927) + Z ** 4 * (-15 * -1634.89996 + 135 * 41.3399927)) / (
                        15 * Z ** 12)) / 190592.559 + (
                          -15.8415841584158 * 41.3399927 * z11 ** 1.0 - 48 * 41.3399927 *
                          sympy.log(0.996688717477339 * z11 ** (1 / 3)) - 0.996688717477339 * (
                                      16 * -1634.89996 - 144 * 41.3399927) * z11 ** (
                                      1 / 3) - 0.0693736006666667 * (-19.8019801980198 * 41.3399927 * z11
                                                                     ** 1.0 - 117.63552592883 * 41.3399927 * z11 ** 2.0
                                                                     - 232.941635502635 * 41.3399927 * z11 **
                                                                     3.0 - 15 * 41.3399927 + 0.967376249734144 *
                                                                     (-120 * -1634.89996 + 1080 * 41.3399927) *
                                                                     z11 ** (10 / 3) + 0.977050012231486 *
                                                                     (-24 * -1634.89996 + 216 * 41.3399927) *
                                                                     z11 ** (7 / 3) + 0.986820512353801 *
                                                                     (-15 * -1634.89996 + 135 * 41.3399927) *
                                                                     z11 ** (4 / 3) + 0.973814723602218 *
                                                                     z11 ** 2.66666666666667 *
                                                                     (-15 * 190592.559 + 90 * -1634.89996 -
                                                                      405 * 41.3399927) + 0.964172993665562
                                                                     * z11 ** 3.66666666666667 *
                                                                     (-60 * 190592.559 + 360 * -1634.89996 -
                                                                      1620 * 41.3399927)) / z11 ** 4.0) / 190592.559)






print('**************************************************')

I1 = sympy.integrate(H1, (Z, 1, x))
print("积分为：", I1)
print('****************')