日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

柯特斯系数计算

發布時間:2025/4/16 编程问答 18 豆豆
生活随笔 收集整理的這篇文章主要介紹了 柯特斯系数计算 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

柯特斯系數理論知識在下面的這個鏈接中(也是我的blog)

牛頓-柯特斯公式


代碼

用上面的blog的公式生成對應的柯特斯系數,運用代碼如下:

from sympy import *def C(n=int(), k=int()):if (n - k) % 2 == 0:ans = 1else:ans = -1for j in range(n+1):if j != k:ans *= (x - j)ans = integrate(ans, (x, 0, n)) / nfor i in range(k):ans /= (i+1)for i in range(n-k):ans /= (i+1)print(ans, end=' ')if __name__ == '__main__':x = symbols('x')for i in range(8):for j in range(i+2):C(i+1, j)print()

輸出

(通過與書上進行對比,發現正確(《數值分析(第五版)》李慶陽 P104))

1/2 1/2 1/6 2/3 1/6 1/8 3/8 3/8 1/8 7/90 16/45 2/15 16/45 7/90 19/288 25/96 25/144 25/144 25/96 19/288 41/840 9/35 9/280 34/105 9/280 9/35 41/840 751/17280 3577/17280 49/640 2989/17280 2989/17280 49/640 3577/17280 751/17280 989/28350 2944/14175 -464/14175 5248/14175 -454/2835 5248/14175 -464/14175 2944/14175 989/28350 《新程序員》:云原生和全面數字化實踐50位技術專家共同創作,文字、視頻、音頻交互閱讀

總結

以上是生活随笔為你收集整理的柯特斯系数计算的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。