関数電卓・グラフ電卓総合 Part3

関数電卓・グラフ電卓 ..

104:名無しさん＠３周年
19/01/29 22:28:49.21 cK9db8kv.net
あの程度の精度は出せるでしょ
出せないとしたら処理系の選択ミス
pi@raspi ~ $ cat simpson.py
from math import *
def f(x):
return x*(sin(x)+cos(x)+1.0)
N = 1000000
a = 0.0
b = 10.0**4.0
h = (b-a)/N
S = (h/3) * sum((f(h*i) + 4*f(h*(i+1)) + f(h*(i+2))) for i in range(0,N-1, 2))
print(S)
pi@raspi ~ $ python simpson.py
50006463.152

次ページ

続きを表示

1を表示