Информационный сервер для программистов: Исходники со всего света. Паскальные исходники со всего света
  Powered by Поисковый сервер Яndex: Найдется ВСЁ!
На Главную Исходники Форум Информер Страны мира
   Математика    >>    solvgj
   
 
 Solve Linear Equations System with the Gauss-Jordan  Unknown 03.04.1985

Четыре варианта программы для решения системы линейных уравнений с использованием метода Гаусса-Жордана. Причем, уравнений может быть больше, чем неизвестных.



5k 
 

program solvec; { -> 119 } { pascal program to perform simultaneous solution by Gauss-Jordan elimination for complex coefficients } const maxr = 8; maxc = 8; type ary = array[1..maxr] of real; arys = array[1..maxc] of real; ary2s = array[1..maxr,1..maxc] of real; aryc2 = array[1..maxr,1..maxc,1..2] of real; aryc = array[1..maxr,1..2] of real; var y : arys; coef : arys; a,b : ary2s; n,m,i,j : integer; error : boolean; external procedure cls; external procedure revon; external procedure revoff; procedure get_data(var a: ary2s; var y: arys; var n,m: integer); { get complex values for n and arrays a,y } var c : aryc2; v : aryc; i,j,k,l : integer; procedure show; { print original data } var i,j,k : integer; begin { show } writeln; for i:=1 to n do begin for j:=1 to m do for k:=1 to 2 do write(c[i,j,k]:7:4,' '); writeln(':',v[i,1]:7:4,':',v[i,2]:7:4) end; n:=2*n; m:=n; writeln; for i:=1 to n do begin for j:=1 to m do write(a[i,j]:7:4,' '); writeln(':',y[i]:9:5) end; writeln end; { show } begin { procedure get_data } writeln; repeat write('How many equations? '); readln(n); m:=n until n<maxr; if n>1 then begin for i:=1 to n do begin writeln('Equation',i:3); k:=0; l:=2*i-1; for j:=1 to n do begin k:=k+1; write('Real',j:3,':'); read(c[i,j,1]); { read real part } a[l,k]:=c[i,j,1]; a[l+1,k+1]:=c[i,j,1]; k:=k+1; write('Imag',j:3,':'); read(c[i,j,2]); { imaginary part } a[l,k]:=-c[i,j,2]; a[l+1,k-1]:=c[i,j,2] end; { j-loop } write('Real const:'); read(v[i,1]); { real constant } y[l]:=v[i,1]; write('Imag const:'); readln(v[i,2]); { imag constant } y[l+1]:=v[i,2] end; { i-loop } show { the original DATA } end { if n>1 } end; { procedure get_data } procedure write_data; { print out