#include "Solver.h"
#include "Header.h"
#include <Eigen/Dense>
#define A 5.
#define B 4.
#define C 0.
#define D 1.

using namespace Eigen;
int main() {
    std::cout << A << "u''" << " + " << B << "u'" << "+ " << C << "u + " << D << " = 0" << std::endl;
    Solver slv(A, B, C, D, 2, 0, 10);
    std::cout << "Linear element:";
    slv.Execute_Linear(0, 5);
    std::cout << "\nCubic element:";
    slv.Execute_Cubic(0, 5);
    return 0;
}