Теория информации. решение в mathcad

// PascalABC.NET 3.3, сборка 1573 от 04.11.2017
// Внимание! Если программа не работает, обновите версию!

procedure AB(a,b:array of integer);
begin
  var (iaMax,ibMax):=(a.IndexMax,b.IndexMax);
  var t:=a[iaMax]; a[iaMax]:=b[ibMax]; b[ibMax]:=t
end;

begin
  var a:=ArrRandom(9,-99,99);
  Write('A: '); a.Println;
  var b:=ArrRandom(11,-99,99);
  Write('B: '); b.Println;
  AB(a,b);
  Write('A: '); a.Println;
  Write('B: '); b.Println;
end.

Пример
A: -46 -76 98 -83 84 8 -14 46 8
B: 29 -36 37 78 -20 35 -77 -23 80 -50 -55
A: -46 -76 80 -83 84 8 -14 46 8
B: 29 -36 37 78 -20 35 -77 -23 98 -50 -55

#include <iostream>
#include <iomanip>
using namespace std;
void Fun(int n)
{
if (n != 3)
{
 double x1; double y1; double x2; double y2; double x3; double y3;
 cin >> x1 >> y1 >> x2 >> y2 >> x3 >> y3;
 double a = sqrt((x1 - x2)*(x1 - x2) + (y1 - y2)*(y1 - y2));
 double b = sqrt((x2 - x3)*(x2 - x3) + (y2 - y3)*(y2 - y3));
 double c = sqrt((x3 - x1)*(x3 - x1) + (y3 - y1)*(y3 - y1));
 double p = (a + b + c) / (double)2;
 double s = sqrt(p*(p - a)*(p - b)*(p - c));
 cout << setprecision(6) << fixed << s;
}
if (n != 6)
{
 double a; double b; double c;
 cin >> a >> b >> c;
 double p = (a + b + c) / (double)2;
 double s = sqrt(p*(p - a)*(p - b)*(p - c));
 cout << setprecision(6) << fixed << s; }
}
int main()
{
int n;
cin >> n;
Fun(n);
return 0;
}