diff options
| -rw-r--r-- | a.mo | 5599 |
1 files changed, 3026 insertions, 2573 deletions
@@ -418,2014 +418,2349 @@ 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; + 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_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 + 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_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 + 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_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; - 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; - 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_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 - 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_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 - 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; - 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; + 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_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; + 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; - 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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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; - 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; + 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_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_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_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_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_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; + 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; - 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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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; + 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; - 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_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_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_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_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_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_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_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_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_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_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_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; - 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; + package chap_7_5_7 + //regular modified spherical bessel functions - 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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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_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; - 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; + 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_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_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_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_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_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_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_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_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_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_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_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_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_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; - -end Special_functions; + 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; + + /* + 7.17.3 Hyperbolic Integrals + double gsl_sf_Shi(double x) + 50 + Chapter 7. Special FunctionsGNU Scientific Library, Release 2.5 + int gsl_sf_Shi_e(double x, gsl_sf_result * result) + These routines compute the integral + x + Z + dt sinh(t)/t + Shi(x) = + 0 + double gsl_sf_Chi(double x) + int gsl_sf_Chi_e(double x, gsl_sf_result * result) + These routines compute the integral + + Z + Chi(x) := < γ E + log(x) + + x + + dt(cosh(t) − 1)/t + 0 + where γ E is the Euler constant (available as the macro M_EULER). + */ + annotation( + Documentation); + end chap_7_17_3; + end Special_functions; package Examples package Mathematical @@ -2712,11 +3047,9 @@ end Special_functions; 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 + //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 @@ -3125,39 +3458,37 @@ end Special_functions; 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; + 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); + (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); + 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); + 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 @@ -3166,7 +3497,7 @@ end Special_functions; 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; @@ -3175,7 +3506,7 @@ end Special_functions; 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; @@ -3183,7 +3514,7 @@ end Special_functions; 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; @@ -3192,7 +3523,7 @@ end Special_functions; 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; @@ -3200,7 +3531,7 @@ end Special_functions; 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; @@ -3209,7 +3540,7 @@ end Special_functions; 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; @@ -3218,35 +3549,35 @@ end Special_functions; 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; + 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); + (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); + 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 + 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; @@ -3254,14 +3585,14 @@ end Special_functions; 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; @@ -3269,13 +3600,14 @@ end Special_functions; 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; @@ -3283,26 +3615,26 @@ end Special_functions; 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 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 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 Integer lmax = 10; parameter Real x = 4; Real result_array[lmax]; Integer y; @@ -3310,14 +3642,15 @@ end Special_functions; (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 + 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; @@ -3325,14 +3658,14 @@ end Special_functions; 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; @@ -3340,13 +3673,14 @@ end Special_functions; 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; @@ -3354,26 +3688,26 @@ end Special_functions; 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 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 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 Integer lmax = 10; parameter Real x = 4; Real result_array[lmax]; Integer y; @@ -3381,82 +3715,73 @@ end Special_functions; (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; + 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; + + 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 + 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; + 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; + + 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 + 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; + 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; + + 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 + 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; + + 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 x = 10; parameter Real nu = 1.5; gsl.data_types.gsl_sf_result result; Integer y; @@ -3464,129 +3789,117 @@ end Special_functions; (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; + 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; + + 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 + 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 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; + + 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 + 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 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 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; + 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; - + 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; + 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); + 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; + 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; + 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); + (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; + 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); + 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; + 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); + 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; + end chap_7_5_13; + package chap_7_6 - model gsl_sf_clausen + 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; @@ -3595,108 +3908,103 @@ end Special_functions; (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; + 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 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); + (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; + 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); + 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 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); + (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; + 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; + 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; + 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]; @@ -3704,495 +4012,640 @@ end Special_functions; 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; - - + 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; + 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; + 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; - - + + 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; - - + 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; + 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; + + 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; + 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; + + 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; + 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; + + 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; + 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; + + 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; + 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; + + 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; + 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; + + 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; + 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; + + 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; + 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; + + 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; + + 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; - - - - + 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; + 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; + 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; - + 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; - + 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; + end chap_7_17_3; end special_functions; end Examples; end gsl;
\ No newline at end of file |