Tengo que hacer un programa que imprima un cuadrado latino de orden N

Ahí van el algoritmo y su implementación en C ya corregidos:
Inicio
para i de 1 a n hacer
para j de 1 a n hacer
si i=1 entonces
matriz[j] <-- j
sino si matriz[i-1][j]=1 entonces
matriz[j] <-- n
Sino
Matriz[j] <-- matriz[i-1][j] - 1
Finsi
Finpara
Finpara
Fin
Main()
{
int **matriz;
int n, i, j;
system("clear");
printf("Introduzca 'n': ");
scanf("%d", &n);
printf("\n");
matriz=(int **)malloc(n*(sizeof(int *)));
for(i=0;i<n;i++)
matriz=(int *)malloc(n*(sizeof(int)));
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
{
if(i==0)
matriz[j]=j+1;
else if(matriz[i-1][j]==1)
matriz[j]=n;
else
matriz[j]=matriz[i-1][j]-1;
printf("%d ", matriz[j]);
}
printf("\n");
}
printf("\n");
}

Las funciones para resolver tu problema son las siguientes
MAX_N es la dimension máxima para la matriz, y N debe ser <= MAX_N y >=1, de lo contrario no la llena
#include <memory.h>
#define MAX_N 10
typedef int tipo;
int Rota(tipo *v, int n, int m) {
tipo *p;
if (m<0 || m>n)
return 1; //no se puede rotar
if (m==n)
return 0; //no hay necesidad de rotar
p = (tipo *)malloc(m*sizeof(tipo));
memcpy(p, &v[n-m], m*sizeof(tipo));
memmove(&v[m], v, (n-m)*sizeof(tipo));
memcpy(v, p, m*sizeof(tipo));
free(p);
return 0;
}
int LlenaCuadroLatino(tipo cl[MAX_N][MAX_N], int n) {
int i;
if (n<1 || n>MAX_N)
return 1;
for (i=0; i<n; i++)
cl[0] = i+1;
for (i=1; i<n; i++) {
memcpy(cl, cl[i-1], n*sizeof(tipo));
Rota(cl, n, 1);
}
return 0;
}