Inhomogeneous Whittaker Equation with Initial and Boundary Conditions

1. Introduction

The homogeneous Whittaker equation, first formulated in 1903, represents a classical second-order linear differential equation and is expressed in the form: [1]:

$${\displaystyle \frac{{\mathrm{d}}^2\mathrm{W}}{\mathrm{d}{z}^2}}+\left(-{\displaystyle \frac{1}{4}}+{\displaystyle \frac{\kappa }{z}}+{\displaystyle \frac{\frac{1}{4}-{\mu }^2}{z^2}}\right)\mathrm{W}=0$$

(1)

wherein $\kappa$ and $\mu$ are parameters and z and $\mathrm{W}$ are variables that could be real or complex. Whittaker [1] introduced the functions $\mathrm{M}\kappa ,\mu \left(z\right)$ and $\mathrm{W}\kappa ,\mu \left(z\right)$ as linearly independent solutions to the homogeneous Whittaker’s equation. The pairs of functions ${\mathrm{M}}_{\kappa ,\mu }\left(z\right)$, ${\mathrm{M}}_{\kappa ,-\mu }\left(z\right)$ and ${\mathrm{W}}_{\kappa ,\mu }\left(z\right),{\mathrm{W}}_{-\kappa ,\mu }\left(z\right)$ discussed below are linearly independent solutions of equation 1. The Wronskian of the Whittaker functions is provided in detail in [2]:

$$\begin{array}{c}\mathcal{W}\left\{{\mathrm{M}}_{\kappa ,\mu }\left(z\right),{\mathrm{M}}_{\kappa ,-\mu }\left(z\right)\right\}=-2\mu ,\\ \mathcal{W}\left\{{\mathrm{M}}_{\kappa ,\mu }\left(z\right),{\mathrm{W}}_{\kappa ,\mu }\left(z\right)\right\}=-{\displaystyle \frac{\Gamma (1+2\mu )}{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}},\\ \mathcal{W}\left\{{\mathrm{M}}_{\kappa ,\mu }\left(z\right),{\mathrm{W}}_{-\kappa ,\mu }\left({\mathrm{e}}^{\pm \pi \mathrm{i}}z\right)\right\}={\displaystyle \frac{\Gamma (1+2\mu )}{\Gamma \left(\frac{1}{2}+\mu +\kappa \right)}}{\mathrm{e}}^{\mp \left({\displaystyle \frac{1}{2}}+\mu \right)\pi \mathrm{i}},\\ \mathcal{W}\left\{{\mathrm{M}}_{\kappa ,-\mu }\left(z\right),{\mathrm{W}}_{\kappa ,\mu }\left(z\right)\right\}=-{\displaystyle \frac{\Gamma (1-2\mu )}{\Gamma \left(\frac{1}{2}-\mu -\kappa \right)}}.\end{array}$$

(2)

Equation 1 represents the reduced form of a degenerate hypergeometric equation and possesses a regular singular point at $z=0$ and an irregular singular point at $z=\infty$.

Whittaker’s equation 1 arises in various areas of physics, particularly in quantum mechanics, and plays a significant role in solving the radial Schrödinger equation in spherical coordinates, and in the solution to the wave equation in parabolic coordinates [3], in addition to describing the behaviour of charged particles in a Coulomb potential [4]. Akbarzadeh [5] utilized the Whittaker differential equation in deriving an exact, analytical solution for convective heat transfer of thermally fully developed laminar nanofluid flow in a circular tube, where the pipe wall is exposed to a constant temperature. Gupta and Bhengra [6] applied Whittaker functions to the derivation of the dispersion equation governing the propagation of torsional surface waves in an anisotropic layer sandwiched between two anisotropic inhomogeneous media. Conway [7] employed the Wronskian of the Whittaker functions to calculate indefinite integrals involving Whittaker’s functions and their products.

In the analysis of processes governed by time fractional diffusion and diffusion-wave equations, the Whittaker functions are important both as special functions and for their broad applications in mathematical physics, as they play a central role in the theory of uniform asymptotic expansions of differential equations with coalescing turning points and simple poles, as discussed in [2,8,9]. Mainardi et al. [10] compared Wright functions of the second kind with Whittaker functions in specific cases of fractional order. Their work in [10] underscores the importance of Whittaker functions in the context of higher transcendental functions. Szmytkowski [11] investigated the orthogonality of Whittaker functions of the second kind, ${\mathrm{W}}_{\kappa ,i\mu }\left(x\right)$, where $\mu \in \mathbb{R}$, with the weight function $w\left(x\right)={\displaystyle \frac{1}{x^2}}$. Chang et al. [12] investigated the asymptotic behaviour of the Whittaker function of the second kind for large values of the parameters and the independent variable z. Dunster [13] derived uniform asymptotic expansions for the Whittaker functions $M_{\kappa ,\mu }\left(z\right)$ and $W_{\kappa ,\mu }\left(z\right)$, as well as the numerically satisfactory companion function $W_{-\kappa ,\mu }\left(z{e}^{-\pi i}\right)$. The expansions are uniformly valid for $\mu \to \infty$ and for a specific ratio of the parameters $\kappa$ and $\mu$, with $0\le arg\left(z\right)\le \pi$. Using appropriate connection and analytic continuation, the expansions are extended to all unbounded non-zero complex values of z.

Izarra et al. [14] applied Pade Approximants in combination with Wynn’s algorithm on a specific asymptotic expansion to achieve precise numerical computations of the Whittaker functions ${\mathrm{W}}_{\kappa ,\mu }\left(z\right)$ for various values of the argument z and the parameters $\kappa$ and $\mu$. Ragab [15] systematically evaluated integrals involving products of Whittaker and Bessel functions.

A general solution to equation 1 can be expressed as a linear combination of these two solutions as:

$$\begin{array}{cc}\hfill & \mathrm{W}\left(z\right)=C_1{\mathrm{M}}_{\kappa ,\mu }\left(z\right)+C_2{\mathrm{W}}_{\kappa ,\mu }\left(z\right),\hfill \\ \hfill & 2\mu \ne -1,-2,\dots \hfill \end{array}$$

(3)

where $\kappa$ and $\mu$ are parameters.

The Whittaker functions can be expressed as [2]:

$$\begin{array}{cc}\hfill & {\mathrm{M}}_{\kappa ,\mu }\left(z\right)={\mathrm{e}}^{-{\displaystyle \frac{1}{2}}z}{z}^{{\displaystyle \frac{1}{2}}+\mu }M\left({\displaystyle \frac{1}{2}}+\mu -\kappa ,1+2\mu ,z\right),\hfill \\ \hfill & {\mathrm{W}}_{\kappa ,\mu }\left(z\right)={\mathrm{e}}^{-{\displaystyle \frac{1}{2}}z}{z}^{{\displaystyle \frac{1}{2}}+\mu }U\left({\displaystyle \frac{1}{2}}+\mu -\kappa ,1+2\mu ,z\right).\hfill \end{array}$$

(4)

Where $M(a,b,z):={}_{1}F_1(a;b;z)$ denotes Kummer’s confluent hypergeometric function of the first kind, and $U(a,b;z)$ denotes the confluent hypergeometric function of the second kind (also known as Tricomi’s function).

