private void rButton3_Click(object sender, EventArgs e)
{
tabControl1.TabPages[1].Enabled = false;
tabControl1.TabPages[0].Enabled = false;
if (sRun)
{
sRun = false;
super.Abort();
rButton3.Text = "запуск";
tabControl1.TabPages[1].Enabled = true;
tabControl1.TabPages[0].Enabled = true;
return;
}
sRun = true;
rButton3.Text = "cтоп";
int selected = comboBox1.SelectedIndex;
super = new Thread(() =>
{
double u0 = Convert.ToDouble(u0TextBox3.Text.Replace('.', ','));
double u01 = Convert.ToDouble(u01TextBox3.Text.Replace('.', ','));
double h = Convert.ToDouble(hTextBox3.Text.Replace('.', ','));
int n = Convert.ToInt32(nTextBox3.Text.Replace('.', ','));
double eps = Convert.ToDouble(epsTextBox3.Text.Replace('.', ','));
double rb = Convert.ToDouble(rbTextBox3.Text.Replace('.', ','));
double a = Convert.ToDouble(aTextBox3.Text.Replace('.', ','));
double b = Convert.ToDouble(bTextBox3.Text.Replace('.', ','));
bool ctrl = !checkBox3.Checked;
SecondMethod m = new SecondMethod(
(x, u1, u2) => (u2),
(x, u1, u2) => (-a * Math.Sqrt(u2 * u2 + 1) - b),
0, u0, u01, h, eps, ctrl);
mainChart.Invoke(new Action(() =>
{
mainChart.Series["Численное решение"].Points.Clear();
mainChart.Series["Точное решение"].Points.Clear();
chart1.Series["h"].Points.Clear();
}));
++nums;
info = new DotForm();
info.label1.Text = "Запуск номер " + nums + "; Метод 3";
minDot = double.MaxValue;
maxDot = double.MinValue;
double minX = double.MaxValue;
double minStep = double.MaxValue;
double mns = 0;
double maxStep = double.MinValue;
double mxs = 0;
MDot p = null;
int count = 0;
double maxOLP = 0;
foreach (var i in Enumerable.Range(0, n))
{
double step = m.Step;
p = m.nextStep(out double contr, out double olp);
if (Math.Abs(p.U1) < 1e-8)
{
p.U2 = 0;
}
if (Math.Abs(p.U1) > 10e+20)
{
break;
}
if (Math.Abs(p.X) < 1e-8)
{
p.X = 0;
}
if (Math.Abs(p.X) > 10e+20)
{
break;
}
Dot add = new Dot();
if (selected == 0)
{
add.X = p.X;
add.Y = p.U1;
}
if (selected == 1)
{
add.X = p.X;
add.Y = p.U2;
}
if (selected == 2)
{
add.X = p.U1;
add.Y = p.U2;
}
mainChart.Invoke(new Action(() => {
mainChart.Series["Численное решение"].Points.AddXY(add.X, add.Y);
chart1.Series["h"].Points.AddXY(i, step);
}));
Console.WriteLine(p.X + " " + p.U1);
mainChart.Invoke(new Action(() =>
{
info.dataGridView1.Rows.Add(i + "", add.Y, contr, olp, step, m.C1, m.C2, "-", "-");
}));
if (minDot > add.Y)
{
minDot = add.Y;
}
if (maxDot < add.Y)
{
maxDot = add.Y;
}
if (minX > add.X)
{
minX = add.X;
}
if (p.X > rb)
{
break;
}
if (maxOLP < Math.Abs(olp))
{
maxOLP = Math.Abs(olp);
}
if (step > maxStep)
{
maxStep = step;
mxs = p.X;
}
if (step < minStep)
{
minStep = step;
mns = p.X;
}
count++;
}
mainChart.Invoke(new Action(() =>
{
correctAxis(rb, minDot - 0.01, maxDot + 0.01);
info.param.Text = $"n = : {count}, \nb-Xn: {rb - p.X}, \nmaxOLP: {maxOLP}, \nC1: {m.C1}, \nC2: {m.C2}, \nmax Hi: {maxStep} -> x: {mxs}, \nmin Hi: {minStep} -> x: {mns}";
info.Show();
sRun = false;
rButton3.Text = "запуск";
tabControl1.TabPages[1].Enabled = true;
tabControl1.TabPages[0].Enabled = true;
if (selected == 0)
{
mainChart.ChartAreas[0].AxisX.Title = "x";
mainChart.ChartAreas[0].AxisY.Title = "u1";
}
if (selected == 1)
{
mainChart.ChartAreas[0].AxisX.Title = "x";
mainChart.ChartAreas[0].AxisY.Title = "u2";
}
if (selected == 2)
{
mainChart.ChartAreas[0].AxisX.Title = "u1";
mainChart.ChartAreas[0].AxisY.Title = "u2";
correctAxis(rb, minDot - 0.01, maxDot + 0.01, minX - 0.01);
}
}));
});
super.Start();
}
private void rButton1_Click(object sender, EventArgs e)
{
++chartCount;
string chartName = "Численное решение №" + chartCount;
var settings = mainChart.Series.Add(chartName);
settings.ChartType = SeriesChartType.Line;
if (sRun)
{
sRun = false;
super.Abort();
rButton1.Text = "запуск";
return;
}
sRun = true;
rButton1.Text = "cтоп";
mainChart.ChartAreas[0].AxisX.Title = "x";
mainChart.ChartAreas[0].AxisY.Title = "u";
super = new Thread(() =>
{
double l = Convert.ToDouble(lTextBox1.Text.Replace('.', ','));
double h = Convert.ToDouble(hTextBox1.Text.Replace('.', ','));
int n = Convert.ToInt32(nTextBox1.Text.Replace('.', ','));
double eps = Convert.ToDouble(epsTextBox1.Text.Replace('.', ','));
double el = Convert.ToDouble(elTextBox1.Text.Replace('.', ','));
double pp = Convert.ToDouble(pTextBox1.Text.Replace('.', ','));
bool ctrl = !checkBox1.Checked;
Method m = new Method(
(x, u1, u2) => ((1 / l) + -x * (1 / (l * l))) * pp * l * l / (el),
(x, u1, u2) => u1,
0, 0, 0, h, eps, ctrl);
mainChart.Invoke(new Action(() => {
mainChart.Series[chartName].Points.Clear();
//chart1.Series["h"].Points.Clear();
}));
++nums;
info = new DotForm();
info.label1.Text = "Запуск номер " + nums + ";";
minDot = double.MaxValue;
maxDot = double.MinValue;
double minStep = double.MaxValue;
double mns = 0;
double maxStep = double.MinValue;
double mxs = 0;
int count = 0;
MDot p = null;
double maxOLP = 0;
double sum = 0;
MDot prev = null;
foreach (var i in Enumerable.Range(0, n))
{
double step = m.Step;
p = m.nextStep(out double contr, out double olp, out double len);
if (Math.Abs(p.U2) < 1e-20)
{
p.U2 = 0;
}
if (Math.Abs(p.U2) > 10e+20)
{
break;
}
if (Math.Abs(p.X) < 1e-20)
{
p.X = 0;
}
if (Math.Abs(p.X) > 10e+20)
{
break;
}
sum += len;
bool needToBreak = false;
if (checkBox2.Checked)
{
if (sum > l)
{
var diff = sum - l;
var proc = diff / len;
var cStep = step * (proc);
var up = (p.U2 - prev.U2) * (proc);
p.X -= cStep;
p.U2 -= up;
needToBreak = true;;
}
}
mainChart.Invoke(new Action(() =>
{
mainChart.Series[chartName].Points.AddXY(p.X, p.U2);
//chart1.Series["h"].Points.AddXY(i, step);
}));
if (!checkBox2.Checked)
{
if (p.X > 5)
{
break;
}
}
else if (needToBreak)
{
break;
}
Console.WriteLine(p.X + " " + p.U2);
mainChart.Invoke(new Action(() =>
{
info.dataGridView1.Rows.Add(i + "", p.X, p.U2, contr, olp, step, m.C1, m.C2);
}));
if (minDot > p.U2)
{
minDot = p.U2;
}
if (maxDot < p.U2)
{
maxDot = p.U2;
}
if (maxOLP < Math.Abs(olp))
{
maxOLP = Math.Abs(olp);
}
if (step > maxStep)
{
maxStep = step;
mxs = p.X;
}
if (step < minStep)
{
minStep = step;
mns = p.X;
}
count++;
prev = p;
}
mainChart.Invoke(new Action(() =>
{
correctAxis(10, minDot - 0.01, maxDot + 0.01);
info.param.Text = $"Номер шага = : {count},\nМаксимальная оценка локальной погрешности (ОЛП): {maxOLP}, \nКоличество делений шага: {m.C1}, \nКоличество удвоений шага: {m.C2}, \nМаксимальный шаг: {maxStep}, достигнут в точке: {mxs}, \nМинимальный шаг: {minStep} , достигнут в точке: {mns}";
info.Show();
sRun = false;
rButton1.Text = "запуск";
}));
});
super.Start();
}
