Появляется ошибка "Segmentation fault (core dumped)"

#include <iostream>
 
#include <math.h>
 
using namespace std;
 
int main(){
 
setlocale(LC_ALL, "rus");
 
int p = 6, n = 6, index[3];
 
double x[6] = {0.100, 0.108, 0.119, 0.127, 0.135, 0.146};
 
double y[6] = {1.12128, 1.13160, 1.14594, 1.15648, 1.16712, 1.18191};
 
double y1[3];
 
double x1[3] = {0.115, 0.130, 0.164};
 
double s_rr = 0, pr_rr, pr_onrr, s_fn;
 
// Линейная интерполяция
 
for (int i = 0; i < 3; i++){
 
for (int j = 1; j < n; j++)
 
if (x1[i] > x[j - 1] && x1[i] < x[j])
 
index[i] = j;
 
}
 
for (int i = 0; i < 3; i++){
 
y1[i] = ((y[index[i]] - y[index[i] - 1]) / (x[index[i]] - x[index[i] - 1])) * x1[i] + (y[index[i] - 1] - (((y[index[i]] - y[index[i] - 1]) / (x[index[i]] - x[index[i] - 1])) * x[index[i] - 1]));
 
}
 
cout << endl << "значения (Линейная интерполяция):" << endl;
 
for (int i = 0; i < 3; i++)
 
cout << "F(" << x1[i] << ") = " << y1[i] << endl;
 
// Интерполяциооный полином Ньютона
 
y1[0] = 0;
 
y1[1] = 0;
 
y1[2] = 0;
 
for (int l = 0; l < 3; l++)
 
{
 
s_fn = 0;
 
for(int k = 1; k <= n; k++)
 
{
 
pr_onrr = 1;
 
s_rr = 0;
 
for (int i = 0; i <= k; i++)
 
{
 
pr_rr = 1;
 
for (int j = 0; j <= k; j++)
 
if(i != j)
 
pr_rr *= x[i] - x[j];
 
s_rr += y[i] / pr_rr;
 
}
 
for (int i = 0; i < k; i++)
 
{
 
pr_onrr *= x1[l] - x[i];
 
}
 
s_fn += pr_onrr * s_rr;
 
}
 
y1[l] += y[0] + s_fn;
 
}
 
cout << endl << "значения (полином Ньютона):" << endl;
 
for (int i = 0; i < 3; i++)
 
cout << "F(" << x1[i] << ") = " << y1[i] << endl;
 
y1[0] = 0;
 
y1[1] = 0;
 
y1[2] = 0;
 
// Сплайн интерполяция
 
double h[6], l[6], D[5], A[5], c[6], b[6], a[6], d[6];
 
for (int k = 1; k < n; k++){
 
h[k] = x[k] - x[k - 1];
 
l[k] = (y[k] - y[k - 1]) / h[k];
 
}
 
D[1] = (-h[2]) / (2 * (h[1] + h[2]));
 
A[1] = (3*(l[2] - l[1])) / (2 * (h[1] + h[2]));
 
for (int k = 3; k < n; k++){
 
D[k - 1] = h[k] / (2 * h[k - 1] + 2 * h[k] + h[k - 1] * D[k - 2]);
 
A[k - 1] = (3 * l[k] - 3 * l[k - 1] - h[k - 1] * A[k - 2]) / (2 * h[k - 1] + 2 * h[k] + h[k - 1] * D[k - 2]);
 
}
 
c[5] = 0;
 
c[0] = 0;
 
for (int k = n - 1; k >= 2; k--)
 
c[k - 1] = D[k - 1] * c[k] + A[k - 1];
 
x1[0] = 0.702;
 
x1[1] = 0.512;
 
x1[2] = 0.608;
 
for (int k = 1; k < n; k++){
 
b[k] = l[k] + (2/3 * c[k] * h[k]) + (1/3 * h[k] * c[k - 1]);
 
}
 
for (int k = 1; k < n; k++)
 
d[k] = (c[k] - c[k - 1]) / 3 * h[k];
 
for (int i = 0; i < 3; i++)
 
{
 
y1[i] = y[index[i]] + b[index[i]] * (x1[i] - x[index[i]]) + c[index[i]] * pow((x1[i] - x[index[i]]), 2) + d[index[i]] * pow((x1[i] - x[index[i]]), 3);
 
}
 
cout << endl << "значения (сплайн интерполяция):" << endl;
 
for (int i = 0; i < 3; i++)
 
cout << "F(" << x1[i] << ") = " << y1[i] << endl;
 
return 0;
}