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