00001 #include <SyFi.h>
00002
00003 int main() {
00004 int order = 2;
00005 int nsd = 2;
00006
00007 ex p = pol(order,nsd, "a");
00008 cout <<"polynom p = "<<p<<endl;
00009 EQUAL_OR_DIE(p,"y^2*a5+x^2*a3+a2*y+y*x*a4+a0+a1*x");
00010
00011 ex dpdx = diff(p,x);
00012 cout <<"dpdx = "<<dpdx<<endl;
00013 EQUAL_OR_DIE(dpdx, "y*a4+a1+2*x*a3");
00014
00015 Triangle triangle(lst(0,0), lst(1,0), lst(0,1));
00016 ex intp = triangle.integrate(p);
00017 cout <<"integral of p over reference triangle = "<<intp<<endl;
00018 EQUAL_OR_DIE(intp, "1/6*a2+1/6*a1+1/12*a5+1/2*a0+1/24*a4+1/12*a3");
00019
00020
00021 }
