C 2차원 배열 동적 해제 질문입니다

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 
void main()
{
	int i, j, k;
	int elementNumber;
	char line_buf[256];
	char char_buf[10];
	char temp[10];
	int sign;
	//int temp;
 
	//Gauss Elimination
	FILE *in1 = NULL;
	in1 = fopen("hw1_in1", "r");
	fscanf(in1, "%d", &elementNumber);
	//elementNumber++;
 
	//입력받은 수로 2차원 배열 동적할당
	int **matrix;
	matrix = (int **)malloc(sizeof(int)*(elementNumber));
	for(i = 0; i <= elementNumber; i++)
		matrix[i] = (int *)malloc(sizeof(int)*(elementNumber));
 
 
	//파일로부터 행렬 받아옴
	i = 0, j = 0;
	for(i; i < elementNumber; i++)	{
		for(j; j <= elementNumber; j++)	{
			fscanf(in1, "%s", &temp);
			if( (temp[0] == 45) || (temp[0] == 43) || (temp[0] == 61) )	{ //계수의 부호처리
				if( temp[0] == 43 ) //계수가 +일 경우
					sign = 1;
				else if( temp[0] == 45 ) //계수가 -일 경우
					sign = -1;
				*temp = NULL;
				fscanf(in1, "%s", &temp);
			}
 
			if( temp[0] == 48 )
				matrix[i][j] = 0;
			else if( (temp[0] != 48) && (sign == -1) )	{
				matrix[i][j] = -atoi(temp);
				if( matrix[i][j] == 0 )
					matrix[i][j] = -1;
			}
			else if( (temp[0] != 48) && (sign == 1) )	{
				matrix[i][j] = atoi(temp);
				if( matrix[i][j] == 0 )
					matrix[i][j] = 1;
			}
			else if( (temp[0] != 48) && (atoi(temp) == 0) )
				matrix[i][j] = 1;
			else
				matrix[i][j] = atoi(temp);
			sign = 0;
			*temp = NULL;
		}
		j = 0;
	}
 
 
	//for debug
	i = 0, j = 0;
	for(i; i < elementNumber; i++)	{
		for(j; j <= elementNumber; j++)	{
			printf("%d\t", matrix[i][j]);
		}
		printf("\n");
		j = 0;
	}
 
	///*
	//매트릭스 동적해제
	for(i = 0; i < elementNumber; i++)	{
		free(matrix[i]);
	}
	free(matrix);
	//*/
 
}