The function $M_{x,\mu }\left(z\right)$ is undefined when $2\mu =-1,-2,-3,\dots$ . Therefore from here on forth, its understood that $2\mu =-1,-2,-3,\dots$, unless otherwise specified. For specific values of these parameters, the Whittaker functions ${\mathrm{M}}_{\kappa ,\mu }\left(x\right)$ and ${\mathrm{W}}_{\kappa ,\mu }\left(x\right)$ can be simplified to various elementary and special functions, including modified Bessel functions, incomplete gamma functions, parabolic cylinder functions, error functions, logarithmic and cosine integrals, as well as generalized Hermite and Laguerre polynomials. For more details, see [16] and the references therein. In view of potential applications to engineering problems, we restrict our analysis of the inhomogeneous Whittaker equation to the case where the independent variable z is real, and hence

$${\mathrm{M}}_{\kappa ,\mu }\left(x\right)=x^{\mu -1/2}{e}^{-x/2}{}_{1}F_1\left(\left.\begin{array}{c}\mu -\kappa +{\displaystyle \frac{1}{2}}\\ 1+2\mu \end{array}\right|x\right).$$

(5)

The Whittaker function ${\mathrm{W}}_{\kappa ,\mu }\left(x\right)$ can be expressed in terms of the Whittaker functions of the first kind ${\mathrm{M}}_{\kappa ,\mu }\left(x\right)$, provided that $2\mu \notin \mathbb{Z}$, as follows [8,9]:

$${\mathrm{W}}_{\kappa ,\mu }\left(x\right)={\displaystyle \frac{\Gamma (-2\mu )}{\Gamma \left(\frac{1}{2}-k-\mu \right)}}{\mathrm{M}}_{\kappa ,\mu }\left(x\right)+{\displaystyle \frac{\Gamma \left(2\mu \right)}{\Gamma \left(\frac{1}{2}-k+\mu \right)}}{\mathrm{M}}_{\kappa ,-\mu }\left(x\right).$$

(6)

It is important to note that both ${\mathrm{M}}_{\kappa ,\mu }\left(x\right)$ and ${\mathrm{M}}_{\kappa ,-\mu }\left(x\right)$ constitute two linearly independent solutions of the homogeneous Whittaker equation.

In this work, we present a method for obtaining the general solution to Whittaker’s inhomogeneous equation. The approach is based on the introduction of an integral function, denoted by ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)$, whose evaluation relies on a generalization of the Whittaker integrals previously discussed by Appleblatt and Santandar[16]. To achieve this objective, the manuscript is organized as follows.

Section 2 introduces the inhomogeneous Whittaker equation and presents the function ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)$ as part of its solution. The section also outlines the properties of ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)$ and explains how it can be evaluated. In Section 3, we evaluate the function ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)$ for specific parameter values and various forms of the function $f\left(x\right)$. Section 4 presents the derivative of ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)$, derived using the known derivatives of the Whittaker functions. In Section 5, we provide solutions to the inhomogeneous Whittaker differential equation for both initial and boundary value problems. Finally, the conclusion summarizes the key results and suggests possible directions for future research.

2. Inhomogeneous Whittaker Equation

In this section, we will investigate the inhomogeneous Whittaker equation, which is given by:

$${\displaystyle \frac{{\mathrm{d}}^{2}W}{\mathrm{d}{x}^2}}+\left(-{\displaystyle \frac{1}{4}}+{\displaystyle \frac{\kappa }{x}}+{\displaystyle \frac{\frac{1}{4}-{\mu }^2}{x^2}}\right)W=f\left(x\right).$$

(7)

A particular solution to the differential equation can be found using the method of variation of parameters. To this end, we introduce the following function:

$${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)=\left[{\mathrm{M}}_{\kappa ,\mu }\left(x\right){\int }_0^{x}f\left(t\right){\mathrm{W}}_{\kappa ,\mu }\left(t\right)dt-{\mathrm{W}}_{\kappa ,\mu }\left(x\right){\int }_0^{x}f\left(t\right){\mathrm{M}}_{\kappa ,\mu }\left(t\right)dt\right].$$

(8)

The Wronskian of the Whittaker functions is given in equation 2 as

$$\mathcal{W}\left\{M_{\kappa ,\mu }\left(z\right),W_{\kappa ,\mu }\left(z\right)\right\}=-{\displaystyle \frac{\Gamma (1+2\mu )}{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}},$$

(9)

Therefore a particular solution of the inhomogeneous Whittaker equation 7 is :

$$\begin{array}{cc}\hfill & {\mathrm{Wp}}_{\mu ,\kappa }\left(z\right)={\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }\left(z\right)\hfill \\ \hfill & 2\mu \ne -1,-2,\dots \hfill \end{array}$$

(10)

Using equation 3 and 10 general solution of the inhomogeneous Whittaker equation 7 is given by :

$$\begin{array}{cc}\hfill & \mathrm{W}\left(z\right)=C_1{\mathrm{M}}_{\kappa ,\mu }\left(z\right)+C_2{\mathrm{W}}_{\kappa ,\mu }\left(z\right)+{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }\left(z\right),\hfill \\ \hfill & 2\mu \ne -1,-2,\dots \hfill \end{array}$$

(11)

The properties of the function ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)$ are influenced by the properties of the Whittaker functions for different values of the parameters $\kappa$ and $\mu$, as well as by the characteristics of the function $f\left(x\right)$ . Before presenting a solution to the inhomogeneous Whittaker equation, we will first investigate the function ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)$ for various parameter values and different forms of $f\left(x\right)$.

In their analysis of integrals involving Whittaker functions, Apelblat et al.[16] define the integral functions ${\mathrm{Mi}}_{\kappa ,\mu }\left(x\right)$ and ${\mathrm{Wi}}_{\kappa ,\mu }\left(x\right)$ that can offer examples of the function ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)$ for various parameter values, especially when $f\left(x\right)={\displaystyle \frac{1}{x}}$. These Whittaker integral functions are define as follows [16]:

$${\mathrm{Mi}}_{\kappa ,\mu }\left(x\right)={\int }_0^x{\displaystyle \frac{{\mathrm{M}}_{\kappa ,\mu }\left(t\right)}{t}}dt$$

(12)

and

$${\mathrm{Wi}}_{\kappa ,\mu }\left(x\right)={\int }_0^x{\displaystyle \frac{{\mathrm{W}}_{\kappa ,\mu }\left(t\right)}{t}}dt$$

(13)

For simplicity, we will generalize the notation from [16] as follows:

$$\begin{array}{c}{}^f\mathrm{Mi}_{\kappa ,\mu }\left(x\right)={\int }_0^{x}f\left(t\right){\mathrm{M}}_{\kappa ,\mu }\left(t\right)dt,\\ {}^f\mathrm{Wi}_{\kappa ,\mu }\left(x\right)={\int }_0^{x}f\left(t\right){\mathrm{W}}_{\kappa ,\mu }\left(t\right)dt.\end{array}.$$

(14)

Hence, in the newly introduced generalized notation, the functions

