Возвести в квадрат

#include "stdafx.h"
#include "iostream"
using namespace std;
#define _USE_MATH_DEFINES
#include "cmath"
 
int main(int argc, char* argv[])
{
    setlocale(LC_ALL,".1251");
    double a[3];
    int i;
    cout << "Введите три действительных числа:" << endl;
    cin >> a[0] >> a[1] >> a[2];
    _flushall();
    for(i = 0;i < 3;i++)
        if (a[i] > 0)
            a[i] = pow(a[i],2);
        else
            a[i] = pow(a[i],4);
    cout << "Три числа после всех манипуляций:" << endl;
    for(i = 0;i < 3;i++)
        cout << a[i] << "  ";
    system("pause");
    return 0;
}

#include "stdafx.h"
#include "iostream"
using namespace std;
 
 
int main(int argc, char* argv[])
{
    setlocale(LC_ALL,".1251");
    int i,n,
        k = 0;
    cout << "Сколько дано натуральных чисел?" << endl;
    cin >> n;
    _flushall();
    for(i = 1;i <= n;i++)
        k += i;
    cout << "Сумма " << n << " натуральных чисел равняется " << k << endl;
    system("pause");
    return 0;
}

#include "stdafx.h"
#include "iostream"
using namespace std;
#include <ctime>
#define n 4
#define k 3
 
int main(int argc, char* argv[])
{
    setlocale(LC_ALL,".1251");
    double a[n][k],p;
    int i,j;
    srand(time(0));
    for(i = 0; i < n;i++){
        for(j = 0;j < k;j++){
            a[i][j] = rand()%401-200;
            cout << a[i][j] << "  ";
        }
        cout << endl;
    }
    p = a[1][1]*a[2][2]*a[3][2]*a[4][3];
    cout << "Произведение элементов главной диагонали равно: " << p << endl;
    system("pause");
    return 0;
}