[R] rgamma function

[Chandler Zuo]

Hi,

Has anyone encountered the problem of rgamma function in C? The 
following simplified program always dies for me, and I wonder if anyone 
can tell me the reason.

#include <Rmath.h>
#include <time.h>
#include <Rinternals.h>

SEXP generateGamma ()
{
     srand(time(NULL));
     return (rgamma(5000,1));
}

Has anyone encountered a similar problem before? Is there another way of 
generating Gamma random variable in C?

P.S. I have no problem compiling and loading this function in R.

Thanks for suggestions in advance!

--Chandler