作业帮 > 综合 > 作业

用mathematica软件编写牛顿-科特斯公式遇到的问题,请高手指教,

来源:学生作业帮 编辑:百度作业网作业帮 分类:综合作业 时间:2024/05/29 20:56:22
用mathematica软件编写牛顿-科特斯公式遇到的问题,请高手指教,
NC[f_, {x_, a_, b_}, n_] :=
 Module[{Z = b - a, C},
  C = Product[(t - j)/(k - j), {j, 0, n, 1}];
  Z*Sum[1/n*Q*(f /. x ->k), {k, 0, n}]]
-------------------------上面函数中C表达式有问题,不知道如何表达-------------------------

 
C式中,分母的j是不能等于k的,但是写成函数语言时却不知道如何表达,求高手指教.
---------PS:我尝试把函数分段写了,也不行,如图--------------

语言如下:
NC[f_, {x_, a_, b_}, n_] :=
 Module[{Z = b - a, W, Q, M},
  Q = Product[(t - j)/(k - j), {j, 0, k - 1, 1}];
  W = Product[(t - j)/(k - j), {j, k + 1, n, 1}];
  M = Integrate[Q*W, {t, 0, n}];
  Z*Sum[1/n*Q*(f /. x-> k), {k, 0, n}]]
用mathematica软件编写牛顿-科特斯公式遇到的问题,请高手指教,
这里的第二页和第三页就是你的公式吧,怎么感觉有点区别啊?你这函数最后要输出啥啊?另,Newton-Cotes算法在mathematica里是有内建函数的:
再问: 我其实是想用牛顿-科特斯公式求一个含有参数的积分方程,方程如下。Q,R,A,L,G,N都是常数,我想积分出含有这些参数的一个方程出来。用NIntegrate貌似只能积不含参数的方程,对于下面这张方程束手无策啊。