Задача поменять столбец с максимальным значением на первый элемент матрицы

#include <iostream>
#include <math.h>
#include <iomanip>
using namespace std;
int main()
{
    setlocale(LC_ALL, "rus");
    float x, A[31], B[31], y, dop, sum, maxsum = -100;
    int j, i, iA = 0, k = 0, iB = 0;
    cout << "\t     X \t\t Y" << endl;
    for (x = -2; x <= 4.2; x += 0.2)
    {
        y = 3 * cos(5 * x) + 2 * sin(x);
        cout << setw(15) << x << setw(15) << y << endl;
        if ((int)y % 3 == 0) { A[iA] = y; iA++; }
        if (y < 0) {
            k++;
            if (k % 2 == 0) {
                B[iB] = y; iB++;
 
            }
        }
    }
    cout << "\n\nМассив А\n";
    for (i = 0; i < iA; i++)
        cout << A[i] << "  ";
    cout << "\n\nМассив B\n";
    for (i = 0; i < iB; i++)
        cout << B[i] << "  ";
    for (i = 0; i < iA; i++)
    {
        for (j = 0; j < iA - 1; j++)
            if (A[j] > A[j + 1])
            {
                dop = A[j];
                A[j] = A[j + 1];
                A[j + 1] = dop;
            }
    }
    cout << "\n\nИзмененный массив А\n";
    for (i = 0; i < iA; i++)
        cout << A[i] << "  ";
    for (i = 0; i < iB; i += 2)
    {
        dop = B[i];
        B[i] = B[i + 1];
        B[i + 1] = dop;
    }
    cout << "\n\nИзмененный массив B\n";
    for (i = 0; i < iB; i++)
        cout << B[i] << "  ";
    float S[31][31], max = S[0][j];;
    int  r;
    if (iA < iB) r = iA;
    else r = iB;
    if (r > 10) r = 10;
    cout << "\n\nМатрица S\n";
    for (i = 0; i < r; i++)
    {
        for (j = 0; j < r; j++)
        {
            S[i][j] = pow(A[i], 3) + fabs(A[j] - B[j]);
            cout << setw(15) << setprecision(5) << S[i][j];
 
        }
        cout << endl;
    }
    cout << "\n\nИзмененая матрица S" << endl;
    for (i = 0; i < r; i++)
    {
        for (j = 0; j< r; j++)
        {
            if (S[i][j] > max) max = S[i][j];
            S[i][j] = pow(A[i], 3) + fabs(A[j] - B[j]);
            S[i][] = S[0][0];
            cout << setw(15) << setprecision(5) << S[i][j];
        }
        cout << endl;
    }
    cout << "\n\nMax =" << max;
}