Не могу исправить ошибку в программе вычисляющей интерполяционные полином Лагранжа

#include <iostream>
 #include <vector>
 #include <functional>
 
 using namespace std;
 
 template <class T>
 T lagrange(vector<T> xs, vector<T> ys, T x) {
 T res = 0;
 T mul;
 size_t len = xs.size();
 for (size_t i=0; i<len; i++) {
 mul = 1;
 for(size_t j=0; j<len; j++) {
 if (i != j) {
 mul *= (x-xs[j])/(xs[i]-xs[j]);
 }
 }
 res += mul * ys[i];
 }
 }
 
 int main()
 {
 vector<double> xs = {.8, .9 , 1.0, 1.1, 1.2}; ///<ошибка>
 vector<double> ys = {-2, -0.1, 0 , 0.1, 0.1}; ///<ошибка>
 
 double x;
 
 while(true) {
 cout << "x ? ";
 cin >> x;
 cout << "lagrange(" << x << ") = " << lagrange(xs, ys, x) << '\n';
 }
 }

#include <iostream>
#include <vector>
#include <functional>
 
using namespace std;
 
template <class T>
T lagrange(vector<T> xs, vector<T> ys, T x) {
    T res = 0;
    T mul;
    size_t len = xs.size();
    for (size_t i=0; i<len; i++) {
        mul = 1;
        for(size_t j=0; j<len; j++) {
            if (i != j) {
                mul *= (x-xs[j])/(xs[i]-xs[j]);
            }
        }
        res += mul * ys[i];
    }
}
 
int main() {
    vector<double> xs = {0.8, 0.9 , 1.0, 1.1, 1.2};  // я так понимаю тут 0.8, а не .8
    vector<double> ys = {-2.0, -0.1, 0.0 , 0.1, 0.1};
 
    double x;
 
    while(true) {
        cout << "x ? ";
        cin >> x;
        cout << "lagrange(" << x << ") = " << lagrange(xs, ys, x) << '\n';
    }
}

vector<double> xs;
xs.push_back(0.8);
xs.push_back(0.9);
xs.push_back(1.0);
...

double _xs_[] = {.8, .9 , 1.0, 1.1, 1.2};
double _ys[] = {-2, -0.1, 0 , 0.1, 0.1};
 
vector<double> xs(_xs, _xs + 5);
vector<double> ys(_ys, _ys + 5);

#include <iostream>
#include <vector>
#include <functional>
 
using namespace std;
 
template <class T>
T lagrange(vector<T> xs, vector<T> ys, T x) {
    T res = 0;
    T mul;
    size_t len = xs.size();
    for (size_t i=0; i<len; i++) {
        mul = 1;
        for(size_t j=0; j<len; j++) {
            if (i != j) {
                mul *= (x-xs[j])/(xs[i]-xs[j]);
            }
        }
        res += mul * ys[i];
    }
}
 
int main()
{
vector<double> xs;
xs.push_back(0.8);
xs.push_back(0.9);
xs.push_back(1.0);
xs.push_back(1.1);
xs.push_back(1.2);
 
vector<double> ys;
xs.push_back(-2);
xs.push_back(-0.1);
xs.push_back(0);
xs.push_back(0.1);
xs.push_back(0.1);
 
 
    double x;
    
    while(true) {
        cout << "x ? ";
        cin >> x;
        cout << "lagrange(" << x << ") = " << lagrange(xs, ys, x) << '\n';
    }
    return 0;
}