我建议你帮个小忙。我需要创建在求积规则中计算定积分的简单应用程序,它是在运行时作为参数给出的。我从添加参数开始,但当我在其中使用括号时,编译器不会编译我的应用程序。然后我陷入了入门,如何将参数作为积分传递到一些计算部分?
谢谢
编辑:
int x;
int i;
float sum=0;
double PI;
printf("Please enter the parts to divide the interval: ");
scanf("%d", &x);
for (i=1; i<x; i+=2){
if (i%4==1)
问题是,我有德拜公式,需要用辛普森法则来写一个函数cv(T)来计算给定温度下的Cv。那么Cv = 9*V_p_k_B*(T/theta_D)3 (积分从0到θ_D/T)x_4*ex/ (ex -1)*2,对于这个积分,如何用辛普森方法建立一个函数来求积分?这个公式的积分是θ_D/T的0。
这是我到目前为止所拥有的
from __future__ import division, print_function
from math import e
import numpy as np
# constants
V = 1000 # cm**3 of solid al
我试图把一种方法放在一起计算一个不规则的体积,但是用其他凸的:它使用三角剖分将多面体分裂成多个子四面体(单纯形),独立计算体积,然后把所有的子体积值汇总起来。
然而,在我的测试中,我得到了以下单元立方体的奇怪结果。有人知道臭虫住在哪里吗?
class Simplex(object):
def __init__(self,coordinates):
if not len(coordinates) == 4:
raise RuntimeError('You must provide only 4 coordinates!')