ALG'6

276__Rozdziału. Algorytmy numeryczne

double simpson_f(double i *f) (double),//wskaźnik do f(x) double a, double b, int N)

// funkcja zwraca całkę znanej w postaci wzoru // funkcji f(x) w przedziale [a,b],

// N - ilość podziałów {

double s=0,h=(b-a)/(double;N; for(int 1=1;i<=N;i++)

s+=h*(f(a+fi-l)*h)+4*f(a-h/2.0 + i-h!+f(a + i*h) ! /6.0; return s;

I

void main()

{

ccut "Wartość całki =" << simpson(f,-5,3) << endl; cout<<"Wartość całki ="«simpson_t (tun,-5,3,8) « endl;

I

11.6.Rozwiązywanie układów równań liniowych metodą Gaussa

Potrzeba rozwiązywania układów równań liniowych zachodzi w wielu dziedzinach, szczególnie technicznych. Biorąc pod uwagę, że w samym rozwiązywaniu układów równań nie ma nic odkrywczego (uczono nas już tego w szkole podstawowej!), cenne wydaje się dysponowanie procedurą komputerową, która wykona za nas tę żmudną pracę.

Aby komputer mógł rozwiązać dany układ równań, musimy go uprzednio zapisać w postaci rozszerzonej, tzn. nie eliminując współczynników równych zero i pisząc zmienne w określonej kolejności. To wszystko ma na celu prawidłowe skonstruowanie macierzy układu.

Układ równań:

5x+z~9 a:-z+y-6 2x-y+z—0

musi zatem zostać przedstawiony jako:

5x+0y+lz=9 lx+Iy-lz=6 2x-Iy+ lz=0