Preprint
Article

This version is not peer-reviewed.

Higher-Order Partial Differentiation Formulas for Extended Srivastava Hypergeometric Functions via a General Kernel-Based Beta Function

Submitted:

23 May 2026

Posted:

25 May 2026

You are already at the latest version

Abstract
In this paper, we investigate the higher-order partial differentiation formulas for the new extended Srivastava hypergeometric functions KA, KB, and KC with respect to the parameters κ, ω, and τ. By employing their triple series representations together with classical identities of the Pochhammer symbol and the beta function, we derive explicit formulas for the r-th order partial derivatives. The obtained results show that repeated differentiation preserves the structural form of these functions, up to multiplicative Pochhammer factors and systematic shifts in the parameters. This reveals an inherent invariance property of the considered extended hypergeometric family under differential operators. The presented formulas provide a unified and consistent framework for all three classes of K-functions and may be useful in the study of fractional differential equations, integral transforms, and related applications in mathematical physics and applied analysis.
Keywords: 
;  ;  ;  ;  

1. Introduction

The theory of special functions continues to play a fundamental role in modern mathematical analysis, particularly in the study of differential equations, integral transforms, and applied mathematical modeling. Classical hypergeometric functions, introduced in the works of Gauss and generalized by Euler, Appell, and Horn, have served as a cornerstone in this development due to their rich analytic structure and remarkable stability properties under various operators, (see, [1,2,3,4,5,6,7,8]).
In recent decades, there has been a growing interest in constructing extended and generalized classes of hypergeometric-type functions in order to describe increasingly complex systems arising in physics, engineering, and applied sciences. These extensions are often motivated by the need to incorporate additional parameters, coupling effects, and multivariable interactions that cannot be adequately captured by classical frameworks. In this direction, Srivastava-type multivariable hypergeometric functions have emerged as a powerful and flexible tool, unifying several known families of special functions within a single analytical structure, (see, [9,10,11,12,13]).
More recently, further extensions involving generalized beta-type kernels and parameter-dependent deformations have been introduced. These constructions not only enrich the algebraic structure of hypergeometric functions but also provide new avenues for analytical manipulation, particularly in the context of operational calculus and fractional analysis. Such generalized functions are typically defined through multi-index series expansions involving Pochhammer symbols, gamma functions, and extended beta-type operators, which allow for a highly structured parameter dependence, (see, [14,15,16,17,18,19,20,21,22,23,24,25]).
Motivated by these developments, we consider in this work the extended Srivastava-type hypergeometric functions H ^ A , H ^ B , and H ^ C defined by a general kernel-based beta function. These functions are represented by triple hypergeometric series involving coupled parameters and a generalized beta-function structure generated through an arbitrary kernel function. This kernel-based formulation provides additional flexibility and unifies various previously known extensions within a common analytical framework. In particular, the parameters κ , ω , and τ play a central role in describing the interaction and deformation structure of the associated series representations.
A central motivation of this paper is to understand how these extended functions behave under repeated differentiation with respect to their deformation parameters. Unlike classical hypergeometric functions, where differentiation typically leads to simple parameter shifts, the presence of extended kernels and coupled Pochhammer structures leads to more intricate transformation rules. Nevertheless, as we demonstrate in this work, a remarkable invariance phenomenon still persists: the family of H ^ -functions remains closed under higher-order partial derivatives, up to explicit multiplicative Pochhammer factors and systematic parameter shifts.
From a theoretical perspective, this invariance property reveals a hidden algebraic stability of the extended Srivastava framework under differential operators. It also suggests that these functions form a natural class of objects for operational methods, particularly in problems involving repeated differentiation, convolution structures, and fractional-type generalizations. From an applied viewpoint, such properties are expected to be useful in the construction of exact solutions to generalized differential and integro-differential equations, especially those arising in mathematical physics and engineering models with memory or multi-scale effects.
The main contributions of this paper are as follows. First, we establish explicit r-th order partial differentiation formulas for H ^ A , H ^ B , and H ^ C with respect to κ , ω , and τ . Second, we show that these operations preserve the structural form of the underlying functions, leading to closed-form transformation rules. Third, we provide a unified treatment that highlights the symmetry and coherence among the three function classes within the extended Srivastava family.
The remainder of this paper is structured as follows. In Section 2, we present the preliminary definitions and classical results related to the gamma and beta functions, the Pochhammer symbol, and Srivastava’s triple hypergeometric functions. Section 3 introduces new extended Srivastava-type hypergeometric functions constructed via a general kernel-based beta function. In Section 4, we derive explicit higher-order partial differentiation formulas for the introduced function classes. Section 5 is devoted to the discussion and interpretation of the obtained analytical results, while concluding remarks and possible directions for future research are provided in Section 6.

2. Preliminaries

This section provides the fundamental mathematical concepts and auxiliary results required throughout the paper. We first recall several classical special functions, including the gamma and beta functions together with the Pochhammer symbol, which frequently arise in the theory of hypergeometric functions. Subsequently, we present the classical Srivastava triple hypergeometric functions H A , H B , and H C , which serve as the basis for the generalized formulations investigated in this work.
Definition 1 
([26]). The gamma function for ( σ ) > 0 is defined by
Γ ( σ ) = 0 ω σ 1 exp ω d ω .
The gamma function satisfies the following properties:
  • Γ ( σ + 1 ) = σ Γ ( σ ) ,
  • Γ ( σ + 1 ) = σ ! ,
  • Γ ( 1 ) = 1 .
Definition 2 
([26]). The beta function for ( σ ) > 0 and ( τ ) > 0 is defined by
B ( σ , τ ) = 0 1 ω σ 1 ( 1 ω ) τ 1 d ω .
The beta function satisfies the following properties:
  • B ( σ , τ ) = B ( τ , σ ) ,
  • B ( σ , τ ) = Γ ( σ ) Γ ( τ ) Γ ( σ + τ ) .
