>>103のプログラム
#include <stdio.h>
#define M_DIM (3)
void sol_gj(double [][M_DIM+1],int);
void sol_g(double [][M_DIM+1],int);
int main(void) {
int n=M_DIM,k,j;
double a[M_DIM][M_DIM+1],aa[M_DIM][M_DIM+1];
for(k=0;k<n;k++){
printf("(%d)ax+by+cz=d a,b,c,d?",k+1);
scanf("%lf,%lf,%lf,%lf",&a[k][0],&a[k][1],&a[k][2],&a[k][3]);
}
for(k=0;k<n;k++)
for(j=0;j<n+1;j++)
aa[k][j]=a[k][j];
sol_gj(a,n);
for(k=0;k<n;k++)
for(j=0;j<n+1;j++)
a[k][j]=aa[k][j];
sol_g(a,n);
return 0;
}