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 ex p2 = lst(a,a+h,a); 00014 00015 Triangle triangle(p0,p1,p2); 00016 00017 00018 ex f = 1; 00019 ex intf = triangle.integrate(f); 00020 cout <<"intf "<<intf<<endl; 00021 EQUAL_OR_DIE(intf, "1/20000"); 00022 00023 return 0; 00024 } 00025