${\mathrm{Mi}}_{\kappa ,\mu }\left(x\right)$ and ${\mathrm{Wi}}_{\kappa ,\mu }\left(x\right)$ are represented as

$$\begin{array}{c}{\mathrm{Mi}}_{\kappa ,\mu }\left(x\right)={}^{{\displaystyle \frac{1}{t}}}\mathrm{Mi}_{\kappa ,\mu }\left(x\right),\\ {\mathrm{Wi}}_{\kappa ,\mu }\left(x\right)={}^{{\displaystyle \frac{1}{t}}}\mathrm{Wi}_{\kappa ,\mu }\left(x\right).\end{array}$$

(15)

Therefore, we can rewrite the function ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(z\right)$ given by equation 8 as follows:

$${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)={\mathrm{M}}_{\kappa ,\mu }\left(x\right){}^f\mathrm{Mi}_{\kappa ,\mu }\left(x\right)-{\mathrm{W}}_{\kappa ,\mu }\left(x\right){}^f\mathrm{Wi}_{\kappa ,\mu }\left(x\right).$$

(16)

It is straight forward that the ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(0\right)=0$. Apelblat et al. [16] used the Mathematica program to express the Whittaker integral functions in terms of elementary special functions and obtained specific cases of the Whittaker integral functions for different values of the parameters $\mu$ and $\kappa$. We used these expressions to derive the following table for the function ${}^{{\displaystyle \frac{1}{x}}}\mathrm{Zi}_{\kappa ,\mu }\left(x\right)$:

Table 1. Example of ${}^{{\displaystyle \frac{1}{x}}}\mathrm{Zi}_{\kappa ,\mu }\left(z\right)$

$\mathit{\kappa }$	$\mathit{\mu }$	${}^{\frac{1}{\mathit{x}}}\mathbf{Zi}_{\mathit{\kappa },\mathit{\mu }}\left(\mathit{x}\right)$
${\displaystyle \frac{3}{2}}$	0	$2\sqrt{x}{e}^{-x/2}\left({\mathbf{W}}_{{\displaystyle \frac{3}{2}},0}\left(x\right)-{\mathbf{M}}_{{\displaystyle \frac{3}{2}},0}\left(x\right)\right)$
2	${\displaystyle \frac{1}{2}}$	$x{e}^{-x/2}\left({\mathbf{W}}_{2,{\displaystyle \frac{1}{2}}}\left(x\right)-\mathbf{2}{M}_{2,{\displaystyle \frac{1}{2}}}\left(x\right)\right)$

Table 1. Example of ${}^{{\displaystyle \frac{1}{x}}}\mathrm{Zi}_{\kappa ,\mu }\left(z\right)$

$\mathit{\kappa }$	$\mathit{\mu }$	${}^{\frac{1}{\mathit{x}}}\mathbf{Zi}_{\mathit{\kappa },\mathit{\mu }}\left(\mathit{x}\right)$
${\displaystyle \frac{3}{2}}$	0	$2\sqrt{x}{e}^{-x/2}\left({\mathbf{W}}_{{\displaystyle \frac{3}{2}},0}\left(x\right)-{\mathbf{M}}_{{\displaystyle \frac{3}{2}},0}\left(x\right)\right)$
2	${\displaystyle \frac{1}{2}}$	$x{e}^{-x/2}\left({\mathbf{W}}_{2,{\displaystyle \frac{1}{2}}}\left(x\right)-\mathbf{2}{M}_{2,{\displaystyle \frac{1}{2}}}\left(x\right)\right)$

Apelblat et al.[16] used the following recurrence relations between the Whittaker functions [8,9]:

$$2\mu \left[{\mathrm{M}}_{\kappa -1/2,\mu -1/2}\left(t\right)-{\mathrm{M}}_{\kappa +1/2,\mu -1/2}\left(t\right)\right]=t^{1/2}{\mathrm{M}}_{\kappa ,\mu }\left(t\right).$$

(17)

$$(\kappa +\mu ){\mathrm{W}}_{\kappa -1/2,\mu }\left(t\right)+{\mathrm{W}}_{\kappa +1/2,\mu }\left(t\right)=t^{1/2}{\mathrm{W}}_{\kappa ,\mu +1/2}\left(t\right).$$

(18)

By rearranging this expression, Equation 18 can thus be reformulated as

$$t^{1/2}{\mathrm{W}}_{\kappa ,\mu }\left(t\right)=\left(\kappa +\mu -{\displaystyle \frac{1}{2}}\right){\mathrm{W}}_{\kappa -1/2,\mu -1/2}\left(t\right)+{\mathrm{W}}_{\kappa +1/2,\mu -1/2}\left(t\right)$$

(19)

Consequently, we express integrals involving Whittaker functions in terms of the Whittaker functions ${\mathrm{Mi}}_{\kappa ,\mu }\left(x\right)$ and ${\mathrm{Wi}}_{\kappa ,\mu }\left(x\right)$ in the following relation

$$\begin{array}{cc}\hfill & {\int }_0^x{\displaystyle \frac{{\mathrm{M}}_{\kappa ,\mu }\left(t\right)}{t^{1/2}}}dt=2\mu \left[{\mathrm{Mi}}_{\kappa -1/2,\mu -1/2}\left(t\right)-{\mathrm{Mi}}_{\kappa +1/2,\mu -1/2}\left(t\right)\right],\hfill \\ \hfill & {\int }_0^x{\displaystyle \frac{{\mathrm{W}}_{\kappa ,\mu }\left(t\right)}{t^{1/2}}}dt=(\kappa +\mu -{\displaystyle \frac{1}{2}}){\mathrm{Wi}}_{\kappa -1/2,\mu -1/2}\left(t\right)+{\mathrm{Wi}}_{\kappa +1/2,\mu -1/2}\left(t\right).\hfill \end{array}$$

(20)

Equation 20 can be used to derive an expression for the function ${}^{t^{{\displaystyle \frac{1}{2}}}}\mathbf{Zi}_{\kappa ,\mu }\left(t\right)$ as follow:

$$\begin{array}{c}\hfill {}^{t^{{\displaystyle \frac{1}{2}}}}\mathbf{Zi}_{\kappa ,\mu }\left(x\right)={\mathrm{M}}_{\kappa ,\mu }\left(x\right)\left[(\kappa +\mu -{\displaystyle \frac{1}{2}}){\mathrm{Wi}}_{\kappa -1/2,\mu -{\displaystyle \frac{1}{2}}}\left(t\right)+{\mathrm{Wi}}_{\kappa +1/2,\mu -{\displaystyle \frac{1}{2}}}\left(t\right)\right]\\ \hfill -2\mu {\mathrm{W}}_{\kappa ,\mu }\left(x\right)\left[{\mathrm{Mi}}_{\kappa -1/2,\mu -1/2}\left(t\right)-{\mathrm{Mi}}_{\kappa +1/2,\mu -1/2}\left(t\right)\right].\end{array}$$

(21)

3. Values of the Function ${}^{\mathit{f}}\mathbf{Zi}_{\mathit{\kappa }\mathbf{,}\mathit{\mu }}\mathbf{\left(}\mathit{x}\mathbf{\right)}$ for Specific Values of $\mathit{f}\mathbf{\left(}\mathit{x}\mathbf{\right)}$ and Parameters $\mathit{\kappa }$ and $\mathit{\mu }$

