main.cpp

研究生期间上《数值计算方法》课的作业

#include <iostream>
#include <stdlib.h>
#include<math.h>
#define n 20

using namespace std;

float fx(float x)     //积被函数  f(x) 
{   float f;
    f=exp(x)*cos(x);
    return(f); 
}

int main(int argc, char *argv[])
{		     
    float s,s1,s2,x[2*n],h,a,b;
    int i;
    printf("输入下限：");    scanf("%f",&a);
    printf("输入上限：");    scanf("%f",&b);
    h=(b-a)/(2*n);
    for(i=0;i<=2*n-1;i++)
        x[i]=a+h*i;
    s1=s2=0;
    for(i=1;i<=n;i++)
        {s1=s1+fx(x[2*i-1]);
         s2=s2+fx(x[2*i-2]);}
    s=(fx(a)+fx(b)+4*s1+2*(s2-x[0]))*h/3;
    printf("利用Simpson公式求得的积分的近似值是：%f\n",s);

    system("PAUSE");	
    return 0;
}