如何用c语言求函数导数

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}。

这是辛普森积分法。

给你写了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=1h=b-a

t1=h*(P(a)+P(b))/2.0

s1=t1

ep=eps+1.0

while (ep>=eps)

{

p=0.0

for (k=0k<=n-1k++)

{

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=t2s1=s2n=n+nh=h/2.0

}

return(s2)

}

void main()

{

double a,b,eps,t

a=0.0b=3.141592653589793238eps=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()

}

---