In this section, we derive closed-form expressions for the function ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)$. The resulting formulas are obtained either by leveraging established relations and known identities of the Whittaker functions or through symbolic computation using Maple to evaluate the relevant integral expressions.

• $\kappa =0$ and $\mu ={\displaystyle \frac{1}{2}}$:
  from [2] We have:

  $${\mathrm{M}}_{0,{\displaystyle \frac{1}{2}}}\left(2\mathrm{x}\right)=2sinhx.$$

  (22)

  $${\mathrm{W}}_{0,{\displaystyle \frac{1}{2}}}\left(\mathrm{x}\right)=e^{{\displaystyle \frac{-x}{2}}}.$$

  (23)
  Hence, the following closed-form expressions for the function ${}^f\mathrm{Zi}_{0,{\displaystyle \frac{1}{2}}}\left(x\right)$can be derived:

  $${}^1\mathrm{Zi}_{0,{\displaystyle \frac{1}{2}}}\left(z\right)=2\left(e^{{\displaystyle \frac{x}{2}}}+e^{-{\displaystyle \frac{x}{2}}}\right)-4$$

  (24)

  $${}^t\mathrm{Zi}_{0,{\displaystyle \frac{1}{2}}}\left(z\right)=4\left(e^{{\displaystyle \frac{x}{2}}}-e^{-{\displaystyle \frac{x}{2}}}\right)-4x$$

  (25)

  $${}^{t^2}\mathrm{Zi}_{0,{\displaystyle \frac{1}{2}}}\left(z\right)=16\left(e^{{\displaystyle \frac{x}{2}}}-e^{-{\displaystyle \frac{x}{2}}}\right)-4x^2-32$$

  (26)
  Additionally, by employing Maple computations, we derive the following closed-form expressions for the function ${}^f\mathrm{Zi}_{0,{\displaystyle \frac{1}{2}}}\left(x\right)$:

  $${}^{e^{-t}}\mathrm{Zi}_{0,{\displaystyle \frac{1}{2}}}\left(z\right)=0.6666666665e^{{\displaystyle \frac{x}{2}}}+1.333333334e^{-x}-2e^{-{\displaystyle \frac{x}{2}}}-2\times {10}^{-10}{e}^{-2x}$$

  (27)

  $${}^{e^t}\mathrm{Zi}_{0,{\displaystyle \frac{1}{2}}}\left(z\right)=-2e^{{\displaystyle \frac{x}{2}}}+0.666666667e^{-{\displaystyle \frac{x}{2}}}+1.33333333333e^x-6\times {10}^{-10}{e}^{-x}+1\times {10}^{-9}$$

  (28)
• $\mu =\kappa -{\displaystyle \frac{1}{2}}$:
  from [2] we have:

  $${\mathrm{M}}_{\kappa ,\kappa -{\displaystyle \frac{1}{2}}}\left(x\right)={\mathrm{W}}_{\kappa ,\kappa -{\displaystyle \frac{1}{2}}}\left(x\right)={\mathrm{W}}_{\kappa ,-\kappa +{\displaystyle \frac{1}{2}}}\left(x\right)={\mathrm{e}}^{-{\displaystyle \frac{1}{2}}x}{x}^{\kappa }.$$

  (29)
  And from [16] we have:

  $${}^{{\displaystyle \frac{1}{t}}}\mathrm{Mi}_{\kappa ,\kappa -1/2}\left(x\right)={}^{{\displaystyle \frac{1}{t}}}\mathrm{Wi}_{\kappa ,\kappa -1/2}\left(x\right)={}^{{\displaystyle \frac{1}{t}}}\mathrm{Wi}_{\kappa ,-\kappa +1/2}\left(x\right)=2^{\kappa }\gamma \left(\kappa ,{\displaystyle \frac{x}{2}}\right).$$

  (30)
  Where $\gamma (a,x)$ is the lower incomplete gamma function defined by

  $$\gamma (a,x)={\int }_0^x{t}^{a-1}{\mathrm{e}}^{-t}\mathrm{d}t.$$

  (31)
  The lower incomplete gamma function can be evaluated from Maple using:

  $$\gamma (a,x)+\Gamma (a,x)=\Gamma \left(a\right).$$

  (32)
  Where $\Gamma (a,x)$ represents the upper incomplete gamma function, defined by:

  $$\Gamma (a,x)={\int }_x^{\infty }{t}^{a-1}{\mathrm{e}}^{-t}\mathrm{d}t.$$

  (33)
  And equation 29 will leads to the following result:

  $${}^f\mathbf{Zi}_{\kappa ,\kappa -1/2}\left(t\right)=0$$

  (34)
• $\kappa =\mu -{\displaystyle \frac{1}{2}}$:
  from [2] we have:

  $$\begin{array}{c}{M}_{\mu -{\displaystyle \frac{1}{2}},\mu }\left(x\right)=2\mu {\mathrm{e}}^{{\displaystyle \frac{1}{2}}x}{x}^{{\displaystyle \frac{1}{2}}-\mu }\gamma (2\mu ,x)\\ W_{\mu -{\displaystyle \frac{1}{2}},\mu }\left(x\right)={\mathrm{e}}^{{\displaystyle \frac{1}{2}}x}{x}^{{\displaystyle \frac{1}{2}}-\mu }\Gamma (2\mu ,x).\end{array}$$

  (35)
  Here, $\gamma (a,x)$ denotes the lower incomplete gamma function as defined in Equation 31, while $\Gamma (a,x)$ represents the Upper incomplete gamma function defiend in Equation 33.
  Using Maple we get the following:

  $$\begin{array}{cc}\hfill {}^1\mathrm{Zi}_{\mu -{\displaystyle \frac{1}{2}},\mu }\left(x\right)=& -x^{-\mu +{\displaystyle \frac{1}{2}}}·e^{0.5x}·[\hfill \\ \hfill & \mu \left({\int }_0^x{e}^{0.5t}{t}^{-\mu +{\displaystyle \frac{1}{2}}}\Gamma (2\mu ,t)\mathrm{d}t\right)\left(-2\Gamma \left(2\mu \right)+2\Gamma (2\mu ,x)\right)\hfill \\ \hfill & +\Gamma (2\mu ,x)\left({\int }_0^x\mu e^{0.5t}{t}^{-\mu +{\displaystyle \frac{1}{2}}}\left(2\Gamma \left(2\mu \right)-2\Gamma (2\mu ,t)\right)\mathrm{d}t\right)]\hfill \end{array}$$

  (36)

  $$\begin{array}{cc}\hfill & {}^x\mathrm{Zi}_{\mu -{\displaystyle \frac{1}{2}},\mu }\left(x\right)=-e^{0.5x}{x}^{-\mu +{\displaystyle \frac{1}{2}}}·[\hfill \\ \hfill & -\left({\int }_0^x{t}^{{\displaystyle \frac{3}{2}}-\mu }{e}^{0.5t}\Gamma (2\mu ,t)\mathrm{d}t\right)\left(2\Gamma \left(2\mu \right)-2\Gamma (2\mu ,x)\right)\mu \hfill \\ \hfill & +\Gamma (2\mu ,x)\left({\int }_0^x\mu t^{{\displaystyle \frac{3}{2}}-\mu }{e}^{0.5t}\left(2\Gamma \left(2\mu \right)-2\Gamma (2\mu ,t)\right)\mathrm{d}t\right)].\hfill \end{array}$$

  (37)
