?? 7.12 求解實函數或復函數方程一個復根的蒙特卡洛法 cmtc.c
字號:
#include "math.h"
#include "3rnd1.c"
void cmtc(x,y,b,m,eps,f)
int m;
double *x,*y,b,eps,(*f)();
{
int k;
double xx,yy,a,r,z,x1,y1,z1;
a=b; k=1; r=1.0; xx=*x; yy=*y;
z=(*f)(xx,yy);
while (a>=eps)
{
x1=-a+2.0*a*rnd1(&r); x1=xx+x1;
y1=-a+2.0*a*rnd1(&r); y1=yy+y1;
z1=(*f)(x1,y1);
k=k+1;
if (z1>=z)
{
if (k>m)
{
k=1; a=a/2.0;
}
}
else
{
k=1; xx=x1; yy=y1; z=z1;
if (z<eps)
{
*x=xx; *y=yy;
return;
}
}
}
*x=xx; *y=yy;
return;
}
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -