请补充函数fun,它的功能：计算并输出n（包括n）以内所有能被3或7整除的自然数的倒数之和.

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