C数列或者C向量以及C矩阵

#include <stdlib.h>
#include <stdio.h>
#define TP double
#define UI unsigned short int
#define SI signed short int
/*1维度n维数C指针（C数列或者C向量）*/
TP *vec(UI n)
{
 TP *s=(TP *)calloc(n,sizeof(TP));
 UI i;
 for(i=0;i<n;i++){s[i]=0;}
 return s;
 free(s);
}
/*C函数vec示例*/
void vector()
{
 TP *s;
 UI n;
 UI i;
 printf("C数列或者C向量的维数n\n");
 scanf("%u",&n);
 s=vec(n);
 printf("C数列或者C向量的所有元素\n");
 for(i=0;i<n;i++){scanf("%lf",&s[i]);}
 printf("C数列或者C向量的所有元素\n");
 for(i=0;i<n;i++){printf("%lf\t",s[i]);}
 printf("\n");
}
/*2维度m行n列C指针（C矩阵）*/
TP **mat(UI m, UI n)
{
 UI i, j;
 TP **A=(TP **)calloc(m,sizeof(TP *));
 for(i=0;i<m;i++)
 {A[i]=(TP *)calloc(n,sizeof(TP));}
 for(i=0;i<m;i++){for(j=0;j<n;j++){A[i][j]=0;}}
 return A;
 free(A);
}
/*C函数mat示例*/
void matrix()
{
 TP **A;
 UI m, n;
 UI i, j;
 printf("C矩阵的行数m，列数n\n");
 scanf("%u%u",&m,&n);
 A=mat(m,n);
 printf("C矩阵的所有元素\n");
 for(i=0;i<m;i++)
 {
  for(j=0;j<n;j++){scanf("%lf",&A[i][j]);}
 }
 printf("C矩阵的所有元素\n");
 for(i=0;i<m;i++)
 {
  for(j=0;j<n;j++){printf("%lf\t",A[i][j]);}
  printf("\n");
 }
}

#define N 2
void main()
{
 TP (*JuZhen)[N];//m行N列的C矩阵
 UI m;
 UI i, j;

 vector();
 system("pause");
 matrix();
 system("pause");

 printf("C矩阵的行数m\n");
 scanf("%u",&m);
 JuZhen=(TP (*)[N])calloc(m*N,sizeof(TP));
 printf("C矩阵的所有元素（2列）\n");
 for(i=0;i<m;i++)
 {
  for(j=0;j<N;j++)
  {
   scanf("%lf",&JuZhen[i][j]);
  }
 }
 printf("C矩阵的所有元素（2列）\n");
 for(i=0;i<m;i++)
 {
  for(j=0;j<N;j++)
  {
   printf("%lf\t",JuZhen[i][j]);
  }
  printf("\n");
 }
 free(JuZhen);
 system("pause");
}