请补充函数fun()，该函数的功能是：计算N×N维矩阵元素的方差，结果由函数返回。维数N在主函数中输入。例如：

注意：部分源程序给出如下。

请勿改动主函数main和其他函数中的任何内容，仅在函数fun()的横线上填入所编写的若干表达式或语句。

试题程序：

#include ＜stdio.h＞

#include ＜conio.h＞

#include ＜stdlib.h＞

#include ＜math.h＞

#define N 20

double fun( 【1】 ，int n)

int i,j；

int k；

double s=0.0；

double f=0.0；

double aver=0.0；

double sd=0.0；

for(i=0；i＜n；i++)

for(j=0；j＜n；j++)

s+=a[i][j]；

aver= 【2】 ;

for(i=0；i＜n；i++)

for(j=0；i＜n；j++)

f+=(a[i][j]-aver)*(a[i][j]-aver)；

f/=(n*n)；

sd= 【3】 ；

return sd；

main()

int a[N][N]；

int n；

int i，j；

double s；

clrscr()；

printf("***+Input the dimension of

array N*****\n")；

scanf("%d"，&n)；

printf("***** The array *****\n")；

for(i=0；i＜n；i++)

for(j=0；j＜n；j++)

a[i][j]=rand()%50；

while(a[i][j]＝0)

a[i][j]=rand()%60；

printf("%4d"，a[i][j])；

printf("\n\n")；

s=fun(a，n)；

printf("******* THE RESULT *******\n")；

printf("%4.3f\n"，s)；