Definition 3 
([26]). The Pochhammer symbol ( σ ) n , where σ C and n N 0 , is defined by the relation
( σ ) n = Γ ( σ + n ) Γ ( σ ) = σ ( σ + 1 ) ( σ + n 1 ) , n N , 1 , n = 0 .
The Pochhammer symbol satisfies the following properties:
  • ( σ ) m + n = ( σ ) m ( σ + m ) n ,
  • ( σ ) n + 1 = σ ( σ + 1 ) n ,
  • ( 1 ) n = n ! .
Definition 4 
([27,28]). The classical Srivastava’s triple hypergeometric functions H A , H B and H C , respectively, are defined by
H A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n + k κ m m ! ω n n ! τ k k ! ,
s < 1 , t < 1 , r < ( 1 s ) ( 1 t ) ,
H B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k κ m m ! ω n n ! τ k k ! ,
r + s + t + 2 r s t < 1 ,
and
H C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) m + n + k κ m m ! ω n n ! τ k k ! ,
r < 1 , s < 1 , t < 1 , r + s + t 2 ( 1 r ) ( 1 s ) ( 1 t ) < 2 ,
where, for simplicity of notation, the coordinates ( r , s , t ) are used in place of ( | κ | , | ω | , | τ | ) .

3. Extended Srivastava Hypergeometric Functions

In this section, we introduce new extensions of the classical Srivastava triple hypergeometric functions via the general kernel-based beta function. This construction provides a unified and flexible framework that includes several known extensions as special cases corresponding to appropriate choices of the kernel function. The proposed approach preserves the analytic structure of the classical Srivastava functions while enriching them through additional deformation parameters embedded in the kernel.
We first recall the general kernel-based beta function.
Definition 5 
([25]). The general kernel-based beta function is defined by
B ^ ( σ , τ ) = 0 1 ω σ 1 ( 1 ω ) τ 1 K ( ω , X ) d ω ,
where X = X ( p , q , α , β ) denotes a multi-parameter variable and K ( ω , X ) denotes a general kernel function. The parameters satisfy
( p ) > 0 , ( q ) > 0 , ( α ) > 0 , ( β ) > 0 ,
and
( σ ) > 0 , ( τ ) > 0 .
The general kernel function K ( ω , X ) admits several important structures, allowing the general kernel-based beta function definition to unify a broad family of extended or generalized beta functions. In this context, the K ( ω , X ) can take the following forms:
(i) Product-type kernels
  • K ( ω , X ) = K p ω α K q ( 1 ω ) β ,
  • K ( ω , X ) = K p ω α K q ( 1 ω ) α ,
  • K ( ω , X ) = K p ω K q 1 ω .
(ii) Additive-type kernels
  • K ( ω , X ) = K p ω α q ( 1 ω ) β ,
  • K ( ω , X ) = K p ω α q ( 1 ω ) α ,
  • K ( ω , X ) = K p ω q 1 ω .
(iii) Coupled-type kernels
  • K ( ω , X ) = K p ω α ( 1 ω ) β ,
  • K ( ω , X ) = K p ω α ( 1 ω ) α ,
  • K ( ω , X ) = K p ω ( 1 ω ) .
The general kernel function K may be selected from a broad class of special functions such as the exponential function, Kummer confluent hypergeometric function, Mittag–Leffler function, Wright function, Fox–Wright function, and the M-series. Consequently, the proposed framework provides a unified representation for many extended or generalized beta-type functions available in the literature.
Definition 6. 
Assume that the coordinates are given by ( r , s , t ) instead of ( | κ | , | ω | , | τ | ) and let the parameters satisfy ( p ) > 0 , ( q ) > 0 , ( α ) > 0 , ( β ) > 0 . Then, using the general kernel-based beta function, the extended forms of the Srivastava hypergeometric functions H A , H B , and H C are defined as follows:
H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = H A , p , q ( α , β ) ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! ω n n ! τ k k ! ,
s < 1 , t < 1 , r < ( 1 s ) ( 1 t ) ,
H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = H B , p , q ( α , β ) ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 + μ 2 ) 2 m + n + k ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + m + n ) B ( μ 1 , μ 2 ) κ m m ! ω n n ! τ k k ! ,
r + s + t + 2 r s t < 1 ,
and
H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = H C , p , q ( α , β ) ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m , n , k = 0 ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) n B ^ ( μ 1 + m + k , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) κ m m ! ω n n ! τ k k ! ,
r < 1 , s < 1 , t < 1 , r + s + t 2 ( 1 r ) ( 1 s ) ( 1 t ) < 2 .
They are called the Srivastava hypergeometric function H ^ A , the Srivastava hypergeometric function H ^ B , and the Srivastava hypergeometric function H ^ C , respectively. These functions, owing to the incorporation of the general kernel-based beta function, provide a unified framework that encompasses both the classical Srivastava functions and various extensions of Srivastava-type hypergeometric functions reported in the literature.

4. Higher-Order Partial DIfferentiation Formulas

