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=dd2goto Lab}。
//多项式求导数
intPolyDeri(list<nodePoly>&polyFunc)
{
list<nodePoly>::iteratoriter
for(iter=polyFunc.begin()iter!=polyFunc.end()++iter)
{
if((*iter).ex>1)
{
(*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
}
其中,多项式的定义是list<nodePoly>,如下:
//多项式节点结构体定义
typedefstructstuPolynomNode
{
doublecoef
intex
}nodePoly
扩展资料
c语言求导数据范围及提示DataSize&Hint
#include<iostream>
#include<cmath>
usingnamespacestd
intmain()
{
intnum=0,i=0
cin>>num
for(i=2i<=sqrt(num)i++)
{
if(num%i==0)
break
}
if(i>sqrt(num)
cout<<num<<"为素数"<<endl
else
cout<<num<<"不是素数"endl
return0
}
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)