氢原子电子云模拟显示
模拟各种氢原子电子云的状态,希望对大家有帮助!!
# include " graphics . h "
# include " math . h "
# include " conio . h "
# include " stdio . h "
# define THETA PI/8
# define PHI PI * 2/5
# define L 800
# define RR 80
# define SQR(X,Y,Z) sqrt(X*X+Y*Y
int x1,y1;
float n1(x,y,z) /*1s*/
float x,y,z;
{
double r;
r=SQR(x,y,z);
return(16.0 * exp(-r));
}
浮点数n2(x,y,z) /*2s*/
浮点数x,y,z;
{
double r;
r=SQR(x,y,z);
return(16 *(2-r)* exp(-r/2));
}
浮点数n3(x,y,z) /*2pz*/
浮点数x,y,z;
{
double r;
r=SQR(x,y,z);
return(16 * z * exp(-r/2));
}
float n4(x,y,z) /*2px*/
float x,y,z;
{
double r;
r=SQR(x,y,z);
return(20 * x * exp(-r/2));
}
浮点n5(x,y,z) /*2py*/
浮点x,y,z;
{
double r;
r=SQR(x,y,z);
return(20 * y * exp(-r/2));
}
浮点n6(x,y,z) /*3s*/
浮点x,y,z;
{
double r;
r=SQR(x,y,z);
return(4 *(27-18 * r+2 * r * r)* exp(-r/3));
}
浮点n7(x,y,z) /*3pz*/
浮点x,y,z;
{
double r;
r=SQR(x,y,z);
return(4 *(6-r)* z * exp(-r/3));
}
0条评论