急求：高斯约当法解线性方程组的通用C程序

#include <stdio.h>
#include <conio.h>
#include <alloc.h>
#include <math.h>

int main()
{
int i,j,n;
double *A,*b,*p,*x,temp;
int gauss_jodan(double *A,double *b,int n,double *x);

n=3; /*3元方程组*/
A=calloc(n*n, sizeof(double));
b=calloc(n,sizeof(double));
x=calloc(n,sizeof(double));
p=&temp;

printf("n=");
scanf("%d",&n);

printf("A=\n");
for (i=0;i<n*n;i++)
scanf("%lf",&A[i]);

printf("b=\n");
for (i=0;i<n;i++)
scanf("%lf",&b[i]);

printf("==================\n");
printf("A=\n");
for (i=0;i<n;i++)
{
for (j=0;j<n;j++)
printf("%lf\t",A[i*n+j]);
printf("\n");
}
printf("b=\n");
for (i=0;i<n;i++)
p