shiyan4.cpp

求根问题

#include <stdio.h>

#include <stdlib.h>



void LDLTDCMP (unsigned int n, double * *a)

{


   //tf ("\n\nCalculate and save result to ASCII file LDLT_KEY.dat ...\n");



    /* Calculate answers. */

    LDLTDCMP (n, a);

    LDLTBKSB (n, a, b);



    /* Show and save answers. */

    printf ("\nAfter calculated:\n");

    fprintf(fp_ldl_out, "\nAfter calculated:\n");

    for (i = 0; i < n; i++)

    {

        for (j = 0; j < n; j++)

        {

            printf ("%8.3g  ", a[i][j]);

            fprintf(fp_ldl_out, "%8.3g  ", a[i][j]);

        }

        printf ("%8.3g\n", b[i]);

        fprintf(fp_ldl_out, "%8.3g\n", b[i]);

    }

    printf ("\nAnswers x[]:\n");

    fprintf(fp_ldl_out, "\nAnswers x[]:\n");

    for (i = 0; i < n; i++)

    {

        printf ("%8.3g", b[i]);

        fprintf(fp_ldl_out, "%8.3g", b[i]);

    }

    fclose (fp_ldl_out);

    printf ("\n\nAll results are saved.\n");

    printf ("Please open ASCII file ../Output/LDLT_KEY.dat to check it.\n");

    for (i = 0; i < n; i++)  

        free (*(a + i));  

    free (a);  

    free (b);  

}