21 unsigned int N=xPoints.size()-1;
26 for (
unsigned int i=1;i<N;i++) {
28 double dxPlus = xPoints[i+1].first -xPoints[i].first;
29 double dyPlus = xPoints[i+1].second-xPoints[i].second;
30 double mPlus = dyPlus/dxPlus;
32 double dx = xPoints[i].first -xPoints[i-1].first;
33 double dy = xPoints[i].second-xPoints[i-1].second;
37 A[i][i] = 2*(1/dxPlus + 1/dx);
40 V[i] = 3*(m/dx+mPlus/dxPlus);
45 double dx = xPoints[1].first -xPoints[0].first;
46 double dy = xPoints[1].second-xPoints[0].second;
54 double dx = xPoints[N].first -xPoints[N-1].first;
55 double dy = xPoints[N].second-xPoints[N-1].second;