Написать двойную функцию doubl(double [], double []) вычисляющая скалярное произведения векторов, заданных как параметры

#include <iostream>
#include <time.h>
#include <stdlib.h>
 
using namespace std;
 
double il(double a[3], double b[3])
{
    int i;
    return a[i]*b[i];
}
 
int main()
{
    double a[3], b[3];
    int k;
    for (int i=0;i<3;i++) {cout<<"Podaj wartosc "<<i<<"-go elementu: "; cin >> a[i]; cout<<endl;}
    for (int i=0;i<3;i++) {cout<<"Podaj wartosc "<<i<<"-go elementu: "; cin >> b[i]; cout<<endl;}
    for (int i=0;i<3;i++)
    k+= il(a[i],b[i]);
    cout << k ;          
 
    return 0;
}

double dot(double* a, double* b){
    return a[0]*b[0] + a[1]*b[1] + a[2]*b[2];
}

//g++  5.4.0
 
#include <iostream>
 
using namespace std;
 
constexpr size_t VECTOR_SIZE = 3; // Чтобы не было магических констант
 
// Делаем более понятное имя функции
double scalar_product(double a[VECTOR_SIZE], double b[VECTOR_SIZE]) 
{
    // Засовываем всё вычисление произведения в функцию
    double sum = 0;
    for (size_t i = 0; i<VECTOR_SIZE; ++i)
        sum += a[i]*b[i];
    return sum;
}
 
int main()
{
    double a[VECTOR_SIZE], b[VECTOR_SIZE];
    double k; // скалярное произведение должно быть того же типа, что и множители 
    for (int i=0;i<3;i++) {cout<<"Podaj wartosc "<<i<<"-go elementu: "; cin >> a[i]; cout<<endl;}
    for (int i=0;i<3;i++) {cout<<"Podaj wartosc "<<i<<"-go elementu: "; cin >> b[i]; cout<<endl;}
    // убрали цикл, так как вычисление скалярного произведения — один вызов функции
    k = scalar_product(a,b);
    cout << k ;          
 
    return 0;
}