すみません、このプログラムの境界条件Tの宣言がうまくいきません。
どこをどう直せばよいのでしょうか？

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#define NX 41
#define NY 41

void SetBoundaryCondition(double T[NY][NX])
{
int x;
int y;
for (x = 0; x < NX; x++)
{T[0][x] = 0.5;}
for (x = 0; x < NX; x++)
{T[NY-1][x] = 0.0;}
for (y = 0; y < NY; y++)
{T[y][0] = 1.0;}
for (y = 0; y < NY; y++)
{T[y][NX-1] = 0.0;}

＊つづく＊