怎樣用c語言求函數(shù)的一個值

void fun(double *in, double *out)

成都創(chuàng)新互聯(lián)公司云計算的互聯(lián)網(wǎng)服務提供商,擁有超過13年的服務器租用、成都棕樹機房、云服務器、網(wǎng)頁空間、網(wǎng)站系統(tǒng)開發(fā)經(jīng)驗，已先后獲得國家工業(yè)和信息化部頒發(fā)的互聯(lián)網(wǎng)數(shù)據(jù)中心業(yè)務許可證。專業(yè)提供云主機、網(wǎng)頁空間、主機域名、VPS主機、云服務器、香港云服務器、免備案服務器等。

{

while (1)

{

*out = (*in);

if (50 == *in)

{

break;

}

*in = *out;

}

}

如何用c語言求函數(shù)導數(shù)

1、首先要有函數(shù)，設置成double類型的參數(shù)和返回值。

2、然后根據(jù)導數(shù)的定義求出導數(shù)，參數(shù)差值要達到精度極限，這是最關鍵的一步。

3、假如函數(shù)是double fun（doube x），那么導數(shù)的輸出應該是（fun（x）-fun（x-e））/e，這里e是設置的無窮小的變量。

4、C由于精度有限，因此需要循環(huán)反復測試，并判斷無窮小e等于0之前，求出上述導數(shù)的值。二級導數(shù)也是一樣，所不同的是要把上述導數(shù)公式按定義再一次求導。這是算法，具體的實現(xiàn)自己嘗試編程。

C語言的數(shù)據(jù)長度和精度都有限，因此用C語言編程求的導數(shù)并不精確，換句話說C語言編程不適合求導和極限。

擴展資料：

舉例說明：

一階導數(shù)，寫一個函數(shù) y = f(x)：

float f(float x){ ...}

設 dx 初值

計算 dy

dy = f(x0) - f(x0+dx);

導數(shù) 初值

dd1=dy/dx;

Lab:；

dx = 0.5 * dx; ?// 減小步長

dy = f(x0) - f(x0+dx);

dd2=dy/dx; ?// 導數(shù) 新值

判斷新舊導數(shù)值之差是否滿足精度，滿足則得結(jié)果，不滿足則返回

if ( ?fabs(dd1-dd2) 1e-06 ) { 得結(jié)果dd2...}

else { dd1=dd2;goto Lab;}。

已知函數(shù)值C語言中用三次樣條插值求某個點值程序

void SPL(int n, double *x, double *y, int ni, double *xi, double *yi)； 是你所要。

已知 n 個點 x,y; x 必須已按順序排好。要插值 ni 點，橫坐標 xi[], 輸出 yi[]。

程序里用double 型，保證計算精度。

SPL調(diào)用現(xiàn)成的程序。

現(xiàn)成的程序很多。端點處理方法不同，結(jié)果會有不同。想同matlab比較，你需 嘗試 調(diào)用 spline（）函數(shù) 時，令 end1 為 1, 設 slope1 的值，令 end2 為 1 設 slope2 的值。

#include

#include

int spline (int n, int end1, int end2,

double slope1, double slope2,

double x[], double y[],

double b[], double c[], double d[],

int *iflag)

