#include <stdio.h>
#include <math.h>  
main()
{
int idum, k;
double sigmae,sigmat;
double *mu, *u, *sigm,*sigmt;
sigmae=0.0;
/*sigm=vector(0,10);
sigmt=vector(0,10);*/
for (k=0; k<10; k++)
  {idum=-10*k*k-2+k;
   sigm[k]=sigmat+ran1(&idum)*(1.5-sigmat);
   u[k]=sigm[k]*gasdev(&idum);
   mu[k]+=u[k];
  }
}

为什么不行？

修改后的程序如下，可还是有问题。麻烦看看。（我用GCC编译）

#include <stdio.h>
#include <math.h>

main()
{
int idum, k,s;
double sigmae,sigmat;
double *mu, *u, *sigm,*sigmt;
double *vector(), ran1(), gasdev();

printf("Enter group number s:");
scanf("%d",&s);
sigmae=0.0;
sigmat=0.0;
sigm=vector(0,s);
sigmt=vector(0,s);

for (k=0; k<s; k++)
  {idum=-10*k*k-2+k;
   mu[k]=0.0; u[k]=0.0; sigm[k]=0.0;
   sigm[k]=sigmat+ran1(&idum)*(1.5-sigmat);
   u[k]=sigm[k]*gasdev(&idum);
   mu[k]+=u[k];
  }
}

谢谢。这个问题已经解决，错在没有给mu 和 u 分配空间。