Thank you very much again. Your code is really helpful and it demontrates the power of R to a new-comer like me.
--- [EMAIL PROTECTED] wrote: > Understood! Though I think you mean > > 1/c = sum(1.0 / pow((double) j, alpha)), > j=1,2,...,N You are right. Weiguang ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html