NUMERICAL EVALUATION OF H-FUNCTION BY CONTINUED FRACTIONS ABSTRACT INTRODUCTION

Similar documents
A Modified ADM for Solving Systems of Linear Fredholm Integral Equations of the Second Kind

New Expansion and Infinite Series

On the Decomposition Method for System of Linear Fredholm Integral Equations of the Second Kind

APPROXIMATE LIMIT CYCLES FOR THE RAYLEIGH MODEL

Application of Exp-Function Method to. a Huxley Equation with Variable Coefficient *

The Riemann-Lebesgue Lemma

Solutions of Klein - Gordan equations, using Finite Fourier Sine Transform

Fredholm Integral Equations of the First Kind Solved by Using the Homotopy Perturbation Method

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı

Im λ + D + C Re λ A B η = λ 1 ln λ. Im η E + F + B + A + C + D + A D. C Re η A C

Some estimates on the Hermite-Hadamard inequality through quasi-convex functions

Definite integral. Mathematics FRDIS MENDELU

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

A basic logarithmic inequality, and the logarithmic mean

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Quadrature Rules for Evaluation of Hyper Singular Integrals

Modification Adomian Decomposition Method for solving Seventh OrderIntegro-Differential Equations

1 The Lagrange interpolation formula

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

Realistic Method for Solving Fully Intuitionistic Fuzzy. Transportation Problems

DISTRIBUTION OF SUB AND SUPER HARMONIC SOLUTION OF MATHIEU EQUATION WITHIN STABLE ZONES

The Solution of Volterra Integral Equation of the Second Kind by Using the Elzaki Transform

A Companion of Ostrowski Type Integral Inequality Using a 5-Step Kernel with Some Applications

arxiv: v1 [math.na] 23 Apr 2018

Czechoslovak Mathematical Journal, 55 (130) (2005), , Abbotsford. 1. Introduction

AN INEQUALITY OF OSTROWSKI TYPE AND ITS APPLICATIONS FOR SIMPSON S RULE AND SPECIAL MEANS. I. Fedotov and S. S. Dragomir

Research Article Numerical Treatment of Singularly Perturbed Two-Point Boundary Value Problems by Using Differential Transformation Method

arxiv: v1 [math.ca] 1 Jan 2019

An iterative method for solving nonlinear functional equations

#A11 INTEGERS 11 (2011) NEW SEQUENCES THAT CONVERGE TO A GENERALIZATION OF EULER S CONSTANT

Arithmetic Mean Derivative Based Midpoint Rule

Undergraduate Research

Jordan Journal of Mathematics and Statistics (JJMS) 11(1), 2018, pp 1-12

Math& 152 Section Integration by Parts

5.2 Volumes: Disks and Washers

A Numerical Method for Solving Nonlinear Integral Equations

Construction and Selection of Single Sampling Quick Switching Variables System for given Control Limits Involving Minimum Sum of Risks

A Brief Note on Quasi Static Thermal Stresses In A Thin Rectangular Plate With Internal Heat Generation

The mth Ratio Convergence Test and Other Unconventional Convergence Tests

Math 31S. Rumbos Fall Solutions to Assignment #16

CLOSED EXPRESSIONS FOR COEFFICIENTS IN WEIGHTED NEWTON-COTES QUADRATURES

Harmonic Mean Derivative - Based Closed Newton Cotes Quadrature

Continuous Random Variables

A Computational Method for Solving Linear Volterra Integral Equations

The Predom module. Predom calculates and plots isothermal 1-, 2- and 3-metal predominance area diagrams. Predom accesses only compound databases.

Section 6.1 INTRO to LAPLACE TRANSFORMS

S. S. Dragomir. 2, we have the inequality. b a

An approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95

Review of Gaussian Quadrature method

Numerical quadrature based on interpolating functions: A MATLAB implementation

AN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

Ordinary differential equations

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.

FUNCTIONS OF α-slow INCREASE

An inverse steady state thermal stresses in a thin clamped circular plate with internal heat generation

arxiv: v1 [math.gm] 30 Dec 2015

FUZZY HOMOTOPY CONTINUATION METHOD FOR SOLVING FUZZY NONLINEAR EQUATIONS

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature

Orthogonal Polynomials

Exact solutions for nonlinear partial fractional differential equations

