1、首先要有函数，设置成double类型的参数和返回值。
成都创新互联电话联系：18980820575，为您提供成都网站建设网页设计及定制高端网站建设服务，成都创新互联网页制作领域10年，包括成都人造雾等多个行业拥有多年的网站推广经验，选择成都创新互联，为网站锦上添花！
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;}。
用差分计算，当自变量趋于0时，前后两次差分收敛到需要精度，计算结束。
例如，一阶导数，写一个函数y=f(x)：
floatf(floatx){...}
设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;gotoLab;};
//多项式求导数
intPolyDeri(listnodePolypolyFunc)
{
listnodePoly::iteratoriter;
for(iter=polyFunc.begin();iter!=polyFunc.end();++iter)
{
if((*iter).ex1)
{
(*iter).coef=((*iter).coef)*((*iter).ex);
(*iter).ex=(*iter).ex-1;
}
elseif(1==(*iter).ex)
{
(*iter).ex=0;
}
elseif(0==(*iter).ex)
{
(*iter).coef=0;
}
}
returnRET_OK;
}
其中，多项式的定义是listnodePoly，如下：
//多项式节点结构体定义
typedefstructstuPolynomNode
{
doublecoef;
intex;
}nodePoly;
扩展资料
c语言求导数据范围及提示DataSizeHint
#includeiostream
#includecmath
usingnamespacestd;
intmain()
{
intnum=0,i=0;
cinnum;
for(i=2;i=sqrt(num);i++)
{
if(num%i==0)
break;
}
if(isqrt(num)
coutnum"为素数"endl;
else
coutnum"不是素数"endl;
return0;
}
求导数有两种，一种是表达式求导，一种是数值求导。
1.
表达式求导：需要对表达式进行词法分析，然后用常见的求导公式进行演算，求得导函数。在这方面，数学软件matrix，maple做得非常好。如果自己用C进行编程，不建议。
2.
数值求导：利用导数的定义，用差分计算，当自变量趋于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) //判断新旧导数值之差是否满足精度，满足则得结果，不满足则返回