とらりもんHOME  Index  Search  Changes  Login

GSL

GSL (Gnu Scietific Library) http://www.gnu.org/software/gsl/manual/html_node/index.html#Top

Ubuntu Linuxで...

sudo apt-get -yV install gsl-bin libgsl0-dbg libgsl0-dev libgsl0ldbl
gcc test.c -lgsl -lgslcblas
// ルジャンドル陪関数
// gcc test.c -lgsl -lgslcblas -lm
// ./a.out 6 2 >t
// gnuplot> splot "t" using ($4*sin($1)*cos($2)):($4*sin($1)*sin($2)):($4*cos($1)):($3) with pm3d

#include <stdio.h>
#include <math.h>
#include <gsl/gsl_sf_legendre.h>

int main (int argc, char** argv)
{int l, m;
 double t, p, x, y;
 l=atoi(argv[1]);
 m=atoi(argv[2]);
 for (t=0; t<M_PI; t+=0.01){
    x=gsl_sf_legendre_sphPlm(l,m,cos(t));
//     x=gsl_sf_legendre_Plm(l,m,cos(t));
    for (p=-M_PI; p<M_PI; p+=0.01){
       y=x*cos(m*p);
       printf("%lf  %lf  %lf  %lf\n", t, p, y, sqrt(y*y));
    }
  printf("\n");
 }
 return 0;
}
Last modified:2014/01/21 20:49:04
Keyword(s):
References:[とらりもんHOME]