legendre.cpp

Go to the documentation of this file.
00001 #include <SyFi.h>
00002 
00003 using namespace GiNaC; 
00004 
00005 int main() {
00006   int order = 2; 
00007   int nsd   = 2; 
00008 
00009   ex polynom_space = legendre(order,nsd, "a"); 
00010   cout <<"polynom_space "<<polynom_space<<endl; 
00011 
00012   ex p = polynom_space.op(0); 
00013   cout <<"polynom p = "<<p<<endl; 
00014 
00015   ex dpdx = diff(p,x);   
00016   cout <<"dpdx = "<<dpdx<<endl; 
00017 
00018   ex p0 = lst(-1,-1);
00019   ex p1 = lst(1,1);
00020 
00021   Rectangle rectangle(p0,p1) ;
00022   ex basis = polynom_space.op(2); 
00023   for (int i=0; i< basis.nops(); i++) {
00024       cout <<"i "<<i<<endl; 
00025       ex integrand = p*basis.op(i); 
00026       ex ai = rectangle.integrate(integrand); 
00027       cout <<"ai "<<ai<<endl; 
00028   }
00029 }
00030 
00031 
00032 

Generated on Tue Jun 13 13:18:38 2006 for SyFi by  doxygen 1.4.4