00001 #include <SyFi.h> 00002 00003 using namespace GiNaC; 00004 00005 int main() { 00006 00007 ex f = x*y*z; 00008 ReferenceBox box; 00009 00010 ex repr = box.repr(); 00011 cout <<"b.repr "<<repr<<endl; 00012 00013 ex intf = box.integrate(f); 00014 cout <<"intf "<<intf<<endl; 00015 00016 00017 00018 ex p0 = lst(-1.0,-1.0,-1.0); 00019 ex p1 = lst( 1.0, 1.0, 1.0); 00020 00021 Box box2(p0,p1); 00022 00023 ex repr2 = box2.repr(); 00024 cout <<"b2.repr "<<repr<<endl; 00025 00026 ex intf2 = box2.integrate(f); 00027 cout <<"intf2 "<<intf2<<endl; 00028 00029 return 0; 00030 } 00031