{

int nm1, ib, i, ascend;

double t;

nm1 = n - 1;

*iflag = 0;

if (n 2)

{ /* no possible interpolation */

*iflag = 1;

goto LeaveSpline;

}

ascend = 1;

for (i = 1; i n; ++i) if (x[i] = x[i-1]) ascend = 0;

if (!ascend)

{

*iflag = 2;

goto LeaveSpline;

}

if (n = 3)

{

d[0] = x[1] - x[0];

c[1] = (y[1] - y[0]) / d[0];

for (i = 1; i nm1; ++i)

{

d[i] = x[i+1] - x[i];

b[i] = 2.0 * (d[i-1] + d[i]);

c[i+1] = (y[i+1] - y[i]) / d[i];

c[i] = c[i+1] - c[i];

}

/* ---- Default End conditions */

b[0] = -d[0];

b[nm1] = -d[n-2];

c[0] = 0.0;

c[nm1] = 0.0;

if (n != 3)

{

c[0] = c[2] / (x[3] - x[1]) - c[1] / (x[2] - x[0]);

c[nm1] = c[n-2] / (x[nm1] - x[n-3]) - c[n-3] / (x[n-2] - x[n-4]);

c[0] = c[0] * d[0] * d[0] / (x[3] - x[0]);

c[nm1] = -c[nm1] * d[n-2] * d[n-2] / (x[nm1] - x[n-4]);

}

/* Alternative end conditions -- known slopes */

if (end1 == 1)

{

b[0] = 2.0 * (x[1] - x[0]);

c[0] = (y[1] - y[0]) / (x[1] - x[0]) - slope1;

}

if (end2 == 1)

{

b[nm1] = 2.0 * (x[nm1] - x[n-2]);

c[nm1] = slope2 - (y[nm1] - y[n-2]) / (x[nm1] - x[n-2]);

}

/* Forward elimination */

for (i = 1; i n; ++i)

{

t = d[i-1] / b[i-1];

b[i] = b[i] - t * d[i-1];

c[i] = c[i] - t * c[i-1];

}

/* Back substitution */

c[nm1] = c[nm1] / b[nm1];

for (ib = 0; ib nm1; ++ib)

{

i = n - ib - 2;

c[i] = (c[i] - d[i] * c[i+1]) / b[i];

}

b[nm1] = (y[nm1] - y[n-2]) / d[n-2] + d[n-2] * (c[n-2] + 2.0 * c[nm1]);

for (i = 0; i nm1; ++i)

{

b[i] = (y[i+1] - y[i]) / d[i] - d[i] * (c[i+1] + 2.0 * c[i]);

d[i] = (c[i+1] - c[i]) / d[i];

c[i] = 3.0 * c[i];

}

c[nm1] = 3.0 * c[nm1];

d[nm1] = d[n-2];

}

else

{

b[0] = (y[1] - y[0]) / (x[1] - x[0]);

c[0] = 0.0;

d[0] = 0.0;

b[1] = b[0];

c[1] = 0.0;

d[1] = 0.0;

}

LeaveSpline:

return 0;

}

double seval (int n, double u,

double x[], double y[],

double b[], double c[], double d[],

int *last)

{

int i, j, k;

double w;

i = *last;

if (i = n-1) i = 0;

if (i 0) i = 0;

if ((x[i] u) || (x[i+1] u))

{

i = 0;

j = n;

do

{

k = (i + j) / 2;

if (u x[k]) j = k;

if (u = x[k]) i = k;

}

while (j i+1);

}

*last = i;

w = u - x[i];

w = y[i] + w * (b[i] + w * (c[i] + w * d[i]));

return (w);

}

void SPL(int n, double *x, double *y, int ni, double *xi, double *yi)

{

double *b, *c, *d;

int iflag,last,i;

b = (double *) malloc(sizeof(double) * n);

c = (double *)malloc(sizeof(double) * n);

d = (double *)malloc(sizeof(double) * n);

if (!d) { printf("no enough memory for b,c,d\n");}

else {

spline (n,0,0,0,0,x,y,b,c,d,iflag);

if (iflag==0) printf("I got coef b,c,d now\n"); else printf("x not in order or other error\n");

for (i=0;ini;i++) yi[i] = seval(ni,xi[i],x,y,b,c,d,last);

free(b);free(c);free(d);

};

}

main(){

double x[6]={0.,1.,2.,3.,4.,5};

double y[6]={0.,0.5,2.0,1.6,0.5,0.0};

double u[8]={0.5,1,1.5,2,2.5,3,3.5,4};

double s[8];

int i;

SPL(6, x,y, 8, u, s);

for (i=0;i8;i++) printf("%lf %lf \n",u[i],s[i]);

return 0;

}

新聞標題：c語言求個點函數(shù)值 c語言編寫函數(shù)求表達式的值
本文鏈接：http://aaarwkj.com/article14/dochode.html

成都網(wǎng)站建設公司_創(chuàng)新互聯(lián)，為您提供軟件開發(fā)、網(wǎng)站內(nèi)鏈、網(wǎng)站建設、微信小程序、網(wǎng)頁設計公司、ChatGPT

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時需注明來源： 創(chuàng)新互聯(lián)