typedef struct lista{
创新互联建站-专业网站定制、快速模板网站建设、高性价比海盐网站开发、企业建站全套包干低至880元,成熟完善的模板库,直接使用。一站式海盐网站制作公司更省心,省钱,快速模板网站建设找我们，业务覆盖海盐地区。费用合理售后完善，十载实体公司更值得信赖。
struct lista *next;
int data;
}list;
void insert(list *h);
void del(list *h);
int main()
{
int flag;
list *head=(list *)malloc(sizeof(list));
head-next=NULL;
while(1)
{
用蒙特卡洛方法计算定积分
计算定积分
利用蒙特卡洛计算方法，核心步骤是求取随机的 g(X1),………,g(Xn),n∈[a,b],由数学期望和大数定理可以近似计算定积分，公式为
原函数:
导函数:
计算导函数在[10,15]上的定积分；
Python
用蒙特卡洛方法计算的定积分：
直接用原函数计算的定积分：
偏差程度为：
有一个相应的公式的话，你只要把找到相应的公式，然后套进去就可以算积分。
在python中,有一个math module,你可以import math,
里面有math.sin(),math.cos(),math.asin()和math.acos()四个函数.
有了这四个函数你就可以求函数值和角度了.
注意:括号里面填的数值,要用弧度制.
科学计算应该要安装numpy，你可以去下Anaconda，全套自动打包好了的。
把代码改了一下：
# -*- coding: utf-8 -*-
import math
import sys
import time
import numpy as np
def area():
x1=1 #起始区间
x2=4.5 #结束区间
dx=0.0001 #步长，步长越小越精确
y1=0
y2=0
y3=0 #结果
for x in np.arange(x1, x2, dx):
y1=y1+dx*x*x #细小矩形取左边为高
y2=y2+ dx*(x+0.5*dx)*(x+0.5*dx) #取中间
y3=y3+dx*(x+dx)*(x+dx)
print (y1, y2, y3)
start = time.clock()
area()
end = time.clock()
print (end-start)
你好，目前python有一个sympy的包。他提供一个积分计算的函数：
integrate(函数，（变量，下限， 上限）)
所以你上面的那个积分计算可以如下：
from sympy import *
x = symbols('x')
print(integrate(abs（x-a）, (x, b1, b2)))
