Алгоритм и программная реализация одномерной линейной интерполяции сеточной функции кубическими сплайнами

NewLine; } chart1.Series[0].Points.Clear(); chart1.Series[1].Points.Clear(); for (int i = 0; i < x.Count; i++) { chart1.Series[1].Points.AddXY(x[i], y[i]); } } button2.Enabled = true; double xMin = x.Min(); double xMax = x.Max(); double dx = (xMax - xMin) / 10; chart1.ChartAreas[0].AxisX.Minimum = xMin - dx; chart1.ChartAreas[0].AxisX.Maximum = xMax + dx; } catch { MessageBox.Show("Неверный файл или его содержимое не соответствует таблице:\nX\tY\nx1\ty1\nx2\ty2\n......", "Ошибка открытия файла", MessageBoxButtons.OK, MessageBoxIcon.Error); } } } void Calc() { int n = x.Count - 1; double[] k = new double[n + 1]; double[] c = new double[n + 2]; k[1] = 0; c[1] = 0; int i, j, m; double a, b, r, q, d; for (i = 2; i <= n; i++) { j = i - 1; m = j - 1; a = x[i] - x[j]; b = x[j] - x[m]; r = 2 * (a + b) - b * c[j]; c[i] = a / r; k[i] = (3.0 * ((y[i] - y[j]) / a - (y[j] - y[m]) / b) - b * k[j]) / r; } c[n + 1] = k[n]; for (i = n; i >= 0; i--) { c[i] = k[i] - c[i] * c[i + 1]; } double xi, yi; int nx = 100; chart1.Series[0].Points.Clear(); for (i = 1; i <= n; i++) { double h = (x[i] - x[i - 1]) / (nx - 1); for (int l = 0; l < nx; l++) { xi = x[i - 1] + l * h; q = x[i] - x[i - 1]; r = xi - x[i - 1]; b = (y[i] - y[i - 1]) / q - (c[i + 1] + 2 * c[i]) * q / 3.0; d = (c[i + 1] - c[i]) / q * r; yi = y[i - 1] + r * (b + r * (c[i] + d / 3.0)); chart1.Series[0].Points.AddXY(xi, yi); } } chart1.Legends[0].Enabled = true; } private void button2_Click(object sender, EventArgs e) { Calc(); } private void button3_Click(object sender, EventArgs e) { Close(); } }}