In this section, we derive a family of r-th order differentiation formulas for the newly introduced extended Srivastava hypergeometric functions H ^ A , H ^ B , and H ^ C . The main objective is to investigate how these functions behave under repeated partial differentiation with respect to the parameters κ , ω , and τ . By working directly from the defining triple series representations, we apply term-by-term differentiation and use classical identities for factorials, the gamma function, and Pochhammer symbols. This procedure reveals that each differentiation preserves the structural form of the original function, while only shifting the involved parameters. As a result, all obtained formulas exhibit a unified pattern: the r-th derivative of each function can be expressed as a scaled version of the same function with appropriately shifted parameters. This demonstrates that the introduced function families are closed under differentiation and possess a stable algebraic structure, which is essential for their applications in fractional calculus and related differential equation models.
Theorem 1. 
For r N 0 , the following partial differentiation formula with respect to κ holds:
r κ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = ( μ 1 ) r ( μ 2 ) r ( μ 4 ) r H ^ A ( μ 1 + r , μ 2 + r , μ 3 ; μ 4 + r , μ 5 ; κ , ω , τ ) .
Proof. 
Starting from the defining series representation
H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! ω n n ! τ k k ! ,
we differentiate term-by-term with respect to κ . Thus,
r κ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) ω n n ! τ k k ! r κ r κ m m ! .
Using
r κ r κ m m ! = κ m r ( m r ) ! , m r ,
we obtain
r κ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m = r n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m r ( m r ) ! ω n n ! τ k k ! .
Now, by shifting the index m m + r , we get
r κ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + r + k ( μ 2 ) m + r + n ( μ 4 ) m + r B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! ω n n ! τ k k ! .
Next, applying the Pochhammer decomposition formula
( a ) m + r = ( a ) r ( a + r ) m ,
we obtain
( μ 1 ) m + r + k = ( μ 1 ) r ( μ 1 + r ) m + k ,
( μ 2 ) m + r + n = ( μ 2 ) r ( μ 2 + r ) m + n ,
and
( μ 4 ) m + r = ( μ 4 ) r ( μ 4 + r ) m .
Substituting these relations into the above expression yields
r κ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = ( μ 1 ) r ( μ 2 ) r ( μ 4 ) r m , n , k = 0 ( μ 1 + r ) m + k ( μ 2 + r ) m + n ( μ 4 + r ) m B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! ω n n ! τ k k ! .
Recognizing the resulting series as
H ^ A ( μ 1 + r , μ 2 + r , μ 3 ; μ 4 + r , μ 5 ; κ , ω , τ ) ,
we finally obtain
r κ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = ( μ 1 ) r ( μ 2 ) r ( μ 4 ) r H ^ A ( μ 1 + r , μ 2 + r , μ 3 ; μ 4 + r , μ 5 ; κ , ω , τ ) .
This completes the proof. □
Theorem 2. 
For r N 0 , the following partial differentiation formula with respect to ω holds:
r ω r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = ( μ 2 ) r ( μ 3 ) r ( μ 5 ) r H ^ A ( μ 1 , μ 2 + r , μ 3 + r ; μ 4 , μ 5 + r ; κ , ω , τ ) .
Proof. 
Starting from the defining series representation
H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! ω n n ! τ k k ! .
Differentiating term-by-term with respect to ω , we obtain
r ω r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! τ k k ! r ω r ω n n ! .
Using
r ω r ω n n ! = ω n r ( n r ) ! , n r ,
we get
r ω r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = n = r m , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! ω n r ( n r ) ! τ k k ! .
Now, by shifting the index n n + r , we obtain
r ω r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n + r ( μ 4 ) m B ^ ( μ 3 + n + k + r , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! ω n n ! τ k k ! .
Applying the Pochhammer decomposition formula
( a ) m + r = ( a ) r ( a + r ) m ,
we have
( μ 2 ) m + n + r = ( μ 2 ) r ( μ 2 + r ) m + n .
Furthermore, using the beta-function identity
B ( μ 3 + r , μ 5 μ 3 ) = ( μ 3 ) r ( μ 5 ) r B ( μ 3 , μ 5 μ 3 ) ,
we obtain
1 B ( μ 3 , μ 5 μ 3 ) = ( μ 3 ) r ( μ 5 ) r 1 B ( μ 3 + r , μ 5 μ 3 ) .
Hence,
r ω r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = ( μ 2 ) r ( μ 3 ) r ( μ 5 ) r m , n , k = 0 ( μ 1 ) m + k ( μ 2 + r ) m + n ( μ 4 ) m B ^ ( μ 3 + r + n + k , μ 5 μ 3 ) B ( μ 3 + r , μ 5 μ 3 ) κ m m ! ω n n ! τ k k ! .
Recognizing the resulting series as
H ^ A ( μ 1 , μ 2 + r , μ 3 + r ; μ 4 , μ 5 + r ; κ , ω , τ ) ,
we finally arrive at
r ω r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = ( μ 2 ) r ( μ 3 ) r ( μ 5 ) r H ^ A ( μ 1 , μ 2 + r , μ 3 + r ; μ 4 , μ 5 + r ; κ , ω , τ ) .
This completes the proof. □
Theorem 3. 
For r N 0 , the following partial differentiation formula with respect to τ holds:
r τ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = ( μ 1 ) r ( μ 3 ) r ( μ 5 ) r H ^ A ( μ 1 + r , μ 2 , μ 3 + r ; μ 4 , μ 5 + r ; κ , ω , τ ) .
Proof. 
Starting from the defining series representation
H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! ω n n ! τ k k ! .
Differentiating term-by-term with respect to τ , we obtain
r τ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! ω n n ! r τ r τ k k ! .
Using
r τ r τ k k ! = τ k r ( k r ) ! , k r ,
we obtain
r τ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = k = r m , n = 0 ( μ 1 ) m + k ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + n + k , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! ω n n ! τ k r ( k r ) ! .
Now, by shifting the index k k + r , we get
r τ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 ) m + k + r ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + n + k + r , μ 5 μ 3 ) B ( μ 3 , μ 5 μ 3 ) κ m m ! ω n n ! τ k k ! .
Applying the Pochhammer decomposition formula
( a ) m + r = ( a ) r ( a + r ) m ,
we obtain
( μ 1 ) m + k + r = ( μ 1 ) r ( μ 1 + r ) m + k .
Furthermore, using the beta-function identity
B ( μ 3 + r , μ 5 μ 3 ) = ( μ 3 ) r ( μ 5 ) r B ( μ 3 , μ 5 μ 3 ) ,
we have
1 B ( μ 3 , μ 5 μ 3 ) = ( μ 3 ) r ( μ 5 ) r 1 B ( μ 3 + r , μ 5 μ 3 ) .
Substituting these relations into the previous expression yields
r τ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = ( μ 1 ) r ( μ 3 ) r ( μ 5 ) r m , n , k = 0 ( μ 1 + r ) m + k ( μ 2 ) m + n ( μ 4 ) m B ^ ( μ 3 + r + n + k , μ 5 μ 3 ) B ( μ 3 + r , μ 5 μ 3 ) κ m m ! ω n n ! τ k k ! .
Recognizing the resulting series as
H ^ A ( μ 1 + r , μ 2 , μ 3 + r ; μ 4 , μ 5 + r ; κ , ω , τ ) ,
we finally obtain
r τ r H ^ A ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 ; κ , ω , τ ) = ( μ 1 ) r ( μ 3 ) r ( μ 5 ) r H ^ A ( μ 1 + r , μ 2 , μ 3 + r ; μ 4 , μ 5 + r ; κ , ω , τ ) .
This completes the proof. □
Theorem 4. 
For r N 0 , the following partial differentiation formula with respect to κ holds:
r κ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = ( μ 1 ) r ( μ 2 ) r ( μ 4 ) r H ^ B ( μ 1 + r , μ 2 + r , μ 3 ; μ 4 + r , μ 5 , μ 6 ; κ , ω , τ ) .
Proof. 
Starting from the defining series representation
H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 + μ 2 ) 2 m + n + k ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + m + n ) B ( μ 1 , μ 2 ) κ m m ! ω n n ! τ k k ! ,
we differentiate term-by-term with respect to κ . Thus,
r κ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 + μ 2 ) 2 m + n + k ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + m + n ) B ( μ 1 , μ 2 ) ω n n ! τ k k ! r κ r κ m m ! .
Using
r κ r κ m m ! = κ m r ( m r ) ! , m r ,
we obtain
r κ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m = r n , k = 0 ( μ 1 + μ 2 ) 2 m + n + k ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + m + n ) B ( μ 1 , μ 2 ) κ m r ( m r ) ! ω n n ! τ k k ! .
Now, by shifting the index m m + r , we get
r κ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 + μ 2 ) 2 m + 2 r + n + k ( μ 3 ) n + k ( μ 4 ) m + r ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + r + m + k , μ 2 + r + m + n ) B ( μ 1 , μ 2 ) κ m m ! ω n n ! τ k k ! .
Next, using the Pochhammer decomposition formula
( a ) m + r = ( a ) r ( a + r ) m ,
we obtain
( μ 4 ) m + r = ( μ 4 ) r ( μ 4 + r ) m .
Furthermore, employing the beta-function identity
B ( a + r , b + r ) = ( a ) r ( b ) r ( a + b ) 2 r B ( a , b ) ,
we deduce that
1 B ( μ 1 , μ 2 ) = ( μ 1 ) r ( μ 2 ) r ( μ 1 + μ 2 ) 2 r 1 B ( μ 1 + r , μ 2 + r ) .
Substituting these relations into the previous expression yields
r κ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = ( μ 1 ) r ( μ 2 ) r ( μ 4 ) r m , n , k = 0 ( μ 1 + μ 2 + 2 r ) 2 m + n + k ( μ 3 ) n + k ( μ 4 + r ) m ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + r + m + k , μ 2 + r + m + n ) B ( μ 1 + r , μ 2 + r ) κ m m ! ω n n ! τ k k ! .
Finally, recognizing the resulting series as
H ^ B ( μ 1 + r , μ 2 + r , μ 3 ; μ 4 + r , μ 5 , μ 6 ; κ , ω , τ ) ,
we arrive at
r κ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = ( μ 1 ) r ( μ 2 ) r ( μ 4 ) r H ^ B ( μ 1 + r , μ 2 + r , μ 3 ; μ 4 + r , μ 5 , μ 6 ; κ , ω , τ ) .
This completes the proof. □
Theorem 5. 
For r N 0 , the following partial differentiation formula with respect to ω holds:
r ω r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = ( μ 2 ) r ( μ 3 ) r ( μ 5 ) r H ^ B ( μ 1 , μ 2 + r , μ 3 + r ; μ 4 , μ 5 + r , μ 6 ; κ , ω , τ ) .
Proof. 
Starting from the defining series representation
H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 + μ 2 ) 2 m + n + k ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + m + n ) B ( μ 1 , μ 2 ) κ m m ! ω n n ! τ k k ! ,
we differentiate term-by-term with respect to ω . Thus,
r ω r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 + μ 2 ) 2 m + n + k ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + m + n ) B ( μ 1 , μ 2 ) κ m m ! τ k k ! r ω r ω n n ! .
Using
r ω r ω n n ! = ω n r ( n r ) ! , n r ,
we obtain
r ω r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = n = r m , k = 0 ( μ 1 + μ 2 ) 2 m + n + k ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + m + n ) B ( μ 1 , μ 2 ) κ m m ! ω n r ( n r ) ! τ k k ! .
Now, by shifting the index n n + r , we get
r ω r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 + μ 2 ) 2 m + n + r + k ( μ 3 ) n + r + k ( μ 4 ) m ( μ 5 ) n + r ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + m + n + r ) B ( μ 1 , μ 2 ) κ m m ! ω n n ! τ k k ! .
Next, applying the Pochhammer decomposition formula
( a ) m + r = ( a ) r ( a + r ) m ,
we obtain
( μ 3 ) n + r + k = ( μ 3 ) r ( μ 3 + r ) n + k ,
and
( μ 5 ) n + r = ( μ 5 ) r ( μ 5 + r ) n .
Furthermore, employing the beta-function identity
B ( a , b + r ) = ( b ) r ( a + b ) r B ( a , b ) ,
we deduce that
1 B ( μ 1 , μ 2 ) = ( μ 2 ) r ( μ 1 + μ 2 ) r 1 B ( μ 1 , μ 2 + r ) ,
and hence the factor ( μ 1 + μ 2 ) r cancels. Therefore,
r ω r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = ( μ 2 ) r ( μ 3 ) r ( μ 5 ) r m , n , k = 0 ( μ 1 + μ 2 + r ) 2 m + n + k ( μ 3 + r ) n + k ( μ 4 ) m ( μ 5 + r ) n ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + r + m + n ) B ( μ 1 , μ 2 + r ) κ m m ! ω n n ! τ k k ! .
Recognizing the resulting series as
H ^ B ( μ 1 , μ 2 + r , μ 3 + r ; μ 4 , μ 5 + r , μ 6 ; κ , ω , τ ) ,
we arrive at
r ω r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = ( μ 2 ) r ( μ 3 ) r ( μ 5 ) r H ^ B ( μ 1 , μ 2 + r , μ 3 + r ; μ 4 , μ 5 + r , μ 6 ; κ , ω , τ ) .
This completes the proof. □
Theorem 6. 
For r N 0 , the following partial differentiation formula with respect to τ holds:
r τ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = ( μ 1 ) r ( μ 3 ) r ( μ 6 ) r H ^ B ( μ 1 + r , μ 2 , μ 3 + r ; μ 4 , μ 5 , μ 6 + r ; κ , ω , τ ) .
Proof. 
Starting from the defining series representation
H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 + μ 2 ) 2 m + n + k ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + m + n ) B ( μ 1 , μ 2 ) κ m m ! ω n n ! τ k k ! ,
we differentiate term-by-term with respect to τ . Thus,
r τ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 + μ 2 ) 2 m + n + k ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + m + n ) B ( μ 1 , μ 2 ) κ m m ! ω n n ! r τ r τ k k ! .
Using
r τ r τ k k ! = τ k r ( k r ) ! , k r ,
we obtain
r τ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = k = r m , n = 0 ( μ 1 + μ 2 ) 2 m + n + k ( μ 3 ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k B ^ ( μ 1 + m + k , μ 2 + m + n ) B ( μ 1 , μ 2 ) κ m m ! ω n n ! τ k r ( k r ) ! .
Now, by shifting the index k k + r , we get
r τ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = m , n , k = 0 ( μ 1 + μ 2 ) 2 m + n + k + r ( μ 3 ) n + k + r ( μ 4 ) m ( μ 5 ) n ( μ 6 ) k + r B ^ ( μ 1 + m + k + r , μ 2 + m + n ) B ( μ 1 , μ 2 ) κ m m ! ω n n ! τ k k ! .
Next, applying the Pochhammer decomposition formula
( a ) m + r = ( a ) r ( a + r ) m ,
we obtain
( μ 3 ) n + k + r = ( μ 3 ) r ( μ 3 + r ) n + k ,
and
( μ 6 ) k + r = ( μ 6 ) r ( μ 6 + r ) k .
Furthermore, using
B ^ ( μ 1 + m + k + r , μ 2 + m + n ) B ( μ 1 , μ 2 ) = ( μ 1 ) r ( μ 1 + μ 2 ) r B ^ ( μ 1 + r + m + k , μ 2 + m + n ) B ( μ 1 + r , μ 2 ) ,
the factor ( μ 1 + μ 2 ) r is cancelled by
( μ 1 + μ 2 ) 2 m + n + k + r = ( μ 1 + μ 2 ) r ( μ 1 + μ 2 + r ) 2 m + n + k .
Hence,
r τ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = ( μ 1 ) r ( μ 3 ) r ( μ 6 ) r m , n , k = 0 ( μ 1 + μ 2 + r ) 2 m + n + k ( μ 3 + r ) n + k ( μ 4 ) m ( μ 5 ) n ( μ 6 + r ) k B ^ ( μ 1 + r + m + k , μ 2 + m + n ) B ( μ 1 + r , μ 2 ) κ m m ! ω n n ! τ k k ! .
Recognizing the resulting series as
H ^ B ( μ 1 + r , μ 2 , μ 3 + r ; μ 4 , μ 5 , μ 6 + r ; κ , ω , τ ) ,
we arrive at
r τ r H ^ B ( μ 1 , μ 2 , μ 3 ; μ 4 , μ 5 , μ 6 ; κ , ω , τ ) = ( μ 1 ) r ( μ 3 ) r ( μ 6 ) r H ^ B ( μ 1 + r , μ 2 , μ 3 + r ; μ 4 , μ 5 , μ 6 + r ; κ , ω , τ ) .
This completes the proof. □
Theorem 7. 
For r N 0 , the following partial differentiation formula with respect to κ holds:
r κ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = ( μ 1 ) r ( μ 2 ) r ( μ 4 ) r H ^ C ( μ 1 + r , μ 2 + r , μ 3 ; μ 4 + r ; κ , ω , τ ) .
Proof. 
Starting from the defining series representation
H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m , n , k = 0 ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) n B ^ ( μ 1 + m + k , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) κ m m ! ω n n ! τ k k ! .
Differentiating term-by-term with respect to κ , we obtain
r κ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m , n , k = 0 ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) n B ^ ( μ 1 + m + k , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) ω n n ! τ k k ! r κ r κ m m ! .
Using
r κ r κ m m ! = κ m r ( m r ) ! , m r ,
we get
r κ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m = r n , k = 0 ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) n B ^ ( μ 1 + m + k , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) κ m r ( m r ) ! ω n n ! τ k k ! .
Now, by shifting the index m m + r , we obtain
r κ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m , n , k = 0 ( μ 2 ) m + n + r ( μ 3 ) n + k ( μ 4 ) n B ^ ( μ 1 + m + k + r , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) κ m m ! ω n n ! τ k k ! .
Applying the Pochhammer decomposition formula
( a ) m + r = ( a ) r ( a + r ) m ,
we have
( μ 2 ) m + n + r = ( μ 2 ) r ( μ 2 + r ) m + n .
Furthermore, using the beta-function identity
B ( a + r , b ) = ( a ) r ( a + b ) r B ( a , b ) ,
we obtain
B ( μ 1 + r , μ 4 + n μ 1 ) = ( μ 1 ) r ( μ 4 + n ) r B ( μ 1 , μ 4 + n μ 1 ) ,
which yields
1 B ( μ 1 , μ 4 + n μ 1 ) = ( μ 1 ) r ( μ 4 + n ) r 1 B ( μ 1 + r , μ 4 + n μ 1 ) .
Since
( μ 4 + n ) r = ( μ 4 ) n + r ( μ 4 ) n ,
we have
1 ( μ 4 ) n 1 ( μ 4 + n ) r = 1 ( μ 4 ) n + r .
Hence,
r κ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = ( μ 1 ) r ( μ 2 ) r ( μ 4 ) r m , n , k = 0 ( μ 2 + r ) m + n ( μ 3 ) n + k ( μ 4 + r ) n B ^ ( μ 1 + r + m + k , μ 4 + n μ 1 ) B ( μ 1 + r , μ 4 + n μ 1 ) κ m m ! ω n n ! τ k k ! .
Recognizing the resulting series as
H ^ C ( μ 1 + r , μ 2 + r , μ 3 ; μ 4 + r ; κ , ω , τ ) ,
we finally obtain
r κ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = ( μ 1 ) r ( μ 2 ) r ( μ 4 ) r H ^ C ( μ 1 + r , μ 2 + r , μ 3 ; μ 4 + r ; κ , ω , τ ) .
This completes the proof. □
Theorem 8. 
For r N 0 , the following partial differentiation formula with respect to ω holds:
r ω r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = ( μ 2 ) r ( μ 3 ) r ( μ 4 ) r H ^ C ( μ 1 , μ 2 + r , μ 3 + r ; μ 4 + r ; κ , ω , τ ) .
Proof. 
Starting from the defining series representation
H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m , n , k = 0 ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) n B ^ ( μ 1 + m + k , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) κ m m ! ω n n ! τ k k ! .
Differentiating term-by-term with respect to ω , we obtain
r ω r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m , n , k = 0 ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) n B ^ ( μ 1 + m + k , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) κ m m ! τ k k ! r ω r ω n n ! .
Using
r ω r ω n n ! = ω n r ( n r ) ! , n r ,
we get
r ω r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = n = r m , k = 0 ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) n B ^ ( μ 1 + m + k , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) κ m m ! ω n r ( n r ) ! τ k k ! .
Now, by shifting the index n n + r , we obtain
r ω r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m , n , k = 0 ( μ 2 ) m + n + r ( μ 3 ) n + k + r ( μ 4 ) n + r B ^ ( μ 1 + m + k , μ 4 + n + r μ 1 ) B ( μ 1 , μ 4 + n + r μ 1 ) κ m m ! ω n n ! τ k k ! .
Applying the Pochhammer decomposition formula
( a ) m + r = ( a ) r ( a + r ) m ,
we have
( μ 2 ) m + n + r = ( μ 2 ) r ( μ 2 + r ) m + n ,
and
( μ 3 ) n + k + r = ( μ 3 ) r ( μ 3 + r ) n + k .
Also,
( μ 4 ) n + r = ( μ 4 ) r ( μ 4 + r ) n .
Hence,
r ω r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = ( μ 2 ) r ( μ 3 ) r ( μ 4 ) r m , n , k = 0 ( μ 2 + r ) m + n ( μ 3 + r ) n + k ( μ 4 + r ) n B ^ ( μ 1 + m + k , μ 4 + r + n μ 1 ) B ( μ 1 , μ 4 + r + n μ 1 ) κ m m ! ω n n ! τ k k ! .
Recognizing the resulting series as
H ^ C ( μ 1 , μ 2 + r , μ 3 + r ; μ 4 + r ; κ , ω , τ ) ,
we finally obtain
r ω r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = ( μ 2 ) r ( μ 3 ) r ( μ 4 ) r H ^ C ( μ 1 , μ 2 + r , μ 3 + r ; μ 4 + r ; κ , ω , τ ) .
This completes the proof. □
Theorem 9. 
For r N 0 , the following partial differentiation formula with respect to τ holds:
r τ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = ( μ 1 ) r ( μ 3 ) r ( μ 4 ) r H ^ C ( μ 1 + r , μ 2 , μ 3 + r ; μ 4 + r ; κ , ω , τ ) .
Proof. 
Starting from the defining series representation
H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m , n , k = 0 ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) n B ^ ( μ 1 + m + k , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) κ m m ! ω n n ! τ k k ! .
Differentiating term-by-term with respect to τ , we obtain
r τ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m , n , k = 0 ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) n B ^ ( μ 1 + m + k , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) κ m m ! ω n n ! r τ r τ k k ! .
Using
r τ r τ k k ! = τ k r ( k r ) ! , k r ,
we obtain
r τ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = k = r m , n = 0 ( μ 2 ) m + n ( μ 3 ) n + k ( μ 4 ) n B ^ ( μ 1 + m + k , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) κ m m ! ω n n ! τ k r ( k r ) ! .
Now, by shifting the index k k + r , we obtain
r τ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = m , n , k = 0 ( μ 2 ) m + n ( μ 3 ) n + k + r ( μ 4 ) n B ^ ( μ 1 + m + k + r , μ 4 + n μ 1 ) B ( μ 1 , μ 4 + n μ 1 ) κ m m ! ω n n ! τ k k ! .
Applying the Pochhammer decomposition formula
( a ) m + r = ( a ) r ( a + r ) m ,
we have
( μ 3 ) n + k + r = ( μ 3 ) r ( μ 3 + r ) n + k .
Furthermore, using the beta-function identity
B ( a + r , b ) = ( a ) r ( a + b ) r B ( a , b ) ,
we obtain
B ( μ 1 + r , μ 4 + n μ 1 ) = ( μ 1 ) r ( μ 4 + n ) r B ( μ 1 , μ 4 + n μ 1 ) .
Hence,
1 B ( μ 1 , μ 4 + n μ 1 ) = ( μ 1 ) r ( μ 4 + n ) r 1 B ( μ 1 + r , μ 4 + n μ 1 ) .
Since
( μ 4 + n ) r = ( μ 4 ) n + r ( μ 4 ) n ,
we get
1 ( μ 4 + n ) r = ( μ 4 ) n ( μ 4 ) n + r .
Substituting these relations into the previous expression yields
r τ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = ( μ 1 ) r ( μ 3 ) r m , n , k = 0 ( μ 2 ) m + n ( μ 3 + r ) n + k ( μ 4 ) n + r B ^ ( μ 1 + r + m + k , μ 4 + n μ 1 ) B ( μ 1 + r , μ 4 + n μ 1 ) κ m m ! ω n n ! τ k k ! .
Also, we have
1 ( μ 4 ) n + r = 1 ( μ 4 ) r ( μ 4 + r ) n .
Then, we get
r τ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = ( μ 1 ) r ( μ 3 ) r ( μ 4 ) r m , n , k = 0 ( μ 2 ) m + n ( μ 3 + r ) n + k ( μ 4 + r ) n B ^ ( μ 1 + r + m + k , μ 4 + n μ 1 ) B ( μ 1 + r , μ 4 + n μ 1 ) κ m m ! ω n n ! τ k k ! .
Recognizing the resulting series as
H ^ C ( μ 1 + r , μ 2 , μ 3 + r ; μ 4 + r ; κ , ω , τ ) ,
we finally obtain
r τ r H ^ C ( μ 1 , μ 2 , μ 3 ; μ 4 ; κ , ω , τ ) = ( μ 1 ) r ( μ 3 ) r ( μ 4 ) r H ^ C ( μ 1 + r , μ 2 , μ 3 + r ; μ 4 + r ; κ , ω , τ ) .
This completes the proof. □