Applicable Analysis and Discrete Mathematics available online at

A new algorithm for generating Pythagorean triples 1

Chapter 5 : Continuous Random Variables

Journal of Inequalities in Pure and Applied Mathematics

37 Kragujevac J. Math. 23 (2001) A NOTE ON DENSITY OF THE ZEROS OF ff-orthogonal POLYNOMIALS Gradimir V. Milovanović a and Miodrag M. Spalević

Frobenius numbers of generalized Fibonacci semigroups

Chapter 1. Basic Concepts

Application Chebyshev Polynomials for Determining the Eigenvalues of Sturm-Liouville Problem

Travelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing

Parametrized inequality of Hermite Hadamard type for functions whose third derivative absolute values are quasi convex

DOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES

Composite Mendeleev s Quadratures for Solving a Linear Fredholm Integral Equation of The Second Kind

ON PREY-PREDATOR MODEL WITH HOLLING-TYPE II AND LESLIE-GOWER SCHEMES AHMED BUSERI ASHINE

On the interval Legendre polynomials

Chapter D05 Integral Equations

RELATIONS ON BI-PERIODIC JACOBSTHAL SEQUENCE

Generalized Techniques in Numerical Integration

Fractional Newton-Raphson Method Accelerated with Aitken s Method

MAC-solutions of the nonexistent solutions of mathematical physics

Efficient Computation of a Class of Singular Oscillatory Integrals by Steepest Descent Method

Research Article Mixed Initial-Boundary Value Problem for Telegraph Equation in Domain with Variable Borders

arxiv: v1 [math.na] 11 Jun 2015

GENERALIZATIONS OF WEIGHTED TRAPEZOIDAL INEQUALITY FOR MONOTONIC MAPPINGS AND ITS APPLICATIONS. (b a)3 [f(a) + f(b)] f x (a,b)

Bailey [1] established a simple but very useful identity: If

