用fun函数实现s=(ln(1)+ln(2)+...+ln(m))^1/2,函数返回s

#include<math.h>
#include<stdio.h>
void NONO();
double fun(int m)
{
double a=0;
for(int i=1;i<=m;i++){
a+=log(i);
//c语言中可调用log(n)函数求ln(n),log()的引用说明：double log(double x)
}
double s=sqrt(a);
//sqrt()为开方函数
return s;
}

main()
{
printf("%f\n",fun(20));
NONO();

getchar();

}
void NONO()
{
FILE *fp,*wf;
int i,n;
double s;

fp=fopen("in.dat","r");
wf=fopen("out.dat","w");
for(i=0;i<10;i++){
	fscanf(fp,"%d",&n);
	s=fun(n);
	fprintf(wf,"%f\n",s);
}
fclose(fp);
fclose(wf);

}