Производная Сципи

У меня есть вопрос о производной функции Scipy. Я использовал его прошлой ночью и получил несколько странных ответов. Сегодня утром я попробовал еще раз с некоторыми простыми функциями и получил несколько правильных ответов, а некоторые — неправильных. Вот мои тесты:

In [1]: def poly1(x):
...:     return x**2

In [3]: derivative(poly1, 0)
Out[3]: 0.0

In [4]: def poly2(x):
...:    return (x-3)**2

In [6]: derivative(poly2, 3)
Out[6]: 0.0

In [8]: def sin1(x):
...:     return sin(x)

In [14]: derivative(sin1, pi/2)
Out[14]: 5.5511151231257827e-17

In [15]: def poly3(x):
....:     return 3*x**4 + 2*x**3 - 10*x**2 + 15*x - 2

In [19]: derivative(poly3, -2)
Out[19]: -39.0

In [20]: derivative(poly3, 2)
Out[20]: 121.0

In [22]: derivative(poly3, 0)
Out[22]: 17.0

Я проверил значения poly3 вручную и -2 = 17, 2 = 95, 0 = 15. Значит, я неправильно использую функцию, или с функцией что-то не так. Спасибо

Используя :Python 2.7.3, IPython 0.12.1, Numpy 1.6.1, Scipy 0.9.0, Linux Mint 13

6
задан Hugh Bothwell 13 July 2012 в 15:05
поделиться