6lab

#define _USE_MATH_DEFINES
#include <iostream>
#include <cmath>

//точное решение задачи
double source_f(double x){
	return pow(M_E, 1 - cos(x));
}

int main(){
	double y_0 = 1;
	double x_0 = 0;
	printf("%8s|%15s|%15s|%15s\n","n","delta_left","delta_mid","delta_right");
	for (int n = 4; n < 1e5; n *= 2){
		double h = M_PI_2 / n;
		double* y - new double[n + 1];
		y[0] = y_0;
		for (int i = 1; i < n; i++){
			y[i] = y[i-1] + h * y[i-1] * sin(x_0 + i * h);
		}
		printf("%8d|%15e|%15e|%15e\n", n, fabs(source_f(M_PI_2) - y[n]));
	}
}