Построить график зависимости количества итераций, необходимых для достижения заданной точности, от параметра релаксации

#include  <iostream>
#include "cmath"
#include "iomanip"
using namespace std;
int sign(double a)
/*variant 9 Mq d=-5 q=-2,85*/
{
 
    if (a > 0) return 1;
 
    else if (a < 0) return -1;
 
}
 
int main()
 
{
 
    int i, j, n, k;
 
    double A[5][5], S[5][5], q, d, o;
 
    cout << " Vvedite q: ";
 
    cin >> q;
 
    cout << " Vvedite d: ";
 
    cin >> d;
 
    cout << " Vvedite n - razmernost' matrici A: ";
 
    cin >> n;
    cout << " Matrix X: ";
 
 
    for (i = 0; i < n; i++) //matrica
 
        for (j = 0; j < n; j++) {
 
            A[i][j] = S[i][j] = 0;
 
            A[0][0] = A[n - 1][n - 1] = q;
 
            for (k = 1; k < n - 1; k++) A[k][k] = -5;
 
        }
 
    for (i = 0; i < n - 1; i++) A[i][i + 1] = A[i + 1][i] = 1;
 
    double B[5], X[5], D[5], Y[5], O[5];
 
    double delta;
 
    for (k = 0; k < n; k++)
 
    {
 
        delta = A[k][k];
 
        if (k != 0) for (i = 0; i <= k - 1; i++) delta = delta - D[i] * S[i][k] * S[i][k];
 
        D[k] = sign(delta);
 
        S[k][k] = sqrt(fabs(delta));
 
        for (j = k + 1; j < n; j++) {
 
            delta = A[k][j];
 
            if (k != 0) for (i = 0; i <= k - 1; i++) delta = delta - D[i] * S[i][k] * S[i][j];
 
            S[k][j] = delta / (S[k][k] * D[k]);
 
        }
 
    }
 
 
 
    B[0] = B[n - 1] = 0;
 
    for (i = 1; i < n - 1; i++) B[i] = d;
 
    Y[0] = B[0] / (S[0][0] * D[0]);
 
    for (i = 1; i < n; i++)
 
    {
 
        delta = B[i];
 
        for (k = 0; k <= i - 1; k++) delta = delta - D[k] * Y[k] * S[k][i];
 
        Y[i] = delta / (S[i][i] * D[i]);
 
    }
 
    X[n - 1] = Y[n - 1] / S[n - 1][n - 1];
 
    for (i = n - 1; i >= 0; i--)
 
    {
 
        delta = Y[i];
 
        for (k = i + 1; k <= n - 1; k++) delta = delta - X[k] * S[i][k];
 
        X[i] = delta / S[i][i];
 
        cout << X[i] << endl;
    }
    return 0;
 
}