5. Results and Discussion

In this section, we summarize the main results obtained in the previous part and discuss their structural and analytical implications for the extended Srivastava-type hypergeometric functions H ^ A , H ^ B , and H ^ C .

Results

We have derived explicit r-th order partial differentiation formulas with respect to the variables κ , ω , and τ for each class of the introduced functions. In all cases, the following common structural properties are observed:
  • The r-th derivative with respect to a given variable preserves the original functional form.
  • Each differentiation introduces multiplicative coefficients expressed in terms of Pochhammer symbols.
  • The parameters of the function are systematically shifted by the order of differentiation r.
More precisely, for each family H ^ A , H ^ B , and H ^ C , we obtain identities of the general form
r x r H ^ ( · ) = ( Pochhammer factors ) × H ^ ( shifted parameters ) ,
where x { κ , ω , τ } .

Discussion

The obtained results show that all three function classes possess a strong closure property under partial differentiation. In particular, repeated differentiation does not alter the functional structure but only modifies the parameter set in a controlled way. This indicates a stable algebraic framework similar to classical hypergeometric functions.
From an analytical point of view, this property is significant for several reasons:
  • It simplifies the treatment of higher-order differential equations involving these functions.
  • It enables direct construction of recursive relations between derivatives.
  • It provides a natural framework for applications in fractional calculus, where such shift-operator behavior is frequently encountered.
