差商表
import numpy as npfrom sympy import *np.set_printoptions(suppress=True)def cs(x, f):n=len(x)res=np.ones([n,n])*np.inffor i in range(0, n):res[i][0] = f[i]for j in range(1, n):for i in range(j, n):res[i][j] = round((res[i][j - 1] - res[i - 1][j - 1]) / (x[i] - x[i - j]),5)return resif __name__=="__main__":x=np.array([0.4,0.55,0.65,0.80,0.90,1.05])f=np.array([0.41075,0.57815,0.69675,0.88811,1.02652,1.25382])print(cs(x,f)) 实现效果:
【python 均差表 差商表 数值分析】[[ 0.41075infinfinfinfinf] [ 0.578151.116infinfinfinf] [ 0.696751.1860.28infinfinf] [ 0.888111.275730.358920.1973infinf] [ 1.026521.38410.433480.213030.03146inf] [ 1.253821.515330.524920.22860.03114 -0.00049]]
- 春季老年人吃什么养肝?土豆、米饭换着吃
- 三八妇女节节日祝福分享 三八妇女节节日语录
- 老人谨慎!选好你的“第三只脚”
- 校方进行了深刻的反思 青岛一大学生坠亡校方整改校规
- 脸皮厚的人长寿!有这特征的老人最长寿
- 长寿秘诀:记住这10大妙招 100%增寿
- 春季老年人心血管病高发 3条保命要诀
- 眼睛花不花要看四十八 老年人怎样延缓老花眼
- 香槟然能防治老年痴呆症? 一天三杯它人到90不痴呆
- 老人手抖的原因 为什么老人手会抖
