public void SetParameters(FunctionVector f, int N, double a, double b, Vector y0)
{
this.f = f;
h = (b - a) / N;
x = new double[N + 1];
for (int i = 0; i <= N; ++i)
{
x[i] = a + i * h;
}
y = new Vector[N + 1];
for (int i = 0; i < N + 1; ++i)
{
y[i] = new Vector(f.Length);
}
for (int i = 1; i <= f.Length; ++i)
{
y[0][i] = y0[i];
}
}
private void button1_Click(object sender, EventArgs e)
{
chart1.Series[0].Points.Clear();
chart1.Series[1].Points.Clear();
chart1.Series[2].Points.Clear();
chart1.Series[3].Points.Clear();
dataGridView1.Rows.Clear();
a = Convert.ToDouble(textBox1.Text);
b = Convert.ToDouble(textBox2.Text);
switch (comboBox1.SelectedIndex)
{
case 0:
{
y1Accurate y1 = new y1Accurate();
y2Accurate y2 = new y2Accurate();
systemOfAccurateFunctions = new FunctionVector(2);
systemOfAccurateFunctions[1] = y1;
systemOfAccurateFunctions[2] = y2;
for (double i = a; i <= b; i += 0.01)
{
chart1.Series[0].Points.AddXY(i, y1.Value(i));
chart1.Series[1].Points.AddXY(i, y2.Value(i));
}
systemOfFunctions = new FunctionVector(2);
systemOfFunctions[1] = new f1();
systemOfFunctions[2] = new f2();
string[] y0String = textBox3.Text.Split(' ');
y0 = new Vector(2);
y0[1] = Convert.ToDouble(y0String[0]);
y0[2] = Convert.ToDouble(y0String[1]);
break;
}
case 1:
{
Custom1y1Accurate y1 = new Custom1y1Accurate();
systemOfAccurateFunctions = new FunctionVector(1);
systemOfAccurateFunctions[1] = y1;
for (double i = a; i <= b; i += 0.01)
{
chart1.Series[0].Points.AddXY(i, y1.Value(i));
}
systemOfFunctions = new FunctionVector(1);
systemOfFunctions[1] = new Custom1f1();
string[] y0String = textBox3.Text.Split(' ');
y0 = new Vector(1);
y0[1] = Convert.ToDouble(y0String[0]);
break;
}
case 2:
{
Custom2y1Accurate y1 = new Custom2y1Accurate();
Custom2y2Accurate y2 = new Custom2y2Accurate();
systemOfAccurateFunctions = new FunctionVector(2);
systemOfAccurateFunctions[1] = y1;
systemOfAccurateFunctions[2] = y2;
for (double i = a; i <= b; i += 0.01)
{
chart1.Series[0].Points.AddXY(i, y1.Value(i));
chart1.Series[1].Points.AddXY(i, y2.Value(i));
}
systemOfFunctions = new FunctionVector(2);
systemOfFunctions[1] = new Custom2f1();
systemOfFunctions[2] = new Custom2f2();
string[] y0String = textBox3.Text.Split(' ');
y0 = new Vector(2);
y0[1] = Convert.ToDouble(y0String[0]);
y0[2] = Convert.ToDouble(y0String[1]);
break;
}
}
N = Convert.ToInt32(numericUpDown1.Value);
RungeKuttaMethodThirdOrder RKMethod = new RungeKuttaMethodThirdOrder();
RKMethod.SetParameters(systemOfFunctions, N, a, b, y0);
RKMethod.Solve();
for (int i = 0; i < RKMethod.x.Length; ++i)
{
for (int j = 1; j <= systemOfAccurateFunctions.Length; ++j)
{
chart1.Series[j + 1].Points.AddXY(RKMethod.x[i], RKMethod.y[i][j]);
}
}
double err, maxErr;
for (int i = 10; i <= 1000000; i *= 10)
{
RKMethod.SetParameters(systemOfFunctions, i, a, b, y0);
RKMethod.Solve();
maxErr = 0;
for (int j = 1; j <= systemOfAccurateFunctions.Length; ++j)
{
err = Math.Abs(systemOfAccurateFunctions[j].Value(RKMethod.x[i]) - RKMethod.y[i][j]);
if (err > maxErr)
{
maxErr = err;
}
}
dataGridView1.Rows.Add(i, maxErr);
}
}
