#include "property_check.h"
#include "smooth_sequence.h"

void check_RC_expansion(double epsilon) {
    double e = std::exp(1);
    double R = e*(1 + epsilon);
    Sequence seq = best_smooth_sequence(R, 0.00000000001);
    std::cout << (e - seq.consistency)/std::sqrt(epsilon) - 2 * std::sqrt(2) * e / (e-1)<< "\n";
}