‘壹’ c语言编程如何编写积分公式
#include<iostream.h>
//定义结构类型
structstudent
{
intnum;
charname[20];
floatgrade;
};
voidmain(void)
{
//声明数组
inti,size;
charstr[]="Thisisastring.";
intint_values[]={51,23,2,44,45,0,11};
floatfloat_values[]={15.1,13.3,22.2,10.4,1.5};
studentst_arr[]={101,"WangLin",92,102,"LiPing",85,103,"ZhaoMin",88};
//显示char类型数组元素及其大小
size=sizeof(str)/sizeof(char);
cout<<"Numberofelementsinstr:";
cout<<size<<endl;
for(i=0;i<size;i++){
cout<<str[i];
}
cout<<endl;
//显示int类型数组元素及其大小
size=sizeof(int_values)/sizeof(int);
cout<<"Numberofelementsinint_values:";
cout<<size<<endl;
for(i=0;i<size;i++){
cout<<int_values[i]<<"";
}
cout<<endl;
//显示float类型数组元素及其大小
size=sizeof(float_values)/sizeof(float);
cout<<"Numberofelementsinfloat_values:";
cout<<size<<endl;
for(i=0;i<size;i++){
cout<<float_values[i]<<"";
}
cout<<endl;
//显示student类型数组元素及其大小
size=sizeof(st_arr)/sizeof(student);
cout<<"Numberofelementsinst_arr:";
cout<<size<<endl;
for(i=0;i<size;i++){
cout<<st_arr[i].num<<"";
cout<<st_arr[i].name<<"";
cout<<st_arr[i].grade<<endl;
}
}
#include<iostream.h>
//add()函数的定义，其有返回值
doubleadd(doublex,doubley)
{
doublez;
z=x+y;
cout<<x<<"+"<<y<<"="<<z<<endl;
return(z);
}
main()
{
doublea=0.5,b=1.0;
//以不同参数形式调用函数add()
cout<<"add(1.5,2.5)="<<add(1.5,2.5)<<endl;
cout<<"add(a,b)="<<add(a,b)<<endl;
cout<<"add(2*a,a+b)="<<add(2*a,a+b)<<endl;
cout<<"----------------------"<<endl;
//以表达式方式调用函数add()
doublec=2*add(a,b);
cout<<"c="<<c<<endl;
cout<<"----------------------"<<endl;
//以语句式方式调用函数add()
add(2*a,b);
cout<<"----------------------"<<endl;
//用其他类型参数调用函数add()
intn=1,m=2;
cout<<"add("<<n<<","<<m<<")="<<add(n,m)<<endl;
}
#include<iostream.h>
//定义符号函数sgn(),其返回值为int类型
intsgn(doublex)
{
if(x>0)return(1);//返回出口1
if(x<0)return(-1);//返回出口2
return(0);//返回出口3
}
//main()函数定义
main()
{
doublex;
inti;
for(i=0;i<=2;i++){
cout<<"x=";
cin>>x;
cout<<"sgn("<<x<<")="<<sgn(x)<<endl;
}
}
‘贰’ 如何用C语言求二重定积分
按二重积分的定义做。
我的这个出错，一起讨论下吧：
#include<stdio.h>
#include<math.h>
#defineN100
voidmain(){
doubleresult,(*p)(double,double);
doublef(double,double);
doubleintegral(double,double,double,double,double(*p)());
p=f;
// printf("请输入abcd的值： ");
// scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
result=integral(2.0,8.0,0,6.0,p);
printf("该二重积分的值为： ");
printf("%f",result);
printf(" ");
}
//求二重积分的函数
doubleintegral(doublea,doubleb,doublec,doubled,double(*p)(double,double)){
inti,j;
doublex,y,h1,h2,volum=0;
h1=(b-a)/N;
h2=(d-c)/N;
x=a;
y=c;
for(i=0;i<N;i++){
for(j=0;j<N;j++){
volum+=(*p)(x,y)*h1*h2;
y+=h2;
}
x+=h1;
}
return(volum);
}
doublef(doublex,doubley){
doubleh;
h=sqrt((2+x*x)*sin(x*y)*sin(x*y))+x*x*y*cos(x)-30*x*y*y;
return(h);
}
运行结果：
‘叁’ 用C语言编写一个求定积分的程序
这是辛普森积分法。
给你写了fun_1( ），fun_2(），请自己添加另外几个被积函数。
调用方法 t=fsimp(a,b,eps,fun_i);
a,b --上下限，eps -- 迭代精度要求。
#include<stdio.h>
#include<stdlib.h>
#include <math.h>
double fun_1(double x)
{
return 1.0 + x ;
}
double fun_2(double x)
{
return 2.0 * x + 3.0 ;
}
double fsimp(double a,double b,double eps, double (*P)(double))
{
int n,k;
double h,t1,t2,s1,s2,ep,p,x;
n=1; h=b-a;
t1=h*(P(a)+P(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+P(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,fun_1);
printf("%g\n",t);
t=fsimp(a,b,eps,fun_2);
printf("%g\n",t);
// ...
printf("\n Press any key to quit...");
getch();
}
‘肆’ C语言编写f(x)=x^2在区间[0,1]定积分
思路就是将区间无限细分
小区间梯形就可以看成矩形，小矩形所有面积相加，就是定积分
#include"stdio.h"
#include"math.h"
doublef(doublex)//被积函数
{
returnx*x;
}
doublefun(doublep,doubleq,intn)//p,q分别表示积分下上限，n为细分的区间数，越大值越精确
{
inti;
doublesum=0,h=(q-p)/n;
for(i=1;i<=n;i++)//f(p+i*h)*h就是小矩形的面积，加到sum上去。
sum+=f(p+i*h)*h;
returnsum;
}
intmain()
{
printf("(x*x)在(0,1)上定积分为%lf ",fun(0,1,1000000));//区间数自己设越大值越精确
}
结果：
数学分析：
f(x)=x^2=x*x;
定积分：x*x*x/3+c(常数)
在区间(0,1)上定积分：1/3=0.333333
结果正确。
‘伍’ 用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
由上面的比较可知，梯形法的精度要高于矩形法。
‘陆’ 怎样编写c语言积分函数
积分分为两种，数值积分，公式积分。
公式积分：部分函数可以直接用公式求得其不定积分函数。C语言中可以直接用积分公式写出其积分函数。
数值积分：按照积分的定义，设置积分范围的步长，用梯形面积累加求得其积分。
以【f(x)=x*sin(x) 从1到2的积分】为例：
#include<math.h>
#include<stdio.h>
doubleintegral(double(*fun)(doublex),doublea,doubleb,int,n){
doubles,h,y;
inti;
s=(fun(a)+fun(b))/2;
h=(b-a)/n;/*积分步长*/
for(i=1;i<n;i++)
s=s+fun(a+i*h);
y=s*h;
returny;/*返回积分值*/
}
doublef(doublex){
return(x*sinx)/*修改此处可以改变被积函数*/
}
intmain(){
doubley;
y=integral(f,1.0,2.0,150);/*修改此处可以改变积分上下限和步数,步长=(上限-下限)/步数*/
printf("y=%f ",y);
return0;
}
‘柒’ 怎么用C语言表示积分
源代码如下：
#include#includefloat f1(float x)
{
return(1.0+x);
}
float f2(float x)
{
return(2.0*x+3.0);
}
float f3(float x)
{
return(exp(x)+1);
}
float f4(float x)
{
return(pow(1+x,2));
}
float f5(float x)
{
return(pow(x,3));
}
float fsimp(float a,float b,float (*p)(float))
{
float c,s;
c=(a+b)/2;
s=(b-a)/6*(p(a)+4*p(c)+p(b));
return s;
}
int main()
{
float a,b;
printf("请输入积分下限a的值：");
scanf("%f",&a);
printf("请输入积分上限b的值：");
scanf("%f",&b);
printf("%f ",fsimp(a,b,f1));
}
(7)c语言求yx2积分扩展阅读
1、对应于一个积分式要有一段程序，可以改变程序的一小部分来改变所要求的积分式。
2、除数不能位0。
3、两个整数相除，结果仍是整数。
4、若被除数其中有一个为浮点数或者两个都为浮点数，则结果为浮点类型。操作数必须为整数，不能是浮点数。
‘捌’ 用C语言求积分
好象要你自己定义函数，库函数中好象是没有直接可以用的积分函数的。
‘玖’ c语言求定积分
问题就是出在数据类型上的选用上，precision=0.0000001时已经超过了float的数据范围，所以导致数据截断后precision=0.000000，从而程序在计算积分时可能陷入死循环，应该采用double型数据类型。其实不推荐楼主用如此多的define语句，程序的可读性和风格应该重于编程员的劳动度。。。
还有楼主对自然对数e的define也已经超过了计算机的可识别范围。。您那样精确的定义e并不会在结果上获得更加精确地结果，其实反倒会起到相反的作用，要知道与其用一个这样可能导致内存出错以及必定会导致数据截断的变量来实现精度的提高远远不如采用一个更精确的积分算法，而且c语言提供了自然数e为底的指数函数~而且貌似您的积分算法是不准确的，梯形积分的定义并非如此，其再两端的函数值应该只取1/2.希望您多加细心~
如果不介意的话，就是你的precision应该改为step~这样会能更加准备的表达了这个变量的作用，在你的程序中precision变量其实是积分步长~在数值计算方法中积分精度的控制往往不是通过细化步长来表达，而是通过后一个积分值-前一个积分值<precision
这样来实现精度控制~呵呵
‘拾’ 用c语言求二重积分~~
这个你用手工积分就行了，当然常数的代入计算还是可以用计算器的。给你的提示是分0<y<a/2和a/2<y<a两次脱绝对值符号，最后都会演变为对sqrt(x^2+u^2)/(x^2+u^2+v^2)^2的积分形式，或者，你只关心数值结果的话，直接采用矩形法或者辛普森积分法，划分的细一点就不会有太大误差了