simple2.cpp

Go to the documentation of this file.

#include <SyFi.h>

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

int main() {

    initSyFi(1); 

    // initialization (f = x^2 + y^2) 
    ex f = x*x + y*y; 
    cout <<"f "<<f<<endl; 

    // differentiation (dfdx = df/dx = 2x) 
    ex dfdx = f.diff(x,1);  
    cout <<"dfdx "<<dfdx<<endl; 

    // integration (intf = 1/3+y^2, integral of f(x,y) from x=0 to x=1) 
    ex intf = integral(x,0,1,f); 
    cout <<"intf "<<intf<<endl; 
    intf = eval_integ(intf); 
    cout <<"intf "<<intf<<endl; 

    return 0; 
}