请补充函数fun,它的功能：计算并输出n（包括n）以内所有能被3或7整除的自然数的倒数之和.
注意：部分源程序给出如下.请勿改动main函数和其它函数中的任何内容,仅在函数fun的横线上填入所编写的若干表达式或语句.
#include
doublefun(int n)
{
intk;
double sum=0.0;
for( k=1; ①;k++)
if( k%3==0② k%7==0)sum+=③ /k;
return sum;
}
voidmain()
{
int n ;
doubles;
printf(“\n input n:”);
scanf(“%d”,&n);
s=fun(n);
printf(“\n\ns=%f\n”,s);
}