• $\kappa =0$:
  from [2] we have:

  $${\mathrm{M}}_{0,\mu }\left(x\right)=2^{2\mu +{\displaystyle \frac{1}{2}}}\Gamma (1+\mu )\sqrt{{\displaystyle \frac{x}{2}}}{I}_{\mu }\left({\displaystyle \frac{x}{2}}\right),$$

  (38)

  $${\mathrm{W}}_{0,\mu }\left(x\right)=\sqrt{{\displaystyle \frac{x}{\pi }}}{K}_{\mu }\left({\displaystyle \frac{x}{2}}\right)$$

  (39)

  and [16] to derive:

  $${}^{{\displaystyle \frac{1}{t}}}\mathrm{Mi}_{0,\mu }\left(x\right)={\displaystyle \frac{x^{\mu +1/2}}{\mu +1/2}}{}_{1}F_2\left(\left.\begin{array}{c}{\displaystyle \frac{2\mu +1}{4}}\\ \mu +1,{\displaystyle \frac{2\mu +5}{4}}\end{array}\right|\begin{array}{c}{\displaystyle \frac{x^2}{16}}\\ \end{array}\right)$$

  (40)

  $${}^{{\displaystyle \frac{1}{t}}}\mathrm{Wi}_{0,\mu }\left(x\right)={\displaystyle \frac{\sqrt{\pi }}{2sin\pi \mu }}\left[{\displaystyle \frac{4^{\mu }{}^{\frac{1}{t}}\mathrm{Mi}_{0,-\mu }\left(x\right)}{\Gamma (1-\mu )}}-{\displaystyle \frac{4^{-\mu }{}^{\frac{1}{t}}\mathrm{Mi}_{0,\mu }\left(x\right)}{\Gamma (1+\mu )}}\right].$$

  (41)
  Therefore,

  $$\begin{array}{cc}\hfill {}^{{\displaystyle \frac{1}{t}}}\mathrm{Zi}_{0,\mu }\left(x\right)=& 2^{2\mu +{\displaystyle \frac{1}{2}}}\Gamma (1+\mu )\sqrt{{\displaystyle \frac{x}{2}}}{I}_{\mu }\left({\displaystyle \frac{x}{2}}\right)·{\displaystyle \frac{\sqrt{\pi }}{2sin\left(\pi \mu \right)}}\hfill \\ \hfill & \times \left[{\displaystyle \frac{4^{\mu }{}^{\frac{1}{t}}\mathrm{Mi}_{0,-\mu }\left(x\right)}{\Gamma (1-\mu )}}-{\displaystyle \frac{4^{-\mu }{}^{\frac{1}{t}}\mathrm{Mi}_{0,\mu }\left(x\right)}{\Gamma (1+\mu )}}\right]\hfill \\ \hfill & -\sqrt{{\displaystyle \frac{x}{\pi }}}{K}_{\mu }\left({\displaystyle \frac{x}{2}}\right)·{\displaystyle \frac{x^{\mu +1/2}}{\mu +1/2}}{}_{1}F_2\left(\left.\begin{array}{c}{\displaystyle \frac{2\mu +1}{4}}\\ \mu +1,{\displaystyle \frac{2\mu +5}{4}}\end{array}\right|{\displaystyle \frac{x^2}{16}}\right),\hfill \end{array}$$

  (42)

  and using Maple we get:

  $$\begin{array}{cc}\hfill & {}^1\mathrm{Zi}_{0,\mu }\left(x\right)=-{\displaystyle \frac{1}{\sqrt{\pi }sin\left(\pi \mu \right)(4{\mu }^2-9)}}[2x^{-\mu +2}·{16}^{\mu }·\Gamma {\left(\mu \right)}^2·\mu ·\hfill \\ & I_{\mu }\left(\frac{x}{2}\right)·\left(\frac{3}{2}+\mu \right)·{}_{2}F_2\left(\left.\begin{array}{c}{\displaystyle \frac{3}{4}}-{\displaystyle \frac{\mu }{2}}\\ 1-\mu ,{\displaystyle \frac{7}{4}}-{\displaystyle \frac{\mu }{2}}\end{array}\right|{\displaystyle \frac{x^2}{16}}\right)+x^{\mu +2}·\pi ·\left(\mu -\frac{3}{2}\right)\hfill \\ & ·I_{\mu }\left(\frac{x}{2}\right)·{}_{2}F_2\left(\left.\begin{array}{c}{\displaystyle \frac{3}{4}}+{\displaystyle \frac{\mu }{2}}\\ 1+\mu ,{\displaystyle \frac{7}{4}}+{\displaystyle \frac{\mu }{2}}\end{array}\right|{\displaystyle \frac{x^2}{16}}\right)+{\displaystyle \frac{\Gamma \left(\mu \right)·sin\left(\pi \mu \right)·\mu ·\sqrt{x}}{2}}·K_{\mu }\left(\frac{x}{2}\right)\hfill \\ \hfill & ·\left({\int }_0^x\sqrt{t}·I_{\mu }\left(\frac{t}{2}\right)dt\right)·\left(4+2{\mu }^3-{94}^{\mu }\right)]\hfill \end{array}$$

  (43)
