Двумерные массивы. Сумма элементов столбца. (Си)

[Sweta]

#include <stdio.h>
#include <math.h>
#include <conio.h>
# define N 20

void Str (double A[][N], int a, int b);
void Stb (double A[][N], int a, int b);
double s(int i, int j);
double Summ(double A[][N], int k, int b);

int main()
{
	double A[N][N], B[N][N];
	int M=20,a,b,k;

	do
	{

		printf ("Введите a: ");
		scanf ("%u", &a);
		printf ("Введите b: ");
		scanf ("%u", &b);
	}
	while (a>M && b>M);
	printf("\n");
	Str (A,a,b);
	Stb (A,a,b);
	do
	{
		printf ("\nВведите k: ");
		scanf ("%u", &k);
	}
	while (0>=k||k>b);
	printf ("\nСумма элементов k-столбца=%-7.4lf", Summ(A,k,a));
	getch();
	return 0;
  }

void Str (double A[][N], int a, int b)
{
	int i,j;
	for (i=0; i<a; i++)
        	for (j=0; j<b; j++)
	     	A[i][j]=s(i,j);
// 	return i,j;
}

void Stb (double A[][N], int a, int b)
{
 int i,j;
 for (i=0; i<a; i++)
  {
   for (j=0; j<b; j++)
   printf ("A[%d][%d]=%-7.4lf", i, j, A[i][j]);
   printf ("\n");
   }
//	return i,j;
}

double s(int i, int j)
{
	if (i<=j)
        return 1.0/((i*i)+2);
        else
	if (i>j)
	return (1.0/(i+j));
}

double Summ(double A[][N], int k, int a)
{
	int i;
	double sum=A[0][k];

	for (i=1; i<=a; i++)
	sum+=A[i][k];
	return sum;
}

У меня ХР глюков нет.