二次関数の接線が水平になってしまう件

Question

pythonで、任意の二次関数と、その接線を描いてみたのですが...思うようになりません。
アドバイス、よろしくお願いいたします。

import numpy as np
import matplotlib.pyplot as plt
plt.rcParams['font.family'] = 'IPAPGothic' 

#x軸の範囲
x=np.arange(0,10,0.1)

#任意の２次関数
y1=0.01*x**2+0.1*x

# 微分の関数
def derivative(f,a,h):
    return  (f(a+h)-f(a-h))/2*h
h = 1e-4 # 0.0001

#  ５におけるf(x)=0.01*x**2+0.1*xの接線
a=5
m=derivative (lambda x: 0.01*x**2 + 0.1*x,a,h)
b=0.01*a**2 + 0.1*a
y2=m*(x-a)+b

plt.plot(x,y1,label='二次曲線')
plt.plot(x,y2,label='x=５における接線')
plt.xlabel('x ')
plt.ylabel('y ')
plt.legend()
plt.show()

Comments

[y_waiwai]

どうなればいいんでしょうか

[okahijiki]

print(m)
1.9999999999908982e-09

print(a)
5

print(b)
0.75

#微分係数1.9999999999908982e-09、点(5, 0.75）を通る
y2=m*(x-a)+b
となります。

[jbpb0]

return (f(a+h)-f(a-h))/2*h
↓ 分母にカッコ付ける
return (f(a+h)-f(a-h))/(2*h)
かな

[okahijiki]

そうですね。微分係数が小さすぎるので、（ほぼ）平行になってしまうのですね。失礼しました。

[okahijiki]

jbpb0さま
解決しました。ありがとうございます！

[ozwk]

...という回答をつけて質問を閉じたほうがいいと思います

Answer 1

    return  (f(a+h)-f(a-h))/2*h

↓ 修正 (分母にカッコを付ける)

    return  (f(a+h)-f(a-h))/(2*h)