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
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.Pow(x, 3)) / (3 + x);
}
public static double primitive(double x)
{
return 9 * x - (3 * x * x) / 2 + (x * x * x) / 3 - 27 * Math.Log(3 + x, Math.E) - 47 / 6 + 27 * Math.Log(4, Math.E);
}
static void Main(string[] args)
{
Integral ret = new Integral(1, 2, 72);
Console.WriteLine("Метод Симпсона\nЗначение интеграла:\t");
Console.WriteLine("{0:f4}\n", ret.summ(Program.functions));
Console.WriteLine("Проверка. Вычислим первообразную в пределах интегрирования\nПервообразная:\t");
Console.WriteLine("{0:f4}\n", Program.primitive(2) - Program.primitive(1));
}
}
}