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(); }