00001 #include <SyFi.h> 00002 00003 using namespace GiNaC; 00004 00005 int main() { 00006 00007 // example that check scaling 00008 00009 numeric h(1,100); // 1.0/100 00010 numeric a(1,2); // 1.0/2 00011 ex p0 = lst(a,a,a); 00012 ex p1 = lst(a+h,a,a); 00013 00014 Line line(p0,p1); 00015 00016 00017 ex f = 1; 00018 ex intf = line.integrate(f); 00019 cout <<"intf "<<intf<<endl; 00020 EQUAL_OR_DIE(intf, "1/100"); 00021 00022 00023 00024 return 0; 00025 00026 } 00027