imslfunction1.qxd
|
|
|
- たしろう つねざき
- 5 years ago
- Views:
Transcription
1 VERSION 5.5 C Numerical Library Function Catalog
2 IMSL TM C Numerical Library Function Catalog IMSL C IMSL Java TM Fortran 8 IMSL C Math / Library IMSL C / Stat / Library
3 IMSL C C IMSL C CNL IMSL C - 3 -
4 IMSL TM C Numerical Library Function Catalog CNL CNL CNL CNL C CNL UNIX/Linux Windows - 4 -
5 IMSL TM C Numerical Library Function Catalog CNL CNL CNL CNL IMSL
6 IMSL TM C Numerical Library Function Catalog IMSL C C/ C Hermite 3 2 Gauss Adams-Gear Runge-Kutta 1 2 Fourier Laplace Bessel 50 CPU
7 IMSL TM C Numerical Library Function Catalog 12 Wilcoxon 2 Shapiro-Wilk - Yates lattice split-plot strip-plot split-split plot strip-split plot Wilcoxon Cochran 2 Kolmogorov/ Smirnov ARMA GARCH (Generalized Autoregressive Conditional Heteroscedasticity) kalman K-means 2 2 Cox Kaplan-Meier 2 2 Poisson Faure - 7 -
8 IMSL TM C Numerical Library Function Catalog IMSL JAVA TM Fortran JMSL TM - JAVA JMSL Java 100 Java Java JMSL IMSL Java JMSL Java JMSL JMSL 100 Java JMSL JMSL Java PC Java system latency IMSL FORTRAN IMSL Fortran IMSL Fortran 5.0 IMSL F90 IMSL Fortran 77 IMSL IMSL IMSL Fortran 30 IMSL IMSL IMSL Fortran IMSL Fortran - 8 -
9 IMSL C/ Math/ Library TM lin_sol_gen Ax = b lin_sol_gen ( ) Ax = b lin_sol_posdef Ax = b lin_sol_ posdef ( ) Hermitian Ax = b lin_sol_gen_band Ax = b lin_sol_gen_band ( ) Ax = b in_sol_posdef _band Ax = b lin_sol_posdef _band ( ) Hermitian Ax = b lin_sol_gen_coordinate Ax = b lin_sol_gen_coordinate ( ) A Hermitian Ax = b lin_sol_posdef_coordinate Ax = b lin_sol_posdef_coordinate ( ) Hermitian Ax = b lin_sol_gen_min_residual (GMRES) Ax = b - 9 -
10 IMSL C/ Math/ Library TM lin_sol_def_cg lin_least_squares_gen lin_lsq_lin_constraints lin_svd_gen 2 Ax = b 2 A SVD A = USV T lin_svd_gen ( ) A SVD A = USV T lin_sol_nonnegdef Ax = b eig_gen A eig_gen ( ) A eig_sym A eig_herm ( ) Hermitian A eig_symgen geneig A B B Ax = Bx A B Ax = Bx geneig ( ) A B Ax = Bx cub_spline_interp_e_cnd cub_spline_interp_shape
11 IMSL C/ Math/ Library TM cub_spline_value cub_spline_integral spline_interp spline_knots spline_2d_interp 2 2 spline_value spline_integral spline_2d_value spline_2d_integral 1 1 user_fcn_least_squares spline_least_squares spline_2d_least_squares cub_spline_smooth spline_lsq_constrained smooth_1d_data scattered_2d_interp
12 IMSL C/ Math/ Library TM radial_scattered_fit n 1 R n radial_evaluate int_fcn_sing int_fcn int_fcn_sing_pts int_fcn_alg_log int_fcn_inf int_fcn_trig int_fcn_fourier int_fcn_cauchy int_fcn_smooth Gauss-Kronrod Gauss-Kronrod - Fourier b f ( x) dx x c Cauchy a int_fcn_2d 2 ( x0,..., x ) dx... dx int_fcn_hyper_rect n 1 n 1 0 int_fcn_qmc Monte-Carlo b0 a0... bn 1 an 1 gauss_quad_rule Gauss Gauss-Radau fcn_derivative
13 IMSL C/ Math/ Library TM RUNGE-KUTTA ode_runge_kutta Runge-Kutta-Verner 5 6 ADAM GEAR ode_adams_gear Adams-Gear pde_method_of_lines ut = f(x, t, u, ux, uxx) bvp_finite_difference 1 POISSON fast_poisson_2d HODIE Poisson 2 Poisson Helmholtz fft_real fft_real_init imsl_f_fft_real fft_complex fft_complex_init imsl_c_fft_complex fft_cosine
14 IMSL C/ Math/ Library TM fft_cosine_init fft_sine fft_sine_init imsl_f_fft_cosine imsl_f_fft_sine fft_2d_complex 2 2 convolution 2 convolution ( ) 2 LAPLACE inverse_laplace Laplace zeros_poly Jenkins -Traub 3 zeros_poly ( ) Jenkins -Traub 3 zeros_fcn Muller zeros_sys_eqn Powell n f (x) =
15 IMSL C/ Math/ Library TM min_uncon min_uncon_deriv min_uncon_multivar nonlin_least_squares 1 f(x) 1 1 f(x) Newton n f(x) Levenberg-Marquardt 2 lin_prog quadratic_prog min_con_gen_lin bounded_least_squares / 2 / Levenberg-Marquardt 2 Constrained_nlp erf erfc erf(x) erf(x) erf_inverse erfce erfe erfc_inverse beta erf(x)-1 erfc(z) erfc-1(x) beta (x, y)
16 IMSL C/ Math/ Library TM log_beta beta_incomplete gamma log_gamma gamma_incomplete beta ln (x, y) beta Ix = x(a,b)/ (a,b) gamma (x) gamma log (x) gamma (a, x) BESSEL bessel_j0 bessel_j1 bessel_jx bessel_y0 bessel_y1 bessel_yx bessel_i0 bessel_exp_i0 bessel_i1 bessel_exp_i1 bessel_ix bessel_k0 bessel_exp_k0 bessel_k1 bessel_exp_k1 bessel_kx 1 0 Bessel J0(x) 1 1 Bessel J1(x) 1 Bessel 2 0 Bessel Y0(x) 2 1 Bessel Y1(x) 2 Bessel 1 0 Bessel I0(x) 1 0 Bessel 1 1 Bessel I1(x) 1 1 Bessel 1 Bessel 3 1 Bessel K0(x) 3 0 Bessel 3 1 Bessel K1(x) 3 1 Bessel 3 Bessel
17 IMSL C/ Math/ Library TM elliptic_integral_k elliptic_integral_e elliptic_integral_rf elliptic_integral_rd elliptic_integral_rj elliptic_integral_rc 1 K(x) 2 E(x) 1 Carlson R F(x, y, z) 2 Carlson R D(x, y, z) 3 Carlson R J(x, y, z) FRESNEL fresnel_integral_c fresnel_integral_s Fresnel Fresnel AIRY airy_ai airy_bi airy_ai_derivative airy_bi_derivative Airy 2 Airy Airy 2 Airy KELVIN kelvin_ber0 kelvin_bei0 kelvin_ker0 kelvin_kei0 kelvin_ber0_derivative kelvin_bei0_derivative kelvin_ker0_derivative kelvin_kei0_derivative 0 1 Kelvin ber 0 1 Kelvin bei 0 2 Kelvin ker 0 2 Kelvin kei 0 1 Kelvin ber 0 1 Kelvin bei 0 2 Kelvin ker 0 2 Kelvin kei
18 IMSL C/ Math/ Library TM normal_cdf normal_inverse_cdf chi_squared_cdf chi_squared_inverse_cdf F_cdf F_inverse_cdf t_cdf t_inverse_cdf gamma_cdf binomial_cdf hypergeometric_cdf poisson_cdf beta_cdf beta_inverse_cdf bivariate_normal_cdf (Gaussian) (Gaussian) 2 2 F F Student t Student t gamma 2 Poisson beta beta 2 cumulative_interest cumulative_principal depreciation_db depreciation_ddb depreciation_sln depreciation_syd depreciation_vdb
19 IMSL C/ Math/ Library TM dollar_decimal dollar_fraction effective_rate future_value future_value_schedule interest_payment interest_rate_annuity internal_rate_of_return internal_rate_schedule modified_internal_rate net_present_value nominal_rate number_of_periods payment present_value present_value_schedule principal_payment 1 accr_interest_maturity accr_inter est_periodic bond_equivalent_yield
20 IMSL C/ Math/ Library TM convexity coupon_days coupon_number days_before_settlement days_to_next_coupon depreciation_amordegrc depreciation_amorlinc discount_price depreciation_amordegrc depreciation_amordegrc 100 discount_rate discount_yield duration interest_rate_security modified_duration next_coupon_date Macauley 100 previous_coupon_date price 100 price_maturity received_maturity 100 treasury_bill_price 100 treasury_bill_yield year_fraction
21 IMSL C/ Math/ Library TM yield_maturity yield_periodic simple_statistics table_oneway chi_squared_test covariances regression poly_regression ranks random_seed_get random_seed_set random_option random_uniform random_normal random_poisson random_gamma random_beta random_exponential faure_next_point IMSL IMSL (0, 1) (0, 1) CDF Poisson gamma beta Faure
22 IMSL C/ Math/ Library TM write_matrix page write_options output_file version ctime date_to_days days_to_date CPU error_options error_code constant machine (integer) machine (float)
23 IMSL C/ Math/ Library TM sort sort (integer) vector_norm. 2 mat_mul_rect 1 3 mat_mul_rect ( ) 1 3 mat_mul_rect_band mat_mul_rect_band ( ) mat_mul_rect_coordinate mat_mul_rect_coordinate ( ) mat_add_band 2 C A B mat_add_band ( ) 2 C A B mat_add_coordinate 2 C A B mat_add_coordinate ( ) 2 C A B matrix_norm matrix_norm_band matrix_norm_coordinate generate_test_band E (n,c) generate_test_band ( ) E (n,c)
24 IMSL C/ Math/ Library TM generate_test_coordinate D(n,c) E(n,c) generate_test_coordinate ( ) D(n,c) E(n,c) c_neg c_add c_sub c_mul c_div c_eq cz_convert zc_convert cf_convert c_conjg c_abs c_arg c_sqrt c_cos c_sin c_exp c_log cf_power cc_power fi_power ii_power
25 IMSL C/ Stat/ Library TM simple_statistics normal_one_sample normal_two_sample table_oneway table_twoway sort_data ranks 1 2 regressors_for_glm regression 2 regression_summary regression_prediction hypothesis_partial
26 IMSL C/ Stat/ Library TM hypothesis_scph hypothesis_test regression_selection regression_stepwise poly_regression poly_prediction nonlinear_regression nonlinear_optimization 2 Powell 2 Lnorm_regression 2 Lp covariances partial_covariances pooled_covariances robust_covariances
27 IMSL C/ Stat/ Library TM anova_oneway anova_factorial anova_nested anova_balanced 1 crd_factorial rcbd_factorial latin_square lattice split_plot split_split_plot strip_plot strip_split_plot split-plot split-split-plot strip-plot strip-split-plot homogeneity multiple_comparisons Yates Bartlett Levene SNK LSD Tukey Duncan Bonferroni Yate contingency_table
28 IMSL C/ Stat/ Library TM exact_enumeration exact_network 2 2 categorical_glm Probit Poisson sign_test wilcoxon_sign_rank noether_cyclical_trend cox_stuart_trends_test tie_statistics Wilcoxon Noether Cox Stuarts 1 wilcoxon_rank_sum kruskal_wallis_test friedmans_test cochran_q_test k_trends_test Wilcoxon Kruskal Wallis Friedman Cochran Q K chi_squared_test normality_test
29 IMSL C/ Stat/ Library TM kolmogorov_one kolmogorov_two multivar_normality_test Kolmogorov -Smirnov 1 Kolmogorov -Smirnov 2 Mardia randomness_test ARIMA arma arma_forecast difference ARMA 2 ARMA: ARMA box_cox_transform autocorrelation crosscorrelation multi_crosscorrelation partial_autocorrelation lack_of_fit Box -Cox 2 2 GARCH garch GARCH (p, q) kalman
30 IMSL C/ Stat/ Library TM dissimilarities cluster_hierarchical cluster_number k-means cluster_k_means K-means principal_components factor_analysis discriminant_analysis aplan_meier_estimates rop_hazards_gen_lin survival_glm survival_estimates Kaplan-Meier Cox nonparam_hazard_rate life_tables
31 IMSL C/ Stat/ Library TM binomial_cdf binomial_pdf hypergeometric_cdf hypergeometric_pdf poisson_cdf poisson_pdf Poisson 2 2 Poisson beta_cdf beta_inverse_cdf bivariate_normal_cdf chi_squared_cdf chi_squared_inverse_cdf non_central_chi_sq non_central_chi_sq_inv F_cdf F_inverse_cdf gamma_cdf normal_cdf normal_inverse_cdf t_cdf t_inverse_cdf F F Gauss Gauss Student t Student t
32 IMSL C/ Stat/ Library TM non_central_t_cdf non_central_t_inv_cdf Student t Student t random_binomial random_geometric random_hypergeometric random_logarithmic random_neg_binomial random_poisson random_uniform_discrete 2 2 Poisson random_general_discrete random_beta random_cauchy random_chi_squared random_exponential random_exponential_mix random_gamma random_lognormal random_normal random_stable Chaucy 2 CDF
33 IMSL C/ Stat/ Library TM random_student_t random_triangular random_uniform random_von_mises random_weibull random_general_continuous continuous_table_setup student-t 0, 1 Von Mises Weibull random_normal_multivariate random_orthogonal_matrix random_mvar_from_data random_multinomial random_sphere random_table_twoway K 2 random_order_normal random_order_uniform (0, 1) random_arma random_npp ARMA random_permutation random_sample_indices random_sample
34 IMSL C/ Stat/ Library TM random_option random_option_get random_seed_get 0, 1 0, 1 IMSL random_substream_seed_get random_seed_set random_table_set random_table_get random_gfsr_table_set random_gfsr_table_get 100,000 IMSL GFSR GFSR faure_next_point Faure write_matrix page write_options output_file version
35 IMSL C/ Stat/ Library TM error_options error_code machine (integer) machine (float) data_sets mat_mul_rect permute_vector permute_matrix binomial_coefficient beta beta_incomplete log_beta gamma gamma_incomplete log_gamma ctime CPU
CNL社内製本用.PDF
VERSION 6.0 C Numerical Library Function Catalog IMSL TM C IMSL C Math / Library IMSL C / Stat / Library 1 IMSL TM C IMSL C IMSL C IMSL C IMSL C IMSL C IMSL C C IMSL C UNIX/Linux Windows 3 IMSL C IMSL
fortranfunction2.qxd
Fortran Numerical Library Function Catalog IMSL TM Fortran Numerical Library Function Catalog IMSL Fortran 4 8 9 10 IMSL Java TM C 12 IMSL Math / Library 13 13 21 23 27 28 29 31 32 34 40 41 IMSL Math /
JMSL 3.0 FC_日_ .doc
V E R S I O N 3. 0 JMSL Numerical Library Function Catalog JMSL Numerical Library Function Catalog Java 4 7 IMSLC C# Fortran 8 com.imsl Package 9 com.imsl.math Package 10 10 10 10 11 11 11 12 12 12 12
1 1 1.1...................................... 1 1.2................................... 5 1.3................................... 7 1.4............................. 9 1.5....................................
18 ( ) I II III A B C(100 ) 1, 2, 3, 5 I II A B (100 ) 1, 2, 3 I II A B (80 ) 6 8 I II III A B C(80 ) 1 n (1 + x) n (1) n C 1 + n C
8 ( ) 8 5 4 I II III A B C( ),,, 5 I II A B ( ),, I II A B (8 ) 6 8 I II III A B C(8 ) n ( + x) n () n C + n C + + n C n = 7 n () 7 9 C : y = x x A(, 6) () A C () C P AP Q () () () 4 A(,, ) B(,, ) C(,,
PowerPoint プレゼンテーション
40 6 Fischer 2 McNemar 2 No Yes No Yes Mann Whitney s U Wilcoxon 1 Kruskal Wallis 2 Friedman No Yes 2 Yes No Yes Bonferroni ANOVA No Scheffe Yes Student s T No Paired T Welch s T 40 6 Fischer 2 McNemar
総合薬学講座 生物統計の基礎
2013 10 22 ( ) 2013 10 22 1 / 40 p.682 1. 2. 3 2 t Mann Whitney U ). 4 χ 2. 5. 6 Dunnett Tukey. 7. 8 Kaplan Meier.. U. ( ) 2013 10 22 2 / 40 1 93 ( 20 ) 230. a t b c χ 2 d 1.0 +1.0 e, b ( ) e ( ) ( ) 2013
: 1g99p038-8
16 17 : 1g99p038-8 1 3 1.1....................................... 4 1................................... 5 1.3.................................. 5 6.1..................................... 7....................................
7 27 7.1........................................ 27 7.2.......................................... 28 1 ( a 3 = 3 = 3 a a > 0(a a a a < 0(a a a -1 1 6
26 11 5 1 ( 2 2 2 3 5 3.1...................................... 5 3.2....................................... 5 3.3....................................... 6 3.4....................................... 7
1 4 1 ( ) ( ) ( ) ( ) () 1 4 2
7 1995, 2017 7 21 1 2 2 3 3 4 4 6 (1).................................... 6 (2)..................................... 6 (3) t................. 9 5 11 (1)......................................... 11 (2)
会津若松市民の運動行動における変化ステージの評価
2 Bandura1977 2000 Bandura 50 36 3 4 MarcusSimkin1993 2002 6 6 6 6 1995 2000.81 1998 1995 2000.91 2000 5 .65 1986.83 Marcus et al1992 2002 Pros.81 Cons.80 Pros Cons Pros Cons 6 7 8 9 10 11 12 13 5 Tukey
-1-1 1 1 1 1 12 31 2 2 3 4
2007 -1-1 1 1 1 1 12 31 2 2 3 4 -2-5 6 CPU 3 Windows98 1 -3-2. 3. -4-4 2 5 1 1 1 -5- 50000 50000 50000 50000 50000 50000 50000 50000 50000 50000-6- -7-1 Windows 2 -8-1 2 3 4 - - 100,000 200,000 500,000
Bessel ( 06/11/21) Bessel 1 ( ) 1.1 0, 1,..., n n J 0 (x), J 1 (x),..., J n (x) I 0 (x), I 1 (x),..., I n (x) Miller (Miller algorithm) Bess
Bessel 5 3 11 ( 6/11/1) Bessel 1 ( ) 1.1, 1,..., n n J (x), J 1 (x),..., J n (x) I (x), I 1 (x),..., I n (x) Miller (Miller algorithm) Bessel (6 ) ( ) [1] n n d j J n (x), d j I n (x) Deuflhard j= j=.1
(pdf) (cdf) Matlab χ ( ) F t
(, ) (univariate) (bivariate) (multi-variate) Matlab Octave Matlab Matlab/Octave --...............3. (pdf) (cdf)...3.4....4.5....4.6....7.7. Matlab...8.7.....9.7.. χ ( )...0.7.3.....7.4. F....7.5. t-...3.8....4.8.....4.8.....5.8.3....6.8.4....8.8.5....8.8.6....8.9....9.9.....9.9.....0.9.3....0.9.4.....9.5.....0....3
1.2 y + P (x)y + Q(x)y = 0 (1) y 1 (x), y 2 (x) y 1 (x), y 2 (x) (1) y(x) c 1, c 2 y(x) = c 1 y 1 (x) + c 2 y 2 (x) 3 y 1 (x) y 1 (x) e R P (x)dx y 2
1 1.1 R(x) = 0 y + P (x)y + Q(x)y = R(x)...(1) y + P (x)y + Q(x)y = 0...(2) 1 2 u(x) v(x) c 1 u(x)+ c 2 v(x) = 0 c 1 = c 2 = 0 c 1 = c 2 = 0 2 0 2 u(x) v(x) u(x) u (x) W (u, v)(x) = v(x) v (x) 0 1 1.2
5 Armitage x 1,, x n y i = 10x i + 3 y i = log x i {x i } {y i } 1.2 n i i x ij i j y ij, z ij i j 2 1 y = a x + b ( cm) x ij (i j )
5 Armitage. x,, x n y i = 0x i + 3 y i = log x i x i y i.2 n i i x ij i j y ij, z ij i j 2 y = a x + b 2 2. ( cm) x ij (i j ) (i) x, x 2 σ 2 x,, σ 2 x,2 σ x,, σ x,2 t t x * (ii) (i) m y ij = x ij /00 y
Introduction to Numerical Analysis of Differential Equations Naoya Enomoto (Kyoto.univ.Dept.Science(math))
Introduction to Numerical Analysis of Differential Equations Naoya Enomoto (Kyoto.univ.Dept.Science(math)) 2001 1 e-mail:s00x0427@ip.media.kyoto-u.ac.jp 1 1 Van der Pol 1 1 2 2 Bergers 2 KdV 2 1 5 1.1........................................
確率論と統計学の資料
5 June 015 ii........................ 1 1 1.1...................... 1 1........................... 3 1.3... 4 6.1........................... 6................... 7 ii ii.3.................. 8.4..........................
1 Edward Waring Lagrange n {(x i, y i )} n i=1 x i p i p i (x j ) = δ ij P (x) = p i p i (x) = n y i p i (x) (1) i=1 n j=1 j i x x j x i x j (2) Runge
Edwrd Wring Lgrnge n {(x i, y i )} n x i p i p i (x j ) = δ ij P (x) = p i p i (x) = n y i p i (x) () n j= x x j x i x j (2) Runge [] [2] = ξ 0 < ξ < ξ n = b [, b] [ξ i, ξ i ] y = /( + 25x 2 ) 5 2,, 0,,
Stata 11 Stata ts (ARMA) ARCH/GARCH whitepaper mwp 3 mwp-083 arch ARCH 11 mwp-051 arch postestimation 27 mwp-056 arima ARMA 35 mwp-003 arima postestim
TS001 Stata 11 Stata ts (ARMA) ARCH/GARCH whitepaper mwp 3 mwp-083 arch ARCH 11 mwp-051 arch postestimation 27 mwp-056 arima ARMA 35 mwp-003 arima postestimation 49 mwp-055 corrgram/ac/pac 56 mwp-009 dfgls
68 A mm 1/10 A. (a) (b) A.: (a) A.3 A.4 1 1
67 A Section A.1 0 1 0 1 Balmer 7 9 1 0.1 0.01 1 9 3 10:09 6 A.1: A.1 1 10 9 68 A 10 9 10 9 1 10 9 10 1 mm 1/10 A. (a) (b) A.: (a) A.3 A.4 1 1 A.1. 69 5 1 10 15 3 40 0 0 ¾ ¾ É f Á ½ j 30 A.3: A.4: 1/10
36 3 D f(z) D z f(z) z Taylor z D C f(z) z C C f (z) C f(z) f (z) f(z) D C D D z C C 3.: f(z) 3. f (z) f 2 (z) D D D D D f (z) f 2 (z) D D f (z) f 2 (
3 3. D f(z) D D D D D D D D f(z) D f (z) f (z) f(z) D (i) (ii) (iii) f(z) = ( ) n z n = z + z 2 z 3 + n= z < z < z > f (z) = e t(+z) dt Re z> Re z> [ ] f (z) = e t(+z) = (Rez> ) +z +z t= z < f(z) Taylor
x h = (b a)/n [x i, x i+1 ] = [a+i h, a+ (i + 1) h] A(x i ) A(x i ) = h 2 {f(x i) + f(x i+1 ) = h {f(a + i h) + f(a + (i + 1) h), (2) 2 a b n A(x i )
![x h = (b a)/n [x i, x i+1 ] = [a+i h, a+ (i + 1) h] A(x i ) A(x i ) = h 2 {f(x i) + f(x i+1 ) = h {f(a + i h) + f(a + (i + 1) h), (2) 2 a b n A(x i )](/thumbs/92/107866678.jpg "x h = (b a)/n [x i, x i+1 ] = [a+i h, a+ (i + 1) h] A(x i ) A(x i ) = h 2 {f(x i) + f(x i+1 ) = h {f(a + i h) + f(a + (i + 1) h), (2) 2 a b n A(x i )")
1 f(x) a b f(x)dx = n A(x i ) (1) ix [a, b] n i A(x i ) x i 1 f(x) [a, b] n h = (b a)/n y h = (b-a)/n y = f (x) h h a a+h a+2h a+(n-1)h b x 1: 1 x h = (b a)/n [x i, x i+1 ] = [a+i h, a+ (i + 1) h] A(x
A
A 2563 15 4 21 1 3 1.1................................................ 3 1.2............................................. 3 2 3 2.1......................................... 3 2.2............................................
( ) C 2 f(x) = 0 (1) x (1) x 3 3x 2 + 9x 8 = 0 (2) 1
2004 8 25 1 1.1 ( ) 4 2 1.2 4 C 2 f(x) = 0 (1) x (1) x 3 3x 2 + 9x 8 = 0 (2) 1 [ ] 1/3 [ 2 1 ( x 1 = 1 2 1 + + 1 + ) ] 1/3 33 33 2 x 2 = 1 + (1 + [ ] 1/3 2 3i) 1 + 1 33 2 (1 [ 1 3i) 2 x 2 = 1 + (1 [ ]
V E R S I O N C# Numerical Library Function Catalog
V E R S I O N 4.0. 2 C# Numerical Library Function Catalog 目次 1 IMSL C# ライブラリ4.02 概要 2 Imsl.Math Package 5 基本関数 5 線形代数 5 固有値解析 5 補間と近似 5 積分 5 微分方程式 6 変換 6 非線形方程式 6 最適化 6 特殊関数 6 その他 7 印刷関数 7 Imsl.Stat Package
() x + y + y + x dy dx = 0 () dy + xy = x dx y + x y ( 5) ( s55906) 0.7. (). 5 (). ( 6) ( s6590) 0.8 m n. 0.9 n n A. ( 6) ( s6590) f A (λ) = det(a λi)
0. A A = 4 IC () det A () A () x + y + z = x y z X Y Z = A x y z ( 5) ( s5590) 0. a + b + c b c () a a + b + c c a b a + b + c 0 a b c () a 0 c b b c 0 a c b a 0 0. A A = 7 5 4 5 0 ( 5) ( s5590) () A ()
( ) kadai4, kadai4.zip.,. 3 cos x [ π, π] Python. ( 100 ), x cos x ( ). (, ). def print cos(): print cos()
![( ) kadai4, kadai4.zip.,. 3 cos x [ π, π] Python. ( 100 ), x cos x ( ). (, ). def print cos(): print cos()](/thumbs/91/105030926.jpg "( ) kadai4, kadai4.zip.,. 3 cos x [ π, π] Python. ( 100 ), x cos x ( ). (, ). def print cos(): print cos()")
4 2010.6 1 :, HP.. HP 4 (, PGM/PPM )., python,,, 2, kadai4,.,,, ( )., ( ) N, exn.py ( 3 ex3.py ). N 3.., ( )., ( ) N, (exn.txt).. 1 ( ) kadai4, kadai4.zip.,. 3 cos x [ π, π] Python. ( 100 ), x cos x (
b3e2003.dvi
15 II 5 5.1 (1) p, q p = (x + 2y, xy, 1), q = (x 2 + 3y 2, xyz, ) (i) p rotq (ii) p gradq D (2) a, b rot(a b) div [11, p.75] (3) (i) f f grad f = 1 2 grad( f 2) (ii) f f gradf 1 2 grad ( f 2) rotf 5.2
一般演題（ポスター）
6 5 13 : 00 14 : 00 A μ 13 : 00 14 : 00 A β β β 13 : 00 14 : 00 A 13 : 00 14 : 00 A 13 : 00 14 : 00 A β 13 : 00 14 : 00 A β 13 : 00 14 : 00 A 13 : 00 14 : 00 A β 13 : 00 14 : 00 A 13 : 00 14 : 00 A
n Y 1 (x),..., Y n (x) 1 W (Y 1 (x),..., Y n (x)) 0 W (Y 1 (x),..., Y n (x)) = Y 1 (x)... Y n (x) Y 1(x)... Y n(x) (x)... Y n (n 1) (x) Y (n 1)
D d dx 1 1.1 n d n y a 0 dx n + a d n 1 y 1 dx n 1 +... + a dy n 1 dx + a ny = f(x)...(1) dk y dx k = y (k) a 0 y (n) + a 1 y (n 1) +... + a n 1 y + a n y = f(x)...(2) (2) (2) f(x) 0 a 0 y (n) + a 1 y
S I. dy fx x fx y fx + C 3 C dy fx 4 x, y dy v C xt y C v e kt k > xt yt gt [ v dt dt v e kt xt v e kt + C k x v + C C k xt v k 3 r r + dr e kt S dt d
S I.. http://ayapin.film.s.dendai.ac.jp/~matuda /TeX/lecture.html PDF PS.................................... 3.3.................... 9.4................5.............. 3 5. Laplace................. 5....
ax 2 + bx + c = n 8 (n ) a n x n + a n 1 x n a 1 x + a 0 = 0 ( a n, a n 1,, a 1, a 0 a n 0) n n ( ) ( ) ax 3 + bx 2 + cx + d = 0 4
20 20.0 ( ) 8 y = ax 2 + bx + c 443 ax 2 + bx + c = 0 20.1 20.1.1 n 8 (n ) a n x n + a n 1 x n 1 + + a 1 x + a 0 = 0 ( a n, a n 1,, a 1, a 0 a n 0) n n ( ) ( ) ax 3 + bx 2 + cx + d = 0 444 ( a, b, c, d
all.dvi
5,, Euclid.,..,... Euclid,.,.,, e i (i =,, ). 6 x a x e e e x.:,,. a,,. a a = a e + a e + a e = {e, e, e } a (.) = a i e i = a i e i (.) i= {a,a,a } T ( T ),.,,,,. (.),.,...,,. a 0 0 a = a 0 + a + a 0
Release Notes for JMP book
8.0.2 JMP, A Business Unit of SAS SAS Campus Drive Cary, NC 27513 SAS Institute Inc. 2009.JMP 8.0.2, Cary, NC: SAS Institute Inc. JMP 8.0.2 Copyright 2009, SAS Institute Inc., Cary, NC, USA All rights
u = u(t, x 1,..., x d ) : R R d C λ i = 1 := x 2 1 x 2 d d Euclid Laplace Schrödinger N := {1, 2, 3,... } Z := {..., 3, 2, 1,, 1, 2, 3
2 2 1 5 5 Schrödinger i u t + u = λ u 2 u. u = u(t, x 1,..., x d ) : R R d C λ i = 1 := 2 + + 2 x 2 1 x 2 d d Euclid Laplace Schrödinger 3 1 1.1 N := {1, 2, 3,... } Z := {..., 3, 2, 1,, 1, 2, 3,... } Q
Microsoft Word - 触ってみよう、Maximaに2.doc
i i e! ( x +1) 2 3 ( 2x + 3)! ( x + 1) 3 ( a + b) 5 2 2 2 2! 3! 5! 7 2 x! 3x! 1 = 0 ",! " >!!! # 2x + 4y = 30 "! x + y = 12 sin x lim x!0 x x n! # $ & 1 lim 1 + ('% " n 1 1 lim lim x!+0 x x"!0 x log x
2.2 ( y = y(x ( (x 0, y 0 y (x 0 (y 0 = y(x 0 y = y(x ( y (x 0 = F (x 0, y(x 0 = F (x 0, y 0 (x 0, y 0 ( (x 0, y 0 F (x 0, y 0 xy (x, y (, F (x, y ( (
(. x y y x f y = f(x y x y = y(x y x y dx = d dx y(x = y (x = f (x y = y(x x ( (differential equation ( + y 2 dx + xy = 0 dx = xy + y 2 2 2 x y 2 F (x, y = xy + y 2 y = y(x x x xy(x = F (x, y(x + y(x 2
I , : ~/math/functional-analysis/functional-analysis-1.tex
I 1 2004 8 16, 2017 4 30 1 : ~/math/functional-analysis/functional-analysis-1.tex 1 3 1.1................................... 3 1.2................................... 3 1.3.....................................
FD 2 ( )..................... 3 FD 14 ( )..................... 15 FD 39 ( ).................... 40................. 54................................
14 15 3 FD 2 ( )..................... 3 FD 14 ( )..................... 15 FD 39 ( ).................... 40................. 54.......................................... 66..........................................
A11 (1993,1994) 29 A12 (1994) 29 A13 Trefethen and Bau Numerical Linear Algebra (1997) 29 A14 (1999) 30 A15 (2003) 30 A16 (2004) 30 A17 (2007) 30 A18
2013 8 29y, 2016 10 29 1 2 2 Jordan 3 21 3 3 Jordan (1) 3 31 Jordan 4 32 Jordan 4 33 Jordan 6 34 Jordan 8 35 9 4 Jordan (2) 10 41 x 11 42 x 12 43 16 44 19 441 19 442 20 443 25 45 25 5 Jordan 26 A 26 A1
C による数値計算法入門 ( 第 2 版 ) 新装版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 新装版 1 刷発行時のものです.
C による数値計算法入門 ( 第 2 版 ) 新装版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009383 このサンプルページの内容は, 新装版 1 刷発行時のものです. i 2 22 2 13 ( ) 2 (1) ANSI (2) 2 (3) Web http://www.morikita.co.jp/books/mid/009383
P1-1 P1-2 P1-3 P1-4 P1-5 P1-6 P3-1 P3-2 P3-3 P3-4 P3-5 P3-6 P5-1 P5-2 P5-3 P5-4 P5-5 P5-6 P7-1 P7-2 P7-3 P7-4 P7-5 P7-6 P9-1 P9-2 P9-3 P9-4 P9-5 P9-6 P11-1 P11-2 P11-3 P11-4 P13-1 P13-2 P13-3 P13-4 P13-5
2S III IV K A4 12:00-13:30 Cafe David 1 2 TA 1 appointment Cafe David K2-2S04-00 : C
2S III IV K200 : April 16, 2004 Version : 1.1 TA M2 TA 1 10 2 n 1 ɛ-δ 5 15 20 20 45 K2-2S04-00 : C 2S III IV K200 60 60 74 75 89 90 1 email 3 4 30 A4 12:00-13:30 Cafe David 1 2 TA 1 email appointment Cafe
漸化式のすべてのパターンを解説しましたー高校数学の達人・河見賢司のサイト
https://www.hmg-gen.com/tuusin.html https://www.hmg-gen.com/tuusin1.html 1 2 OK 3 4 {a n } (1) a 1 = 1, a n+1 a n = 2 (2) a 1 = 3, a n+1 a n = 2n a n a n+1 a n = ( ) a n+1 a n = ( ) a n+1 a n {a n } 1,
−g”U›ß™ö‡Æ…X…y…N…g…‰
1 / 74 ( ) 2019 3 8 URL: http://www.math.kyoto-u.ac.jp/ ichiro/ 2 / 74 Contents 1 Pearson 2 3 Doob h- 4 (I) 5 (II) 6 (III-1) - 7 (III-2-a) 8 (III-2-b) - 9 (III-3) Pearson 3 / 74 Pearson Definition 1 ρ
A S- hara/lectures/lectures-j.html r A = A 5 : 5 = max{ A, } A A A A B A, B A A A %
A S- http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html r A S- 3.4.5. 9 phone: 9-8-444, e-mail: hara@math.kyushu-u.ac.jp, http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html Office
solutionJIS.dvi
May 0, 006 6 morimune@econ.kyoto-u.ac.jp /9/005 (7 0/5/006 1 1.1 (a) (b) (c) c + c + + c = nc (x 1 x)+(x x)+ +(x n x) =(x 1 + x + + x n ) nx = nx nx =0 c(x 1 x)+c(x x)+ + c(x n x) =c (x i x) =0 y i (x
I ( ) 1 de Broglie 1 (de Broglie) p λ k h Planck ( Js) p = h λ = k (1) h 2π : Dirac k B Boltzmann ( J/K) T U = 3 2 k BT
I (008 4 0 de Broglie (de Broglie p λ k h Planck ( 6.63 0 34 Js p = h λ = k ( h π : Dirac k B Boltzmann (.38 0 3 J/K T U = 3 k BT ( = λ m k B T h m = 0.067m 0 m 0 = 9. 0 3 kg GaAs( a T = 300 K 3 fg 07345
I L01( Wed) : Time-stamp: Wed 07:38 JST hig e, ( ) L01 I(2017) 1 / 19
I L01(2017-09-20 Wed) : Time-stamp: 2017-09-20 Wed 07:38 JST hig e, http://hig3.net ( ) L01 I(2017) 1 / 19 ? 1? 2? ( ) L01 I(2017) 2 / 19 ?,,.,., 1..,. 1,2,.,.,. ( ) L01 I(2017) 3 / 19 ? I. M (3 ) II,
2
2 3 Page 4 5 6 A-1B-1 C0 D0 E0 F0 G0 A0 B0 C1 D1 E1 F1 G1 A1 B1 C2 D2 E2 F2 G2 A2 B2 C3 D3 E3 7 F3 G3 A3 B3 C4 D4 E4 F4 G4 A4 B4 C5 D5 E5 F5 G5 A5 B5 C6 D6 E6 F6 G6 A6 B6 C7 8 Page 9 1 2 3 1 2 10 1 11
1 (bit ) ( ) PC WS CPU IEEE754 standard ( 24bit) ( 53bit)
GNU MP BNCpack tkouya@cs.sist.ac.jp 2002 9 20 ( ) Linux Conference 2002 1 1 (bit ) ( ) PC WS CPU IEEE754 standard ( 24bit) ( 53bit) 10 2 2 3 4 5768:9:; = %? @BADCEGFH-I:JLKNMNOQP R )TSVU!" # %$ & " #
takei.dvi
0 Newton Leibniz ( ) α1 ( ) αn (1) a α1,...,α n (x) u(x) = f(x) x 1 x n α 1 + +α n m 1957 Hans Lewy Lewy 1970 1 1.1 Example 1.1. (2) d 2 u dx 2 Q(x)u = f(x), u(0) = a, 1 du (0) = b. dx Q(x), f(x) x = 0
( ) p.2/70
p.1/70 ( ) p.2/70 ( p.3/70 (1)20 1970 cf.box and Jenkins(1970) (2)1980 p.4/70 Kolmogorov(1940),(1941) Granger(1966) Hurst(1951) Hurst p.5/70 Beveridge Wheat Prices Index p.6/70 Nile River Water Level p.7/70
1 filename=mathformula tex 1 ax 2 + bx + c = 0, x = b ± b 2 4ac, (1.1) 2a x 1 + x 2 = b a, x 1x 2 = c a, (1.2) ax 2 + 2b x + c = 0, x = b ± b 2
filename=mathformula58.tex ax + bx + c =, x = b ± b 4ac, (.) a x + x = b a, x x = c a, (.) ax + b x + c =, x = b ± b ac. a (.3). sin(a ± B) = sin A cos B ± cos A sin B, (.) cos(a ± B) = cos A cos B sin
6. Euler x
...............................................................................3......................................... 4.4................................... 5.5......................................
p12.dvi
301 12 (2) : 1 (1) dx dt = f(x,t) ( (t 0,t 1,...,t N ) ) h k = t k+1 t k. h k k h. x(t k ) x k. : 2 (2) :1. step. 1 : explicit( ) : ξ k+1 = ξ k +h k Ψ(t k,ξ k,h k ) implicit( ) : ξ k+1 = ξ k +h k Ψ(t k,t
> > <., vs. > x 2 x y = ax 2 + bx + c y = 0 2 ax 2 + bx + c = 0 y = 0 x ( x ) y = ax 2 + bx + c D = b 2 4ac (1) D > 0 x (2) D = 0 x (3
13 2 13.0 2 ( ) ( ) 2 13.1 ( ) ax 2 + bx + c > 0 ( a, b, c ) ( ) 275 > > 2 2 13.3 x 2 x y = ax 2 + bx + c y = 0 2 ax 2 + bx + c = 0 y = 0 x ( x ) y = ax 2 + bx + c D = b 2 4ac (1) D >
FX ) 2
(FX) 1 1 2009 12 12 13 2009 1 FX ) 2 1 (FX) 2 1 2 1 2 3 2010 8 FX 1998 1 FX FX 4 1 1 (FX) () () 1998 4 1 100 120 1 100 120 120 100 20 FX 100 100 100 1 100 100 100 1 100 1 100 100 1 100 101 101 100 100
time2print4.dvi
iii 2 P. J. Brockwell and R.A. Davis, Introduction to Time Series and Forecasting, 2nd edition (Springer, 2002) 1 2 ITSM2000(version 7) 6 10 7 2 2 2 ITSM2000 student version professional version CD-ROM
情報理論 第5回 情報量とエントロピー
5 () ( ) ( ) ( ) p(a) a I(a) p(a) p(a) I(a) p(a) I(a) (2) (self information) p(a) = I(a) = 0 I(a) = 0 I(a) a I(a) = log 2 p(a) = log 2 p(a) bit 2 (log 2 ) (3) I(a) 7 6 5 4 3 2 0 0.5 p(a) p(a) = /2 I(a)
6 6.1 B A: Γ d Q S(B) S(A) = S (6.1) T (e) Γ (6.2) : Γ B A R (reversible) 6-1
6 6.1 B A: Γ d Q S(B) S(A) = S (6.1) (e) Γ (6.2) : Γ B A R (reversible) 6-1 (e) = Clausius 0 = B A: Γ B A: Γ d Q A + d Q (e) B: R d Q + S(A) S(B) (6.3) (e) // 6.2 B A: Γ d Q S(B) S(A) = S (6.4) (e) Γ (6.5)
Microsoft Word - 03-数値計算の基礎.docx
δx f x 0 + δ x n=0 a n = f ( n) ( x 0 ) n δx n f x x=0 sin x = x x3 3 + x5 5 x7 7 +... x ( ) = a n δ x n ( ) = sin x ak = (-mod(k,2))**(k/2) / fact_k 10 11 I = f x dx a ΔS = f ( x)h I = f a h I = h b (
(1) (2) (3) (4) HB B ( ) (5) (6) (7) 40 (8) (9) (10)
2017 12 9 4 1 30 4 10 3 1 30 3 30 2 1 30 2 50 1 1 30 2 10 (1) (2) (3) (4) HB B ( ) (5) (6) (7) 40 (8) (9) (10) (1) i 23 c 23 0 1 2 3 4 5 6 7 8 9 a b d e f g h i (2) 23 23 (3) 23 ( 23 ) 23 x 1 x 2 23 x
A B C A B C X Y Z
Argonauta 8: 27-37 (2003) 1 ANOVA, analysis of variance t 1980 90 ANOVA ANOVA Underwood, 1997 ANOVA ANOVA ANOVA 1975 Sokal & Rohlf 1981 1990 Underwood 1997 1997 Zar 1999 1-way ANOVA 2-way, 3-way, ANOVA,
卓球の試合への興味度に関する確率論的分析
17 i 1 1 1.1..................................... 1 1.2....................................... 1 1.3..................................... 2 2 5 2.1................................ 5 2.2 (1).........................
II
II 16 16.0 2 1 15 x α 16 x n 1 17 (x α) 2 16.1 16.1.1 2 x P (x) P (x) = 3x 3 4x + 4 369 Q(x) = x 4 ax + b ( ) 1 P (x) x Q(x) x P (x) x P (x) x = a P (a) P (x) = x 3 7x + 4 P (2) = 2 3 7 2 + 4 = 8 14 +
A
A05-132 2010 2 11 1 1 3 1.1.......................................... 3 1.2..................................... 3 1.3..................................... 3 2 4 2.1............................... 4 2.2
2 p T, Q
270 C, 6000 C, 2 p T, Q p: : p = N/ m 2 N/ m 2 Pa : pdv p S F Q 1 g 1 1 g 1 14.5 C 15.5 1 1 cal = 4.1855 J du = Q pdv U ( ) Q pdv 2 : z = f(x, y). z = f(x, y) (x 0, y 0 ) y y = y 0 z = f(x, y 0 ) x x =
4 2016 3 8 2.,. 2. Arakawa Jacobin., 2 Adams-Bashforth. Re = 80, 90, 100.. h l, h/l, Kármán, h/l 0.28,, h/l.., (2010), 46.2., t = 100 t = 2000 46.2 < Re 46.5. 1 1 4 2 6 2.1............................
d dt A B C = A B C d dt x = Ax, A 0 B 0 C 0 = mm 0 mm 0 mm AP = PΛ P AP = Λ P A = ΛP P d dt x = P Ax d dt (P x) = Λ(P x) d dt P x =
3 MATLAB Runge-Kutta Butcher 3. Taylor Taylor y(x 0 + h) = y(x 0 ) + h y (x 0 ) + h! y (x 0 ) + Taylor 3. Euler, Runge-Kutta Adams Implicit Euler, Implicit Runge-Kutta Gear y n+ y n (n+ ) y n+ y n+ y n+
June 2016 i (statistics) F Excel Numbers, OpenOffice/LibreOffice Calc ii *1 VAR STDEV 1 SPSS SAS R *2 R R R R *1 Excel, Numbers, Microsoft Office, Apple iwork, *2 R GNU GNU R iii URL http://ruby.kyoto-wu.ac.jp/statistics/training/
Design of highly accurate formulas for numerical integration in weighted Hardy spaces with the aid of potential theory 1 Ken ichiro Tanaka 1 Ω R m F I = F (t) dt (1.1) Ω m m 1 m = 1 1 Newton-Cotes Gauss
( 30 ) 30 4 5 1 4 1.1............................................... 4 1.............................................. 4 1..1.................................. 4 1.......................................
ランダムウォークの境界条件・偏微分方程式の数値計算
B L06(2018-05-22 Tue) : Time-stamp: 2018-05-22 Tue 21:53 JST hig,, 2, multiply transf http://hig3.net L06 B(2018) 1 / 38 L05-Q1 Quiz : 1 M λ 1 = 1 u 1 ( ). M u 1 = u 1, u 1 = ( 3 4 ) s (s 0)., u 1 = 1
2/50 Auction: Theory and Practice 3 / 50 (WTO) 10 SDR ,600 Auction: Theory and Practice 4 / 50 2
stakagi@econ.hokudai.ac.jp June 24, 2011 2.... 3... 4... 7 8... 9.... 10... 11... 12 IPV 13 SPSB... 15 SPSB.... 17 SPSB.... 19 FPSB... 20 FPSB.... 22 FPSB.... 23... 24 Low Price Auction.... 27 APV 29...
1 1 [1] ( 2,625 [2] ( 2, ( ) /
![1 1 [1] ( 2,625 [2] ( 2, ( ) /](/thumbs/95/122470987.jpg "1 1 [1] ( 2,625 [2] ( 2, ( ) /")
[] (,65 [] (,3 ( ) 67 84 76 7 8 6 7 65 68 7 75 73 68 7 73 7 7 59 67 68 65 75 56 6 58 /=45 /=45 6 65 63 3 4 3/=36 4/=8 66 7 68 7 7/=38 /=5 7 75 73 8 9 8/=364 9/=864 76 8 78 /=45 /=99 8 85 83 /=9 /= ( )
webkaitou.dvi
( c Akir KANEKO) ).. m. l s = lθ m d s dt = mg sin θ d θ dt = g l sinθ θ l θ mg. d s dt xy t ( d x dt, d y dt ) t ( mg sin θ cos θ, sin θ sin θ). (.) m t ( d x dt, d y dt ) = t ( mg sin θ cos θ, mg sin
/02/18
3 09/0/8 i III,,,, III,?,,,,,,,,,,,,,,,,,,,,?,?,,,,,,,,,,,,,,!!!,? 3,,,, ii,,,!,,,, OK! :!,,,, :!,,,,,, 3:!,, 4:!,,,, 5:!,,! 7:!,,,,, 8:!,! 9:!,,,,,,,,, ( ),, :, ( ), ( ), 6:!,,, :... : 3 ( )... iii,,
2 1 Introduction (1.1.2) Logistic ث Malthus (1.1.3) (( ) ث)( ) α = ( ) ( + ) [Verhulst 1845] 0 ( ) ( + ) lim ( ) = 0 t (1.1.4) (( ) ث)( ) α = ( ) Logi
![2 1 Introduction (1.1.2) Logistic ث Malthus (1.1.3) (( ) ث)( ) α = ( ) ( + ) [Verhulst 1845] 0 ( ) ( + ) lim ( ) = 0 t (1.1.4) (( ) ث)( ) α = ( ) Logi](/thumbs/81/82814896.jpg "2 1 Introduction (1.1.2) Logistic ث Malthus (1.1.3) (( ) ث)( ) α = ( ) ( + ) [Verhulst 1845] 0 ( ) ( + ) lim ( ) = 0 t (1.1.4) (( ) ث)( ) α = ( ) Logi")
1 1 Introduction 1.1 ( ) ( ) Malthus ( ) (1.1.1) ( ) α = ( ) ( + ) 0 ( ) ( + ) lim ( ) = 0 t ( ) α = ( ) αt م = ( ) ( ) = ( 0 ) 0 (1.1.2) αt م 0 = ( ) (1.1.1) ( ) ( α + (1 = ( + ) = 0 1 ( ( ) 2... ) geometric
seminar0220a.dvi
1 Hi-Stat 2 16 2 20 16:30-18:00 2 2 217 1 COE 4 COE RA E-MAIL: ged0104@srv.cc.hit-u.ac.jp 2004 2 25 S-PLUS S-PLUS S-PLUS S-code 2 [8] [8] [8] 1 2 ARFIMA(p, d, q) FI(d) φ(l)(1 L) d x t = θ(l)ε t ({ε t }