Moreover, the obtained parameter-shift structures appearing in the differentiation formulas indicate that the function families H ^ A , H ^ B , and H ^ C possess a common analytical framework. This demonstrates that the proposed kernel-based construction provides a unified framework for generalized Srivastava-type hypergeometric functions.
Overall, the results demonstrate that these functions are not only well-defined analytically but also highly structured under differentiation, making them suitable candidates for further study in integral transforms, fractional differential equations, and applied mathematical modeling.

6. Conclusions

In this work, we have established a systematic investigation of the differential properties of the extended Srivastava-type hypergeometric functions H ^ A , H ^ B , and H ^ C with respect to the variables κ , ω , and τ . By employing their series representations together with classical identities of the Pochhammer symbol and the beta function, we derived explicit r-th order differentiation formulas.
The obtained results show that differentiation with respect to each variable preserves the structural form of the underlying functions, while inducing systematic parameter shifts and multiplicative Pochhammer factors. This invariance-type property highlights the algebraic stability of the considered extended hypergeometric family under repeated differentiation.
Moreover, the presented formulas provide a unified framework that connects the different members of the H ^ -family, revealing consistent patterns in parameter evolution across H ^ A , H ^ B , and H ^ C . These results may be useful in further analytical studies, particularly in solving fractional differential equations, constructing integral transforms, and investigating applied models in mathematical physics and engineering.
Future research may focus on extending these results to fractional-order derivatives, integral representations, and generating function approaches, as well as exploring potential applications in special function theory and fractional calculus.

