1、首先要有函数，设置成double类型的参数和返回值。
创新互联拥有网站维护技术和项目管理团队，建立的售前、实施和售后服务体系，为客户提供定制化的成都网站制作、做网站、网站维护、成都二枢机房解决方案。为客户网站安全和日常运维提供整体管家式外包优质服务。我们的网站维护服务覆盖集团企业、上市公司、外企网站、成都商城网站开发、政府网站等各类型客户群体，为全球上千家企业提供全方位网站维护、服务器维护解决方案。
2、然后根据导数的定义求出导数，参数差值要达到精度极限，这是最关键的一步。
3、假如函数是double fun（doube x），那么导数的输出应该是（fun（x）-fun（x-e））/e，这里e是设置的无穷小的变量。
4、C由于精度有限，因此需要循环反复测试，并判断无穷小e等于0之前，求出上述导数的值。二级导数也是一样，所不同的是要把上述导数公式按定义再一次求导。这是算法，具体的实现自己尝试编程。
C语言的数据长度和精度都有限，因此用C语言编程求的导数并不精确，换句话说C语言编程不适合求导和极限。
扩展资料：
举例说明：
一阶导数，写一个函数 y = f(x)：
float f(float x){ ...}
设 dx 初值
计算 dy
dy = f(x0) - f(x0+dx);
导数 初值
dd1=dy/dx;
Lab:；
dx = 0.5 * dx; // 减小步长
dy = f(x0) - f(x0+dx);
dd2=dy/dx; // 导数 新值
判断新旧导数值之差是否满足精度，满足则得结果，不满足则返回
if ( fabs(dd1-dd2) 1e-06 ) { 得结果dd2...}
else { dd1=dd2;goto Lab;}。
求导数有两种，一种是表达式求导，一种是数值求导。
表达式求导：需要对表达式进行词法分析，然后用常见的求导公式进行演算，求得导函数。在这方面，数学软件matrix，maple做得非常好。如果自己用C进行编程，不建议。
数值求导：利用导数的定义，用差分计算，当自变量趋于0时，前后两次差分收敛到需要精度，计算结束。这种方法可以求得某一点的导数。
例如：
求一阶导数，原函数 y = f(x)， 程序中是float f(float x){ ...}
dx=0.01; //设 dx 初值
do{
dd1=(f(x0) - f(x0+dx))/dx; //计算导数dd1
dx = 0.5 * dx; // 减小步长
dd2=(f(x0) - f(x0+dx))/dx; //计算导数dd2
}while (fabs(dd1-dd2) = 1e-06) //判断新旧导数值之差是否满足精度，满足则得结果，不满足则返回
这是辛普森积分法。
给你写了fun_1( ），fun_2(），请自己添加另外几个被积函数。
调用方法 t=fsimp(a,b,eps,fun_i);
a,b --上下限，eps -- 迭代精度要求。
#includestdio.h
#includestdlib.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();
}
对于一重定积分来说其求解可以使用梯形法进行求解，计算公式如下所示：
其中，f(x)为被积函数，为横坐标的两点间的间隔，越小，则计算出的结果越精确。
对于求解此类问题可以使用C语言中的回调函数编写通用的计算函数，代码如下：
#include stdio.h
#include stdlib.h
#includemath.h
//功能：返回f(x)在积分区间[a,b]的值
//参数：FunCallBack 指向用于计算f(x)的函数
// a 积分区间的起始值
// b 积分区间的结束值
// dx 横坐标的间隔数，越小计算结果越准确
double Calculate(double (*FunCallBack)(double x),
double a,double b,double dx)
{
double doui;
double total = 0; //保存最后的计算结果
for (doui = a; doui = b; doui += dx)
{
total += FunCallBack(doui)*dx;
}
return total;
}
double f2(double x)
{
return x*x;
}
double f(double x)
{
return x;
}
double f3(double x)
{
return x*x*x ;
}
int main()
{
double total;
total = (Calculate(f, 2, 3, 0.000001));
printf("total = %lf\n", total);
total = (Calculate(f2, 2, 3, 0.000001));
printf("total = %lf\n", total);
total = (Calculate(f3, 2, 3, 0.000001));
printf("total = %lf\n", total);
return 0 ;
}
其中，函数f，f2，f3为自行编写的关于x的被积函数。
运行结果：
total = 2.500000
total = 6.333331
total = 16.249991