Python 代码如下:
import numpy as np
from math import sqrt
import sympy
from matplotlib import pyplot as plt...):
f_np = sympy.lambdify((x,y_x), f_xy, 'numpy')
x_vec = np.linspace(x_lim[0], x_lim[1], n)...y_vec = np.linspace(y_lim[0], y_lim[1], n)
dx = x_vec[1] - x_vec[0]
dy = y_vec[1] - y_vec[0]...(sympy.Eq(y(x).diff(x), f_xy))),fontsize=16)
return ax
if __name__ == '__main__':
x = sympy.symbols...('x')
y = sympy.Function('y')
fig, axes = plt.subplots(1, 3, figsize=(24, 8))
plot_directtion_field