• $\kappa =-{\displaystyle \frac{1}{4}}$ and $\mu ={\displaystyle \frac{1}{4}}$: from [2] we have:

  $$\begin{array}{cc}\hfill & {\mathrm{M}}_{-{\displaystyle \frac{1}{4}},{\displaystyle \frac{1}{4}}}\left(x^2\right)={\displaystyle \frac{1}{2}}{\mathrm{e}}^{{\displaystyle \frac{1}{2}}{x}^2}\sqrt{\pi x}erf\left(x\right)\hfill \\ \hfill & {\mathrm{W}}_{-{\displaystyle \frac{1}{4}},\pm {\displaystyle \frac{1}{4}}}\left(x^2\right)={\mathrm{e}}^{{\displaystyle \frac{1}{2}}{x}^2}\sqrt{\pi x}erfc\left(x\right)\hfill \end{array}$$

  (44)
  Furthermore, using Maple symbolic computation, we derive the following closed-form expressions for the function ${}^f\mathrm{Zi}_{-{\displaystyle \frac{1}{4}},{\displaystyle \frac{1}{4}}}\left(x\right)$:

  $$\begin{array}{cc}\hfill {}^1\mathrm{Zi}_{-{\displaystyle \frac{1}{4}},{\displaystyle \frac{1}{4}}}\left(x\right)=& -{\displaystyle \frac{1}{2}}{\mathrm{e}}^{{\displaystyle \frac{x}{2}}}\sqrt{\pi x}[erf\left(\sqrt{x}\right)\left({\int }_0^x{\mathrm{e}}^{{\displaystyle \frac{t}{2}}}\sqrt{\pi t}(-1+erf\left(\sqrt{t}\right))\mathrm{d}t\right)\hfill \\ \hfill & -2\left({\int }_0^x{\mathrm{e}}^{{\displaystyle \frac{t}{2}}}\sqrt{\pi t}erf\left(\sqrt{t}\right)\mathrm{d}t\right)erf\left(\sqrt{x}\right)+2\left({\int }_0^x{\mathrm{e}}^{{\displaystyle \frac{t}{2}}}\sqrt{\pi t}erf\left(\sqrt{t}\right)\mathrm{d}t\right)]\hfill \end{array}$$

  (45)

  $$\begin{array}{cc}\hfill {}^{{\displaystyle \frac{1}{x}}}\mathrm{Zi}_{-{\displaystyle \frac{1}{4}},{\displaystyle \frac{1}{4}}}\left(x\right)=-{\displaystyle \frac{1}{2}}{e}^{{\displaystyle \frac{x}{2}}}\pi x^{1/4}[& erf\left(\sqrt{x}\right){\int }_0^x{\displaystyle \frac{(erf\left(\sqrt{t}\right)-1)e^{\frac{t}{2}}}{t^{3/4}}}dt\hfill \\ \hfill & -erf\left(\sqrt{x}\right){\int }_0^x{\displaystyle \frac{e^{\frac{t}{2}}erf\left(\sqrt{t}\right)}{t^{3/4}}}dt+{\int }_0^x{\displaystyle \frac{e^{\frac{t}{2}}erf\left(\sqrt{t}\right)}{t^{3/4}}}dt]\hfill \end{array}$$

  (46)

  $$\begin{array}{cc}\hfill {}^{{\displaystyle \frac{1}{\sqrt{x}}}}\mathrm{Zi}_{-{\displaystyle \frac{1}{4}},{\displaystyle \frac{1}{4}}}\left(x\right)=-{\displaystyle \frac{1}{2}}{e}^{{\displaystyle \frac{x}{2}}}\pi x^{1/4}[& erf\left(\sqrt{x}\right){\int }_0^x{\displaystyle \frac{(erf\left(\sqrt{t}\right)-1)e^{\frac{t}{2}}}{t^{1/4}}}dt\hfill \\ \hfill & -erf\left(\sqrt{x}\right){\int }_0^x{\displaystyle \frac{e^{\frac{t}{2}}erf\left(\sqrt{t}\right)}{t^{1/4}}}dt+{\int }_0^x{\displaystyle \frac{e^{\frac{t}{2}}erf\left(\sqrt{t}\right)}{t^{1/4}}}dt]\hfill \end{array}$$

  (47)

  $$\begin{array}{cc}\hfill {}^{e^{-x/2}}\mathrm{Zi}_{-{\displaystyle \frac{1}{4}},{\displaystyle \frac{1}{4}}}\left(x\right)={\displaystyle \frac{4}{7}}[& -{\displaystyle \frac{7\pi x^{1/4}erf\left(\sqrt{x}\right)}{8}}{\int }_0^x{t}^{1/4}\left(erf\left(\sqrt{t}\right)-1\right)dt\hfill \\ \hfill & +\sqrt{\pi }{x^2}_2{F}_2\left(\left.\begin{array}{c}{\displaystyle \frac{1}{2}},{\displaystyle \frac{7}{4}}\\ {\displaystyle \frac{3}{2}},{\displaystyle \frac{11}{4}}\end{array}\right|-x\right)\left(erf\left(\sqrt{x}\right)-1\right)]e^{x/2}\hfill \end{array}$$

  (48)
