❶ c語言求定積分的問題。
根據梯形法求積分的原理,設間隔h= (b-a)/n,則積分近似計算公式為:
s = h/2 *[f(a)+f(a+h)] +h/2 *[f(a+h)+f(a+2h)] +...+h/2 *[f(b-h)+f(b)]
=h/2 *[f(a)+f(b)] + h*[ f(a+h) + f(a+2h) +f(a+3h) + ... +f(b -h)]
令積分s初始值為h/2 *[f(a)+f(b)] ,後面令i=1,...,n-1來迭代s = s+h*f(a+ih)。
❷ C語言利用梯形法求定積分,小梯形數量200
你不會微積分的話俺寫了也白寫
=================================================================
/* ∫(0.8-0)cosxdx=sin(0.8)-sin(0)=sin(0.8)
當切割的梯形個數趨於無窮時,答案會等於sin(0.8)
不過計算機位數有限,能切的分數也很有限*/
main()
{
double t=(0.8-0)/200;//每次增量為0.8至0的1/200,這是一小段的長
double x=0.0,s=0.0;//x為坐標圖上的x值,用來求高
short i;
for(i=0;i<200;i++)
{
s+=t*cos(x);//用一小段的長*一小段的高 =一小段面積,一小段面積累加出定積分的值
x+=t;//下一個x值
}
printf("%lf\n",s);
}
❸ C語言:用數值積分的矩形法和梯形法計算(sin x)/x在區間[1,2]上的定積分
#include <stdio.h>
#include <math.h>
double integrate(double (*f)(double),double a,double b)
{
double i=1.0/1000000;
double sum=0,x;
for(x=a;x<=b;)
{
sum+=f(x)*i;
x+=i;
}
return sum;
}
double f(double x)
{
return sin(x)/x;
}
int main()
{
int a=1,b=2;
printf("%lf\n",integrate(f, a, b));
return 0;
}
❹ c語言用梯形法求積分
a區間起點,b積分區間中點, n是這個區間劃分數。每一個劃分的長度為h=(b-a)/n。每一個劃分利用梯形演算法來求積分,即:上底+下底乘以高再除以2。第一個小梯形面積為 (f(a) + f(a+h))*h/2,第二個為 (f(a+h)+f(a+2h))*h/2,依次類推。將這些加和加起來有:
h*(f(a)/2 + f(a+h) + ... +f(a+(n-1)h) + f(b))/2)。預先將f(a)/2 和f(b)/2加起來,其它數據即f(a+h) + ... f(a+(n-1)n)用循環處理,然後將所有結果乘以h即可。
❺ C語言,用梯形法編程求定積分x^3+x/2+1的值
這個其實很簡單,給你做了一個很簡單的通用的積分函數,
已經做過優化,自己一調用就可以了,效率很高。
#include <stdio.h>
double f(double x) //返回函數的值
{ double y;
y=x*x+x*2+1.0;
return y;
}
double Integrate(double (*f)(double),double a,double b,int n) //計算函數的定積分
{ double dx,y,sy,x;
if(n<=1)return 0.0;
dx=(b-a)/n;
x=a;
sy=0.0;
for(int i=0;i<=n;i++,x+=dx)
sy+=f(x);
sy-=f(a)/2.0;
sy-=f(b)/2.0;
return sy*dx;
}
void main()
{ double s;
s=Integrate(f, 1, 3, 100000);
printf("s=%lf\n",s);
}
❻ 有關c語言用梯形法求定積分的一個程序,請幫忙修改
#include<stdio.h>
#include<math.h>
doublef(doublex)
{
returnsqrt(4-x*x);
}
main()
{
doublea,b,s1,s2,s,h,m=1E-6;
intn,i;
printf("請輸入積分下限a和積分上限b ");
scanf("%lf%lf",&a,&b);
s2=0;
n=10000;
// do
{
s1=s2;
n=2*n;
h=(b-a)/n;
i=1;
s2=0;
do
{
s=(f(a+i*h)+f(a+(i-1)*h))*h/2;
s2=s2+s;
i=i+1;
}
while(i<=n);
}
// while((fabs(s1-s2))>(m*(fabs(s2)-fabs(s1))));
printf("f(x)在[a,b]上的定積分為%.10lf",s2);
return0;
}
❼ c語言課設用累計梯形求定積分
#include
<stdio.h>
#include<math.h>
void
main()
{int
i,n;
float
t,v1,v2,s=0,d,s=0;
printf("請輸入分段個數");
scanf("%d",&n);
printf("請輸入右端點的值");
scanf("%f",&t);
v1=cos(0);
v2=cos(t);
d=t/n;
for(i=1;i<n;i++)
s+=cos(i*d)
s=s*d+(v1+v2)/2;
printf("cosx從0到%f的定積分的近似值(梯形法)是:%f",t,s);
}
❽ C語言實慣用梯形法或辛普森法求解定積分的值
//梯形法求定積分
#include<stdio.h>
#include<math.h>
//定義被積函數
double func(double x){
return sin(x)*cos(x);
}
void main(){
double a,b,h,x,sum;
int i,n;
printf("Input a b and n: ");
scanf("%lf%lf%d",&a,&b,&n);
h=(b-a)/n;
x=a;
sum=(func(a)+func(b))/2;
for(i=1; i<n; i++){
x += h;
sum += func(x);
}
sum *= h;
printf("sum=%.4lf\n",sum);
}
❾ 用C語言求定積分
實際問題描述:
求定積分近似值
程序代碼如下:
#include
#include
void main()
{
int i,n=1000;
float a,b,h,t1,t2,s1,s2,x;
printf("請輸入積分限a,b:");
scanf("%f,%f",&a,&b);
h=(b-a)/n;
for(s1=0,s2=0,i=1;i<=n;i++)
{
x=a+(i-1)*h;
t1=(float)exp(-x*x/2);t2(float)=exp(-(x+h)*(x+h)/2);
s1=s1+t1*h; /*矩形面積累加*/
s2=s2+(t1+t2)*h/2; /*梯形面積累加*/
}
printf("矩形法算得積分值:%f.
",s1);
printf("梯形法算得積分值:%f.
",s2);
}
程序運行結果如下:
矩形法算得積分值:0.855821
梯形法算得積分值:0.855624
由上面的比較可知,梯形法的精度要高於矩形法。