Lab10/Lab_10/Source.c - master - aytsan_ns/Num_meth

Lab10/Lab_10/Source.c

82 строки2 KB

aytsan-ns

Изменила код в 10 работе так, чтобы не перевычислять старые значения функции, а добавлять лишь новые, также объединила две функции в одну с проверкой на число итераций (если дано не нулевое, значит, выполнение с фиксированным n, иначе вычисление n в ходе вычисления интеграла)

22 сен 2024, 19:50

22 сен 2024, 19:50a9cc7ea

#define _CRT_SECURE_NO_WARNINGS
#include <math.h>
#include <stdlib.h>
#include <stdio.h>
#include <time.h>
#define PI 3.14159265358979323846
#define I 8.00


double function(double x) {
	return (pow(x, 5) - 9.2 * pow(x, 3) + 5.5 * x * x - 7 * x);
}

double ComputeIntegralTrapez(double a, double b, double epsilon, int* n) {
    if (*n != 0) {
        double h = (b - a) / *n;
        double sum = (function(a) + function(b)) / 2;
        for (int i = 1; i < *n; i++) {
            sum += function(a + i * h);
        }
        return h * sum;
    }

    double I1, I2;
    *n = 1;
    double h = (b - a) / (*n);
    double sum = (function(a) + function(b)) / 2;

    I1 = h * sum;

    do {
        *n *= 2;
        h = (b - a) / (*n);

        for (int i = 1; i < *n; i+=2) {
            sum += function(a + i * h);
        }

        I2 = h * sum;

        if (fabs(I2 - I1) / 3.0 < epsilon) break;

        I1 = I2;
    } while (1);

    return I2;
}

void Task3_4() {

    double a = -2, b = 3;
    int n;
    FILE* error = fopen("error.txt", "w");
    FILE* error_n = fopen("error_n.txt", "w");

    for (double epsilon = 0.1; epsilon >= 1e-10; epsilon /= 10) {
        n = 0;
        double actual_error = fabs(I - ComputeIntegralTrapez(a, b, epsilon, &n));
        fprintf(error, "%e\n", actual_error);
        fprintf(error_n, "%i\n", n);
    }

    fclose(error);
}

void Task5() {
    double a = -2, b = 3;
    FILE* error_h = fopen("error_h.txt", "w");

    for (int n = 2; n <= 1024; n *= 2) {
        double h = (b - a) / n;
        double integral = ComputeIntegralTrapez(a, b, 0, &n);
        fprintf(error_h, "%.10lf %.10lf\n", h, fabs(I - integral));
    }

    fclose(error_h);
}

void main() {
    Task3_4();
    Task5();
}