00001 #include <SyFi.h>
00002
00003 using namespace GiNaC;
00004
00005 int main() {
00006 int order = 3;
00007
00008 int nsd = 1;
00009 cout <<"third order homogenous polynomial in 1D"<<endl;
00010 ex polh = homogenous_pol(order, nsd, "a");
00011 cout <<"polh "<<polh<<endl;
00012
00013
00014 nsd = 2;
00015 cout <<"third order homogenous polynomial in 2D"<<endl;
00016 polh = homogenous_pol(order, nsd, "a");
00017 cout <<"polh "<<polh<<endl;
00018
00019 nsd = 3;
00020 cout <<"third order homogenous polynomial in 3D"<<endl;
00021 polh = homogenous_pol(order, nsd, "a");
00022 cout <<"polh "<<polh<<endl;
00023
00024 }