#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include "interpolation.h"

using namespace alglib;


int main(int argc, char **argv)
{
   real_1d_array x = "[-4.71, -3.77, -2.83, -1.89, -0.25]";
  real_1d_array y = "[0, 1, 2, 3, 4]";
  ae_int_t m = 3;
  double t = 0;
  ae_int_t info;
  barycentricinterpolant p;
  polynomialfitreport rep;
  double v;
  polynomialfit(x, y, m, info, p, rep);
  v = barycentriccalc(p, t);
  printf("%.2f\n", double(v));
  real_1d_array a_eqdist;
  polynomialbar2pow(p, a_eqdist);
  printf("%s\n", a_eqdist.tostring(4).c_str());
  return 0;
}