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
53
54
55
56
57
58
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
namespace cw_12
{
class Integral
{
public double a;// нижняя граница интегрирования
public double 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)
{
return Math.Sqrt(Math.Pow(Math.E, x) - 1);
}
public static double primitive(double x)
{
double retVal = 2 * Math.Sqrt(Math.Pow(Math.E, x) - 1) - 2 * Math.Atan(Math.Sqrt(Math.Pow(Math.E, x) - 1));
return retVal;
}
static void Main(string[] args)
{
Integral ret = new Integral(0, Math.Log(2, Math.E), 104);
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.Log(2, Math.E)) - Program.primitive(0));
}
}
}