/*
*
* 用变步长梯形求积公式求定积分
* 1
* ∫ (x/(4+x^2))dx
* 0
*/
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <conio.h>
#include <stdlib.h>
#define epsilon 0.00001
float f(float x)
{
return(x/(4+x*x));
}
float computeT(float a,float b)
{
float T=0,h=(b-a)/2;
T=h*(f(a)+2*T+f(b))/2;
return(T);
}
main()
{
int i;
float a,b,T1,T2,S,x,h;
printf("\n Input the begin :");
scanf("%f",&a);
printf("\n Input the end :");
scanf("%f",&b);
h=b-a;
T1=computeT(a,b);
while(1)
{
S=0;
x=a+h/2;
do
{
S+=f(x);
x+=h;
}while(x<b);
T2=(T1+h*S)/2;
if(fabs(T2-T1)>=epsilon)
{
T1=T2;
h=h/2;
continue;
}
else
break;
}
printf("solve is : %f",T2);
getch();
return(S);
}
分享到:
相关推荐
变步长复合梯形公式算法公式 变步长复合梯形公式算法公式 变步长复合梯形公式算法公式
龙贝格 变步长 梯形求积法求积分的程序代码,有简单程序说明
变步长梯形求积分法,图解及代码表示 利用梯形公式计算积分值
分别写出变步长梯形法与 Romberge法计算定积分的算法,编写程序上机调试出结果, 要求所编程序适用于任何类型的定积分,即能解决这一类问题,而不是某一个问题
这样会将原来的求积分问题带到另一个求积分问题:如何求n次多项式的积分,而且当次数变高时,会出现龙悲歌现象,误差反而可能会增大,并且高次的插值求积公式有可能会变得不稳定:详细原因不赘述。 牛顿-科特斯公式...
数值分析与算法_变步长梯形求积法计算定积分.doc
基于MATLAB的变步长梯形数值积分法的研究与实验.pdf
变步长梯形法求积分的C语言实验 复化梯形法或者叫复合梯形法 代码摘要 double Multy(double a, double b) { double h = b - a, T[2]; T[0] = h * (f(b) + f(a)) / 2.0; // 求第一个Tn,用T[0]存储 int m ...
用变步长辛普生求积公式和龙贝格积分法计算习题,完整报告
根据计算方法中阐述的求积分方法编程:变步长梯形求积分 变步长抛物线求积分 复合抛物线求积分
[计算方法作业]利用python中matplotlib实现绘制龙贝格公式和变步长梯形法图像,利用龙贝格公式及变步长积分两种算法计算定积分,用python中的matplotlib库绘制图像
1. 用不同数值方法计算积分 (1) 取不同的步长h. 分别用复合梯形及复合辛普森求积计算积分...附录1 复合梯形求积MATLAB程序 附录2复合辛普森求积MATLAB程序 附录3龙贝格求积MATLAB程序 附录4 自适应辛普森求积MATLAB程序
变步长梯形公式 复化辛普森公式 龙贝格公式 复化梯形公式
定步长复合梯形公式算法公式定步长复合梯形公式算法公式定步长复合梯形公式算法公式
变步长复合梯形公式代.rar变步长复合梯形公式代.rar变步长复合梯形公式代.rar变步长复合梯形公式代.rar变步长复合梯形公式代.rar变步长复合梯形公式代.rar
数值计算 变步长梯形公式-实验6 详细过程加代码
MATLAB数值分析实验二(复合梯形、辛普森和龙贝格求积-以及二重积分计算等)
验证定步长复合梯形公式、变步长复合梯形公式;进一步加深对数值积分的理解。 1、编写定步长复合梯形程序, 求 的近似值。
包含代码和文档 采用复化梯形公式和复化辛普森公式求积分,并与精确值进行比较得下表。 采用复化梯形公式和复化辛普森公式求积分,并与精确值进行比较得下表。
高斯求积公式的PPT课件。计算方法。 熟悉复合梯形公式、复合抛物线公式及其余项; 熟悉Newton-Cotes公式; 熟悉代数精度法构造求积公式的思想;...了解变步长梯形公式和Romberg算法; 了解Guass求积公式的特征。