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

　　# include "graphics.h"
　　# include "math.h"
　　# include "conio.h"
　　# include "stdio.h"
　　# define PI 3.14159
　　# 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+Z*Z)

　　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));
　　}

　　float n2(x,y,z) /*2s*/
　　float x,y,z;
　　{
　　double r;
　　r=SQR(x,y,z);
　　return(16*(2-r)*exp(-r/2));
　　}

　　float n3(x,y,z) /*2pz*/
　　float 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));
　　}

　　float n5(x,y,z) /*2py*/
　　float x,y,z;
　　{
　　double r;
　　r=SQR(x,y,z);
　　return(20*y*exp(-r/2));
　　}

　　float n6(x,y,z) /*3s*/
　　float x,y,z;
　　{
　　double r;
　　r=SQR(x,y,z);
　　return(4*(27-18*r+2*r*r)*exp(-r/3));
　　}

　　float n7(x,y,z) /*3pz*/
　　float x,y,z;
　　{
　　double r;
　　r=SQR(x,y,z);
　　return(4*(6-r)*z*exp(-r/3));
　　}