擅长:python、mysql、java
<p>{{55y>在这种情况下,{55y}可以用来定义cd1}中的临界点。在</p>
<p>数组是构造的,但是如果在<code>numpy.transpose</code>步骤中绘制的常量值有问题,则可能不具有相同的长度。因此,在第一个<code>matplotlib.pyplot</code>之前传递一个条件。在</p>
<pre><code>from sympy import *
import numpy as np
import matplotlib.pyplot as plt
x = symbols('x')
equation = x**2 + 2 #OK = x*exp(-x**2/10) OR x**2 + 2
firstDeriv = equation.diff(x)
secondDeriv = firstDeriv.diff(x)
criticalPoints = list(solveset(firstDeriv, x))
criticalPointsY = [equation.subs(x, a) for a in criticalPoints]
xx = np.linspace(-10, 10, 100)
lam_f= lambdify(x, [equation, firstDeriv, secondDeriv])
yy=[elem if type(elem) == np.ndarray else np.full(len(xx), elem) for elem in lam_f(xx)]
plt.plot(xx, np.transpose(yy))
plt.plot(criticalPoints, criticalPointsY, 'k*')
plt.show()
</code></pre>
<p><a href="https://i.stack.imgur.com/EgJUg.png" rel="nofollow noreferrer">enter image description here</a></p>