Реферат Рішення систем лінійних рівнянь за методом Гаусса

nclude ObratniHod (int n, double ** a, double * b, double * x);

# endif

// obratni_hod.cpp

// реалізація модуля "obratni_hod"

# include "obratni_hod.h"

ObratniHod (int n, double ** a, double * b, double * x)

{s = 0; [n - 1] = 1.0 * b [n - 1]/a [n - 1] [n - 1]; (int i = n - 2, j; 0 <= i; i -)

{= 0; (j = i + 1; j

{= s + a [i] [j] * x [j];

} [i] = 1.0 * (b [i] - s)/a [i] [i];

}

}

7.3. Модуль В«pryamoi_hodВ»

// pryamoi_hod.h

// інтерфейс модуля "pryamoi_hod.h"

# ifndef pryamoi_hod_h

# define pryamoi_hod_h

# include // i/oShowVector (int n, double * vec);

# endif

// pryamoi_hod.cpp

// реалізація модуля "pryamoi_hod"

# include "pryamoi_hod.h"

PryamoiHod (int n, double ** a, double * b)

{v; (int k = 0, i, j, im; k

{= k; (i = k + 1; i

{(fabs (a [im] [k])

{= i;

}

} (im! = k)

{(j = 0; j

{= a [im] [j]; [im] [j] = a [k] [j]; [k] [j] = v;

} = b [im]; [im] = b [k]; [k] = v;

} (i = k + 1; i

{= 1.0 * a [i] [k]/a [k] [k]; [i] [k] = 0; [i] = b [i] - v * b [ k]; (v! = 0) (j = k + 1; j

{[i] [j] = a [i] [j] - v * a [k] [j];

}

}

}

}

.4. Модуль В«show_vectorВ»

// show_vector.h

// інтерфейс модуля "show_vector"

# ifndef show_vector_h

# define show_vector_h

# include

# include

PryamoiHod (int n, double ** a, double * b);

# endif

// show_vector.cpp

// реалізація модуля "show_vector"

# include "show_vector.h"

ShowVector (int n, double * vec)

{(int i = 0; i

}