void main()
{
setup_oscillator(OSC_NORMAL);
while (true)
{
output_A(0xff);
output_B(0xff);
output_C(0xff);
output_D(0xff);
output_E(0xff);
output_high(pin_C1);
output_high(pin_C0);
output_high(pin_A4);
Delay_ms(500);
output_A(0x0);
output_B(0x0);
output_C(0x0);
output_D(0x0);
output_E(0x0);
output_low(pin_C1);
output_low(pin_C0);
output_low(pin_A4);
Delay_ms(500);
}
}
void afficheur(int valeur){
Set_tris_A(0b00011110);
if(valeur<10 && valeur >=0)
{
output_high(PIN_A1);
valeur = valeur << 1;
while(flag){
output_A(valeur);
}
}
}
int main(int argc, char **argv)
{
int n;
double **pFile, **e, **s;
if (argc == 1) {
printf("Please enter the file name which contain Matrix data!\n\nUsage: %s <filename>\n", argv[0]);
exit(EXIT_FAILURE);
}
pFile = readMatrix(argv[1], &n, 1);
#ifdef DEBUG
output_A(pFile, n, n);
#endif
e = eliminate((const double **)pFile, n, n+1);
if (!e)
exit(EXIT_FAILURE);
printf("RESULT...\n");
output_A(e, n, n+1);
s = solve(e, n);
printf("方程的解为:\n");
output_A(s, 1, n);
// 求逆阵
double **T = readMatrix("gauss_inverse", &n, 0);
double **M_I = matrix_I(n);
output_A(T, n, n);
#ifdef DEBUG_INVERSE
printf("输出单位矩阵:\n");
output_A(M_I, n, n);
#endif
double **II = invers(T, M_I, n);
revers(II, n, n);
printf("输出逆阵:\n");
output_A(II, n, n);
printf("测试:\n");
output_A(multi(T, II, n), n, n);
return 0;
}
