Python Sympy导子

def dyson(x, B_0, delta_B, alpha):
  return (((delta_B + alpha*(x - B_0))/(4*(x - B_0)**2 + delta_B**2)) 
  + ((delta_B - alpha*(x + B_0))/(4*(x + B_0)**2 + delta_B**2)))

x = Symbol('x')
delta_B = Symbol('delta_B')
B_0 = Symbol('B_0')
alpha = Symbol('alpha')
y = (((delta_B + alpha*(x - B_0))/(4*(x - B_0)**2 + delta_B**2)) 
  + ((delta_B - alpha*(x + B_0))/(4*(x + B_0)**2 + delta_B**2)))
der_dyson = lambdify((x, delta_B, B_0, alpha), y.diff(x))

def der_dyson_(x, B_0, delta_B, alpha, h=1):
  return (dyson(x + h, B_0, delta_B, alpha) - dyson(x - h, B_0, delta_B, alpha))/(2*h)
fig_test, ax_test = plt.subplots()
x_test = np.linspace(-600, 600, 10000)
ax_test.plot(x_test, der_dyson_(x_test, 200, 250, 0.5, h=0.12), label='finite diff')
ax_test.plot(x_test, der_dyson(x_test, 200, 250, 0.5), label='Sympy')
ax_test.legend()

der_dyson = lambdify((x, delta_B, B_0, alpha), y.diff(x))

def dyson(x, B_0, delta_B, alpha):
    ...

def get_dyson_deriv():
    x = Symbol('x')
    delta_B = Symbol('delta_B')
    B_0 = Symbol('B_0')
    alpha = Symbol('alpha')
    y = dyson(x, B_0, delta_B, alpha)
    # NOTE order of parameters in the tuple here:
    return lambdify((x, B_0, delta_B, alpha), y.diff(x))

der_dyson = get_dyson_deriv()

Python相关问答推荐

如何避免Chained when/then分配中的Mypy不兼容类型警告？

运行Python脚本时，用作命令行参数的SON文本

发生异常：TclMessage命令名称无效.！listbox"

PMMLPipeline._ fit()需要2到3个位置参数，但给出了4个位置参数

如何在Raspberry Pi上检测USB并使用Python访问它？

如何创建一个缓冲区周围的一行与manim？

Scrapy和Great Expectations(great_expectations)—不合作

什么是最好的方法来切割一个相框到一个面具的第一个实例？

在www.example.com中使用`package_data`包含不包含__init__. py的非Python文件

需要帮助重新调整python fill_between与数据点

为什么numpy. vectorize调用vectorized函数的次数比vector中的元素要多？

Flash只从html表单中获取一个值

在极中解析带有数字和SI前缀的字符串

python panda ExcelWriter切换动态公式到数组公式

在Python中使用yaml渲染(多行字符串)

(Python/Pandas)基于列中非缺失值的子集DataFrame

使用__json__的 pyramid 在客户端返回意外格式

时长超过24小时如何从Excel导入时长数据

Stats.ttest_ind：提取df值

Groupby并在组内比较单独行上的两个时间戳