Quantcast
Channel: 円周率近似値の日に生まれて理系じゃないわけないだろ! - knifeのblog
Viewing all articles
Browse latest Browse all 5376

tan xのマクローリン展開をやらないの?

$
0
0

sin xとcos xのマクローリン展開のアニメGIFの記事が、なぜか人気がある。

なんで、tan xのマクローリン展開をやらないの?と自分へと問いかけてみる。

出来なくはないんだけど、高次方程式の各項の係数の表現が、sin xやcos xのように簡潔に表せないからである。

そこを簡潔に表せたら、やってもいいかな。

係数の分母は階乗になるんだけど分子が面倒なんです。

f(x) = tan(x)

f'(x) = (tan(x))' = 1/cos2(x) = 1 + tan2(x)
f'(0) = 1

f''(x) = (tan(x))'' = (1 + tan2(x))' = (1)' + (tan2(x))' = (0 + 2tan(x))(1 + tan2(x)) = 2tan(x) + 2tan3(x)
f''(0) = 0

f'''(x) = (2tan(x) + 2tan3(x))'(1 + tan2(x)) = (2 + 6tan(x))(1 + tan2(x)) = 2 + 8tan2(x) + 6tan4(x)
f'''(0) = 2

f''''(x) = (2 + 8tan2(x) + 6tan4(x))'(1 + tan2(x)) = (16tan(x) + 24tan3(x))(1 + tan2(x)) = 16tan(x) + 40tan3(x) + 24tan5(x)
f''''(0) = 0

f'''''(x) = (16tan(x) + 40tan3(x) + 24tan5(x))'(1 + tan2(x)) = (16 + 120tan2(x) + 120tan4(x))(1 + tan2(x)) = 16 + 136tan2(x) + 240tan4(x) + 120tan6(x)
f'''''(0) = 16

f''''''(x) = (16 + 136tan2(x) + 240tan4(x) + 120tan6(x))'(1 + tan2(x)) = (272tan(x) + 960tan3(x) + 720tan5(x))(1 + tan2(x)) = 272tan(x) + 1232tan3(x) + 1680tan5(x) + 720tan7(x)
f''''''(0) = 0

f'''''''(x) = (272tan(x) + 1232tan3(x) + 1680tan5(x) + 720tan7(x))'(1 + tan2(x)) = (272 + 3696tan2(x) + 8400tan4(x) + 5040tan6(x))(1 + tan2(x)) = 272 + 3968tan2(x) + 12096tan4(x) + 13440tan6(x) + 5040tan8(x)
f'''''''(0) = 272

7階微分まで手でやってみたが、面倒だな。
まぁ、プログラムでやるならば問題はないけど。

偶数階微分にx=0を代入すると0になるので、奇数階微分のx=0の値が分子に来て、分母は奇数階乗です。

tan x = (1/1!)x + (2/3!)x3 + (16/5!)x5 + (272/7!)x7 + …

まぁ、分子の計算式を書かないならば、どうにか出来そうな気がしてきたな。



こんなような計算になるだろうか。


ではでは


Viewing all articles
Browse latest Browse all 5376

Trending Articles