?? rofunc.c
字號:
#include <math.h>#define NRANSI#include "nrutil.h"#define EPS 1.0e-7extern int ndatat;extern float *xt,*yt,aa,abdevt;float rofunc(float b){ float select(unsigned long k, unsigned long n, float arr[]); int j; float *arr,d,sum=0.0; arr=vector(1,ndatat); for (j=1;j<=ndatat;j++) arr[j]=yt[j]-b*xt[j]; if (ndatat & 1) { aa=select((ndatat+1)>>1,ndatat,arr); } else { j=ndatat >> 1; aa=0.5*(select(j,ndatat,arr)+select(j+1,ndatat,arr)); } abdevt=0.0; for (j=1;j<=ndatat;j++) { d=yt[j]-(b*xt[j]+aa); abdevt += fabs(d); if (yt[j] != 0.0) d /= fabs(yt[j]); if (fabs(d) > EPS) sum += (d >= 0.0 ? xt[j] : -xt[j]); } free_vector(arr,1,ndatat); return sum;}#undef EPS#undef NRANSI
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -