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
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
namespace Trapezoidal_rule {
class Program {
// значение подынтегральной функции в точке x
public static double getFunction(double x) {
double retValue = (Math.Pow(Math.E, x) * (1.0 + Math.Sin(x))) / (1.0 + Math.Cos(x));
return retValue;
}
// значение первообразной функции в точке x
public static double getAntiderivative(double x) {
double retValue = Math.Pow(Math.E, x) * Math.Tan(x / 2.0);
return retValue;
}
// вычисление интеграла
public static double calcIntegral(double a, double b, int n) {
double retValue = 0.0,
h = (b - a) / n; // вычисляем высоту трапеции
for (double x = a; x <= b; x += h)
retValue += h * (getFunction(x) + getFunction(x + h)) / 2;
return retValue;
}
static void Main(string[] args) {
double leftRange = 0.0, // левая граница
rightRange = 1.5, // правая граница
midValue = (leftRange + rightRange) / 2; // середина отрезка
int n = 700; // возьмем большое разбиение для большей точности
double J1 = 0.0,
J2 = 0.0;
//метод трапеций
J1 = calcIntegral(leftRange, rightRange, n);
//считаем интеграл по методу Ньютона-Лейбница (Проверка)
J2 = getAntiderivative(rightRange) - getAntiderivative(leftRange);
Console.WriteLine("Значение интеграла: {0:f5}", J1);
Console.WriteLine("Проверка: F(b) - F(a) = {0:f5}", J2);
Console.WriteLine("Относительная точность приблеженного значения интеграла: {0}%",
Math.Abs(J1 - J2) / J2 * 100);
Console.ReadKey();
}
}
}