Матричные операции

DECLARE FUNCTION Y# (Mas() AS DOUBLE)
 e = 2.718281828459#
 CLS
 INPUT "Vvedite n : ", n
DIM A(n, n) AS DOUBLE, X(n) AS DOUBLE
   ' Расчет значений матрицы
   FOR j = 1 TO n: FOR i = 1 TO n
     A(i, j) = (2 ^ i + e ^ j) / 3 ^ j
   NEXT: NEXT
   ' Расчет значений вектора
   FOR i = 1 TO n
     FOR k = 1 TO i
     X(i) = X(i) + A(i, k)
   NEXT: NEXT
   ' Упорядочивание элементов
   FOR j = 1 TO n - 1
     FOR i = 1 TO n - j
     IF X(i) < X(i + 1) THEN SWAP X(i), X(i + 1)
   NEXT: NEXT
   ' Вывод результата
   FOR j = 1 TO n: FOR i = 1 TO n
     C$ = STR$(A(i, j))
     C$ = LEFT$(C$, INSTR(C$, ".") - 1) + MID$(C$, INSTR(C$, "."), 3) + " "
     LOCATE i + 1, j * 5 - 4: PRINT C$
       NEXT: C$ = STR$(X(j))
     C$ = LEFT$(C$, INSTR(C$, ".") - 1) + MID$(C$, INSTR(C$, "."), 5) + " "
     LOCATE n + 3, (j * 7 - 6) MOD 80: PRINT C$
   NEXT
   '  Вывод функции
   PRINT " Y ="; Y#(X())
 
FUNCTION Y# (Mas() AS DOUBLE)
  FOR ii = 1 TO UBOUND(Mas)
   S = S + Mas(ii) + Mas(n)
  NEXT
  Y# = S / 2
END FUNCTION