summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorbansodanurag2019-05-06 10:50:23 +0530
committerGitHub2019-05-06 10:50:23 +0530
commit5003ab5014dad982927b985fbe67e585a29a5ded (patch)
tree803a73680a90b28092e038b01ecf5e6177c5fc30
parent89ce1e127b022226ccf0e824dbc94fa722ef7c87 (diff)
downloadOpenmodelica-gsl-library-5003ab5014dad982927b985fbe67e585a29a5ded.tar.gz
Openmodelica-gsl-library-5003ab5014dad982927b985fbe67e585a29a5ded.tar.bz2
Openmodelica-gsl-library-5003ab5014dad982927b985fbe67e585a29a5ded.zip
Add files via upload
6/05/2019
Diffstat
-rw-r--r--External_Functions/a.mo4391
1 files changed, 4362 insertions, 29 deletions
diff --git a/External_Functions/a.mo b/External_Functions/a.mo
index 7c67be0..c1b03c1 100644
--- a/External_Functions/a.mo
+++ b/External_Functions/a.mo
@@ -23,29 +23,16 @@ package gsl
end destructor;
end gsl_permutation;
- /* model gsl_complex
- extends ExternalObject;
-
- function constructor
- input Integer N;
- output gsl_complex k;
-
- external "C" annotation(
- Library = "gsl",
- Library = "gslcblas",
- Include = "#include<gsl/gsl_complex.h>");
- end constructor;
+ record gsl_sf_result
+ Real val;
+ Real err;
+ end gsl_sf_result;
- function destructor
- input gsl_complex k;
-
- external "C" annotation(
- Library = "gsl",
- Library = "gslcblas",
- Include = "#include<gsl/gsl_complex.h>");
- end destructor;
- end gsl_complex;
- */
+ record gsl_sf_result_e10
+ Real val;
+ Real err;
+ Integer e10;
+ end gsl_sf_result_e10;
end data_types;
package mathematical
@@ -382,9 +369,9 @@ package gsl
package COMPLEX
function gsl_complex_rect
- input Real x = 4;
- input Real y = 5;
- output Complex z;
+ input Real x;
+ input Real y;
+ output Real z[2];
external "C" z = gsl_complex_rect(x, y) annotation(
Library = "gsl",
@@ -431,11 +418,2543 @@ package gsl
/*gsl_ran_shuffle(r, p.data, N, sizeof(size_t))*/
end Permutation;
+ package Special_functions
+ package chap_7_5_1
+ function gsl_sf_bessel_J0
+ // The function takes a real number x as the input and calculates the zeroth order bessel function
+ input Real x;
+ //input to the bessel function
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_J0(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_J0;
+
+ function gsl_sf_bessel_J0_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_J0_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_J0_e;
+
+ function gsl_sf_bessel_J1
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_J1(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_J1;
+
+ function gsl_sf_bessel_J1_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_J1_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_J1_e;
+
+ function gsl_sf_bessel_Jn
+ input Real x;
+ input Integer n;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Jn(n, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Jn;
+
+ function gsl_sf_bessel_Jn_e
+ input Real x;
+ input Integer n;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Jn_e(n, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Jn_e;
+
+ function gsl_sf_bessel_Jn_array
+ input Integer nmin;
+ input Integer nmax;
+ input Real x;
+ output Real result_array[nmax - nmin];
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Jn_array(nmin, nmax, x, result_array) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Jn_array;
+ end chap_7_5_1;
+
+ package chap_7_5_2
+ function gsl_sf_bessel_Y0
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Y0(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Y0;
+
+ function gsl_sf_bessel_Y0_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Y0_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Y0_e;
+
+ function gsl_sf_bessel_Y1
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Y1(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Y1;
+
+ function gsl_sf_bessel_Y1_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Y1_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Y1_e;
+
+ function gsl_sf_bessel_Yn
+ input Real x;
+ input Integer n;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Yn(n, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Yn;
+
+ function gsl_sf_bessel_Yn_e
+ input Real x;
+ input Integer n;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Yn_e(n, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Yn_e;
+
+ function gsl_sf_bessel_Yn_array
+ input Integer nmin;
+ input Integer nmax;
+ input Real x;
+ output Real result_array[nmax - nmin];
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Yn_array(nmin, nmax, x, result_array) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Yn_array;
+ end chap_7_5_2;
+
+ package chap_7_5_3
+ function gsl_sf_bessel_I0
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_I0(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_I0;
+
+ function gsl_sf_bessel_I0_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_I0_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_I0_e;
+
+ function gsl_sf_bessel_I1
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_I1(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_I1;
+
+ function gsl_sf_bessel_I1_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_I1_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_I1_e;
+
+ function gsl_sf_bessel_In
+ input Real x;
+ input Integer n;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_In(n, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_In;
+
+ function gsl_sf_bessel_In_e
+ input Real x;
+ input Integer n;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_In_e(n, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_In_e;
+
+ function gsl_sf_bessel_In_array
+ input Integer nmin;
+ input Integer nmax;
+ input Real x;
+ output Real result_array[nmax - nmin];
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_In_array(nmin, nmax, x, result_array) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_In_array;
+
+ function gsl_sf_bessel_I0_scaled
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_I0_scaled(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_I0_scaled;
+
+ function gsl_sf_bessel_I0_scaled_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_I0_scaled_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_I0_scaled_e;
+
+ function gsl_sf_bessel_I1_scaled
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_I1_scaled(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_I1_scaled;
+
+ function gsl_sf_bessel_I1_scaled_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_I1_scaled_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_I1_scaled_e;
+
+ function gsl_sf_bessel_In_scaled
+ input Real x;
+ input Integer n;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_In_scaled(n, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_In_scaled;
+
+ function gsl_sf_bessel_In_scaled_e
+ input Real x;
+ input Integer n;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_In_scaled_e(n, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_In_scaled_e;
+
+ function gsl_sf_bessel_In_scaled_array
+ input Integer nmin;
+ input Integer nmax;
+ input Real x;
+ output Real result_array[nmax - nmin];
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_In_scaled_array(nmin, nmax, x, result_array) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_In_scaled_array;
+ end chap_7_5_3;
+
+ package chap_7_5_4
+ function gsl_sf_bessel_K0
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_K0(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_K0;
+
+ function gsl_sf_bessel_K0_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_K0_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_K0_e;
+
+ function gsl_sf_bessel_K1
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_K1(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_K1;
+
+ function gsl_sf_bessel_K1_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_K1_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_K1_e;
+
+ function gsl_sf_bessel_Kn
+ input Real x;
+ input Integer n;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Kn(n, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Kn;
+
+ function gsl_sf_bessel_Kn_e
+ input Real x;
+ input Integer n;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Kn_e(n, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Kn_e;
+
+ function gsl_sf_bessel_Kn_array
+ input Integer nmin;
+ input Integer nmax;
+ input Real x;
+ output Real result_array[nmax - nmin];
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Kn_array(nmin, nmax, x, result_array) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Kn_array;
+
+ function gsl_sf_bessel_K0_scaled
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_K0_scaled(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_K0_scaled;
+
+ function gsl_sf_bessel_K0_scaled_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_K0_scaled_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_K0_scaled_e;
+
+ function gsl_sf_bessel_K1_scaled
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_K1_scaled(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_K1_scaled;
+
+ function gsl_sf_bessel_K1_scaled_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_K1_scaled_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_K1_scaled_e;
+
+ function gsl_sf_bessel_Kn_scaled
+ input Real x;
+ input Integer n;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Kn_scaled(n, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Kn_scaled;
+
+ function gsl_sf_bessel_Kn_scaled_e
+ input Real x;
+ input Integer n;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Kn_scaled_e(n, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Kn_scaled_e;
+
+ function gsl_sf_bessel_Kn_scaled_array
+ input Integer nmin;
+ input Integer nmax;
+ input Real x;
+ output Real result_array[nmax - nmin];
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Kn_scaled_array(nmin, nmax, x, result_array) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Kn_scaled_array;
+ end chap_7_5_4;
+
+ package chap_7_5_5
+ function gsl_sf_bessel_j0
+ //this function computes the regular spherical bessels function j0(x)=sin(x)/x
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_j0(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_j0;
+
+ function gsl_sf_bessel_j0_e
+ // This function computes spherical bessels zero order function j(x)=sin(x)/x
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_j0_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_j0_e;
+
+ function gsl_sf_bessel_j1
+ //This function computes the 1st order spherical bessels function j1(x)=(sin(x)/x-cos(x))/x
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_j1(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_j1;
+
+ function gsl_sf_bessel_j1_e
+ //This function computes the 1st order spherical bessels function j1(x)=(sin(x)/x-cos(x))/x
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_j1_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_j1_e;
+
+ function gsl_sf_bessel_j2
+ //This function computes the 2nd order regular spherical bessels function j 2 (x) = ((3/x 2 − 1) sin(x) −3 cos(x)/x)/x.
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_j2(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_j2;
+
+ function gsl_sf_bessel_j2_e
+ //This function computes the 2nd order regular spherical bessels function j 2 (x) = ((3/x 2 − 1) sin(x) −3 cos(x)/x)/x.
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_j2_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_j2_e;
+
+ function gsl_sf_bessel_jl
+ input Integer l;
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_jl(l, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_jl;
+
+ function gsl_sf_bessel_jl_e
+ input Integer l;
+ //order of the bessel function
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_jl_e(l, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_jl_e;
+
+ function gsl_sf_bessel_jl_array
+ // input Integer nmin;
+ input Integer lmax;
+ input Real x;
+ output Real result_array[lmax];
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_jl_array(lmax, x, result_array) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_jl_array;
+
+ function gsl_sf_bessel_jl_steed_array
+ // input Integer nmin;
+ input Integer lmax;
+ input Real x;
+ output Real result_array[lmax];
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_jl_steed_array(lmax, x, result_array) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_jl_steed_array;
+ end chap_7_5_5;
+
+ package chap_7_5_6
+ function gsl_sf_bessel_y0
+ //this function computes the regular spherical bessels function y0(x)=sin(x)/x
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_y0(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_y0;
+
+ function gsl_sf_bessel_y0_e
+ // This function computes spherical bessels zero order function j(x)=sin(x)/x
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_y0_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_y0_e;
+
+ function gsl_sf_bessel_y1
+ //This function computes the 1st order spherical bessels function y1(x)=(sin(x)/x-cos(x))/x
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_y1(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_y1;
+
+ function gsl_sf_bessel_y1_e
+ //This function computes the 1st order spherical bessels function y1(x)=(sin(x)/x-cos(x))/x
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_y1_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_y1_e;
+
+ function gsl_sf_bessel_y2
+ //This function computes the 2nd order regular spherical bessels function j 2 (x) = ((3/x 2 − 1) sin(x) −3 cos(x)/x)/x.
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_y2(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_y2;
+
+ function gsl_sf_bessel_y2_e
+ //This function computes the 2nd order regular spherical bessels function j 2 (x) = ((3/x 2 − 1) sin(x) −3 cos(x)/x)/x.
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_y2_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_y2_e;
+
+ function gsl_sf_bessel_yl
+ input Integer l;
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_yl(l, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_yl;
+
+ function gsl_sf_bessel_yl_e
+ input Integer l;
+ //order of the bessel function
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_yl_e(l, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_yl_e;
+
+ function gsl_sf_bessel_yl_array
+ // input Integer nmin;
+ input Integer lmax;
+ input Real x;
+ output Real result_array[lmax];
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_yl_array(lmax, x, result_array) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_yl_array;
+ end chap_7_5_6;
+
+ package chap_7_5_7
+ //regular modified spherical bessel functions
+
+ function gsl_sf_bessel_i0_scaled
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_i0_scaled(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_i0_scaled;
+
+ function gsl_sf_bessel_i0_scaled_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_i0_scaled_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_i0_scaled_e;
+
+ function gsl_sf_bessel_i1_scaled
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_i1_scaled(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_i1_scaled;
+
+ function gsl_sf_bessel_i1_scaled_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_i1_scaled_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_i1_scaled_e;
+
+ function gsl_sf_bessel_i2_scaled
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_i2_scaled(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_i2_scaled;
+
+ function gsl_sf_bessel_i2_scaled_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_i2_scaled_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_i2_scaled_e;
+
+ function gsl_sf_bessel_il_scaled
+ input Real x;
+ input Integer l;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_il_scaled(l, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_il_scaled;
+
+ function gsl_sf_bessel_il_scaled_e
+ input Real x;
+ input Integer n;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_il_scaled_e(n, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_il_scaled_e;
+
+ function gsl_sf_bessel_il_scaled_array
+ input Integer lmax;
+ input Real x;
+ output Real result_array[lmax];
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_il_scaled_array(lmax, x, result_array) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_il_scaled_array;
+ end chap_7_5_7;
+
+ package chap_7_5_8
+ function gsl_sf_bessel_k0_scaled
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_k0_scaled(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_k0_scaled;
+
+ function gsl_sf_bessel_k0_scaled_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_k0_scaled_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_k0_scaled_e;
+
+ function gsl_sf_bessel_k1_scaled
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_k1_scaled(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_k1_scaled;
+
+ function gsl_sf_bessel_k1_scaled_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_k1_scaled_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_k1_scaled_e;
+
+ function gsl_sf_bessel_k2_scaled
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_k2_scaled(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_k2_scaled;
+
+ function gsl_sf_bessel_k2_scaled_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_k2_scaled_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_k2_scaled_e;
+
+ function gsl_sf_bessel_kl_scaled
+ input Real x;
+ input Integer l;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_kl_scaled(l, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_kl_scaled;
+
+ function gsl_sf_bessel_kl_scaled_e
+ input Real x;
+ input Integer n;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_kl_scaled_e(n, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_kl_scaled_e;
+
+ function gsl_sf_bessel_kl_scaled_array
+ input Integer lmax;
+ input Real x;
+ output Real result_array[lmax];
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_kl_scaled_array(lmax, x, result_array) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_kl_scaled_array;
+ end chap_7_5_8;
+
+ package chap_7_5_9
+ function gsl_sf_bessel_Jnu
+ input Real x;
+ input Real nu;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Jnu(nu, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Jnu;
+
+ function gsl_sf_bessel_Jnu_e
+ input Real x;
+ input Real nu;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Jnu_e(nu, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Jnu_e;
+ end chap_7_5_9;
+
+ package chap_7_5_10
+ function gsl_sf_bessel_Ynu
+ input Real x;
+ input Real nu;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Ynu(nu, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Ynu;
+
+ function gsl_sf_bessel_Ynu_e
+ input Real x;
+ input Real nu;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Ynu_e(nu, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Ynu_e;
+ end chap_7_5_10;
+
+ package chap_7_5_11
+ function gsl_sf_bessel_Inu
+ input Real x;
+ input Real nu;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Inu(nu, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Inu;
+
+ function gsl_sf_bessel_Inu_e
+ input Real x;
+ input Real nu;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Inu_e(nu, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Inu_e;
+
+ function gsl_sf_bessel_Inu_scaled
+ input Real x;
+ input Real nu;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Inu_scaled(nu, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Inu_scaled;
+
+ function gsl_sf_bessel_Inu_scaled_e
+ input Real x;
+ input Real nu;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Inu_scaled_e(nu, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Inu_scaled_e;
+ end chap_7_5_11;
+
+ package chap_7_5_12
+ //Irregular modified bessels function
+
+ function gsl_sf_bessel_Knu
+ input Real x;
+ input Real nu;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Knu(nu, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Knu;
+
+ function gsl_sf_bessel_Knu_e
+ input Real x;
+ input Real nu;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Knu_e(nu, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Knu_e;
+
+ //Irregular modified bessels function
+
+ function gsl_sf_bessel_lnKnu
+ input Real x;
+ input Real nu;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_lnKnu(nu, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_lnKnu;
+
+ function gsl_sf_bessel_lnKnu_e
+ input Real x;
+ input Real nu;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_lnKnu_e(nu, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_lnKnu_e;
+
+ function gsl_sf_bessel_Knu_scaled
+ input Real x;
+ input Real nu;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_Knu_scaled(nu, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Knu_scaled;
+
+ function gsl_sf_bessel_Knu_scaled_e
+ input Real x;
+ input Real nu;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_Knu_scaled_e(nu, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_Knu_scaled_e;
+ end chap_7_5_12;
+
+ package chap_7_5_13
+ function gsl_sf_bessel_zero_J0
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_zero_J0(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_zero_J0;
+
+ function gsl_sf_bessel_zero_J0_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_zero_J0_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_zero_J0_e;
+
+ function gsl_sf_bessel_zero_J1
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_zero_J1(x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_zero_J1;
+
+ function gsl_sf_bessel_zero_J1_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_zero_J1_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_zero_J1_e;
+
+ function gsl_sf_bessel_zero_Jnu
+ input Real x;
+ input Real nu;
+ output Real y;
+
+ external "C" y = gsl_sf_bessel_zero_Jnu(nu, x) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_zero_Jnu;
+
+ function gsl_sf_bessel_zero_Jnu_e
+ input Real x;
+ input Real nu;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_bessel_zero_Jnu_e(nu, x, result) annotation(
+ Include = "#include <gsl/gsl_sf_bessel.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_bessel_zero_Jnu_e;
+ end chap_7_5_13;
+
+ package chap_7_6
+ //clausen function
+
+ function gsl_sf_clausen
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_clausen(x) annotation(
+ Library = "gsl",
+ Library = "gslcblas",
+ Include = "#include<gsl/gsl_sf_clausen.h>");
+ end gsl_sf_clausen;
+
+ function gsl_sf_clausen_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_clausen_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_clausen.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_clausen_e;
+ end chap_7_6;
+
+ package chap_7_7_1
+ function gsl_sf_hydrogenicR_1
+ input Real Z;
+ //atomic number of the atom
+ input Real r;
+ //radus of the atom
+ output Real R;
+ //lowest-order normalized hydrogenic bound state radial wavefunction
+
+ external "C" R = gsl_sf_hydrogenicR_1(Z, r) annotation(
+ Include = "#include <gsl/gsl_sf_coulomb.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_hydrogenicR_1;
+
+ function gsl_sf_hydrogenicR_1_e
+ input Real Z;
+ //atomic number of the atom
+ input Real r;
+ //radus of the atom
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ //lowest-order normalized hydrogenic bound state radial wavefunction
+
+ external "C" y = gsl_sf_hydrogenicR_1_e(Z, r, result) annotation(
+ Include = "#include <gsl/gsl_sf_coulomb.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_hydrogenicR_1_e;
+
+ function gsl_sf_hydrogenicR
+ input Integer n;
+ //order of the wavefunction
+ input Integer l;
+ input Real Z;
+ //atomic number of the atom
+ input Real r;
+ //radus of the atom
+ output Real R;
+ //lowest-order normalized hydrogenic bound state radial wavefunction
+
+ external "C" R = gsl_sf_hydrogenicR(n, l, Z, r) annotation(
+ Include = "#include <gsl/gsl_sf_coulomb.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_hydrogenicR;
+
+ function gsl_sf_hydrogenicR_e
+ input Integer n;
+ input Integer l;
+ input Real Z;
+ //atomic number of the atom
+ input Real r;
+ //radus of the atom
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ //lowest-order normalized hydrogenic bound state radial wavefunction
+
+ external "C" y = gsl_sf_hydrogenicR_e(n, l, Z, r, result) annotation(
+ Include = "#include <gsl/gsl_sf_coulomb.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_hydrogenicR_e;
+ end chap_7_7_1;
+
+ package chap_7_7_2
+ function gsl_sf_coulomb_wave_FG_e
+ input Real eta;
+ // eta=Z/k,where Z is atomic number
+ input Real x;
+ input Real L_F;
+ input Integer k;
+ //L_F-k>-1/2
+ output gsl.data_types.gsl_sf_result F;
+ //coulomb wave function
+ output gsl.data_types.gsl_sf_result Fp;
+ //derivative of wave function
+ output gsl.data_types.gsl_sf_result G;
+ //coulomb wave function
+ output gsl.data_types.gsl_sf_result Gp;
+ //derivative of wave function
+ output Real exp_F;
+ //out of bound error exponenet stored here
+ output Real exp_G;
+ //out of bound error exponent stored here
+ output Integer y;
+
+ external "C" y = gsl_sf_coulomb_wave_FG_e(eta, x, L_F, k, F, Fp, G, Gp, exp_F, exp_G) annotation(
+ Include = "#include<gsl/gsl_sf_coulomb.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coulomb_wave_FG_e;
+
+ /*
+ int gsl_sf_coulomb_wave_F_array(double L_min, int kmax, double eta, double x, double fc_array[],
+ double * F_exponent)
+ */
+
+ function gsl_sf_coulomb_wave_F_array
+ input Real L_min;
+ input Integer kmax;
+ input Real eta;
+ input Real x;
+ output Real fc_array[kmax];
+ output Real F_exponent;
+ output Integer y;
+
+ external "C" y = gsl_sf_coulomb_wave_F_array(L_min, kmax, eta, x, fc_array, F_exponent) annotation(
+ Include = "#include<gsl/gsl_sf_coulomb.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coulomb_wave_F_array;
+
+ /*
+
+ int gsl_sf_coulomb_wave_FG_array(double L_min, int kmax, double eta, double x, double fc_array[],
+ double gc_array[], double * F_exponent, double * G_exponent)*/
+
+ function gsl_sf_coulomb_wave_FG_array
+ input Real L_min;
+ input Integer kmax;
+ input Real eta;
+ input Real x;
+ output Real fc_array[kmax];
+ output Real gc_array[kmax];
+ output Real F_exponent;
+ output Real G_exponent;
+ output Integer y;
+
+ external "C" y = gsl_sf_coulomb_wave_FG_array(L_min, kmax, eta, x, fc_array, gc_array, F_exponent, G_exponent) annotation(
+ Include = "#include<gsl/gsl_sf_coulomb.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coulomb_wave_FG_array;
+
+ function gsl_sf_coulomb_wave_FGp_array
+ input Real L_min;
+ input Integer kmax;
+ input Real eta;
+ input Real x;
+ output Real fc_array[kmax];
+ output Real fcp_array[kmax];
+ output Real gc_array[kmax];
+ output Real gcp_array[kmax];
+ output Real F_exponent;
+ output Real G_exponent;
+ output Integer y;
+
+ external "C" y = gsl_sf_coulomb_wave_FGp_array(L_min, kmax, eta, x, fc_array, fcp_array, gc_array, gcp_array, F_exponent, G_exponent) annotation(
+ Include = "#include<gsl/gsl_sf_coulomb.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coulomb_wave_FGp_array;
+
+ function gsl_sf_coulomb_wave_sphF_array
+ input Real L_min;
+ input Integer kmax;
+ input Real eta;
+ input Real x;
+ output Real fc_array[kmax];
+ output Real F_exponent;
+ output Integer y;
+
+ external "C" y = gsl_sf_coulomb_wave_sphF_array(L_min, kmax, eta, x, fc_array, F_exponent) annotation(
+ Include = "#include<gsl/gsl_sf_coulomb.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coulomb_wave_sphF_array;
+ end chap_7_7_2;
+
+ package chap_7_7_3
+ function gsl_sf_coulomb_CL_e
+ input Real L;
+ input Real eta;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_coulomb_CL_e(L, eta, result) annotation(
+ Include = "#include<gsl/gsl_sf_coulomb.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coulomb_CL_e;
+
+ function gsl_sf_coulomb_CL_array
+ input Real Lmin;
+ input Integer kmax;
+ input Real eta;
+ output Real cl[kmax];
+ output Integer y;
+
+ external "C" y = gsl_sf_coulomb_CL_array(Lmin, kmax, eta, cl) annotation(
+ Include = "#include<gsl/gsl_sf_coulomb.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coulomb_CL_array;
+ end chap_7_7_3;
+
+ package chap_7_8_1
+ /*
+
+ double
+ int gsl_sf_coupling_3j_e(int two_ja, int two_jb, int two_jc, int two_ma, int two_mb, int two_mc,
+ gsl_sf_result * result)
+ These routines compute the Wigner 3-j coefficient,*/
+
+ function gsl_sf_coupling_3j
+ input Integer two_ja;
+ input Integer two_jb;
+ input Integer two_jc;
+ input Integer two_ma;
+ input Integer two_mb;
+ input Integer two_mc;
+ output Real y;
+
+ external "C" y = gsl_sf_coupling_3j(two_ja, two_jb, two_jc, two_ma, two_mb, two_mc) annotation(
+ Include = "#include<gsl/gsl_sf_coupling.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coupling_3j;
+
+ function gsl_sf_coupling_3j_e
+ input Integer two_ja;
+ input Integer two_jb;
+ input Integer two_jc;
+ input Integer two_ma;
+ input Integer two_mb;
+ input Integer two_mc;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_coupling_3j_e(two_ja, two_jb, two_jc, two_ma, two_mb, two_mc, result) annotation(
+ Include = "#include<gsl/gsl_sf_coupling.h>",
+ Include = "#include<gsl/gsl_sf_result>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coupling_3j_e;
+ end chap_7_8_1;
+
+ package chap_7_8_2
+ function gsl_sf_coupling_6j
+ input Integer two_ja;
+ input Integer two_jb;
+ input Integer two_jc;
+ input Integer two_ma;
+ input Integer two_mb;
+ input Integer two_mc;
+ output Real y;
+
+ external "C" y = gsl_sf_coupling_6j(two_ja, two_jb, two_jc, two_ma, two_mb, two_mc) annotation(
+ Include = "#include<gsl/gsl_sf_coupling.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coupling_6j;
+
+ function gsl_sf_coupling_6j_e
+ input Integer two_ja;
+ input Integer two_jb;
+ input Integer two_jc;
+ input Integer two_ma;
+ input Integer two_mb;
+ input Integer two_mc;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_coupling_6j_e(two_ja, two_jb, two_jc, two_ma, two_mb, two_mc, result) annotation(
+ Include = "#include<gsl/gsl_sf_coupling.h>",
+ Include = "#include<gsl/gsl_sf_result>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coupling_6j_e;
+ end chap_7_8_2;
+
+ /*
+
+ double gsl_sf_coupling_9j(int two_ja, int two_jb, int two_jc, int two_jd, int two_je, int two_jf,
+ int two_jg, int two_jh, int two_ji)
+ int gsl_sf_coupling_9j_e(int two_ja, int two_jb, int two_jc, int two_jd, int two_je, int two_jf, int two_jg,
+ int two_jh, int two_ji, gsl_sf_result * result)
+ These routines compute the Wigner 9-j coeffici
+
+ */
+
+ package chap_7_8_3
+ function gsl_sf_coupling_9j
+ input Integer two_ja;
+ input Integer two_jb;
+ input Integer two_jc;
+ input Integer two_jd;
+ input Integer two_je;
+ input Integer two_jf;
+ input Integer two_jg;
+ input Integer two_jh;
+ input Integer two_ji;
+ output Real y;
+
+ external "C" y = gsl_sf_coupling_9j(two_ja, two_jb, two_jc, two_jd, two_je, two_jf, two_jg, two_jh, two_ji) annotation(
+ Include = "#include<gsl/gsl_sf_coupling.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coupling_9j;
+
+ function gsl_sf_coupling_9j_e
+ input Integer two_ja;
+ input Integer two_jb;
+ input Integer two_jc;
+ input Integer two_jd;
+ input Integer two_je;
+ input Integer two_jf;
+ input Integer two_jg;
+ input Integer two_jh;
+ input Integer two_ji;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_coupling_9j_e(two_ja, two_jb, two_jc, two_jd, two_je, two_jf, two_jg, two_jh, two_ji, result) annotation(
+ Include = "#include<gsl/gsl_sf_coupling.h>",
+ Include = "#include<gsl/gsl_sf_result>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_coupling_9j_e;
+ end chap_7_8_3;
+
+ package chap_7_9
+ function gsl_sf_dawson
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_dawson(x) annotation(
+ Include = "#include <gsl/gsl_sf_dawson.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_dawson;
+
+ function gsl_sf_dawson_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_dawson_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_dawson.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_dawson_e;
+ end chap_7_9;
+
+ package chap_7_10
+ // double gsl_sf_debye_1(double x)
+
+ function gsl_sf_debye_1
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_debye_1(x) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_1;
+
+ function gsl_sf_debye_1_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_debye_1_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_1_e;
+
+ function gsl_sf_debye_2
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_debye_2(x) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_2;
+
+ function gsl_sf_debye_2_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_debye_2_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_2_e;
+
+ function gsl_sf_debye_3
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_debye_3(x) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_3;
+
+ function gsl_sf_debye_3_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_debye_3_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_3_e;
+
+ function gsl_sf_debye_4
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_debye_4(x) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_4;
+
+ function gsl_sf_debye_4_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_debye_4_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_4_e;
+
+ function gsl_sf_debye_5
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_debye_5(x) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_5;
+
+ function gsl_sf_debye_5_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_debye_5_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_5_e;
+
+ function gsl_sf_debye_6
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_debye_6(x) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_6;
+
+ function gsl_sf_debye_6_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_debye_6_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_debye.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_debye_6_e;
+ end chap_7_10;
+
+ package chap_7_11_1
+ function gsl_sf_dilog
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_dilog(x) annotation(
+ Include = "#include <gsl/gsl_sf_dilog.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_dilog;
+
+ function gsl_sf_dilog_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_dilog_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_dilog.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_dilog_e;
+ end chap_7_11_1;
+
+ package chap_7_11_2
+ function gsl_sf_complex_dilog_e
+ input Real r;
+ input Real theta;
+ output gsl.data_types.gsl_sf_result result_re;
+ output gsl.data_types.gsl_sf_result result_im;
+ output Integer y;
+
+ external "C" y = gsl_sf_complex_dilog_e(r, theta, result_re, result_im) annotation(
+ Include = "#include<gsl/gsl_sf_dilog.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_complex_dilog_e;
+
+ /* int gsl_sf_complex_dilog_e(double r, double theta, gsl_sf_result * result_re, gsl_sf_result * re-
+ sult_im)*/
+ end chap_7_11_2;
+
+ package chap_7_12
+ function gsl_sf_multiply
+ input Real x;
+ input Real y;
+ output Real z;
+
+ external "C" z = gsl_sf_multiply(x, y) annotation(
+ Include = "#include<gsl/gsl_sf_elementary.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_multiply;
+
+ function gsl_sf_multiply_e
+ input Real x;
+ input Real y;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer z;
+
+ external "C" z = gsl_sf_multiply_e(x, y, result) annotation(
+ Include = "#include<gsl/gsl_sf_elementary.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_multiply_e;
+
+ function gsl_sf_multiply_err_e
+ input Real x;
+ input Real dx;
+ input Real y;
+ input Real dy;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer z;
+
+ external "C" z = gsl_sf_multiply_err_e(x, dx, y, dy, result) annotation(
+ Include = "#include<gsl/gsl_sf_elementary.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_multiply_err_e;
+ end chap_7_12;
+
+ package chap_7_15_1
+ //double gsl_sf_erf(double x)
+ //int gsl_sf_erf_e(double x, gsl_sf_result * result)
+
+ function gsl_sf_erf
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_erf(x) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_erf;
+
+ function gsl_sf_erf_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_erf_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Include = "#include<gsl/gsl.sf.result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_erf_e;
+ end chap_7_15_1;
+
+ package chap_7_15_2
+ function gsl_sf_erfc
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_erfc(x) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_erfc;
+
+ function gsl_sf_erfc_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_erfc_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Include = "#include<gsl/gsl.sf.result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_erfc_e;
+ end chap_7_15_2;
+
+ package chap_7_15_3
+ function gsl_sf_log_erfc
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_log_erfc(x) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_log_erfc;
+
+ function gsl_sf_log_erfc_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_log_erfc_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Include = "#include<gsl/gsl.sf.result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_log_erfc_e;
+ end chap_7_15_3;
+
+ package chap_7_15_4
+ /*
+ double gsl_sf_hazard(double x)
+ int gsl_sf_erf_Z_e(double x, gsl_sf_result * result)
+ double gsl_sf_erf_Q(double x)
+ int gsl_sf_erf_Q_e(double x, gsl_sf_result * result)*/
+
+ function gsl_sf_erf_Z
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_erf_Z(x) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_erf_Z;
+
+ function gsl_sf_erf_Z_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_erf_Z_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Include = "#include<gsl/gsl.sf.result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_erf_Z_e;
+
+ function gsl_sf_erf_Q
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_erf_Q(x) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_erf_Q;
+
+ function gsl_sf_erf_Q_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_erf_Q_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Include = "#include<gsl/gsl.sf.result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_erf_Q_e;
+
+ /*double gsl_sf_hazard(double x)
+ int gsl_sf_hazard_e(double x, gsl_sf_result * result)*/
+
+ function gsl_sf_hazard
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_hazard(x) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_hazard;
+
+ function gsl_sf_hazard_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_hazard_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_erf.h>",
+ Include = "#include<gsl/gsl.sf.result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_hazard_e;
+ end chap_7_15_4;
+
+ package chap_7_16_1
+ /*double gsl_sf_exp(double x)
+ int gsl_sf_exp_e(double x, gsl_sf_result * result)
+ These routines provide an exponential function exp(x) using GSL semantics and error checking.
+ int gsl_sf_exp_e10_e(double x, gsl_sf_result_e10 * result)
+ This function computes the exponential exp(x) using the gsl_sf_result_e10 type to return a result with
+ extended range. This function may be useful if the value of exp(x) would overflow the numeric range of
+ double.
+ double gsl_sf_exp_mult(double x, double y)
+
+ These routines exponentiate x and multiply by the factor y to return the product y exp(x).
+ int gsl_sf_exp_mult_e10_e(const double x, const double y, gsl_sf_result_e10 * result)*/
+
+ function gsl_sf_exp
+ input Real x;
+ //input to the exp function
+ output Real y;
+
+ external "C" y = gsl_sf_exp(x) annotation(
+ Include = "#include <gsl/gsl_sf_exp.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exp;
+
+ function gsl_sf_exp_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_exp_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_exp.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exp_e;
+
+ function gsl_sf_exp_e10_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result_e10 result;
+ output Integer y;
+
+ external "C" y = gsl_sf_exp_e10_e(x, result) annotation(
+ Include = "#include <gsl/gsl_sf_exp.h>",
+ Include = "#include <gsl/gsl_sf_result_e10.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exp_e10_e;
+
+ function gsl_sf_exp_mult
+ input Real x;
+ input Real y;
+ output Real z;
+
+ external "C" z = gsl_sf_exp_mult(x, y) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Library = "gsl",
+ Library = "gslcbla)");
+ end gsl_sf_exp_mult;
+
+ function gsl_sf_exp_mult_e
+ input Real x;
+ input Real y;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer z;
+
+ external "C" z = gsl_sf_exp_mult_e(x, y, result) annotation(
+ Include = "#include <gsl/gsl_sf_exp.h>",
+ Include = "#include <gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exp_mult_e;
+
+ function gsl_sf_exp_mult_e10_e
+ input Real x;
+ input Real y;
+ output gsl.data_types.gsl_sf_result_e10 result;
+ output Integer z;
+
+ external "C" z = gsl_sf_exp_mult_e10_e(x, y, result) annotation(
+ Include = "#include <gsl/gsl_sf_exp.h>",
+ Include = "#include <gsl/gsl_sf_result_e10.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exp_mult_e10_e;
+ end chap_7_16_1;
+
+ package chap_7_16_2
+ function gsl_sf_expm1
+ //This routines compute the quantity exp(x) − 1 using an algorithm that is accurate for small x.
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_expm1(x) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_expm1;
+
+ function gsl_sf_expm1_e
+ //This routines compute the quantity exp(x) − 1 using an algorithm that is accurate for small x.
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_expm1_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_expm1_e;
+
+ function gsl_sf_exprel
+ //These routines compute the quantity (exp(x) − 1)/x using an algorithm that is accurate for small x
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_exprel(x) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exprel;
+
+ function gsl_sf_exprel_e
+ //These routines compute the quantity (exp(x) − 1)/x using an algorithm that is accurate for small x
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_exprel_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exprel_e;
+
+ function gsl_sf_exprel_2
+ //This routines compute the quantity 2(exp(x) − 1 − x)/x^2 using an algorithm that is accurate for small x
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_exprel_2(x) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exprel_2;
+
+ function gsl_sf_exprel_2_e
+ //This routines compute the quantity 2(exp(x) − 1 − x)/x^2 using an algorithm that is accurate for small x
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_exprel_2_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exprel_2_e;
+
+ function gsl_sf_exprel_n
+ input Integer n;
+ input Real x;
+ output Real z;
+
+ external "C" z = gsl_sf_exprel_n(n, x) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exprel_n;
+
+ function gsl_sf_exprel_n_e
+ input Integer n;
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_exprel_n_e(n, x, result) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exprel_n_e;
+ end chap_7_16_2;
+
+ package chap_7_16_3
+ function gsl_sf_exp_err_e
+ //This function exponentiates x with an associated absolute error dx
+ input Real x;
+ input Real dx;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_exp_err_e(x, dx, result) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exp_err_e;
+
+ function gsl_sf_exp_err_e10_e
+ //This function exponentiates a quantity x with an associated absolute error dx using the gsl_sf_result_e10
+ input Real x;
+ input Real dx;
+ output gsl.data_types.gsl_sf_result_e10 result;
+ output Integer y;
+
+ external "C" y = gsl_sf_exp_err_e10_e(x, dx, result) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Include = "#include<gsl/gsl_sf_result_e10.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exp_err_e10_e;
+
+ function gsl_sf_exp_mult_err_e
+ //This routine computes the product y exp(x) for the quantities x, y with associated absolute errors dx, dy.
+ input Real x;
+ input Real dx;
+ input Real y;
+ input Real dy;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer z;
+
+ external "C" z = gsl_sf_exp_mult_err_e(x, dx, y, dy, result) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exp_mult_err_e;
+
+ function gsl_sf_exp_mult_err_e10_e
+ //This routine computes the product y exp(x) for the quantities x, y with associated absolute errors dx, dy usingthe gsl_sf_result_e10 type to return a result with extended range
+ input Real x;
+ input Real dx;
+ input Real y;
+ input Real dy;
+ output gsl.data_types.gsl_sf_result_e10 result;
+ output Integer z;
+
+ external "C" z = gsl_sf_exp_mult_err_e10_e(x, dx, y, dy, result) annotation(
+ Include = "#include<gsl/gsl_sf_exp.h>",
+ Include = "#include<gsl/gsl_sf_result_e10.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_exp_mult_err_e10_e;
+ end chap_7_16_3;
+
+ package chap_7_17_1
+ function gsl_sf_expint_E1
+ //These routines compute the exponential integral E 1 (x),
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_expint_E1(x) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_expint_E1;
+
+ function gsl_sf_expint_E1_e
+ //These routines compute the exponential integral E 1 (x),
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_expint_E1_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_expint_E1_e;
+
+ function gsl_sf_expint_E2
+ //These routines compute the exponential integral E 2 (x),
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_expint_E2(x) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_expint_E2;
+
+ function gsl_sf_expint_E2_e
+ //These routines compute the exponential integral E 2 (x),
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_expint_E2_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_expint_E2_e;
+
+ function gsl_sf_expint_En
+ //These routines compute the exponential integral E n (x),
+ input Integer n;
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_expint_En(n, x) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_expint_En;
+
+ function gsl_sf_expint_En_e
+ //These routines compute the exponential integral E 2 (x),
+ input Integer n;
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_expint_En_e(n, x, result) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_expint_En_e;
+ end chap_7_17_1;
+
+ package chap_7_17_2
+ function gsl_sf_expint_Ei
+ //These routines compute the exponential integral Ei(x)
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_expint_Ei(x) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_expint_Ei;
+
+ function gsl_sf_expint_Ei_e
+ //These routines compute the exponential integral Ei(x)
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_expint_Ei_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_expint_Ei_e;
+ end chap_7_17_2;
+
+ package chap_7_17_3
+ function gsl_sf_Shi
+ //hyperbolic integrals
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_Shi(x) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_Shi;
+
+ function gsl_sf_Shi_e
+ //hyperbolic integral
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_Shi_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_Shi_e;
+
+
+
+ function gsl_sf_Chi
+ //hyperbolic integrals
+ input Real x;
+ output Real y;
+
+ external "C" y = gsl_sf_Chi(x) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_Chi;
+
+ function gsl_sf_Chi_e
+ //hyperbolic integral
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+
+ external "C" y = gsl_sf_Chi_e(x, result) annotation(
+ Include = "#include<gsl/gsl_sf_expint.h>",
+ Include = "#include<gsl/gsl_sf_result.h>",
+ Library = "gsl",
+ Library = "gslcblas");
+ end gsl_sf_Chi_e;
+ end chap_7_17_3;
+
+ package chap_7_17_4
+ function gsl_sf_expint_3
+ input Real x;
+ output Real y;
+ external "C" y=gsl_sf_expint_3(x)annotation(Include="#include<gsl/gsl_sf_expint.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_expint_3;
+
+ function gsl_sf_expint_3_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ external "C" y=gsl_sf_expint_3_e(x,result)annotation(Include="#include<gsl/gsl_sf_expint.h>",Include="#include<gsl/gsl_sf_result.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_expint_3_e;
+ end chap_7_17_4;
+
+ package chap_7_17_5
+ function gsl_sf_Si
+ input Real x;
+ output Real y;
+ external "C" y=gsl_sf_Si(x)annotation(Include="#include<gsl/gsl_sf_expint.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_Si;
+
+ function gsl_sf_Si_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ external "C" y=gsl_sf_Si_e(x,result)annotation(Include="#include<gsl/gsl_sf_expint.h>",Include="#include<gsl/gsl_sf_result>",Library="gsl",Library="gslcbas");
+ end gsl_sf_Si_e;
+
+ function gsl_sf_Ci
+ input Real x;
+ output Real y;
+ external "C" y=gsl_sf_Ci(x)annotation(Include="#include<gsl/gsl_sf_expint.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_Ci;
+
+ function gsl_sf_Ci_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ external "C" y=gsl_sf_Ci_e(x,result)annotation(Include="#include<gsl/gsl_sf_expint.h>",Include="#include<gsl/gsl_sf_result.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_Ci_e;
+ end chap_7_17_5;
+
+ package chap_7_18_1
+
+
+ function gsl_sf_fermi_dirac_m1
+ input Real x;
+ output Real y;
+ external "C" y=gsl_sf_fermi_dirac_m1(x)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_m1;
+
+
+ //int gsl_sf_fermi_dirac_m1_e(double x, gsl_sf_result * result)
+ function gsl_sf_fermi_dirac_m1_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ external "C" y=gsl_sf_fermi_dirac_m1_e(x,result)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Library="gsl",Library="gslcblas",Include="#include<gsl/gsl_sf_result.h>");
+ end gsl_sf_fermi_dirac_m1_e;
+
+ //double gsl_sf_fermi_dirac_0(double x)
+ function gsl_sf_fermi_dirac_0
+ input Real x;
+ output Real y;
+ external "C" y=gsl_sf_fermi_dirac_0(x)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_0;
+
+ function gsl_sf_fermi_dirac_0_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ external "C" y=gsl_sf_fermi_dirac_0_e(x,result)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Include="#include<gsl/gsl_sf_result.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_0_e;
+
+
+ //double gsl_sf_fermi_dirac_1(double x)
+ function gsl_sf_fermi_dirac_1
+ input Real x;
+ output Real y;
+ external "C" y=gsl_sf_fermi_dirac_1(x)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_1;
+
+ //int gsl_sf_fermi_dirac_1_e(double x, gsl_sf_result * result)
+ function gsl_sf_fermi_dirac_1_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ external "C" y=gsl_sf_fermi_dirac_1_e(x,result)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Include="#include<gsl/gsl_sf_result.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_1_e;
+
+
+
+ //double gsl_sf_fermi_dirac_2(double x)
+ function gsl_sf_fermi_dirac_2
+ input Real x;
+ output Real y;
+ external "C" y=gsl_sf_fermi_dirac_2(x)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_2;
+
+ //int gsl_sf_fermi_dirac_2_e(double x, gsl_sf_result * result)
+ function gsl_sf_fermi_dirac_2_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ external "C" y=gsl_sf_fermi_dirac_2_e(x,result)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Include="#include<gsl/gsl_sf_result.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_2_e;
+
+
+
+
+
+
+ //double gsl_sf_fermi_dirac_int(int j, double x)
+ function gsl_sf_fermi_dirac_int
+ input Integer j;
+ input Real x;
+ output Real y;
+ external "C" y=gsl_sf_fermi_dirac_int(j,x)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_int;
+
+ function gsl_sf_fermi_dirac_int_e
+ input Integer j;
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Real y;
+ external "C" y=gsl_sf_fermi_dirac_int_e(j,x,result)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Include="#include<gsl/gsl_sf_result.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_int_e;
+
+
+
+
+
+
+
+
+
+
+
+
+
+ function gsl_sf_fermi_dirac_mhalf
+ input Real x;
+ output Real y;
+ external "C" y=gsl_sf_fermi_dirac_mhalf(x)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_mhalf;
+
+
+
+ function gsl_sf_fermi_dirac_mhalf_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ external "C" y=gsl_sf_fermi_dirac_mhalf_e(x,result)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Library="gsl",Library="gslcblas",Include="#include<gsl/gsl_sf_result.h>");
+ end gsl_sf_fermi_dirac_mhalf_e;
+
+
+ function gsl_sf_fermi_dirac_half
+ input Real x;
+ output Real y;
+ external "C" y=gsl_sf_fermi_dirac_half(x)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_half;
+
+ function gsl_sf_fermi_dirac_half_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ external "C" y=gsl_sf_fermi_dirac_half_e(x,result)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Include="#include<gsl/gsl_sf_result.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_half_e;
+
+
+
+ function gsl_sf_fermi_dirac_3half
+ input Real x;
+ output Real y;
+ external "C" y=gsl_sf_fermi_dirac_3half(x)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_3half;
+
+
+ function gsl_sf_fermi_dirac_3half_e
+ input Real x;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ external "C" y=gsl_sf_fermi_dirac_3half_e(x,result)annotation(Include="#include<gsl/gsl_sf_fermi_dirac.h>",Include="#include<gsl/gsl_sf_result.h>",Library="gsl",Library="gslcblas");
+ end gsl_sf_fermi_dirac_3half_e;
+
+
+ end chap_7_18_1;
+ end Special_functions;
+
package Examples
package Mathematical
model gsl_log1p
//This model computes the value of log(1 + x) in a way that is accurate for small x by calling the function gsl_log1p(x)
- parameter Real x = 0.01;
+ parameter Real x = -2;
Real y;
algorithm
y := gsl.mathematical.gsl_log1p(x);
@@ -515,8 +3034,8 @@ package gsl
model gsl_pow_int
//this function computes x^n
- parameter Real x = 2.0;
- parameter Integer n = 2;
+ parameter Real x = 0;
+ parameter Integer n = 0;
Real y;
algorithm
y := gsl.mathematical.gsl_pow_int(x, n);
@@ -648,7 +3167,7 @@ package gsl
model GSL_MIN_INT
// This model gives the minimum of two integers by calling the function GSL_MIN_INT
- parameter Real a = 2;
+ parameter Real a = 2.5;
parameter Real b = 3;
Real c;
algorithm
@@ -705,6 +3224,1820 @@ package gsl
end Permutation;
package COMPLEX
+ model gsl_complex_rect
+ /*this model calls the function gsl_complex_rect*/
+ parameter Real a = 10;
+ parameter Real b = 4;
+ Real z[2];
+ algorithm
+ z := gsl.COMPLEX.gsl_complex_rect(a, b);
+ end gsl_complex_rect;
end COMPLEX;
+
+ package special_functions
+ package chap_7_5_1
+ model gsl_sf_bessel_J0
+ //this model calls the function gsl_sf_bessel_J0 by taking value of x as input and returning y as output
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_1.gsl_sf_bessel_J0(x);
+ end gsl_sf_bessel_J0;
+
+ model gsl_sf_bessel_J0_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_1.gsl_sf_bessel_J0_e(x);
+ end gsl_sf_bessel_J0_e;
+
+ model gsl_sf_bessel_J1
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_1.gsl_sf_bessel_J1(x);
+ end gsl_sf_bessel_J1;
+
+ model gsl_sf_bessel_J1_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_1.gsl_sf_bessel_J1_e(x);
+ end gsl_sf_bessel_J1_e;
+
+ model gsl_sf_bessel_Jn
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_1.gsl_sf_bessel_Jn(x, n);
+ end gsl_sf_bessel_Jn;
+
+ model gsl_sf_bessel_Jn_e
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_1.gsl_sf_bessel_Jn_e(x, n);
+ end gsl_sf_bessel_Jn_e;
+
+ model gsl_sf_bessel_Jn_array
+ parameter Integer nmin = 2;
+ parameter Integer nmax = 10;
+ parameter Real x = 4;
+ Real result_array[nmax - nmin];
+ Integer y;
+ algorithm
+ (result_array, y) := gsl.Special_functions.chap_7_5_1.gsl_sf_bessel_Jn_array(nmin, nmax, x);
+ end gsl_sf_bessel_Jn_array;
+ end chap_7_5_1;
+
+ package chap_7_5_2
+ model gsl_sf_bessel_Y0
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_2.gsl_sf_bessel_Y0(x);
+ end gsl_sf_bessel_Y0;
+
+ model gsl_sf_bessel_Y0_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_2.gsl_sf_bessel_Y0_e(x);
+ end gsl_sf_bessel_Y0_e;
+
+ model gsl_sf_bessel_Y1
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_2.gsl_sf_bessel_Y1(x);
+ end gsl_sf_bessel_Y1;
+
+ model gsl_sf_bessel_Y1_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_2.gsl_sf_bessel_Y1_e(x);
+ end gsl_sf_bessel_Y1_e;
+
+ model gsl_sf_bessel_Yn
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_2.gsl_sf_bessel_Yn(x, n);
+ end gsl_sf_bessel_Yn;
+
+ model gsl_sf_bessel_Yn_e
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_2.gsl_sf_bessel_Yn_e(x, n);
+ end gsl_sf_bessel_Yn_e;
+
+ model gsl_sf_bessel_Yn_array
+ parameter Integer nmin = 2;
+ parameter Integer nmax = 10;
+ parameter Real x = 4;
+ Real result_array[nmax - nmin];
+ Integer y;
+ algorithm
+ (result_array, y) := gsl.Special_functions.chap_7_5_2.gsl_sf_bessel_Yn_array(nmin, nmax, x);
+ end gsl_sf_bessel_Yn_array;
+ end chap_7_5_2;
+
+ package chap_7_5_3
+ model gsl_sf_bessel_I0
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_I0(x);
+ end gsl_sf_bessel_I0;
+
+ model gsl_sf_bessel_I0_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_I0_e(x);
+ end gsl_sf_bessel_I0_e;
+
+ model gsl_sf_bessel_I1
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_I1(x);
+ end gsl_sf_bessel_I1;
+
+ model gsl_sf_bessel_I1_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_I1_e(x);
+ end gsl_sf_bessel_I1_e;
+
+ model gsl_sf_bessel_In
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_In(x, n);
+ end gsl_sf_bessel_In;
+
+ model gsl_sf_bessel_In_e
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_In_e(x, n);
+ end gsl_sf_bessel_In_e;
+
+ model gsl_sf_bessel_In_array
+ parameter Integer nmin = 2;
+ parameter Integer nmax = 10;
+ parameter Real x = 4;
+ Real result_array[nmax - nmin];
+ Integer y;
+ algorithm
+ (result_array, y) := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_In_array(nmin, nmax, x);
+ end gsl_sf_bessel_In_array;
+
+ model gsl_sf_bessel_I0_scaled
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_I0_scaled(x);
+ end gsl_sf_bessel_I0_scaled;
+
+ model gsl_sf_bessel_I0_scaled_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_I0_scaled_e(x);
+ end gsl_sf_bessel_I0_scaled_e;
+
+ model gsl_sf_bessel_I1_scaled
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_I1_scaled(x);
+ end gsl_sf_bessel_I1_scaled;
+
+ model gsl_sf_bessel_I1_scaled_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_I1_scaled_e(x);
+ end gsl_sf_bessel_I1_scaled_e;
+
+ model gsl_sf_bessel_In_scaled
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_In_scaled(x, n);
+ end gsl_sf_bessel_In_scaled;
+
+ model gsl_sf_bessel_In_scaled_e
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_In_scaled_e(x, n);
+ end gsl_sf_bessel_In_scaled_e;
+
+ model gsl_sf_bessel_In_scaled_array
+ parameter Integer nmin = 2;
+ parameter Integer nmax = 10;
+ parameter Real x = 4;
+ Real result_array[nmax - nmin];
+ Integer y;
+ algorithm
+ (result_array, y) := gsl.Special_functions.chap_7_5_3.gsl_sf_bessel_In_scaled_array(nmin, nmax, x);
+ end gsl_sf_bessel_In_scaled_array;
+ end chap_7_5_3;
+
+ package chap_7_5_4
+ model gsl_sf_bessel_K0
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_K0(x);
+ end gsl_sf_bessel_K0;
+
+ model gsl_sf_bessel_K0_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_K0_e(x);
+ end gsl_sf_bessel_K0_e;
+
+ model gsl_sf_bessel_K1
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_K1(x);
+ end gsl_sf_bessel_K1;
+
+ model gsl_sf_bessel_K1_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_K1_e(x);
+ end gsl_sf_bessel_K1_e;
+
+ model gsl_sf_bessel_Kn
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_Kn(x, n);
+ end gsl_sf_bessel_Kn;
+
+ model gsl_sf_bessel_Kn_e
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_Kn_e(x, n);
+ end gsl_sf_bessel_Kn_e;
+
+ model gsl_sf_bessel_Kn_array
+ parameter Integer nmin = 2;
+ parameter Integer nmax = 10;
+ parameter Real x = 4;
+ Real result_array[nmax - nmin];
+ Integer y;
+ algorithm
+ (result_array, y) := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_Kn_array(nmin, nmax, x);
+ end gsl_sf_bessel_Kn_array;
+
+ model gsl_sf_bessel_K0_scaled
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_K0_scaled(x);
+ end gsl_sf_bessel_K0_scaled;
+
+ model gsl_sf_bessel_K0_scaled_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_K0_scaled_e(x);
+ end gsl_sf_bessel_K0_scaled_e;
+
+ model gsl_sf_bessel_K1_scaled
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_K1_scaled(x);
+ end gsl_sf_bessel_K1_scaled;
+
+ model gsl_sf_bessel_K1_scaled_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_K1_scaled_e(x);
+ end gsl_sf_bessel_K1_scaled_e;
+
+ model gsl_sf_bessel_Kn_scaled
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_Kn_scaled(x, n);
+ end gsl_sf_bessel_Kn_scaled;
+
+ model gsl_sf_bessel_Kn_scaled_e
+ parameter Real x = 10;
+ parameter Integer n = 1;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_Kn_scaled_e(x, n);
+ end gsl_sf_bessel_Kn_scaled_e;
+
+ model gsl_sf_bessel_Kn_scaled_array
+ parameter Integer nmin = 2;
+ parameter Integer nmax = 10;
+ parameter Real x = 4;
+ Real result_array[nmax - nmin];
+ Integer y;
+ algorithm
+ (result_array, y) := gsl.Special_functions.chap_7_5_4.gsl_sf_bessel_Kn_scaled_array(nmin, nmax, x);
+ end gsl_sf_bessel_Kn_scaled_array;
+ end chap_7_5_4;
+
+ package chap_7_5_5
+ model gsl_sf_bessel_j0
+ //This model calls the function gsl_sf_bessel_j0 to calculate the value of spherical bessel function j0(x)=sin(x)/x
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j0(x);
+ end gsl_sf_bessel_j0;
+
+ model gsl_sf_bessel_j0_e
+ //This model calls the function gsl_sf_bessel_j0 to calculate the value of spherical bessel function j0(x)=sin(x)/x
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j0_e(x);
+ end gsl_sf_bessel_j0_e;
+
+ model gsl_sf_bessel_j1
+ //This model calls the function gsl_sf_bessel_j1 to calculate the value of spherical bessel function j1(x)=(sin(x)/x - cos(x))/x
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j1(x);
+ end gsl_sf_bessel_j1;
+
+ model gsl_sf_bessel_j1_e
+ //This model calls the function gsl_sf_bessel_j1 to calculate the value of spherical bessel function j1(x)=(sin(x)/x -cos(x))/x
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j1_e(x);
+ end gsl_sf_bessel_j1_e;
+
+ model gsl_sf_bessel_j2
+ //These routines compute the regular spherical Bessel function of second order, j 2 (x) = ((3/x 2 − 1) sin(x) −3 cos(x)/x)/x
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j2(x);
+ end gsl_sf_bessel_j2;
+
+ model gsl_sf_bessel_j2_e
+ //These routines compute the regular spherical Bessel function of second order, j 2 (x) = ((3/x 2 − 1) sin(x) −3 cos(x)/x)/x
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j2_e(x);
+ end gsl_sf_bessel_j2_e;
+
+ model gsl_sf_bessel_jl
+ //These routines compute the regular spherical Bessel function of order l, j1(x), for l ≥ 0 and x ≥ 0.
+ parameter Integer l = 10;
+ parameter Real x = 10;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_jl(l, x);
+ end gsl_sf_bessel_jl;
+
+ model gsl_sf_bessel_jl_e
+ parameter Integer l = 10;
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_jl_e(l, x);
+ end gsl_sf_bessel_jl_e;
+
+ model gsl_sf_bessel_jl_array
+ // parameter Integer nmin;
+ parameter Integer lmax = 10;
+ parameter Real x = 10;
+ Real result_array[lmax];
+ Integer y;
+ algorithm
+ (result_array, y) := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_jl_array(lmax, x);
+ end gsl_sf_bessel_jl_array;
+
+ model gsl_sf_bessel_jl_steed_array
+ // parameter Integer nmin;
+ parameter Integer lmax = 10;
+ parameter Real x = 10;
+ Real result_array[lmax];
+ Integer y;
+ algorithm
+ (result_array, y) := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_jl_steed_array(lmax, x);
+ end gsl_sf_bessel_jl_steed_array;
+ end chap_7_5_5;
+
+ package chap_7_5_6
+ model gsl_sf_bessel_j0
+ //This model calls the function gsl_sf_bessel_j0 to calculate the value of irregular spherical irregular spherical bessel function j0(x)=sin(x)/x
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j0(x);
+ end gsl_sf_bessel_j0;
+
+ model gsl_sf_bessel_j0_e
+ //This model calls the function gsl_sf_bessel_j0 to calculate the value of irregular spherical bessel function j0(x)=sin(x)/x
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j0_e(x);
+ end gsl_sf_bessel_j0_e;
+
+ model gsl_sf_bessel_j1
+ //This model calls the function gsl_sf_bessel_j1 to calculate the value of irregular spherical bessel function j1(x)=(sin(x)/x - cos(x))/x
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j1(x);
+ end gsl_sf_bessel_j1;
+
+ model gsl_sf_bessel_j1_e
+ //This model calls the function gsl_sf_bessel_j1 to calculate the value of irregular spherical bessel function j1(x)=(sin(x)/x -cos(x))/x
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j1_e(x);
+ end gsl_sf_bessel_j1_e;
+
+ model gsl_sf_bessel_j2
+ //These routines compute the regular irregular spherical Bessel function of second order, j 2 (x) = ((3/x 2 − 1) sin(x) −3 cos(x)/x)/x
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j2(x);
+ end gsl_sf_bessel_j2;
+
+ model gsl_sf_bessel_j2_e
+ //These routines compute the irregular spherical Bessel function of second order, j 2 (x) = ((3/x 2 − 1) sin(x) −3 cos(x)/x)/x
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_j2_e(x);
+ end gsl_sf_bessel_j2_e;
+
+ model gsl_sf_bessel_jl
+ //These routines compute the irregular spherical Bessel function of order l, j1(x), for l ≥ 0 and x ≥ 0.
+ parameter Integer l = 10;
+ parameter Real x = 10;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_jl(l, x);
+ end gsl_sf_bessel_jl;
+
+ model gsl_sf_bessel_jl_e
+ parameter Integer l = 10;
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_jl_e(l, x);
+ end gsl_sf_bessel_jl_e;
+
+ model gsl_sf_bessel_jl_array
+ // parameter Integer nmin;
+ parameter Integer lmax = 10;
+ parameter Real x = 10;
+ Real result_array[lmax];
+ Integer y;
+ algorithm
+ (result_array, y) := gsl.Special_functions.chap_7_5_5.gsl_sf_bessel_jl_array(lmax, x);
+ end gsl_sf_bessel_jl_array;
+ end chap_7_5_6;
+
+ package chap_7_5_7
+ model gsl_sf_bessel_i0_scaled
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_7.gsl_sf_bessel_i0_scaled(x);
+ end gsl_sf_bessel_i0_scaled;
+
+ model gsl_sf_bessel_i0_scaled_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_7.gsl_sf_bessel_i0_scaled_e(x);
+ end gsl_sf_bessel_i0_scaled_e;
+
+ model gsl_sf_bessel_i1_scaled
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_7.gsl_sf_bessel_i1_scaled(x);
+ end gsl_sf_bessel_i1_scaled;
+
+ model gsl_sf_bessel_i1_scaled_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_7.gsl_sf_bessel_i1_scaled_e(x);
+ end gsl_sf_bessel_i1_scaled_e;
+
+ model gsl_sf_bessel_i2_scaled
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_7.gsl_sf_bessel_i2_scaled(x);
+ end gsl_sf_bessel_i2_scaled;
+
+ model gsl_sf_bessel_i2_scaled_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_7.gsl_sf_bessel_i2_scaled_e(x);
+ end gsl_sf_bessel_i2_scaled_e;
+
+ model gsl_sf_bessel_il_scaled
+ parameter Real x = 10;
+ parameter Integer l = 1;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_7.gsl_sf_bessel_il_scaled(x, l);
+ end gsl_sf_bessel_il_scaled;
+
+ model gsl_sf_bessel_il_scaled_e
+ parameter Real x = 10;
+ parameter Integer l = 1;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_7.gsl_sf_bessel_il_scaled_e(x, l);
+ end gsl_sf_bessel_il_scaled_e;
+
+ model gsl_sf_bessel_il_scaled_array
+ parameter Integer lmax = 10;
+ parameter Real x = 4;
+ Real result_array[lmax];
+ Integer y;
+ algorithm
+ (result_array, y) := gsl.Special_functions.chap_7_5_7.gsl_sf_bessel_il_scaled_array(lmax, x);
+ end gsl_sf_bessel_il_scaled_array;
+ end chap_7_5_7;
+
+ package chap_7_5_8
+ model gsl_sf_bessel_k0_scaled
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_8.gsl_sf_bessel_k0_scaled(x);
+ end gsl_sf_bessel_k0_scaled;
+
+ model gsl_sf_bessel_k0_scaled_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_8.gsl_sf_bessel_k0_scaled_e(x);
+ end gsl_sf_bessel_k0_scaled_e;
+
+ model gsl_sf_bessel_k1_scaled
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_8.gsl_sf_bessel_k1_scaled(x);
+ end gsl_sf_bessel_k1_scaled;
+
+ model gsl_sf_bessel_k1_scaled_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_8.gsl_sf_bessel_k1_scaled_e(x);
+ end gsl_sf_bessel_k1_scaled_e;
+
+ model gsl_sf_bessel_k2_scaled
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_8.gsl_sf_bessel_k2_scaled(x);
+ end gsl_sf_bessel_k2_scaled;
+
+ model gsl_sf_bessel_k2_scaled_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_8.gsl_sf_bessel_k2_scaled_e(x);
+ end gsl_sf_bessel_k2_scaled_e;
+
+ model gsl_sf_bessel_kl_scaled
+ parameter Real x = 10;
+ parameter Integer l = 1;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_8.gsl_sf_bessel_kl_scaled(x, l);
+ end gsl_sf_bessel_kl_scaled;
+
+ model gsl_sf_bessel_kl_scaled_e
+ parameter Real x = 10;
+ parameter Integer l = 1;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_8.gsl_sf_bessel_kl_scaled_e(x, l);
+ end gsl_sf_bessel_kl_scaled_e;
+
+ model gsl_sf_bessel_kl_scaled_array
+ parameter Integer lmax = 10;
+ parameter Real x = 4;
+ Real result_array[lmax];
+ Integer y;
+ algorithm
+ (result_array, y) := gsl.Special_functions.chap_7_5_8.gsl_sf_bessel_kl_scaled_array(lmax, x);
+ end gsl_sf_bessel_kl_scaled_array;
+ end chap_7_5_8;
+
+ package chap_7_5_9
+ model gsl_sf_bessel_Jnu
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_5_9.gsl_sf_bessel_Jnu(nu, x);
+ end gsl_sf_bessel_Jnu;
+
+ model gsl_sf_bessel_Jnu_e
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_9.gsl_sf_bessel_Jnu_e(x, nu);
+ end gsl_sf_bessel_Jnu_e;
+ end chap_7_5_9;
+
+ package chap_7_5_10
+ model gsl_sf_bessel_Ynu
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_5_10.gsl_sf_bessel_Ynu(nu, x);
+ end gsl_sf_bessel_Ynu;
+
+ model gsl_sf_bessel_Ynu_e
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_10.gsl_sf_bessel_Ynu_e(x, nu);
+ end gsl_sf_bessel_Ynu_e;
+ end chap_7_5_10;
+
+ package chap_7_5_11
+ model gsl_sf_bessel_Inu
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_5_11.gsl_sf_bessel_Inu(nu, x);
+ end gsl_sf_bessel_Inu;
+
+ model gsl_sf_bessel_Inu_e
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_11.gsl_sf_bessel_Inu_e(x, nu);
+ end gsl_sf_bessel_Inu_e;
+
+ model gsl_sf_bessel_Inu_scaled
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_11.gsl_sf_bessel_Inu_scaled(x, nu);
+ end gsl_sf_bessel_Inu_scaled;
+
+ model gsl_sf_bessel_Inu_scaled_e
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_11.gsl_sf_bessel_Inu_scaled_e(x, nu);
+ end gsl_sf_bessel_Inu_scaled_e;
+ end chap_7_5_11;
+
+ package chap_7_5_12
+ model gsl_sf_bessel_Knu
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_5_12.gsl_sf_bessel_Knu(nu, x);
+ end gsl_sf_bessel_Knu;
+
+ model gsl_sf_bessel_Knu_e
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_12.gsl_sf_bessel_Knu_e(x, nu);
+ end gsl_sf_bessel_Knu_e;
+
+ model gsl_sf_bessel_lnKnu
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_5_12.gsl_sf_bessel_lnKnu(nu, x);
+ end gsl_sf_bessel_lnKnu;
+
+ model gsl_sf_bessel_lnKnu_e
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_12.gsl_sf_bessel_lnKnu_e(x, nu);
+ end gsl_sf_bessel_lnKnu_e;
+
+ model gsl_sf_bessel_Knu_scaled
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_5_12.gsl_sf_bessel_Knu_scaled(x, nu);
+ end gsl_sf_bessel_Knu_scaled;
+
+ model gsl_sf_bessel_Knu_scaled_e
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_12.gsl_sf_bessel_Knu_scaled_e(x, nu);
+ end gsl_sf_bessel_Knu_scaled_e;
+ end chap_7_5_12;
+
+ package chap_7_5_13
+ model gsl_sf_bessel_zero_J0
+ parameter Real x = 10;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_5_13.gsl_sf_bessel_zero_J0(x);
+ end gsl_sf_bessel_zero_J0;
+
+ model gsl_sf_bessel_zero_J0_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_13.gsl_sf_bessel_zero_J0_e(x);
+ end gsl_sf_bessel_zero_J0_e;
+
+ model gsl_sf_bessel_zero_J1
+ parameter Real x = 10;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_5_13.gsl_sf_bessel_zero_J1(x);
+ end gsl_sf_bessel_zero_J1;
+
+ model gsl_sf_bessel_zero_J1_e
+ parameter Real x = 10;
+ output gsl.data_types.gsl_sf_result result;
+ output Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_13.gsl_sf_bessel_zero_J1_e(x);
+ end gsl_sf_bessel_zero_J1_e;
+
+ model gsl_sf_bessel_zero_Jnu
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ output Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_5_13.gsl_sf_bessel_zero_Jnu(nu, x);
+ end gsl_sf_bessel_zero_Jnu;
+
+ model gsl_sf_bessel_zero_Jnu_e
+ parameter Real x = 10;
+ parameter Real nu = 1.5;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_5_13.gsl_sf_bessel_zero_Jnu_e(nu, x);
+ end gsl_sf_bessel_zero_Jnu_e;
+ end chap_7_5_13;
+
+ package chap_7_6
+ model gsl_sf_clausen
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_6.gsl_sf_clausen(x);
+ end gsl_sf_clausen;
+
+ model gsl_sf_clausen_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_6.gsl_sf_clausen_e(x);
+ end gsl_sf_clausen_e;
+ end chap_7_6;
+
+ package chap_7_7_1
+ model gsl_sf_hydrogenicR_1
+ parameter Real Z = 10;
+ parameter Real r = 10;
+ Real R;
+ algorithm
+ R := gsl.Special_functions.chap_7_7_1.gsl_sf_hydrogenicR_1(Z, r);
+ end gsl_sf_hydrogenicR_1;
+
+ model gsl_sf_hydrogenicR_1_e
+ parameter Real Z = 10;
+ parameter Real r = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_7_1.gsl_sf_hydrogenicR_1_e(Z, r);
+ end gsl_sf_hydrogenicR_1_e;
+
+ model gsl_sf_hydrogenicR
+ parameter Integer n = 1;
+ parameter Integer l = 0;
+ parameter Real Z = 10;
+ parameter Real r = 0.1;
+ Real R;
+ algorithm
+ R := gsl.Special_functions.chap_7_7_1.gsl_sf_hydrogenicR(n, l, Z, r);
+ end gsl_sf_hydrogenicR;
+
+ model gsl_sf_hydrogenicR_e
+ parameter Integer n = 1;
+ parameter Integer l = 0;
+ parameter Real Z = 10;
+ parameter Real r = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_7_1.gsl_sf_hydrogenicR_e(n, l, Z, r);
+ end gsl_sf_hydrogenicR_e;
+ end chap_7_7_1;
+
+ package chap_7_7_2
+ model gsl_sf_coulomb_wave_FG_e
+ parameter Real eta = 1;
+ // eta=Z/k,where Z is atomic number
+ parameter Real x = 0.1;
+ parameter Real L_F = 1;
+ parameter Integer k = 1;
+ //L_F-k>-1/2
+ gsl.data_types.gsl_sf_result F;
+ //coulomb wave function
+ gsl.data_types.gsl_sf_result Fp;
+ //derivative of wave function
+ gsl.data_types.gsl_sf_result G;
+ //coulomb wave function
+ gsl.data_types.gsl_sf_result Gp;
+ //derivative of wave function
+ Real exp_F;
+ //out of bound error exponenet stored here
+ Real exp_G;
+ //out of bound error exponent stored here
+ Integer y;
+ algorithm
+ (F, Fp, G, Gp, exp_F, exp_G, y) := gsl.Special_functions.chap_7_7_2.gsl_sf_coulomb_wave_FG_e(eta, x, L_F, k);
+ end gsl_sf_coulomb_wave_FG_e;
+
+ model gsl_sf_coulomb_wave_F_array
+ parameter Real L_min = 1;
+ parameter Integer kmax = 2;
+ parameter Real eta = 1;
+ parameter Real x = 1;
+ Real fc_array[kmax];
+ Real F_exponent;
+ Integer y;
+ algorithm
+ (fc_array, F_exponent, y) := gsl.Special_functions.chap_7_7_2.gsl_sf_coulomb_wave_F_array(L_min, kmax, eta, x);
+ end gsl_sf_coulomb_wave_F_array;
+
+ model gsl_sf_coulomb_wave_FG_array
+ parameter Real L_min = 1;
+ parameter Integer kmax = 2;
+ parameter Real eta = 1;
+ parameter Real x = 1.5;
+ Real fc_array[kmax];
+ Real gc_array[kmax];
+ Real F_exponent;
+ Real G_exponent;
+ Integer y;
+ algorithm
+ (fc_array, gc_array, F_exponent, G_exponent, y) := gsl.Special_functions.chap_7_7_2.gsl_sf_coulomb_wave_FG_array(L_min, kmax, eta, x);
+ end gsl_sf_coulomb_wave_FG_array;
+
+ model gsl_sf_coulomb_wave_FGp_array
+ parameter Real L_min = 1;
+ parameter Integer kmax = 2;
+ parameter Real eta = 1;
+ parameter Real x = 1.5;
+ Real fc_array[kmax];
+ Real fcp_array[kmax];
+ Real gc_array[kmax];
+ Real gcp_array[kmax];
+ Real F_exponent;
+ Real G_exponent;
+ Integer y;
+ algorithm
+ (fc_array, fcp_array, gc_array, gcp_array, F_exponent, G_exponent, y) := gsl.Special_functions.chap_7_7_2.gsl_sf_coulomb_wave_FGp_array(L_min, kmax, eta, x);
+ end gsl_sf_coulomb_wave_FGp_array;
+
+ model gsl_sf_coulomb_wave_sphF_array
+ parameter Real L_min = 1;
+ parameter Integer kmax = 2;
+ parameter Real eta = 1;
+ parameter Real x = 1.5;
+ Real fc_array[kmax];
+ Real F_exponent;
+ Integer y;
+ algorithm
+ (fc_array, F_exponent, y) := gsl.Special_functions.chap_7_7_2.gsl_sf_coulomb_wave_sphF_array(L_min, kmax, eta, x);
+ end gsl_sf_coulomb_wave_sphF_array;
+ end chap_7_7_2;
+
+ package chap_7_7_3
+ model gsl_sf_coulomb_CL_e
+ parameter Real L = 1;
+ parameter Real eta = 1.5;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_7_3.gsl_sf_coulomb_CL_e(L, eta);
+ end gsl_sf_coulomb_CL_e;
+
+ model gsl_sf_coulomb_CL_array
+ parameter Real Lmin = 1;
+ parameter Integer kmax = 2;
+ parameter Real eta = 1.5;
+ Real cl[kmax];
+ Integer y;
+ algorithm
+ (cl, y) := gsl.Special_functions.chap_7_7_3.gsl_sf_coulomb_CL_array(Lmin, kmax, eta);
+ end gsl_sf_coulomb_CL_array;
+ end chap_7_7_3;
+
+ package chap_7_8_1
+ model gsl_sf_coupling_3j
+ parameter Integer two_ja = 0;
+ parameter Integer two_jb = 0;
+ parameter Integer two_jc = 0;
+ parameter Integer two_ma = 0;
+ parameter Integer two_mb = 0;
+ parameter Integer two_mc = 0;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_8_1.gsl_sf_coupling_3j(two_ja, two_jb, two_jc, two_ma, two_mb, two_mc);
+ end gsl_sf_coupling_3j;
+
+ model gsl_sf_coupling_3j_e
+ parameter Integer two_ja = 0;
+ parameter Integer two_jb = 0;
+ parameter Integer two_jc = 0;
+ parameter Integer two_ma = 0;
+ parameter Integer two_mb = 0;
+ parameter Integer two_mc = 0;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_8_1.gsl_sf_coupling_3j_e(two_ja, two_jb, two_jc, two_ma, two_mb, two_mc);
+ end gsl_sf_coupling_3j_e;
+ end chap_7_8_1;
+
+ package chap_7_8_2
+ model gsl_sf_coupling_6j
+ parameter Integer two_ja = 0;
+ parameter Integer two_jb = 0;
+ parameter Integer two_jc = 0;
+ parameter Integer two_ma = 0;
+ parameter Integer two_mb = 0;
+ parameter Integer two_mc = 0;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_8_2.gsl_sf_coupling_6j(two_ja, two_jb, two_jc, two_ma, two_mb, two_mc);
+ end gsl_sf_coupling_6j;
+
+ model gsl_sf_coupling_6j_e
+ parameter Integer two_ja = 0;
+ parameter Integer two_jb = 0;
+ parameter Integer two_jc = 0;
+ parameter Integer two_ma = 0;
+ parameter Integer two_mb = 0;
+ parameter Integer two_mc = 0;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_8_2.gsl_sf_coupling_6j_e(two_ja, two_jb, two_jc, two_ma, two_mb, two_mc);
+ end gsl_sf_coupling_6j_e;
+ end chap_7_8_2;
+
+ package chap_7_8_3
+ model gsl_sf_coupling_9j
+ parameter Integer two_ja = 0;
+ parameter Integer two_jb = 0;
+ parameter Integer two_jc = 0;
+ parameter Integer two_jd = 0;
+ parameter Integer two_je = 0;
+ parameter Integer two_jf = 0;
+ parameter Integer two_jg = 0;
+ parameter Integer two_jh = 0;
+ parameter Integer two_ji = 0;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_8_3.gsl_sf_coupling_9j(two_ja, two_jb, two_jc, two_jd, two_je, two_jf, two_jg, two_jh, two_ji);
+ end gsl_sf_coupling_9j;
+
+ model gsl_sf_coupling_9j_e
+ parameter Integer two_ja = 0;
+ parameter Integer two_jb = 0;
+ parameter Integer two_jc = 0;
+ parameter Integer two_jd = 0;
+ parameter Integer two_je = 0;
+ parameter Integer two_jf = 0;
+ parameter Integer two_jg = 0;
+ parameter Integer two_jh = 0;
+ parameter Integer two_ji = 0;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_8_3.gsl_sf_coupling_9j_e(two_ja, two_jb, two_jc, two_jd, two_je, two_jf, two_jg, two_jh, two_ji);
+ end gsl_sf_coupling_9j_e;
+ end chap_7_8_3;
+
+ package chap_7_9
+ model gsl_sf_dawson
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_9.gsl_sf_dawson(x);
+ end gsl_sf_dawson;
+
+ model gsl_sf_dawson_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_9.gsl_sf_dawson_e(x);
+ end gsl_sf_dawson_e;
+ end chap_7_9;
+
+ package chap_7_10
+ model gsl_sf_debye_1
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_10.gsl_sf_debye_1(x);
+ end gsl_sf_debye_1;
+
+ model gsl_sf_debye_1_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_10.gsl_sf_debye_1_e(x);
+ end gsl_sf_debye_1_e;
+
+ model gsl_sf_debye_2
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_10.gsl_sf_debye_2(x);
+ end gsl_sf_debye_2;
+
+ model gsl_sf_debye_2_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_10.gsl_sf_debye_2_e(x);
+ end gsl_sf_debye_2_e;
+
+ model gsl_sf_debye_3
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_10.gsl_sf_debye_3(x);
+ end gsl_sf_debye_3;
+
+ model gsl_sf_debye_3_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_10.gsl_sf_debye_3_e(x);
+ end gsl_sf_debye_3_e;
+
+ model gsl_sf_debye_4
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_10.gsl_sf_debye_4(x);
+ end gsl_sf_debye_4;
+
+ model gsl_sf_debye_4_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_10.gsl_sf_debye_4_e(x);
+ end gsl_sf_debye_4_e;
+
+ model gsl_sf_debye_5
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_10.gsl_sf_debye_5(x);
+ end gsl_sf_debye_5;
+
+ model gsl_sf_debye_5_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_10.gsl_sf_debye_5_e(x);
+ end gsl_sf_debye_5_e;
+
+ model gsl_sf_debye_6
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_10.gsl_sf_debye_6(x);
+ end gsl_sf_debye_6;
+
+ model gsl_sf_debye_6_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_10.gsl_sf_debye_6_e(x);
+ end gsl_sf_debye_6_e;
+ end chap_7_10;
+
+ package chap_7_11_1
+ model gsl_sf_dilog
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_11_1.gsl_sf_dilog(x);
+ end gsl_sf_dilog;
+
+ model gsl_sf_dilog_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_11_1.gsl_sf_dilog_e(x);
+ end gsl_sf_dilog_e;
+ end chap_7_11_1;
+
+ package chap_7_11_2
+ model gsl_sf_complex_dilog_e
+ parameter Real r = 10;
+ parameter Real theta = 1;
+ gsl.data_types.gsl_sf_result result_re;
+ gsl.data_types.gsl_sf_result result_im;
+ Integer y;
+ algorithm
+ (result_re, result_im, y) := gsl.Special_functions.chap_7_11_2.gsl_sf_complex_dilog_e(r, theta);
+ end gsl_sf_complex_dilog_e;
+ end chap_7_11_2;
+
+ package chap_7_12
+ model gsl_sf_multiply
+ parameter Real x = 2.5;
+ parameter Real y = 2.5;
+ Real z;
+ algorithm
+ z := gsl.Special_functions.chap_7_12.gsl_sf_multiply(x, y);
+ end gsl_sf_multiply;
+
+ model gsl_sf_multiply_e
+ parameter Real x = 2.5;
+ parameter Real y = 2.5;
+ gsl.data_types.gsl_sf_result result;
+ Integer z;
+ algorithm
+ (result, z) := gsl.Special_functions.chap_7_12.gsl_sf_multiply_e(x, y);
+ end gsl_sf_multiply_e;
+
+ model gsl_sf_multiply_err_e
+ parameter Real x = 2;
+ parameter Real dx = 0.5;
+ parameter Real y = 2;
+ parameter Real dy = 0.2;
+ gsl.data_types.gsl_sf_result result;
+ Integer z;
+ algorithm
+ (result, z) := gsl.Special_functions.chap_7_12.gsl_sf_multiply_err_e(x, dx, y, dy);
+ end gsl_sf_multiply_err_e;
+ end chap_7_12;
+
+ package chap_7_15_1
+ model gsl_sf_erf
+ parameter Real x = 2;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_15_1.gsl_sf_erf(x);
+ end gsl_sf_erf;
+
+ model gsl_sf_erf_e
+ parameter Real x = 2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_15_1.gsl_sf_erf_e(x);
+ end gsl_sf_erf_e;
+ end chap_7_15_1;
+
+ package chap_7_15_2
+ model gsl_sf_erfc
+ parameter Real x = 2;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_15_2.gsl_sf_erfc(x);
+ end gsl_sf_erfc;
+
+ model gsl_sf_erfc_e
+ parameter Real x = 2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_15_2.gsl_sf_erfc_e(x);
+ end gsl_sf_erfc_e;
+ end chap_7_15_2;
+
+ package chap_7_15_3
+ model gsl_sf_log_erfc
+ parameter Real x = 2;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_15_3.gsl_sf_log_erfc(x);
+ end gsl_sf_log_erfc;
+
+ model gsl_sf_log_erfc_e
+ parameter Real x = 2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_15_3.gsl_sf_log_erfc_e(x);
+ end gsl_sf_log_erfc_e;
+ end chap_7_15_3;
+
+ package chap_7_15_4
+ model gsl_sf_erf_Z
+ parameter Real x = 2;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_15_4.gsl_sf_erf_Z(x);
+ end gsl_sf_erf_Z;
+
+ model gsl_sf_erf_Z_e
+ parameter Real x = 2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_15_4.gsl_sf_erf_Z_e(x);
+ end gsl_sf_erf_Z_e;
+
+ model gsl_sf_erf_Q
+ parameter Real x = 2;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_15_4.gsl_sf_erf_Q(x);
+ end gsl_sf_erf_Q;
+
+ model gsl_sf_erf_Q_e
+ parameter Real x = 2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_15_4.gsl_sf_erf_Q_e(x);
+ end gsl_sf_erf_Q_e;
+
+ model gsl_sf_hazard
+ parameter Real x = 2;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_15_4.gsl_sf_hazard(x);
+ end gsl_sf_hazard;
+
+ model gsl_sf_hazard_e
+ parameter Real x = 2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_15_4.gsl_sf_hazard_e(x);
+ end gsl_sf_hazard_e;
+ end chap_7_15_4;
+
+ package chap_7_16_1
+ model gsl_sf_exp
+ //this model calls the function gsl_sf_exp by taking value of x as input and returning y as output
+ parameter Real x = 10;
+ Real y(start = 1);
+ algorithm
+ y := gsl.Special_functions.chap_7_16_1.gsl_sf_exp(x);
+ end gsl_sf_exp;
+
+ model gsl_sf_exp_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_16_1.gsl_sf_exp_e(x);
+ end gsl_sf_exp_e;
+
+ model gsl_sf_exp_e10_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result_e10 result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_16_1.gsl_sf_exp_e10_e(x);
+ end gsl_sf_exp_e10_e;
+
+ model gsl_sf_exp_mult
+ //this model calls the function gsl_sf_exp by taking value of x as input and returning y as output
+ parameter Real x = 10;
+ parameter Real y = 10;
+ Real z;
+ algorithm
+ z := gsl.Special_functions.chap_7_16_1.gsl_sf_exp_mult(x, y);
+ end gsl_sf_exp_mult;
+
+ model gsl_sf_exp_mult_e
+ parameter Real x = 10;
+ parameter Real y = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer z;
+ algorithm
+ (result, z) := gsl.Special_functions.chap_7_16_1.gsl_sf_exp_mult_e(x, y);
+ end gsl_sf_exp_mult_e;
+
+ model gsl_sf_exp_mult_e10_e
+ parameter Real x = 10;
+ parameter Real y = 10;
+ gsl.data_types.gsl_sf_result_e10 result;
+ Integer z;
+ algorithm
+ (result, z) := gsl.Special_functions.chap_7_16_1.gsl_sf_exp_mult_e10_e(x, y);
+ end gsl_sf_exp_mult_e10_e;
+ end chap_7_16_1;
+
+ package chap_7_16_2
+ model gsl_sf_expm1
+ parameter Real x = 0.00034;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_16_2.gsl_sf_expm1(x);
+ end gsl_sf_expm1;
+
+ model gsl_sf_expm1_e
+ parameter Real x = 0.000002;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_16_2.gsl_sf_expm1_e(x);
+ end gsl_sf_expm1_e;
+
+ model gsl_sf_exprel
+ parameter Real x = 0.000002;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_16_2.gsl_sf_exprel(x);
+ end gsl_sf_exprel;
+
+ model gsl_sf_exprel_e
+ parameter Real x = 0.002;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_16_2.gsl_sf_exprel_e(x);
+ end gsl_sf_exprel_e;
+
+ model gsl_sf_exprel_2
+ parameter Real x = 0.002;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_16_2.gsl_sf_exprel_2(x);
+ end gsl_sf_exprel_2;
+
+ model gsl_sf_exprel_2_e
+ parameter Real x = 0.002;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_16_2.gsl_sf_exprel_2_e(x);
+ end gsl_sf_exprel_2_e;
+
+ model gsl_sf_exprel_n
+ parameter Integer n = 2;
+ parameter Real x = 0.002;
+ Real z;
+ algorithm
+ z := gsl.Special_functions.chap_7_16_2.gsl_sf_exprel_n(n, x);
+ end gsl_sf_exprel_n;
+
+ model gsl_sf_exprel_n_e
+ parameter Integer n = 2;
+ parameter Real x = 0.2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_16_2.gsl_sf_exprel_n_e(n, x);
+ end gsl_sf_exprel_n_e;
+ end chap_7_16_2;
+
+ package chap_7_16_3
+ model gsl_sf_exp_err_e
+ //This model calls the function gsl_sf_exp_err_e which exponentiates x with an associated absolute error dx
+ parameter Real x = 0.15;
+ parameter Real dx = 0.00015;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_16_3.gsl_sf_exp_err_e(x, dx);
+ end gsl_sf_exp_err_e;
+
+ model gsl_sf_exp_err_e10_e
+ //this model calls the function gsl_sf_exp_err_e10_e which exponentiates x with an associated absolute error dx
+ parameter Real x = 0.15;
+ parameter Real dx = 0.000015;
+ gsl.data_types.gsl_sf_result_e10 result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_16_3.gsl_sf_exp_err_e10_e(x, dx);
+ end gsl_sf_exp_err_e10_e;
+
+ model gsl_sf_exp_mult_err_e
+ parameter Real x = 0.15;
+ parameter Real dx = 0.0000000015;
+ parameter Real y = 1;
+ parameter Real dy = 0.00000001;
+ gsl.data_types.gsl_sf_result result;
+ Integer z;
+ algorithm
+ (result, z) := gsl.Special_functions.chap_7_16_3.gsl_sf_exp_mult_err_e(x, dx, y, dy);
+ end gsl_sf_exp_mult_err_e;
+
+ model gsl_sf_exp_mult_err_e10_e
+ parameter Real x = 0.15;
+ parameter Real dx = 0.0000000015;
+ parameter Real y = 1;
+ parameter Real dy = 0.000000001;
+ gsl.data_types.gsl_sf_result_e10 result;
+ Integer z;
+ algorithm
+ (result, z) := gsl.Special_functions.chap_7_16_3.gsl_sf_exp_mult_err_e10_e(x, dx, y, dy);
+ end gsl_sf_exp_mult_err_e10_e;
+ end chap_7_16_3;
+
+ package chap_7_17_1
+ model gsl_sf_expint_E1
+ parameter Real x = 10;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_17_1.gsl_sf_expint_E1(x);
+ end gsl_sf_expint_E1;
+
+ model gsl_sf_expint_E1_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_17_1.gsl_sf_expint_E1_e(x);
+ end gsl_sf_expint_E1_e;
+
+ model gsl_sf_expint_E2
+ parameter Real x = 10;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_17_1.gsl_sf_expint_E2(x);
+ end gsl_sf_expint_E2;
+
+ model gsl_sf_expint_E2_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_17_1.gsl_sf_expint_E2_e(x);
+ end gsl_sf_expint_E2_e;
+
+ model gsl_sf_expint_En
+ parameter Integer n = 2;
+ parameter Real x = 10;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_17_1.gsl_sf_expint_En(n, x);
+ end gsl_sf_expint_En;
+
+ model gsl_sf_expint_En_e
+ parameter Integer n = 4;
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_17_1.gsl_sf_expint_En_e(n, x);
+ end gsl_sf_expint_En_e;
+ end chap_7_17_1;
+
+ package chap_7_17_2
+ model gsl_sf_expint_Ei
+ parameter Real x = 10;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_17_2.gsl_sf_expint_Ei(x);
+ end gsl_sf_expint_Ei;
+
+ model gsl_sf_expint_Ei_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_17_2.gsl_sf_expint_Ei_e(x);
+ end gsl_sf_expint_Ei_e;
+ end chap_7_17_2;
+
+ package chap_7_17_3
+ model gsl_sf_Shi
+ parameter Real x = 10;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_17_3.gsl_sf_Shi(x);
+ end gsl_sf_Shi;
+
+ model gsl_sf_Shi_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_17_3.gsl_sf_Shi_e(x);
+ end gsl_sf_Shi_e;
+
+
+ model gsl_sf_Chi
+ parameter Real x = 10;
+ Real y;
+ algorithm
+ y := gsl.Special_functions.chap_7_17_3.gsl_sf_Chi(x);
+ end gsl_sf_Chi;
+
+ model gsl_sf_Chi_e
+ parameter Real x = 10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result, y) := gsl.Special_functions.chap_7_17_3.gsl_sf_Chi_e(x);
+ end gsl_sf_Chi_e;
+ end chap_7_17_3;
+
+ package chap_7_17_4
+ model gsl_sf_expint_3
+ parameter Real x=10;
+ Real y;
+ algorithm
+ y:=gsl.Special_functions.chap_7_17_4.gsl_sf_expint_3(x);
+ end gsl_sf_expint_3;
+
+ model gsl_sf_expint_3_e
+ parameter Real x=10;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result,y):=gsl.Special_functions.chap_7_17_4.gsl_sf_expint_3_e(x);
+ end gsl_sf_expint_3_e;
+ end chap_7_17_4;
+
+ package chap_7_17_5
+
+ model gsl_sf_Si
+ parameter Real x=10;
+ Real y;
+ algorithm
+ y:=gsl.Special_functions.chap_7_17_5.gsl_sf_Si(x);
+ end gsl_sf_Si;
+
+ model gsl_sf_Si_e
+ parameter Real x=2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result,y):=gsl.Special_functions.chap_7_17_5.gsl_sf_Si_e(x);
+ end gsl_sf_Si_e;
+
+ model gsl_sf_Ci
+ parameter Real x=2;
+ Real y;
+ algorithm
+ y:=gsl.Special_functions.chap_7_17_5.gsl_sf_Ci(x);
+ end gsl_sf_Ci;
+
+ model gsl_sf_Ci_e
+ parameter Real x=2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result,y):=gsl.Special_functions.chap_7_17_5.gsl_sf_Ci_e(x);
+ end gsl_sf_Ci_e;
+ end chap_7_17_5;
+
+ package chap_7_18_1
+
+
+ model gsl_sf_fermi_dirac_m1
+ parameter Real x=2;
+ Real y;
+ algorithm
+ y:=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_m1(x);
+ end gsl_sf_fermi_dirac_m1;
+
+
+ model gsl_sf_fermi_dirac_m1_e
+ parameter Real x=2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result,y):=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_m1_e(x);
+ end gsl_sf_fermi_dirac_m1_e;
+
+ model gsl_sf_fermi_dirac_0
+ parameter Real x=2;
+ Real y;
+ algorithm
+ y:=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_0(x);
+ end gsl_sf_fermi_dirac_0;
+
+ model gsl_sf_fermi_dirac_0_e
+ parameter Real x=2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result,y):=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_0_e(x);
+ end gsl_sf_fermi_dirac_0_e;
+
+ model gsl_sf_fermi_dirac_1
+ parameter Real x=2;
+ Real y;
+ algorithm
+ y:=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_1(x);
+ end gsl_sf_fermi_dirac_1;
+
+ model gsl_sf_fermi_dirac_1_e
+ parameter Real x=2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result,y):=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_1_e(x);
+ end gsl_sf_fermi_dirac_1_e;
+
+
+
+ model gsl_sf_fermi_dirac_2
+ parameter Real x=2;
+ Real y;
+ algorithm
+ y:=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_2(x);
+ end gsl_sf_fermi_dirac_2;
+
+ model gsl_sf_fermi_dirac_2_e
+ parameter Real x=2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result,y):=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_2_e(x);
+ end gsl_sf_fermi_dirac_2_e;
+
+
+ model gsl_sf_fermi_dirac_int
+ parameter Integer j=2;
+ parameter Real x=2;
+ Real y;
+ algorithm
+ y:=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_int(j,x);
+ end gsl_sf_fermi_dirac_int;
+
+ model gsl_sf_fermi_dirac_int_e
+ parameter Integer j=2;
+ parameter Real x=3;
+ gsl.data_types.gsl_sf_result result;
+ Real y;
+ algorithm
+ (result,y):=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_int_e(j,x);
+ end gsl_sf_fermi_dirac_int_e;
+
+
+ model gsl_sf_fermi_dirac_mhalf
+ parameter Real x=2;
+ Real y;
+ algorithm
+ y:=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_mhalf(x);
+ end gsl_sf_fermi_dirac_mhalf;
+
+
+ model gsl_sf_fermi_dirac_mhalf_e
+ parameter Real x=2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result,y):=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_mhalf_e(x);
+ end gsl_sf_fermi_dirac_mhalf_e;
+
+ model gsl_sf_fermi_dirac_half
+ parameter Real x=2;
+ Real y;
+ algorithm
+ y:=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_half(x);
+ end gsl_sf_fermi_dirac_half;
+
+ model gsl_sf_fermi_dirac_half_e
+ parameter Real x=2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result,y):=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_half_e(x);
+ end gsl_sf_fermi_dirac_half_e;
+
+ model gsl_sf_fermi_dirac_3half
+ parameter Real x=2;
+ Real y;
+ algorithm
+ y:=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_3half(x);
+ end gsl_sf_fermi_dirac_3half;
+
+ model gsl_sf_fermi_dirac_3half_e
+ parameter Real x=2;
+ gsl.data_types.gsl_sf_result result;
+ Integer y;
+ algorithm
+ (result,y):=gsl.Special_functions.chap_7_18_1.gsl_sf_fermi_dirac_3half_e(x);
+ end gsl_sf_fermi_dirac_3half_e;
+
+
+
+
+ end chap_7_18_1;
+
+
+ end special_functions;
end Examples;
end gsl; \ No newline at end of file
generated by cgit (git 2.34.1) at 2024-12-19 10:10:33 +0000