Author Contributions

E.A.-Conceptualization, Methodology, Data Curation, Writing Original Draft, Visualization, Investigation, Supervision, Formal Analysis, Validation. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All data supporting the findings of this study are included in the article.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Goswami, Y.; Sharma, V. Adjoint Appell-Euler polynomials and approximation. Filomat 2026, 40(7), 2375–2392. [Google Scholar]
  2. Usman, T.; Khan, N.; Ghayasuddin, M. An analysis of generalized Appell’s and Lauricella’s hypergeometric functions and their attributes. Moroc. J. Pure Appl. Anal. 2025, 11(3), 266–277. [Google Scholar]
  3. Ahmad, N.; Khan, W.A. A new generalization of q-Laguerre-based Appell polynomials and quasi-monomiality. Symmetry 2025, 17(3), 439. [Google Scholar] [CrossRef]
  4. Husain, S. A study of generalized Appell’s functions and RL fractional derivative operator. Honam Math. J. 2025, 47(3), 367–387. [Google Scholar]
  5. Kim, I.; Rathie, A.K. A note on certain general transformation formulas for the Appell and the Horn functions. Symmetry 2023, 15(3), 696. [Google Scholar] [CrossRef]
  6. Ata, E. M-Lauricella hypergeometric functions: integral representations and solutions of fractional differential equations. Commun. Fac. Sci. Univ. Ank. Ser. A1 Math. Stat. 2023, 72(2), 512–529. [Google Scholar] [CrossRef]
  7. Dmytryshyn, R.I.; Lutsiv, I.A. Three- and four-term recurrence relations for Horn’s hypergeometric function H4. Res. Math. 2022, 30(1), 21–29. [Google Scholar] [CrossRef]
  8. Issenova, A.A.; Tasmambetov, Z.N.; Talipova, M.Z. Construction of solutions hypergeometric system of Horn type in the form of Laguerre polynomials. Lobachevskii J. Math. 2022, 43(11), 3167–3173. [Google Scholar] [CrossRef]
  9. Ata, E. M-Srivastava hypergeometric functions: integral representations and solution of fractional differential equations. Miskolc Math. Notes 2024, 25(1), 93–108. [Google Scholar] [CrossRef]
  10. Srivastava, H.M.; Chandel, R.S.; Kumar, H. Some general Hurwitz-Lerch type zeta functions associated with the Srivastava-Daoust multiple hypergeometric functions. J. Nonlinear Var. Anal. 2022, 6(4), 299–315. [Google Scholar]
  11. Srivastava, H.M.; Gupta, B.; Qureshi, M.I.; Baboo, M.S. Some summation theorems and transformations for hypergeometric functions of Kampé de Fériet and Srivastava. Georgian Math. J. 2024, 31(5), 885–897. [Google Scholar] [CrossRef]
  12. Dar, S.A.; Kamarujjama, M.; Daud, M. Some properties of generalized Srivastava’s triple hypergeometric function HC,p,q(·). Int. J. Appl. Comput. Math. 2022, 8(4), 168. [Google Scholar] [CrossRef]
  13. Verma, A. Some results on Srivastava’s triple hypergeometric matrix functions. Asian-Eur. J. Math. 2021, 14(04), 2150056. [Google Scholar] [CrossRef]
  14. Ata, E. A general kernel-based beta function approach to the forward expansion of generalized M-series with applications to fractional operators and integral transforms. Boletín De La Soc. Matemática Mex. 2026, 32(1), 41. [Google Scholar] [CrossRef]
  15. Parida, S. On Mittag-Leffler-Kummer extended beta function. Honam Math. J. 2025, 47(1), 24–33. [Google Scholar]
  16. Ata, E.; Kiymaz, I.O. Modified special functions: properties, integral transforms and applications to fractional differential equations. Bol. Da Soc. Parana. De Matemática 2024, 42, 1–22. [Google Scholar] [CrossRef]
  17. Ali, S.; Regar, N.K.; Parida, S. On generalized extended beta and hypergeometric functions. Honam Math. J. 2024, 46(2), 313–334. [Google Scholar]
  18. Ata, E.; Kiymaz, I.O. Generalized gamma, beta and hypergeometric functions defined by Wright function and applications to fractional differential equations. Cumhur. Sci. J. 2022, 43(4), 684–695. [Google Scholar] [CrossRef]
  19. Khudhuir, L.T.; Al-Janaby, H.F.; Ghanim, F.; Lupas, A.A. A Wright-based generalization of the Euler beta function with statistical applications. Mathematics 2026, 14(6), 1069. [Google Scholar] [CrossRef]
  20. Ata, E.; Kiymaz, I.O. New generalized special functions with two generalized M-series at their kernels and solution of fractional PDEs via double Laplace transform. Comput. Methods Differ. Equ. 2024, 12(1), 31–43. [Google Scholar]
  21. Mubeen, S.; Aslam, I.; Khammash, G.S.; Nadarajah, S.; Shehata, A. New modified gamma and beta functions with an application. Contemp. Math. 2025, 6(1), 850–861. [Google Scholar] [CrossRef]
  22. Ata, E.; Kiymaz, I.O. A study on certain properties of generalized special functions defined by Fox-Wright function. Appl. Math. Nonlinear Sci. 2020, 5(1), 147–162. [Google Scholar] [CrossRef]
  23. Ata, E. Generalized beta function defined by Wright function. arXiv 2018, arXiv:1803.03121. [Google Scholar]
  24. Ata, E. Modified special functions defined by generalized M-series and their properties. arXiv 2022, arXiv:2201.00867. [Google Scholar] [CrossRef]
  25. Ata, E.; Kiymaz, I.O. Special functions with general kernel: properties and applications to fractional partial differential equations. Int. J. Math. Comput. Eng. 2025, 3(2), 153–170. [Google Scholar] [CrossRef]
  26. Andrews, G.E.; Askey, R.; Roy, R. Special Functions; Cambridge University Press: Cambridge, 1999. [Google Scholar]
  27. Srivastava, H.M.; Karlsson, P.W. Multiple Gaussian Hypergeometric Series; Halsted Press: New York, 1985. [Google Scholar]
  28. Srivastava, H.M.; Manocha, H.L. A Treatise On Generating Functions; Halsted Press: New York, 1984. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2026 MDPI (Basel, Switzerland) unless otherwise stated