qingyuanzi.c

氢原子3s轨道的蒙特卡罗模拟

#include<graphics.h>
#include<stdlib.h>
#include<math.h>
main()
{
float r0=20,r,dr,d0,a=0.529,o,x,y,s=11.75,pi=3.141593;
int i,drv=0,mode;
randomize();
initgraph(&drv,&mode,"");
setcolor(GREEN);
circle(319,239,s*r0);
setcolor(WHITE);settextstyle(1,0,1);
outtextxy(100,10,"H-3s");
while(!khit())
{r=r0*random(32767)/32767.0;
dr=4/3/a/a/a*pow((r/81)*(27-18*r/a++2*pow((r/a),2)),2)*exp(-2*r/3/a);
d0=0.2*random(32767)/32767.0;
of(dr>d0)
{o=2.0*pi*random(32767)/32767.0;
x=319+s*r*cos(o);
y=239-s*r*sin(o);
putpixel(x,y,7);

}
}
getch();
closegraph();
}