C語言用辛普森公式求sinx在0到π上的定積分的源程序

用梯形法估算，再用辛普森法。

成都創(chuàng)新互聯(lián)公司是一家集網(wǎng)站建設(shè),蒙陰企業(yè)網(wǎng)站建設(shè),蒙陰品牌網(wǎng)站建設(shè),網(wǎng)站定制,蒙陰網(wǎng)站建設(shè)報(bào)價(jià),網(wǎng)絡(luò)營(yíng)銷,網(wǎng)絡(luò)優(yōu)化,蒙陰網(wǎng)站推廣為一體的創(chuàng)新建站企業(yè)，幫助傳統(tǒng)企業(yè)提升企業(yè)形象加強(qiáng)企業(yè)競(jìng)爭(zhēng)力。可充分滿足這一群體相比中小企業(yè)更為豐富、高端、多元的互聯(lián)網(wǎng)需求。同時(shí)我們時(shí)刻保持專業(yè)、時(shí)尚、前沿，時(shí)刻以成就客戶成長(zhǎng)自我，堅(jiān)持不斷學(xué)習(xí)、思考、沉淀、凈化自己，讓我們?yōu)楦嗟钠髽I(yè)打造出實(shí)用型網(wǎng)站。

fsimpf 積分函數(shù)

a,b 積分下上限，eps 精度。

#includestdlib.h

#include math.h

double fsimpf(double x)

{

return sin(x) ;

}

double fsimp(double a,double b,double eps)

{

int n,k;

double h,t1,t2,s1,s2,ep,p,x;

n=1; h=b-a;

t1=h*(fsimpf(a)+fsimpf(b))/2.0;

s1=t1;

ep=eps+1.0;

while (ep=eps)

{

p=0.0;

for (k=0;k=n-1;k++)

{

x=a+(k+0.5)*h;

p=p+fsimpf(x);

}

t2=(t1+h*p)/2.0;

s2=(4.0*t2-t1)/3.0;

ep=fabs(s2-s1);

t1=t2; s1=s2; n=n+n; h=h/2.0;

}

return(s2);

}

void main()

{

double a,b,eps,t;

a=0.0; b=3.141592653589793238; eps=0.0000001;

a definite integral by Simpson Method.

t=fsimp(a,b,eps);

printf("%e\n",t);

printf("\n Press any key to quit...");

getch();

}

用c語言求定積分

#include

#include

double integral(double(*fun)(double x),double a,double b,int,n){

double s,h,y;

int i;

s=(fun(a)+fun(b))/2;

h=(b-a)/n; /*積分步長(zhǎng)*/

for(i=1;in;i++)

s=s+fun(a+i*h);

y=s*h;

return y;/*返回積分值*/

}

double f(double x){

return(x*sinx) /*修改此處可以改變被積函數(shù)*/

}

int main(){

double y;

y=integral(f,1.0,2.0,150);/*修改此處可以改變積分上下限和步數(shù),步長(zhǎng)=(上限-下限)/步數(shù)*/

printf("y=%f\n",y);

return 0;

}

int main()

C語言編程求函數(shù)y=sinx在（0,1）區(qū)間的定積分？

#includestdio.h

#includemath.h

#define N 10000000

void main()

{

double sum=0;

int i;

for(i=0; iN; i++)

{

sum += sin((double)(i)/N)/N;

}

printf("%lf\n%lf",sum,1-cos(1));

}

N后面的0有點(diǎn)多了，不過這個(gè)數(shù)剛好能精確到小數(shù)點(diǎn)后6位。

當(dāng)然，這個(gè)程序的算法實(shí)在是太垃圾了，效率非常的低下……

怎么用c語言編正弦函數(shù)計(jì)算

1、C語言中要編寫sin函數(shù)，實(shí)質(zhì)上要利用sin的泰勒公式，然后根據(jù)泰勒公式，將其中的每一項(xiàng)進(jìn)行分解，最后用循環(huán)，累加計(jì)算出最終結(jié)果。

2、下面用for循環(huán)實(shí)現(xiàn)sin的算法，程序代碼如下：

#includestdio.h

#includemath.h

void?main()

{

int??i;

float??x,sum,a,b;??//sum代表和，a為分子，b為分母

char?s;

printf("please?input?x");

scanf("%f",x);

s=1;

sum=0;

a=x;?????//分母賦初值

b=1;?????//分子賦初值

for（i=1;a/b=1e-6;i++）

{

sum=sum+s*a/b;????//累加一項(xiàng)

a=a*x*x;?????//求下一項(xiàng)分子

b=b*2*i*(2*i+1);???//求下一項(xiàng)分母

s*=-1;

}

printf("sum=%f\n",sum);

}

3、? 關(guān)于上述程序的幾點(diǎn)說明：上述程序的計(jì)算結(jié)果精確到小數(shù)點(diǎn)后六位；上述程序運(yùn)用了sin的泰勒展開式 sin x=x-x^3/3!+x^5/5! ...... ，程序中將sin泰勒公式中的每一項(xiàng)拆成了分子，分母以及每一項(xiàng)前的符號(hào)這三項(xiàng)，以便于每一項(xiàng)的累加。