using System using System Collections Generic using System Linq using

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
namespace Метод_Симпсона {
class Integral {
public double a, // нижняя граница интегрирования
b; // верхняя граница интегрирования
public int n; // количество частей разбиений
public Integral(double a, double b, int n) {
this.a = a;
this.b = b;
this.n = n;
}
public delegate double D(double x);
public double summ(D F) {
double h = 0,
s = 0,
x = 0;
h = ((b - a) / (3 * n));
s = F(a) + F(b);
int m = 3 * n - 1;
for (int i = 1; i <= m; i++) {
x = a + h * i;
if (i % 3 == 0)
s = s + 2 * F(x);
else
s = s + 3 * F(x);
}
s = s * 3 / 8 * h;
return s;
}
}
class Program {
public static double functions(double x) {
double retValue = Math.Pow(Math.Tan(x), 2) + Math.Pow(Math.Cos(x) / Math.Sin(x), 2);
return retValue;
}
public static double primitive(double x) {
double retValue = Math.Tan(x) - 2.0 * x - Math.Cos(x) / Math.Sin(x);
return retValue;
}
static void Main(string[] args) {
Integral ret = new Integral(Math.PI / 6, Math.PI / 3, 100);
Console.WriteLine("Метод Симпсона\nЗначение интеграла:\t");
Console.WriteLine("{0:f4}\n", ret.summ(Program.functions));
Console.WriteLine("Проверка. Вычислим первообразную в пределах интегрирования\nПервообразная:\t");
Console.WriteLine("{0:f4}\n", Program.primitive(Math.PI / 3) - Program.primitive(Math.PI / 6));
}
}
}