ガウス ジョルダン法（連立方程式）(十進basic)

INPUT  m
INPUT  n
DIM  a(m,n)
PRINT  m×n;"個";"の数字を入力してください。"
FOR  x=1  TO   m
   FOR  y=1  TO  n  
      INPUT  w
      LET a(x,y)=w
   NEXT Y
NEXT X
REM  データー例　m=3  n=4     w=2,3,1,4,4,1,-3,-2,-1,2,2,2  解  x=2  y=-1  z=3

DIM  d(m)
REM--------------------------------

FOR   k=1  TO  m
   LET pivot=a(k,k)
   FOR  j=1   TO  n
    
      LET a(k,j)=a(k,j)/pivot
   NEXT  j
   FOR  i=1  TO m
      IF  i<>k  THEN  LET d(i)=a(i,k)
   NEXT  i
   FOR  ii=1  TO  m
      FOR  jj=1  TO  n
       
       
         IF  ii<>k  THEN      LET a(ii,jj)=a(ii,jj)-d(ii)*a(k,jj)  
      NEXT  jj
   NEXT  ii   
    
    
NEXT  k
REM--------------------------------------
FOR  i=1  TO   m
   FOR j=1  TO  n
      PRINT  a(i,j);
   NEXT J
   PRINT
NEXT I


END