氢原子电子云模拟显示

模拟各种氢原子电子云的状态，希望对大家有帮助！！

# 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))；
}