• $\kappa =-{\displaystyle \frac{1}{2}}$ and $\mu ={\displaystyle \frac{1}{4}}$:
  Using Maple we get the following:

  $$\begin{array}{cc}\hfill & {}^1\mathrm{Zi}_{-1/2,1/4}\left(x\right)={\displaystyle \frac{x^{9/4}}{9\sqrt{\pi }}}\left[-{\displaystyle \frac{24x^{3/2}}{11}}\left(K_{1/4}\left({\displaystyle \frac{x}{2}}\right)-K_{3/4}\left({\displaystyle \frac{x}{2}}\right)\right)·{}_{2}F_2\left(\left.\begin{array}{c}{\displaystyle \frac{11}{8}}\\ {\displaystyle \frac{7}{4}},{\displaystyle \frac{19}{8}}\end{array}\right|\begin{array}{c}{\displaystyle \frac{x^2}{16}}\\ \mid \end{array}\right)\right.\hfill \\ \hfill & +{\displaystyle \frac{9\Gamma {\left(\frac{3}{4}\right)}^2}{5}}\left(I_{3/4}\left({\displaystyle \frac{x}{2}}\right)+I_{-1/4}\left({\displaystyle \frac{x}{2}}\right)\right)·{}_{2}F_2\left(\left.\begin{array}{c}{\displaystyle \frac{5}{8}}\\ {\displaystyle \frac{1}{4}},{\displaystyle \frac{19}{8}}\end{array}\right|\begin{array}{c}{\displaystyle \frac{x^2}{16}}\\ \end{array}\right)\hfill \\ \hfill & -{\displaystyle \frac{9\sqrt{x}}{7}}\left(K_{1/4}\left({\displaystyle \frac{x}{2}}\right)-K_{3/4}\left({\displaystyle \frac{x}{2}}\right)\right)·{}_{2}F_2\left(\left.\begin{array}{c}{\displaystyle \frac{7}{8}}\\ {\displaystyle \frac{3}{4}},{\displaystyle \frac{15}{8}}\end{array}\right|\begin{array}{c}{\displaystyle \frac{x^2}{16}}\\ 16\end{array}\right)\hfill \\ \hfill & +\left\{x\Gamma {\left({\displaystyle \frac{3}{4}}\right)}^2·{}_{2}F_2\left(\left.\begin{array}{c}{\displaystyle \frac{9}{8}}\\ 2,{\displaystyle \frac{17}{8}}\end{array}\right|\begin{array}{c}{\displaystyle \frac{x^2}{16}}\\ \end{array}\right)\right.-{\displaystyle \frac{3\pi \sqrt{2}{x}^{3/2}}{22}}·{}_{2}F_2\left(\left.\begin{array}{c}{\displaystyle \frac{11}{8}}\\ 2,{\displaystyle \frac{19}{8}}\end{array}\right|\begin{array}{c}{\displaystyle \frac{x^2}{16}}\\ \end{array}\right)\hfill \\ \hfill & \left.-{\displaystyle \frac{9\pi \sqrt{2}\sqrt{x}}{14}}·{}_{2}F_2\left(\left.\begin{array}{c}{\displaystyle \frac{7}{8}}\\ {\displaystyle \frac{3}{4}},{\displaystyle \frac{17}{8}}\end{array}\right|\begin{array}{c}{x}^2\\ {\displaystyle \frac{16}{16}}\end{array}\right)·\left(I_{3/4}\left({\displaystyle \frac{x}{2}}\right)+I_{-1/4}\left({\displaystyle \frac{x}{2}}\right)\right)\right]\hfill \end{array}$$

  (49)

Figure 1 and Figure 2 illustrate the graph of the function ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(x\right)$ for various values of the function $f\left(x\right)$ .

4. Derivatives of the Function ${}^{\mathit{f}}\mathbf{Zi}_{\mathit{\kappa },\mathit{\mu }}\left(\mathit{x}\right)$

The derivatives of the Whittaker functions are given in [2]

$${\mathrm{M}}_{\kappa ,\mu }^{\prime }\left(x\right)=\left({\displaystyle \frac{1}{2}}-{\displaystyle \frac{\kappa }{x}}\right){\mathrm{M}}_{\kappa ,\mu }\left(x\right)+\left({\displaystyle \frac{1}{2}}+\mu +\kappa \right){\displaystyle \frac{{\mathrm{M}}_{\kappa +1,\mu }\left(x\right)}{x}}.$$

(50)

$${\mathrm{W}}_{\kappa ,\mu }^{\prime }\left(x\right)=\left({\displaystyle \frac{1}{2}}-{\displaystyle \frac{\kappa }{x}}\right){\mathrm{W}}_{\kappa ,\mu }\left(x\right)-{\displaystyle \frac{{\mathrm{W}}_{\kappa +1,\mu }\left(x\right)}{x}}.$$

(51)

the derivative of the function ${}^{f}Z_{\kappa ,\mu }\left(x\right)$ is equal to

$${}^{f}Z_{\kappa ,\mu }^{\prime }\left(z\right)=M_{\kappa ,\mu }^{\prime }\left(z\right){}^f\mathrm{Wi}_{\kappa ,\mu }\left(z\right)-W_{\kappa ,\mu }{\left(z\right)}^{\prime }{}^f\mathrm{Mi}_{\kappa ,\mu }\left(z\right)$$

(52)

Using equations 50 and 51 in 52 we reach the following form of the derivative of the function ${}^{f}Z_{\kappa ,\mu }\left(z\right)$

$$\begin{array}{cc}\hfill {}^f\mathrm{Z}_{\kappa ,\mu }^{\prime }\left(x\right)& =\left[\left({\displaystyle \frac{1}{2}}-{\displaystyle \frac{\kappa }{x}}\right){\mathrm{M}}_{\kappa ,\mu }\left(x\right)+\left({\displaystyle \frac{1}{2}}+\mu +\kappa \right){\displaystyle \frac{{\mathrm{M}}_{\kappa +1,\mu }\left(x\right)}{x}}\right]{}^f\mathrm{Wi}_{\kappa ,\mu }\left(x\right)\hfill \\ \hfill & -\left[\left({\displaystyle \frac{1}{2}}-{\displaystyle \frac{\kappa }{x}}\right){\mathrm{W}}_{\kappa ,\mu }\left(x\right)-{\displaystyle \frac{{\mathrm{W}}_{\kappa +1,\mu }\left(x\right)}{x}}\right]{}^f\mathrm{Mi}_{\kappa ,\mu }\left(x\right).\hfill \end{array}$$

(53)

5. Initial and Boundary Value Problems

In this section, we explore examples of initial and boundary value problems related to the inhomogeneous Whittaker equation. To compute the solutions efficiently, we will employ Maple for the necessary calculations and provide the corresponding graphical representations. By doing so, we aim to illustrate the behaviour of the solutions under various conditions and gain insights into the underlying mathematical structures of the equation.

5.1. Initial Value Problems

Consider Inhomogeneous Whittaker equation 7 with the initial conditions:

$$W\left(a\right)=\alpha ;W^{\prime }\left(a\right)=\beta ,$$

(54)

where $\alpha$ and $\beta$ are known constants. the solution of equation 7 is given by 11.

By applying the initial conditions, we obtain:

$$\begin{array}{cc}\hfill & C_1{\mathrm{M}}_{\kappa ,\mu }\left(a\right)+C_2{\mathrm{W}}_{\kappa ,\mu }\left(a\right)+{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }\left(a\right)=\alpha ,\hfill \\ \hfill & C_1{\mathrm{M}}_{\kappa ,\mu }^{\prime }\left(a\right)+C_2{\mathrm{W}}_{\kappa ,\mu }^{\prime }\left(a\right)+{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }^{\prime }\left(a\right)=\beta ,\hfill \end{array}$$

(55)

which can be rewritten in matrix form as:

$$\left[\begin{array}{cc}{\mathrm{M}}_{k,\mu }\left(a\right)\hfill & {\mathrm{W}}_{k,\mu }\left(a\right)\hfill \\ {\mathrm{M}}_{k,\mu }^{\prime }\left(a\right)\hfill & {\mathrm{W}}_{k,\mu }^{\prime }\left(a\right)\hfill \end{array}\right]\left[\begin{array}{c}{C}_1\hfill \\ C_2\hfill \end{array}\right]=\left[\begin{array}{c}\alpha -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }\left(a\right)\hfill \\ \beta -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }^{\prime }\left(a\right)\hfill \end{array}\right].$$

(56)

Upon solving the system we get:

$$\begin{array}{cc}\hfill C_1={\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}[& {\mathrm{W}}_{\kappa ,\mu }\left(a\right)\left(\beta -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }^{\prime }\left(a\right)\right)\hfill \\ \hfill -& {\mathrm{W}}_{\kappa ,\mu }^{\prime }\left(a\right)\left(\alpha -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }\left(a\right)\right)]\hfill \end{array}$$

(57)

$$\begin{array}{cc}\hfill C_2={\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}[& {\mathrm{M}}_{\kappa ,\mu }^{\prime }\left(a\right)\left(\alpha -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }\left(a\right)\right)\hfill \\ \hfill -& {\mathrm{M}}_{\kappa ,\mu }\left(a\right)\left(\beta -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }^{\prime }\left(a\right)\right)]\hfill \end{array}$$

(58)

Since ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(0\right)={}^f\mathrm{Zi}_{\kappa ,\mu }^{\prime }\left(0\right)=0$, the value of the constants $C_1$ and $C_2$ in case $a=0$ are:

$$C_1={\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}\left[\beta {\mathrm{W}}_{k,\mu }\left(0\right)-\alpha {\mathrm{W}}_{k,\mu }^{\prime }\left(0\right)\right]$$

(59)

$$C_2={\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}\left[\alpha {\mathrm{M}}_{k,\mu }^{\prime }\left(0\right)-\beta {\mathrm{M}}_{k,\mu }\left(0\right)\right]$$

(60)

The limiting behaviour of the Whittaker functions for $x\to 0$ and real parameters [2]:

$$\begin{array}{cc}\hfill & {\mathrm{M}}_{\kappa ,\mu }\left(x\right)=x^{\mu +{\displaystyle \frac{1}{2}}}(1+O\left(x\right))\hfill \\ \hfill & 2\mu \ne -1,-2,-3,\dots \hfill \end{array}$$

(61)

