C语言如何求导函数

用计算，当趋于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;};