#include <imsl.h>

main()
{
    int         n = 3;
    float       *x;
    static float        a[] = { 1.0, 3.0, 3.0,
                                1.0, 3.0, 4.0,
                                1.0, 4.0, 3.0 };
    static float        b[] = { 1.0, 4.0, -1.0 };

    /*
     * Verify the version of the library we are running by
     * retrieving the version number via imsl_version().
     * Verify correct installation of the error message file
     * by retrieving the customer number via imsl_version().
     */
    char        *library_version = imsl_version(IMSL_LIBRARY_VERSION);
    char        *customer_number = imsl_version(IMSL_LICENSE_NUMBER);

    printf("Library version:  %s\n", library_version);
    printf("Customer number:  %s\n", customer_number);

                                /* Solve Ax = b for x */
    x = imsl_f_lin_sol_gen(n, a, b, 0);
                                /* Print x */
    imsl_f_write_matrix("Solution, x of Ax = b", 1, n, x, 0);
}