$$\begin{array}{cc}{\mathrm{W}}_{\kappa ,\mu }\left(x\right)={\displaystyle \frac{\Gamma \left(2\mu \right)}{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}}{x}^{{\displaystyle \frac{1}{2}}-\mu }+O\left(x^{{\displaystyle \frac{3}{2}}-\mu }\right),& \mu \ge {\displaystyle \frac{1}{2}},\mu \ne {\displaystyle \frac{1}{2}},\\ {\mathrm{W}}_{\kappa ,{\displaystyle \frac{1}{2}}}\left(x\right)={\displaystyle \frac{1}{\Gamma (1-\kappa )}}+O(xlnx),& \\ {\mathrm{W}}_{\kappa ,\mu }\left(x\right)={\displaystyle \frac{\Gamma \left(2\mu \right)}{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}}{x}^{{\displaystyle \frac{1}{2}}-\mu }+{\displaystyle \frac{\Gamma (-2\mu )}{\Gamma \left(\frac{1}{2}-\mu -\kappa \right)}}{x}^{{\displaystyle \frac{1}{2}}+\mu }+O\left(x^{{\displaystyle \frac{3}{2}}-\mu }\right),& 0\le \mu <{\displaystyle \frac{1}{2}},\mu \ne 0,\\ {\mathrm{W}}_{\kappa ,0}\left(x\right)=-{\displaystyle \frac{\sqrt{x}}{\Gamma \left(\frac{1}{2}-\kappa \right)}}\left(lnx+\psi \left({\displaystyle \frac{1}{2}}-\kappa \right)+2\gamma \right)+O\left(x^{3/2}lnx\right),\end{array}$$

(62)

Where, $\psi \left(z\right)$ is given by

$$\begin{array}{cc}\hfill & \psi \left(x\right)={\Gamma }^{\prime }\left(x\right)/\Gamma \left(x\right),\hfill \\ \hfill & x\ne 0,-1,-2,\dots \hfill \end{array}$$

(63)

Solution of the Initial value problem is shown in in Figure 3 for various $f\left(x\right)$.

5.2. Boundary Value Problems

In the previous section, we derived the solution to the inhomogeneous Whittaker equation under a set of prescribed initial conditions. This formulation allowed us to explore the time-evolution and dynamic behaviour of the system starting from known initial states. The focus was primarily on the temporal aspect of the solution, assuming the spatial domain and boundary influences were either trivial or not the primary concern.

In the present section, we shift our attention to a more general and practically significant case: the solution of the inhomogeneous Whittaker equation subject to boundary conditions. This framework is especially pertinent to physical systems in which spatial boundary behaviour significantly influences the overall dynamics, such as quantum mechanical potential problems, heat conduction processes, fluid flow through porous media, and wave propagation in confined or bounded domains. To this end, we consider the inhomogeneous Whittaker equation 8 subject to the following boundary conditions:

$$W\left(a\right)=\alpha ;W\left(b\right)=\beta ,$$

(64)

where $\alpha$ and $\beta$ are known constants. Applying the boundary conditions we get:

$$\begin{array}{cc}\hfill & C_1{\mathrm{M}}_{\kappa ,\mu }\left(a\right)+C_2{\mathrm{W}}_{\kappa ,\mu }\left(a\right)+{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }\left(a\right)=\alpha ,\hfill \\ \hfill & C_1{\mathrm{M}}_{\kappa ,\mu }\left(b\right)+C_2{\mathrm{W}}_{\kappa ,\mu }\left(b\right)+{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }\left(b\right)=\beta .\hfill \end{array}$$

(65)

Equations 65 can be rewritten in matrix form as:

$$\left[\begin{array}{cc}{\mathrm{M}}_{k,\mu }\left(a\right)\hfill & {\mathrm{W}}_{k,\mu }\left(a\right)\hfill \\ {\mathrm{M}}_{k,\mu }\left(b\right)\hfill & {\mathrm{W}}_{k,\mu }\left(b\right)\hfill \end{array}\right]\left[\begin{array}{c}{C}_1\hfill \\ C_2\hfill \end{array}\right]=\left[\begin{array}{c}\alpha -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }\left(a\right)\hfill \\ \beta -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{}^f\mathrm{Zi}_{\kappa ,\mu }\left(b\right)\hfill \end{array}\right].$$

(66)

Upon solving the system we get:

$$\begin{array}{cc}\hfill C_1={\displaystyle \frac{1}{{\mathrm{M}}_{\kappa ,\mu }\left(a\right){\mathrm{W}}_{\kappa ,\mu }\left(b\right)-{\mathrm{M}}_{\kappa ,\mu }\left(b\right){\mathrm{W}}_{\kappa ,\mu }\left(a\right)}}[& {\mathrm{W}}_{\kappa ,\mu }\left(b\right)\left(\alpha -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{\mathrm{Zi}}_{\kappa ,\mu }\left(a\right)\right)\hfill \\ \hfill -& {\mathrm{W}}_{\kappa ,\mu }\left(a\right)\left(\beta -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{\mathrm{Zi}}_{\kappa ,\mu }\left(b\right)\right)],\hfill \end{array}$$

(67)

$$\begin{array}{cc}\hfill C_2={\displaystyle \frac{1}{{\mathrm{M}}_{\kappa ,\mu }\left(a\right){\mathrm{W}}_{\kappa ,\mu }\left(b\right)-{\mathrm{M}}_{\kappa ,\mu }\left(b\right){\mathrm{W}}_{\kappa ,\mu }\left(a\right)}}[& {\mathrm{M}}_{\kappa ,\mu }\left(a\right)\left(\beta -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{\mathrm{Zi}}_{\kappa ,\mu }\left(b\right)\right)\hfill \\ \hfill -& {\mathrm{M}}_{\kappa ,\mu }\left(b\right)\left(\alpha -{\displaystyle \frac{\Gamma \left(\frac{1}{2}+\mu -\kappa \right)}{\Gamma (1+2\mu )}}{\mathrm{Zi}}_{\kappa ,\mu }\left(a\right)\right)].\hfill \end{array}$$

(68)

Figure 4 shows the solution of the inhomogeneous Whittaker equation for $\kappa =-{\displaystyle \frac{1}{4}}$ and $\mu ={\displaystyle \frac{1}{4}}$ with boundary conditions.

6. Conclusions

The primary objective of this work has been to develop a method for solving the inhomogeneous Whittaker differential equation. To achieve this, we introduced a novel class of special integral functions, denoted by ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(a\right)$ , which facilitate the construction of particular solutions to the inhomogeneous form of the Whittaker equation. A detailed analytical study of the function ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(a\right)$ was carried out, including the investigation of its fundamental properties and the derivation of expressions for its derivatives. Computational formulas were established for certain representative cases, especially when the forcing function $f\left(x\right)$ and the associated parameters $\kappa$ and $\mu$ take specific, analytically tractable forms.

To complement the theoretical analysis, we provided graphical representations of the function ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(a\right)$ , which reveal its behaviour under various parameter regimes. These visualizations not only illustrate the properties of the ${}^f\mathrm{Zi}_{\kappa ,\mu }\left(a\right)$ function but also highlight the structure of the corresponding solutions to initial and boundary value problems associated with the Whittaker inhomogeneous differential equation.

The results presented in this work contribute to the broader understanding of special function theory in the context of linear differential equations with inhomogeneous terms. The approach developed here opens up further avenues for the study of related differential systems and may serve as a foundation for applications in mathematical physics and engineering where such equations naturally arise.