WHEN IS A FUNCTION NOT FLAT? 1. Introduction. {e 1 0, x = 0. f(x) =

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

The Hadamard s inequality for quasi-convex functions via fractional integrals

On Error Sum Functions Formed by Convergents of Real Numbers

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Application of Kudryashov method for the Ito equations

Predict Global Earth Temperature using Linier Regression

Adomian Decomposition Method with Green s. Function for Solving Twelfth-Order Boundary. Value Problems

Method of stationary phase

A Nonclassical Collocation Method For Solving Two-Point Boundary Value Problems Over Infinite Intervals

Numerical Analysis: Trapezoidal and Simpson s Rule

Lecture 21: Order statistics

KRASNOSEL SKII TYPE FIXED POINT THEOREM FOR NONLINEAR EXPANSION

8 Laplace s Method and Local Limit Theorems

Transcription:

NMERICAL EVALATION OF H-FNCTION BY CONTINED FRACTIONS B.S. Rn & H.S. Dhmi Deptt. Of Mthemtics niversity of Kumun S.S.J. Cmpus, Almor (ttrnchl INDIA- 6360 ABSTRACT In the present pper n ttempt hs been mde to evlute for different vlues of prmeters m, n, p, q nd the vrible in the rnge 0. to 0.0 by the ppliction of continued frctions INTRODCTION Fox H-function is one of the most generlized function mong ll specil functions. It covers wider rnge nd gives deeper, more generl nd more useful results directly pplicble in vrious problems rising in physicl nd biologicl sciences, engineering nd sttistics. With n im to provide ssistnce to those mthemticins, engineers, scientists nd sttisticins who discover tht they need to generte numericl vlue of the specil function in the course of solving their problems, Lozier nd Olver [3] hve prepred report which consist reviews of pcges, librries nd systems usble in the evlution of clssicl functions but not for more generlized functions lie E-function, G-function, H-function etc. Lozier [] hs developed softwre test device for use in

testing the ccurcy or numericl precision of mthemticl softwre for elementry specil functions. Specil functions cn be evluted by continued frctions using qd scheme. Such type of wor ws initited long bc in 95. Detiled ccount of which cn be seen to the boo of P. Henrici [6]. Historicl survey of continued frctions nd Pde pproximtion hs been compiled in the wor of Brezinsi []. A code of evlute modified Bessel s functions bsed on the continued frction method hs been investigted by Segur et. l [] nd for higher order function, the device hs been developed by Rtis nd Crdob [8]. Rtios of contiguous hypergeometric functions of the type 3 F hve been represented in the form of continued frctions by Singh [7]. A set of fst codes to clculte Bessel functions of integrl nd frctionl orders, modified Bessel functions etc., bsed on the continued frctions method hs been presented by Monrel et. l [5]. We re ming n ttempt to evlute more generlized function from mong specil functions, the H- function by the ppliction of continued frction for vlues of prmeters in the rnges :- m vrying from to n hs been fixed s 0 p vrying from 0 to q vrying from to nd vrible vrying in the rnges 0. to 0.0

3. GENERAL FORMLATION Appliction of initil conditions to the qd schemes ssocited with confluent hypergeoemtric function, Guss hypergeometric function nd Bessel s function hve yielded results, which hs been compiled in the boo of Henirici [6] but does not throw ny light on the evlution of generlized function lie G-function, H-function etc. We hve generted following results under qd scheme for the H- function q e q ( n ( n ( n ( n γ ( n n 03,,,, ( n( n γ ( n ( n (,,3--------- ( γ ( n 3( n ( 3,, --------- V W (. q q e ( n ( n ( n ( n( n δ ( γ n n 03,,,, ( n ( γ δ ( γ n 3( γ n ( 3,, --------- ( γ n δ ( γ n ( γ n (,,3--------- V W (.

Corresponding contined frction for the function cn be obtined with the help of following results:- q 0 e 0 q 0 The division lgorithm for the H-function bsed hs been produced s Jν( Jν ( ( ν ( ν ( ν ( v 6 ( v 8 V W (. 3

5. NMERICAL EVALATION OF THE FNCTION In the present section we hve defined vrious H- functions by RIT method nd then qd - scheme hs been utilized in order to tbulte the function for vlue of, nd. Defining the function L NM H 3 [ /(,0 0, 3 O QP,] V W by RIT method ---------(. nd then ting qd Schemes ( ( where,,3,----------------- V W ----------(. we cn tbulte the functions by ming the substitutings 0., 0., 0.3, 0.,0.5,0.6,0.7,0.8,0.9,.0 0., 0.3,0.,0.5,.0 in the from of tble (Appendix-I

6 The function H cn be defined with the help of RIT 0 0,, (,,(, - method in the following form:- e e e 6 5 3 3 5 3 3/ π π (.3 we cn tbulted the function by ming the substitutings 8 V W b gb g b g ( (,.3 - - - --, -------(. 0.,0.,0.3,.0.,0.5,0.6,0.7,0.8,0.9,.0,.0,3.0,.0,5.0,,,3,, in the form of tble (Appendix-II

7 The function H (,,(,,(, (, 0, 0, the help of RIT - method in the following form cn be defined with bπg L NM e 3 L NM 5 3 5 3 5 bπg e F e I H K 3 5 6 O QP O P Q 5 P (. nd then ting qd schemes b 8 gb g 6 ( 6 b ( g where 3,,, V W ( 6. we cn tbulte the functions by ming the substituting 0.,0.,0.3,0.,0.5,0.6,0.7,0.8,0.9,.0,, In the form of tble (Appendix-III

8 The function H,, cn be defined with the help of RIT method in the following form Γ( 0,, 3 5 6 7 Γ( Γ( 5 6 V W ( 7. nd then ting qd schemes s ( ( ( ( ( 3,, V W ( 8. we cn tbulte the function by ming the substitutings ccording condition z< 0.,0.,0.3,0.,0.5,0.6,0.7,0.8,0.9,3/, /, in the form of tble (Appendix-IV

9 The function H 0,, method in the following form (, (, (, cn be defined with the help of RIT π o L NM cos( π( 3 5 6 Γ( cos( π πγ( 3 5 πγ( t O QP V W ( 9. nd then ting qd schemes s ( ( ( 3,, V W ( 0. we cn tbulte the function by ming the substituting 0.,0.,0.6,0.8,.0,.0,3.0,.0,5.0,6.0,7.0,8.0,9.0,0.0 0,, in the form of tble (Appendix-V

0 0, (, - The function H, (, (, cn be defined with the help of RIT method in the following form e L NM 3 5 6 e O QP V W e 3 5 6 (. nd then ting qd schemes s F HG I K JF HG ( ( ( 3,, we cn tbulte the function by ming the substituting 0.,0.,0.6,0.8,.0,.0,3.0,.0,5.0,6.0,7.0,8.0,9.0,0.0 0,, in the form of tble (Appendix-VI I K J V W (.

Appliction of division lgorithm to the expnsion of H 0 0,, (,,(, by RIT method yields. 3 ( ( ( ( ( ( ( ( ( 6 (( (( (( (( ( 6 (( (( (( (( ( J J J ----(.3 The tble cn be constructed by ming the substitutings 0.,0.,0.3,0.,0.5,0.6,0.7,0.8,0.9,.0,.0,3.0,.0,5.0,,3, s compiled in (ppendix - VII

Above defined proceedure, when pplied to the H (,,(,,(, (, 0, 0, produces following expression :- J ( ( L NM l q l q l q l q ( ( ( ( 6 O QP J ( ( ( ( ( ( 3 l q l q l q l q --(. we cn tbulted the function by ming the substitutings 0.,0.,0.3,.0.,0.5,0.6,0.7,0.8,0.9,.0,, in the form of tble (Appendix-VIII

3 Appedix-I 0. 0.3 0. 0.5.0 0. 0.5709 0.99 0.360 0.86 0.090 0. 0.5933 0.787 0.300 0.366 0.637 0.3 0.58 0.959 0.576 0.057 0. 0. 0.5580 0.938 0.66 0.39 0.68 0.5 0.580 0.8 0.596 0.88 0.303 0.6 0.955 0.68 0.73 0.5 0.39 0.7 0.63 0.09 0.3970 0.5 0.376 0.8 0.97 0.7 0.09 0.08 0.359 0.9 0.3980 0.395 0.3897 0.3857 0.3659.0 0.3678 0.3678 0.3678 0.3678 0.3678

Appedix-II 0. 0. 0.3 0. 0.5 0.6 0.7 0.8 0.9.0.0 3.0.0 5.0 3 0.555 0.083 0.0 0. 0.6 0.5 0.7 0.573 0.56 0.978 0.9 0.700 0.58 0.98 0.3337 0.79 0.36 0.0 0.3509 0.8 0.565 0.578 0.0 0.90 0.63 0.88 0.789.99 0.67 0.996 0.57.9 0.66 0.38 0.565.70 0.77 0.797 0.5075.97 0.75 0.890 0.66.876 0.37 0.597 0.6599.0 0.089 0.997 0.5568.96 0.07 0.77 0.5375.978

5 Appedix-III,,,, 0. (0 0. (- 0.78 (- 0. 0. (0 0.39 (0 0.87 (- 0.78 0.3 (0 0.7 (0 0.50 (- 0.3 0. (0 0.556 (0 0.75 (- 0.6 0.5 (0 0.56 (0 0.85 (- 0.7730 0.6 (0 0.78 (0 0.887 (- 0.8868 0.7 (0 0.538 (0 0.098 (0 0.076 0.8 (0 0.57 (0 0.63 (0 0.56 0.9 (0 0.580 (0 0.386 (0 0..0 (00.0 (0 0.87 (0 0.0

6 Appedix-IV,, 3, 3,, 0. 0.880 0.36 0.059 0.000 0. 0.80 0.7 0.6 0.000 0.3 0.3693 0.577 0.03 0.3000 0. 0.606 0.63 0.93 0.000 0.5 0.56 0.707 0.3989 0.5000 0.6 0.6909 0.775 0.535 0.6000 0.7 0.868 0.8366 0.9708 0.7000 0.8.83 0.89.009 0.8000 0.9.695 0.986.6057 0.9000

7 Appedix-V 0 0. 0. 0.6 0.8.0.0 3.0.0 5.0 6.0 7.0 8.0 9.0 0.0 (0 0.69 (- -0.98 (-3 0.787 (0 0.6958 (0-0.895 (- 0.338 (0 0.7786 (- -0.803 (- 0.886 (0 0.875 (- -0.7 (- 0.637 (0 0.9890 (- -0.9978 (- 0.89 (0.9 (0-0.3603 (0 0.998 (0.68 (0 -.036 (0 0.800 (0 9.50 (0 -.7586 (0.757 ( 0.05 (0-6.6068 (0 8. ( 0.559 ( 0.863 ( 0. (3 0.378 ( 0.8 ( 0.687 (3 0.380 (3 0.50 (3 0.899 (3 0.8876 (3 0.39 (3 0.588 ( 0.80 (3 0.87 ( 0.07

8 Appedix-VI 0 0. 0. 0.6 0.8.0.0 3.0.0 5.0 6.0 7.0 8.0 9.0 0.0 (0.387 (0 5.093 (3 0.08 (0 0.809 (0.056 ( 0.86 (0 0.5735 (0.770 (0 9.0870 (0 0. (0 0.859 (0.550 (0 0.36 (0 0.5667 (0.583 (- 0.8736 (0 0.9 (0 0.337 (- 0.69 (- 0.39 (-0.736 (- 0.869 (- 0.067 (-0.937 (- 0.886 (- 0.3 (-0.563 (-3 0.9756 (- 0.7 (-0.77 (-3 0.3338 (-3 0.3788 (-30.5503 (-3 0.5 (-3 0.89 (-30.797 (- 0.3953 (-3 0.3 (-0.5980 (- 0.0 (-0.537 (-0.057

9 Appedix-VII,,, 3, 0. 0. 0.3 0. 0.5 0.6 0.7 0.8 0.9.0.0 3.0.0 5.0 0.09 0.0090 0.0093 0.000 0.69 0.039 0.0338 0.0030 0.95 0.073 0.0688 0.003 0.68 0.60 0.073 0.05 0.83 0.96 0.7 0.058 0.3070 0.593 0.8 0.0835 0.387 0.367 0.3 0.59 0.383 0.3 0.57 0.803 0.3059 0.965 0..753 0.79 0.38 0.5767 0.38 0.358-0.3700.588-0.70.907 -.098 0.93-3.786.75 -.5885-0.583-6.353 5.860 -.7 -.670-8.556 6.5

0 Appedix-VIII,,,, 0. (- 0.03 (- 0.757 (- 0.03 0. (- -0.936 (- 0. (- -0.387 0.3 (- -0.679 (- 0.668 (- - 0.08 0. (0-0.68 (- 0.76 (- -0.673 0.5 (0-0.7 (- 0.780 (- -0.8607 0.6 (0-0.80 (- 0.638 (0-0.368 0.7 (0-0.7 (- 0.559 (0-0.99 0.8 (0-0.393 (- 0.385 (0-0.555 0.9 (0-0.363 (- 0.09 (0-0.360.0 (0-0.397 (- -0.660 (0-0.397

REFRENCES [] C. Brezinsi, History of continued frctions nd pde pproximnts Springer series in computtionl mthemtics, Vo.. Springer Verlg Berlin, 99 [] D.W. Lozier (997 A proposed softwre test service for specil functions in the boo of R. Boisvert (ed.. The qulity of Numericl Softwre : Assessment nd Enhncement Chpmn nd Hll. Londn. [3] D.W. Lozier nd F.W.J. Olver (Dec. 000 Numericl evlution of specil function in the boo of Wlter Gurtchi (ed. Mthemtics of computtion 93-993: A Hlf century of computtionl mthemtics, proceedings of symposi in Applied mthemtics, Americn mthemticl society province Rhode Islnd 090. []. J. Segur P Fermndez de Cord b nd Yu. C. rtis (997 : A code evlute modified Bessel functions bsed on the continued frction method Compute. Phy. Comm. 05 No. -3, 63-7. [5]. L. Monrel, P. Frenndez de cord J. A. Monsoriu nd S. Abrhm Ibrhim (003 specil function evlution using continued frctions methods nd Engineering AIP Conference Proceedings, 66 (I PP. 67-86 [6]. P. Henrici (977 Specil function, Intregl trnsform symptotics nd continuted frctions Applied nd computtionl Complex Anlysis-, Wiley New Yor. [7]. R.Y. Denis, S.N. Singh (000 Representtion of rtios of contiguous hypergeometric functions in the form of continued frctions, Fr. Est J. Mth Sci. (FJMS (000 03, 385-00. [8]. Yu. L. Rtis nd P. Fremndez de cordobe, A code to clculte (high order Bessel functions bsed on the continued frctions method, Compute Phys. Comm. 76 (993, 38-388. ****