[交流] 【求助】建立体心晶胞结构代码，运行不出结果！ 已有2人参与

用VC++6.0编的一个建立体心立方晶胞结构的代码，编译成功了，但是运行不出结果，找不到错误所在，请高手帮忙一下啊！十分感谢！ void CBCC_CrystalDlg::OnCalculate() {         // TODO: Add your control notification handler code here         FILE * fileFe_xyz = fopen("Fe_xyz.txt", "w+";         fprintf(fileFe_xyz, "i  xbas  ybas  zbas\n";                 a=2.8664;  //Fe的晶格常数         xunit[0]=0.0;  //以下为单胞原子坐标         yunit[0]=0.0;         zunit[0]=0.0;         xunit[1]=0.50;         yunit[1]=0.50;         zunit[1]=0.50;           ncellx=M;   //ncellx,ncelly,ncellz是在xyz方向的重复次数，其中M任意整数        ncelly=M;         ncellz=M;         iatoms=0;         natoms=2*ncellx*ncelly*ncellz;         for(int i=0;i         {            xbas=0;                ybas=0;            zbas=0;         }         for(ix=1;i<=ncellx;ix++)         {            for(iy=1;iy<=ncelly;iy++)            {               for(iz=1;iz<=ncellz;iz++)               {                  for(ibas=0;ibas<2;ibas++)                  {                 xbas[iatoms]=(float(ix-1)+xunit[ibas])*a;                 ybas[iatoms]=(float(iy-1)+yunit[ibas])*a;                 zbas[iatoms]=(float(iz-1)+zunit[ibas])*a;                 iatoms=iatoms+1;                  }               }            }         }                 for(i=0;i         { fprintf(fileFe_xyz, "%6d    %10.5f    %10.5f    %10.5f\n", i,xbas,ybas,zbas);         }         fclose(fileFe_xyz);         } 这是头文件中定义的变量及函数： #define M 5 public:         int i;         int ix;         int iy;         int iz;         int ncellx;         int ncelly;         int ncellz;         int iatoms;         int natoms;         int ibas;         double xunit[2];         double yunit[2];         double zunit[2];         double a;         double xbas[4000];         double ybas[4000];         double zbas[4000];         void OnCalculate(); [ Last edited by _yixiao on 2010-6-3 at 20:01 ]

回复此楼