With Ada.Text_IO with _Real_Arrays procedure Gaussian_Eliminations is type Real is new Float package Real_Arrays is new _Real_Arrays (Real ) use Real_Arrays function Gaussian_Elimination ( A : in Real_Matrix B : in Real_Vector ) return Real_Vector is procedure Swap_Row ( A : in out Real_Matrix B : in out Real_Vector R_1, R_2 : in Integer ) is Temp : Real begin if R_1 = R_2 then return end if - Swal matrix row for Col in A ' Range ( 1 ) loop Temp := A ( R_1, Col ) A ( R_1, Col ) := A ( R_2, Col ) A ( R_2, Col ) := Temp end loop - Swap vector row Temp := B ( R_1 ) B ( R_1 ) := B ( R_2 ) B ( R_2 ) := Temp end Swap_Row AC : Real_Matrix := A BC : Real_Vector := B X : Real_Vector ( A ' Range ( 1 )) := BC Max, Tmp : Real Max_Row : Integer begin if A ' Length ( 1 ) /= A ' Length ( 2 ) or A ' Length ( 1 ) /= B ' Length then raise Constraint_Error with 'Dimensions do not match' end if if A ' First ( 1 ) /= A ' First ( 2 ) or A ' First ( 1 ) /= B ' First then raise Constraint_Error with 'First index must be same' end if for Dia in Ac ' Range ( 1 ) loop Max_Row := Dia Max := Ac ( Dia, Dia ) for Row in Dia + 1. L R13,4(0,R13) restore previous savearea pointer

Gauss seidel method code python