integral_ex1.cpp

Go to the documentation of this file.

#include <SyFi.h>
#include <ginac/ginac.h>

using namespace std; 
using namespace GiNaC; 
using namespace SyFi; 

void check_integral(ex& f){
}

int main() {

    initSyFi(1); 

    ReferenceLine line; 
    Lagrange fe; 
    fe.set_polygon(line);
    fe.set_order(3); 
    fe.compute_basis_functions(); 
    ex integrand = fe.N(0)*fe.N(1);  
    cout <<integrand<<endl; 
    ex integral1 = line.integrate(integrand); 
    cout <<"Integral " <<integral1<<endl; 
    EQUAL_OR_DIE(integral1, "33/560"); 


}