Abstract and Applied Analysis

Volume 2013 (2013), Article ID 304739, 17 pages

http://dx.doi.org/10.1155/2013/304739

## Numerical Solution of Fuzzy Fractional Pharmacokinetics Model Arising from Drug Assimilation into the Bloodstream

^{1}Institute for Mathematical Research, Universiti Putra Malaysia (UPM), 43400 Serdang, Selangor, Malaysia^{2}Mathematics Department, Science Faculty, Universiti Putra Malaysia (UPM), 43400 Serdang, Selangor, Malaysia^{3}Institute of Microengineering and Nanoelectronics (IMEN), Universiti Kebangsaan Malaysia (UKM), 43600 Bangi, Selangor, Malaysia^{4}Department of Computer Engineering, Mashhad Branch, Islamic Azad University, Mashhad, Iran

Received 22 August 2013; Revised 8 October 2013; Accepted 9 October 2013

Academic Editor: Ali H. Bhrawy

Copyright © 2013 Ali Ahmadian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We propose a Jacobi tau method for solving a fuzzy fractional pharmacokinetics. This problem can model the concentration of the drug in the blood as time increases. The proposed approach is based on the Jacobi tau (JT) method. To illustrate the reliability of the method, some special cases of the equations are solved as test examples. The method reduces the solution of the problem to the solution of a system of algebraic equations. Error analysis included the fractional derivative error estimation, and the upper bound of the absolute errors is introduced for this method.

#### 1. Introduction

Pharmacokinetics is defined as the study of the time course of drug absorption, distribution, metabolism, and excretion [1]. Pharmacokinetics determines bioavailability, volume of distribution, and clearance. Bioavailability is the fraction of a drug absorbed into the systemic circulation, while volume of distribution and clearance are measure of apparent space in the body available to contain the drug and measure of the body’s ability to eliminate the drug, respectively. In pharmacokinetics, as a basic hypothesis, always a relationship exists between the pharmacologic or toxic response to a drug and the concentration of the drug in the blood (plasma). However, for some drugs there is no straightforward relationship between concentration in plasma and pharmacologic effects [2]. Generally, the concentration of drug in the systemic circulation is related to the concentration of drug at its sites of action. Absorption, drug distribution (target tissues), drug concentration in circulation system, rate of elimination/plasma concentration, and elimination (metabolism, excretion) are some parameters that directly modify the pharmacokinetics parameters.

Mathematics is widely used for the quantitative description of drug absorption, distribution, metabolism, and excretion (ADME). Some parameters in pharmacokinetics are defined by mathematic equations which can be obtained directly by measurement or through calculation using experimental data based on developed mathematical equations. For example, the first and simplest model used for drug absorption from the gastrointestinal (GI) tract assumes a “pseudosteady state” and utilizes the physicochemical properties of the compound in conjunction with the “pH-partition” hypothesis to predict the fraction of dose absorbed [3–5]. Tissue distribution is another important determinant of the pharmacokinetics profile of a drug. Hence, in drug development, the prediction of tissue distribution would help predict the in vivo pharmacokinetics of a compound prior to any experiments in animals or man. Two methods are available to predict tissue distribution. These methods predict either tissue: plasma ratios or the volume of distribution at steady state (V). From molecular descriptors, properties such as lipophilicity/hydrophobicity are estimated [6]. In some approaches, plasma protein binding is also taken into consideration [7]. Physiological information on tissue composition, the blood composition, and blood flow to the tissue is utilized to develop a partitioning model [8–10]. Möhler et al. [11] proposed a model that describes cell growth and detailed investigations on the metabolism and the kinetics of the influenza infection cycle which allow for the optimization of influenza virus vaccine production. In addition to these simple models, more complex mathematical manipulations called mathematical models have also been used to describe pharmacokinetics [12]. There are many classifications of pharmacokinetics models. In a specific classification for pharmacokinetics modeling, two general approaches can be considered: compartment-based modeling [13] and noncompartment-based modeling [14]. Both types of modeling take advantage of the quantitative structure-pharmacokinetics relationships that are described by empirical mathematical algorithms. They can be used to estimate the activity of a compound based on its chemical structure in a numeric format [15].

In this study, our focus is to find the approximate solution of the fuzzy fractional model of the compartment model for the flow of antihistamine in the blood. To this end, we use the operational matrix of the Caputo fractional derivative of the Jacobi tau approximation based on the Jacobi polynomials to derive the fuzzy fractional approximate solution under Hukuhara differentiability (H-differentiability). To the best of our knowledge, there are no results on the JT method for solving pharmacokinetics equations arising in mathematical physics. This partially motivated our interest in such model and solution method.

##### 1.1. Compartment Model

The compartment model framework is an extremely natural and valuable means with which to formulate models for processes which have inputs and outputs over time [16]. Compartment model itself can be classified in one-compartment model and two-compartment model. In the former one, the body is depicted as a kinetically homogeneous unit, while in the latter one the body resolves into a central compartment and a peripheral compartment [17]. According to the one-compartment model, the drug instantaneously distributes throughout the body, and at the same time it equilibrates between tissues. Thus, the drug concentration time profile shows a monophasic response. On the other hand, in two-compartment model the drug does not achieve instantaneous equilibration between the two compartments. Although, in this model, two compartments have no anatomical or physiological meaning, it is usually surmised that the central compartment is composed of tissues that are highly perfused such as heart, lungs, kidneys, liver, and brain while the peripheral compartment comprises less well-perfused tissues such as muscle, fat, and skin [17].

##### 1.2. Methods Have Been Proposed to Solve Fractional Pharmacokinetics Models

While an increasing number of fractional order integrals and differential equations applications have been reported in the physics [18, 19], signal processing [20], engineering, and bioengineering literatures [21, 22], little attention has been paid to this class of models in the pharmacokinetics-pharmacodynamic (PKPD) literature. The reasons for the lack of application to pharmacodynamics are mainly two: first and foremost, PKPD models incorporating fractional calculus have not been proposed; second is computational: while the analytical solution of fractional differential equations is available in special cases, it turns out that even the simplest PKPD models that can be constructed using fractional calculus do not allow analytical solutions [23, 24].

In this paper, we propose new category of PKPD models incorporating fuzzy fractional calculus and investigate their behavior using purposely written algorithm. The main purpose of this paper is to attract the attention of the field into these possibly interesting families of PKPD models. In particular, we will focus on compartment models, direct concentration-response relationships describing models relating drug concentration (in plasma or biophase) to a pharmacodynamic effect.

The concept of fractional or noninteger order derivation and integration can be traced back to the genesis of integer order calculus itself [25–27]. Due to its tremendous scope and applications in several disciplines, a considerable attention has been given to exact and numerical solutions of fractional differential equations which is an extremely difficult task. Moreover, the solution techniques and their reliability are rather more important aspects. Several methods have been proposed by many researchers to solve the fractional order differential equations such as Taylor series [28], variational iteration method [29–31], adomian decomposition method [32, 33], fractional differential transform method (FDTM) [34, 35], and homotopy analysis method [36, 37].

As it is known, the spectral method is one of the flexible methods of discretization for most types of differential equations [38–40]. Historically, spectral method has been relegated to fractional calculus, but in few years, it has been successfully applied for the fractional equation models based on the different types of orthogonal polynomials such as Block pulse functions [41, 42], Legendre polynomials [43–46], Chebyshev polynomials [47–49], Laguerre polynomials [50–53], and Bernstein polynomials [54–56]. Doha et al. [57] introduced the shifted Jacobi operational matrix of fractional derivative which is based on Jacobi tau method for solving numerically linear multiterm fractional differential equations with initial or boundary conditions. Subsequently, Kazem [58] generalized Jacobi integral operational matrix to fractional calculus. Afterwards, Doha et al. [59] proposed a direct solution techniques for solving the linear multiorder fractional initial value problems with constant and variable coefficients using shifted Jacobi tau method and quadrature shifted Jacobi tau method, respectively. In this paper, we intend to extend the application of the Jacobi polynomials to solve fuzzy fractional pharmacokinetics model of order.

On the other hand, a considerable attention has been made to the fractional differential equations in the sense of fuzzy setting theory. Agarwal et al. [60] is among the pioneers who presented the concept of solutions for fractional differential equations with uncertainty. Thereafter, the existence and uniqueness of the solution of the fuzzy fractional differential equations (FFDEs) under the Riemann-Liouville and Caputo’s fuzzy fractional differentiability were investigated in the literature [61–64]. However, the application of the numerical methods for solving FFDEs is unknown and traceless in the literature of the fuzzy fractional calculus, and only a few number of researches have been reported for the approximate solution of FFDEs [64–69].

In the present paper, we intend to introduce new families of PKPD models based on the application of fractional calculus to PKPD models. The aim of the paper is not to claim the superiority of fractional dynamics models with respect to standard ones, but it is simply to define the new families and provide some insights into their qualitative behavior. The main purpose is to apply the fuzzy logic in differential equations of fractional order which has been used as an effective tool for considering uncertainty in modeling the processes; FFDEs can also offer a more comprehensive account of the process or phenomenon, specifically for analyzing the behavior of the PKPD models. Furthermore, we suggest possible applications and stimulate further research, which might, or might not, demonstrate the applicability and importance of spectral methods by using of the orthogonal polynomials for finding the approximate solutions of the PKPD models based on the fuzzy fractional calculus.

The structure of this paper is as follows. In the next section, we briefly recall the mathematical foundations of fractional calculus, required definitions of fuzzy setting theory, and summarize the properties of Jacobi polynomials. In Section 3, we then provide illustrations of the governing fraction equation. In Section 4, the proposed method is explained for numerical solution of the derived FFDE. Section 5 is devoted to the numerical solution of the problem with different dose of the drug, and the error analysis is provided to demonstrate the applicability and validity of the method. A final discussion comments on the results are presented in Section 6.

#### 2. Preliminaries and Notations

In this section, we are going to state the definition and preliminaries of fuzzy mathematics [63, 70, 71], fractional calculus [25, 27], and some properties of shifted Jacobi polynomials.

##### 2.1. The Fuzzy Settings Definitions

*Definition 1. *Letbe a fuzzy set in.is called a fuzzy number if:(i)is normal: there existssuch that;(ii)is convex: for alland, it holds that
(iii)is upper semicontinuous: for any, it holds that
(iv)is a compact subset of.

In this paper, the set of all fuzzy numbers is denoted by.

*Definition 2. *Letand. The-cut ofis the crisp setthat contains all elements with membership degree ingreater than or equal to; that is
For a fuzzy number, its-cuts are closed intervals in, and we denote them by

According to Zadeh’s extension principle, the operation of addition onis defined as follows: and scalar multiplication of a fuzzy number is given by where.

*Definition 3 (see [70]). *The distancebetween two fuzzy numbersandis defined as
where
is the Hausdorff distance betweenand.

It is easy to see thatis a metric inand has the following properties (see [70, 72]):(i), (ii), (iii), (iv), (v)is a complete metric space.

*Definition 4 (see [73]). *Letandbe the two fuzzy-number-valued functions on the interval,that is,. The uniform distance between fuzzy-number-valued functions is defined by

*Remark 5 (see [73]). *Letbe fuzzy continuous. Then from property (iv) of Hausdorff distance, we can define

*Definition 6 (see [74]). *Let. If there existssuch that, thenis called the H-difference ofand, and it is denoted by.

In this paper, the sign “” always stands for H-difference, and note that. Also throughout the paper, it is assumed that the Hukuhara difference and generalized Hukuhara differentiability exist.

Theorem 7 (see [75]). *Letbe a function and denote, for each. Then,*(1)*ifis (1)-differentiable, thenandare differentiable functions and
*(2)*ifis (2)-differentiable, thenandare differentiable functions and
*

*Definition 8 (see [76]). *Consider thelinear system of the following equations:
The matrix form of the above equations is
where the coefficient matrix,is a crispmatrix and. This system is called a fuzzy linear system (FLS).

*Definition 9 (see [76]). *A fuzzy number vectorgiven byis called a solution of the fuzzy linear system (2) if
If for a particular,, we simply get
To solve fuzzy linear systems, see [77].

In this part, we firstly give some basic definitions and some properties of fractional calculus [27]. Afterwards, the extension of the fractional differentiability in the sense of fuzzy concept is provided, and some relevant properties which are used in the rest of the paper are given [61, 63, 64].

Letbe the smallest integer that exceeds, then Caputo’s fractional derivative operator of order is defined as where For the Caputo derivative, we have Caputo’s fractional differentiation is a linear operation, namely, whereandare constants.

Theorem 10 (generalized Taylor formula [78]). *Suppose thatforwhere; then
**
with
**
where
*

Here,shows fractional derivative operator in the Caputo sense, and is called fractional Taylor expansion with degreeof, andis reminder term in fractional Taylor expansion of.

Now, we present some fuzzy fractional notations which are used later in the paper.(i)is the set of all fuzzy-valued measurable functionsonwhere.(ii)is a space of fuzzy-valued functions which are continuous on.(iii)indicates the set of all fuzzy-valued functions which are continuous up to order.(iv)denotes the set of all fuzzy-valued functions which are absolutely continuous.

Note that one can easily find these notations in the crisp context in [25, 27] and references therein.

*Definition 11 (see [63]). *Let. The Riemann-Liouville integral of fuzzy-valued functionis defined as

*Definition 12 (see [63]). *Let . Then, is said to be Caputo’s H-differentiable atwhen

*Definition 13 (see [63]). *Letandand . We say that is fuzzy Caputo fractional differentiable of order at, if there exists an elementsuch that for all,,
or
or
or

For the sake of simplicity, we say that the fuzzy-valued functionis-differentiable if it is differentiable as in Definition 13 case (i), andis-differentiable if it is differentiable as in Definition 13 case (ii), and so on for the other cases.

Theorem 14 (see [63]). *Letand; then the fuzzy Caputo fractional derivative exists almost everywhere on, and for allone has
**
whenis (1)-differentiable, and
**
when f is (2)-differentiable.*

Theorem 15 (fuzzy generalized Taylor’s formula [72]). *Let, and suppose that forwhere,and. Then, one has
**
where .*

##### 2.2. Jacobi Polynomials

The well-known Jacobi polynomials associated with the parameters () (see, e.g., Luke [79] and Szegö [80]), are a sequence of polynomials, each, respectively, of degree. For using these polynomials on, we present the shifted Jacobi polynomials by implementing the change of variable. Let the shifted Jacobi polynomials be denoted by, satisfying the orthogonality relation whereand The shifted Jacobi polynomialof degreehas the form where Also we can state the shifted Jacobi polynomial by the following recurrence relation: whereand.

We notice that a function, square integrable in, can be expanded in terms of shifted Jacobi polynomials as where the coefficientsare For a given particular problem, only the first-terms shifted Jacobi polynomials are used. Therefore, we have

Theorem 16 (see [59]). *The Caputo fractional derivative of orderof the shifted Jacobi polynomials of degreeis obtained from
**
where
*

#### 3. Pharmacokinetics Model Equation

##### 3.1. Drug Assimilation into the Blood

The drug dissolves in the gastrointestinal tract (GI), and each ingredient is diffused into the bloodstream. They are carried to the locations in which they act and are removed from the blood by the kidneys and liver. Generally, the problem of drug assimilation into the body can be considered as a two-compartment model: GI-tract and the bloodstream [16]. Different compartments and the input/output of the model are depicted in Figure 1.

For each compartment by applying the balance law we can obtain
In this study, we consider a case of a *single cold pill*. Also, there is no ingestion of the drug except that which occurs initially.

##### 3.2. Case of a Single Cold Pill

Let us considerto be the amount of drug in the GI-tract at timeand, the amount in the bloodstream at time. In the GI-tract we suppose that the pill is to be swallowed, and so after this event (over subsequent time) we have nothing more entering the GI-tract. The pill dissolves and diffuses into the bloodstream from the GI-tract. So, the GI-tract is only an output term. Assuming that the output rate is proportional to the amount of drug in the bloodstream, which is in fact proportional to the GI-tract drug concentration, then [16], whereis the amount of a drug in the pill andis a positive coefficient. We suppose that, which means that the initial amount of the drug in the bloodstream is zero. As the drug diffuses from the GI-tract, the level increases, and as the kidneys and liver remove it, the level of drug gradually decreases. Thus, withanother positive constant. Decongestant and an antihistamine are the component of the cold pill, and the coefficient of proportionality (and) is different for the different component drugs in the pill.

In order to obtain the growth and decay of antihistamine levels in the GI tract and bloodstream, rearrange the first rate equation, multiply by the integrating factor, integrate, and then use the initial data to obtain Insert this formula forinto the second rate equation, which becomes In this paper, we are concerned with fractional time derivatives regarding (17) and (18). Having defined , we can now turn to fractional differential equations and systems of fractional differential equations, which will be used to specify PKPD models and will need to be solved over an interval, in accordance with appropriate initial conditions.

A typical feature of differential equations (both classical and fractional) is the need to specify additional conditions in order to produce a unique solution. For the case of Caputo fractional differential equations, these additional conditions are just the initial conditions which are similarly required by classical ODEs. In contrast, for Riemann-Liouville fractional differential equations, these additional conditions constitute certain fractional derivatives (and/or integrals) of the unknown solution at the initial point[81], which are functions of. These initial conditions are not physical; furthermore, it is not clear how such quantities are to be measured from experiment, say, so that they can be appropriately assigned in an analysis [82]. If for no other reason, the need to solve fractional differential equations is justification enough for choosing Caputo’s definition for fractional differentiation over the more commonly used (at least in mathematical analysis) definition of Liouville and Riemann, and this is the operator that we choose to use in the following.

We start by representing drug concentration in the effect compartment by the (Caputo) fractional differential equation:

In the standard direct action model, the effect at time,, is expressed by an arbitrary (memory-less) function of drug concentration in the effect site at time,; however to generate a wider class of relationships, we assume that the effect at timeis related to the fuzzy Caputo fractional derivative of. So, we have in whichis a continuous fuzzy-valued function anddenotes the fuzzy Caputo fractional derivative of order.

*Remark 17. *In this paper, the drug concentration is modeled by an oscillation-relaxation fuzzy fractional differential equation. So the right hand side of (49) is determined based on the model parameters. It should mention here that the proposed method can be extended easily for solving other types of FFDEs with more complicated right hand side functions. For more details, one can refer to (33) in [65] and [57–59].

We have assumed that, an assumption that is justified by the pharmaceutical data. For the “average” person, a pharmaceutical company estimates that the values of the rate constants for the antihistamine in the cold pills areand. It can be observed from (46) that level of antihistamine in the bloodstream increases as the time increases and saturate in a maximum value of antihistamine; however, (47) can conclude that as the time increases the amount of antihistamine in GI-tract declines and reaches a minimum value. We now consider two different cases of patients who are not average. We often defineas clearance coefficient of medication from the blood. The value ofoften varies from old and sick patients than young and healthy cases. This means that the level of medication in the blood may become and then remain excessively high with a standard dosage for the cases thatis much lower, which is normally observed in old and sick cases.

We investigate sensitivity of the medication over a 24-hour period by keeping the value offixed at some value such as 1.386 but setting. Different values ofcorrespond to people of different ages and states of health. For all cases as it is expected that, for the first few hours, the amount of antihistamine in the bloodstream increases and then declines gradually. However, for the cases with lower value of clearance coefficient (old and sick), the maximum value of the medication in bloodstream was much higher than the cases with high, and it did not decline for remaining times. This means that the level of medication in the bloodstream stays at high level for a long time, and it could not be absorbed from the blood. The same trend is observed for the case thatis constant () and the value ofvaries (, and. In this case, theis a constant for the GI tract which is analogous to thefor the bloodstream.

#### 4. Description of the Method

In this study, by developing the Jacobi polynomial approximation [57, 58, 65] with the help of the matrix operations, the tau method, and the fuzzy Caputo fractional derivative, we obtain an approximate solution of the problem (49) for different values ofand. As it is known, the existence and uniqueness of the FFDEs have been introduced in [60, 61, 63]. Here, we assume that the shifted Jacobi polynomials are defined on, so we have [65] where the fuzzy coefficientsare gained by and,, andis as the same as the shifted Jacobi polynomials presented in Section 2.2, andmeans addition with respect toin.

*Remark 18. *In the remainder of paper, for more simplicity, we considerinstead ofandinstead of.

*Remark 19. *Practically, only the first-terms shifted Jacobi polynomials are taken into consideration. So, we have
where the fuzzy shifted Jacobi coefficient vectorand shifted Jacobi polynomials vectorare presented by
We can state the fuzzy approximate function (52) in the parametric form as follows.

*Definition 20 (see [65]). *Let; the approximation of fuzzy-valued functionin the parametric form is

We call thatforms a complete orthogonal system in. Hence, we set The shifted Jacobi tau method to (49) is to obtainsuch that whereanddenotes the fuzzy inner product inandis the fuzzy initial value of the problem.

Let us define the following notations: Hence, the variational formulation of (56) according to Relation (14) in [65], by means of a typical tau method like in the crisp context [29] and (52), is equivalent to Denoting we investigate that (58) is equivalent to the matrix system The elements of the matrices mentioned above are determined in the next theorem.

Theorem 21 (see [59]). *Let us denote and ,; then the nonzero elements ofandare given as
*

##### 4.1. Error Analysis

In this section, error analysis of the method will be presented for the FFDEs. Firstly, an upper bound of the absolute errors will be given for the technique by using generalized Taylor formula. Secondly, an error bound will be introduced for the approximation of the fuzzy fractional Caputo’s derivative using Jacobi polynomials.

Lemma 22. *Letandbe the Jacobi approximate solution (40) and the exact solution of (48), respectively. If, then
**whereandare defined according to Theorem 10.*

*Proof. *Sincecan be expanded to the fractional Taylor series
and its reminder term is
Therefore,
From (65) and triangle inequality, we obtain
Therefore, an upper bound of the absolute errors is obtained for the method in the crisp cases.

Now, we provide an upper bound for the absolute errors of the fuzzy approximate solution by using the proposed method.

Theorem 23. *Letbe the fuzzy exact solution of (49), andis the best fuzzy Jacobi approximate function (52), and suppose thatforwhereand. Then, one has
**
where regarding Theorem 15, one has
*

*Proof. *It is straightforward from Lemma 22 and the fuzzy function definition.

Also in the following theorem, according to the Relation (14) in [65], an upper bound for the absolute errors of the approximate function of fuzzy fractional Caputo’s derivative is provided.

Theorem 24 (see [65]). *Assume that the error function of fuzzy Caputo fractional derivative operator by using the shifted Jacobi polynomials,, is continuously fuzzy differentiable for,. Additionally,and; then the error bound is given by
*

#### 5. Numerical Results

In this section, the fuzzy fractional PKPD model (49) is solved for different values ofandby using the JT method presented in Section 4. We have performed all numerical computations with a computer program written in MATLAB. Also, absolute errors between fuzzy approximate solutionand the corresponding exact solutions , that is, , are considered.

Now, we recall the FFDE (49) as follows: in whichis a continuous fuzzy-valued function anddenotes the fuzzy Caputo fractional derivative of order.

Two cases are considered for solving by the proposed technique. In the first case, we assume thatis unchangeable and different values ofare substituted in (49) to get the fuzzy approximate solution. Conversely, in the second case we try to solve the problem such thathas an invariable value andis varied. The problem is analyzed for both cases in details.

*Case I. *Let us consider from Section 3 that, but letvary (e.g.,, and). So (49), with the assumption, namely,and Dose, is as follows:

By using Theorem 7 for the above equation under differentiability, we have the following systems:

Solving (72) leads to determining the exact solution of (71) as follows:

We seek the fuzzy approximate solutions by applying the technique described in Section 4 withand as Here, we have where. Regarding (54), we have

By solving the fuzzy linear system (60) in the parametric form, one can get fuzzy unknown coefficients in (76) as where we assume thatin (77).

The comparison between absolute errors of differentobtained by our method is shown in Table 1. In Figure 2 logarithmic plot of absolute error with different values ofandwas obtained, and different number of Jacobi functions was experienced for solving Case I by using the proposed method in Figure 3. Finally, the approximate fuzzy solutions are shown for different values ofwith the fractional order in Figure 4.

The absolute error value for-cut varied from 0 to 1 for different value ofwhich is calculated in Table 1. As it can be observed at a constant-cut by increasing the value ofwhich is the ratio of variation of the drug in the bloodstream to the amount of drug in the GI-tract the value of absolute error increases. This is analogous to the previous reports for various values of the[12]. It should be noted here that the value ofis directly proportional to the variation of exact solutions and corresponding fuzzy approximate solution. As a result, increase ofwith increasingis also expected.

*Case II. *In this case, we assume thatand is not variable, but takes different values (i.e.,, and 0.1386) and Dose. So (70) alters with, as follows:
with the exact solution as:

By applying the JT method in Section 4, we can get the fuzzy unknown coefficientsby using (58) which is in the parametric form as follows: whereis given as

withand. By substituting the above matric in (80) we can reach the fuzzy linear algebraic system (60) that can solve easily to determine the fuzzy coefficients as shown in Table 2.

We compared JT method results for different values ofand the outcomes are tabulated in Table 3. The absolute errors of the proposed method for this case are exhibited in Figure 5 with four choices ofand . Clearly, the best approximation is achieved when the values ofand. Additionally, we compared the approximate solutions obtained by the present method at, andwhich is shown in Figure 6. From Figure 6, one can conclude that with increasing the number of Jacobi polynomials, the absolute errors are decreasing dramatically. Finally, Figure 7 shows the numerical results for different values ofatIt can be seen that for all the values ofthe approximate solution is the fuzzy number.

In Table 3, the value of absolute error for-cut varied from 0 to 1 for different value ofis calculated. At a constant-cut value, a descending trend is observed with increasing the value of. This trend is opposed to the trend observed in Table 1 for different value of.is clearance constant and it is a coefficient forwhich is the amount of drug in the bloodstream in (44). We already explained that by decreasing the value ofwhich happens in the case of old and sick patient, the drug absorption from the bloodstream by kidney and liver also decreases and the drug stays in the bloodstream for longer time. Since the constant value ofappears with negative sign in the main differential equation as we expected, the value ofwhich is directly proportional to the exact solutions is decreased with increasing the value of.

#### 6. Conclusion

Fuzzy theory provides a suitable way to objectively account for parameter uncertainty in models. Fuzzy logic approaches appear promising in preclinical applications and might be useful in drug discovery and design. Considerable progress has been made in the last few years in the development of computational approaches for prediction of drug absorption, distribution, metabolism, and excretion. Whilst several approaches have been developed in pharmacokinetics, most of these approaches have not yet been adequately used in the complex process such as prediction of metabolism, and they require further improvement.

In summary in this research, a tau method based on the Jacobi operational matrix was utilized to numerically solve the PKPD equation, arising from drug assimilation into the bloodstream. The comparison of the results shows that the present method is a powerful mathematical tool for finding the numerical solutions of a generalized linear fuzzy fractional PKPD equation.

Although we concentrated on applying our algorithm to solve fuzzy fractional PKPD equation, we show that such algorithm can be applied to solve other types of fractional equations models in science and engineering fields. Our algorithm for the fuzzy fractional PKPD equations is efficient and numerically stable. Numerical results are presented which exhibit the high accuracy of the proposed algorithm.

By using the obtained operational matrix, the following objectives were achieved.(i)The given fuzzy fractional differential equation was converted into a fuzzy algebraic system of equations, simplifying the solution procedure.(ii)The method is computer oriented; thus, solving fuzzy fractional differential equations of different fractional ordersbecomes a matter of changingonly.

#### References

- B. G. Katzung, S. B. Masters, and A. J. Trevor,
*Basic & Clinical Pharmacology*, McGraw-Hill, 2004. - H. Peng and B. Cheung, “A review on pharmacokinetic modeling and the effects of environmental stressors on pharmacokinetics for operational medicine: operational pharmacokinetics,” Tech. Rep. ADA509469, DTIC, 2009. View at Google Scholar
- A. Suzuki, W. I. Higuchi, and N. F. Ho, “Theoretical model studies of drug absorption and transport in the gastrointestinal tract. I,”
*Journal of Pharmaceutical Sciences*, vol. 59, no. 5, pp. 644–651, 1970. View at Publisher · View at Google Scholar · View at Scopus - H. Boxenbaum, “Absorption potential and its variants,”
*Pharmaceutical Research*, vol. 16, no. 12, p. 1893, 1999. View at Publisher · View at Google Scholar · View at Scopus - A. Boobis, U. Gundert-Remy, P. Kremers, P. Macheras, and O. Pelkonen, “In silico prediction of ADME and pharmacokinetics: report of an expert meeting organised by COST B15,”
*European Journal of Pharmaceutical Sciences*, vol. 17, no. 4-5, pp. 183–193, 2002. View at Publisher · View at Google Scholar · View at Scopus - P. Poulin, K. Schoenlein, and F. P. Theil, “Prediction of adipose tissue: plasma partition coefcients for structurally unrelated drugs,”
*Journal of Pharmaceutical Sciences*, vol. 90, pp. 436–447, 2001. View at Google Scholar - F. Lombardo, R. S. Obach, M. Y. Shalaeva, and F. Gao, “Prediction of volume of distribution values in humans for neutral and basic drugs using physicochemical measurements and plasma protein binding data,”
*Journal of Medicinal Chemistry*, vol. 45, no. 13, pp. 2867–2876, 2002. View at Publisher · View at Google Scholar · View at Scopus - P. Poulin and F. P. Theil, “A priori prediction of tissue: plasma partition coefcients of drugs to facilitate the use of physiologically-based pharmacokinetic models in drug discovery,”
*Journal of Pharmaceutical Sciences*, vol. 89, pp. 16–35, 1999. View at Google Scholar - P. Poulin and F.-P. Theil, “Prediction of pharmacokinetics prior to in vivo studies. 1. Mechanism-based prediction of volume of distribution,”
*Journal of Pharmaceutical Sciences*, vol. 91, no. 1, pp. 129–156, 2002. View at Publisher · View at Google Scholar · View at Scopus - P. Poulin and F.-P. Theil, “Prediction of pharmacokinetics prior to in vivo studies. II. Generic physiologically based pharmacokinetic models of drug disposition,”
*Journal of Pharmaceutical Sciences*, vol. 91, no. 5, pp. 1358–1370, 2002. View at Publisher · View at Google Scholar · View at Scopus - L. Möhler, D. Flockerzi, H. Sann, and U. Reichl, “Mathematical model of influenza a virus production in large-scale microcarrier culture,”
*Biotechnology and Bioengineering*, vol. 90, no. 1, pp. 46–58, 2005. View at Publisher · View at Google Scholar · View at Scopus - P. L. Bonate,
*Pharmacokinetic-Pharmacodynamic Modeling and Simulation*, Springer, 2011. - M. Holz and A. Fahr, “Compartment modeling,”
*Advanced Drug Delivery Reviews*, vol. 48, no. 2-3, pp. 249–264, 2001. View at Publisher · View at Google Scholar · View at Scopus - P. Veng-Pedersen, “Noncompartmentally-based pharmacokinetic modeling,”
*Advanced Drug Delivery Reviews*, vol. 48, no. 2-3, pp. 265–300, 2001. View at Publisher · View at Google Scholar · View at Scopus - D. E. Mager, “Quantitative structure-pharmacokinetic/pharmacodynamic relationships,”
*Advanced Drug Delivery Reviews*, vol. 58, no. 12-13, pp. 1326–1356, 2006. View at Publisher · View at Google Scholar · View at Scopus - B. Barnes and G. R. Fulford,
*Mathematical Modelling with Case Studies: A Differential Equation Approach Using Maple and MATLAB*, CRC Press, 2011. - S. Dhillon and K. Gill, “Basic pharmacokinetics,” in
*Clinical Pharmacokinetics*, pp. 1–44, Pharmaceutical Press, 2006. View at Google Scholar - L. Gaul, P. Klein, and S. Kemple, “Damping description involving fractional operators,”
*Mechanical Systems and Signal Processing*, vol. 5, no. 2, pp. 81–88, 1991. View at Google Scholar · View at Scopus - W. G. Glockle and T. F. Nonnenmacher, “A fractional calculus approach to self-similar protein dynamics,”
*Biophysical Journal*, vol. 68, no. 1, pp. 46–53, 1995. View at Google Scholar · View at Scopus - R. J. Marks II and M. W. Hall, “Differintegral interpolation from a bandlimited signal's samples,”
*IEEE Transactions on Acoustics, Speech, and Signal Processing*, vol. 29, no. 4, pp. 872–877, 1981. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - A. Atangana and D. Baleanu, “Nonlinear fractional Jaulent-Miodek and Whitham-Broer-Kaup equations within Sumudu transform,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 160681, 8 pages, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Mehdinejadiani, A. A. Naseri, H. Jafari, A. Ghanbarzadeh, and D. Baleanu, “A mathematical model for simulation of a water table profile between two parallel subsurface drains using fractional derivatives,”
*Computers & Mathematics with Applications*, vol. 66, no. 5, pp. 785–794, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - C. Csajka and D. Verotta, “Pharmacokinetic-pharmacodynamic modelling: history and perspectives,”
*Journal of Pharmacokinetics and Pharmacodynamics*, vol. 33, no. 3, pp. 227–279, 2006. View at Publisher · View at Google Scholar · View at Scopus - D. Verotta, “Fractional dynamics pharmacokinetics-pharmacodynamic models,”
*Journal of Pharmacokinetics and Pharmacodynamics*, vol. 37, no. 3, pp. 257–276, 2010. View at Publisher · View at Google Scholar · View at Scopus - D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo,
*Fractional Calculus: Models and Numerical Methods*, vol. 3 of*Series on Complexity, Nonlinearity and Chaos*, World Scientific Publishing, Hackensack, NJ, USA, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - D. Baleanu, Z. B. G. Güvenç, and J. A. Tenreiro Machado,
*New Trends in Nanotechnology and Fractional Calculus Applications*, Springer, New York, NY, USA, 2010. View at Publisher · View at Google Scholar · View at MathSciNet - I. Podlubny,
*Fractional Differential Equations*, vol. 198 of*Mathematics in Science and Engineering*, Academic Press, San Diego, Calif, USA, 1999. View at MathSciNet - M. Gülsu, Y. Öztürk, and A. Anapalı, “Numerical approach for solving fractional relaxation-oscillation equation,”
*Applied Mathematical Modelling*, vol. 37, no. 8, pp. 5927–5937, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - M. Dehghan and A. Saadatmandi, “A tau method for the one-dimensional parabolic inverse problem subject to temperature overspecification,”
*Computers and Mathematics with Applications*, vol. 52, no. 6-7, pp. 933–940, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - H. Jafari, H. Tajadodi, and D. Baleanu, “A modified variational iteration method for solving fractional Riccati differential equation by Adomian polynomials,”
*Fractional Calculus and Applied Analysis*, vol. 16, no. 1, pp. 109–122, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - G.-C. Wu and D. Baleanu, “Variational iteration method for fractional calculus—a universal approach by Laplace transform,”
*Advances in Difference Equations*, vol. 2013, article 18, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - Y. Hu, Y. Luo, and Z. Lu, “Analytical solution of the linear fractional differential equation by Adomian decomposition method,”
*Journal of Computational and Applied Mathematics*, vol. 215, no. 1, pp. 220–229, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - S. Saha Ray and R. K. Bera, “Analytical solution of the Bagley Torvik equation by Adomian decomposition method,”
*Applied Mathematics and Computation*, vol. 168, no. 1, pp. 398–410, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - A. Arikoglu and I. Ozkol, “Solution of fractional differential equations by using differential transform method,”
*Chaos, Solitons & Fractals*, vol. 34, no. 5, pp. 1473–1481, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - P. Darania and A. Ebadian, “A method for the numerical solution of the integro-differential equations,”
*Applied Mathematics and Computation*, vol. 188, no. 1, pp. 657–668, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - I. Hashim, O. Abdulaziz, and S. Momani, “Homotopy analysis method for fractional IVPs,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 14, no. 3, pp. 674–684, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - H. M. Jaradat, F. Awawdeh, and E. A. Rawashdeh, “An analytical scheme for multi-order fractional differential equations,”
*Tamsui Oxford Journal of Mathematical Sciences*, vol. 26, no. 3, pp. 305–320, 2010. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang,
*Spectral Methods in Fluid Dynamics*, Springer Series in Computational Physics, Springer, New York, NY, USA, 1988. View at MathSciNet - B. Fornberg,
*A Practical Guide to Pseudospectral Methods*, Cambridge University Press, Cambridge, Mass, USA, 1998. - E. L. Ortiz and H. Samara, “Numerical solution of differential eigenvalue problems with an operational approach to the Tau method,”
*Computing*, vol. 31, no. 2, pp. 95–103, 1983. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - Y. Li and N. Sun, “Numerical solution of fractional differential equations using the generalized block pulse operational matrix,”
*Computers and Mathematics with Applications*, vol. 62, no. 3, pp. 1046–1054, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - M. Yi, J. Huang, and J. Wei, “Block pulse operational matrix method for solving fractional partial differential equation,”
*Applied Mathematics and Computation*, vol. 221, pp. 121–131, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - A. Saadatmandi and M. Dehghan, “A new operational matrix for solving fractional-order differential equations,”
*Computers and Mathematics with Applications*, vol. 59, no. 3, pp. 1326–1336, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - S. Kazem, S. Abbasbandy, and S. Kumar, “Fractional-order Legendre functions for solving fractional-order differential equations,”
*Applied Mathematical Modelling*, vol. 37, no. 7, pp. 5498–5510, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - A. H. Bhrawy and M. M. Al-Shomrani, “A shifted Legendre spectral method for fractional-order multipoint boundary value problems,”
*Advances in Difference Equations*, vol. 2012, article 8, 2012. View at Publisher · View at Google Scholar - A. H. Bhrawy, A. S. Alofi, and S. S. Ezz-Eldien, “A quadrature tau method for fractional differential equations with variable coefficients,”
*Applied Mathematics Letters*, vol. 24, no. 12, pp. 2146–2152, 2011. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “Efficient Chebyshev spectral methods for solving multi-term fractional orders differential equations,”
*Applied Mathematical Modelling*, vol. 35, no. 12, pp. 5662–5672, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “A Chebyshev spectral method based on operational matrix for initial and boundary value problems of fractional order,”
*Computers and Mathematics with Applications*, vol. 62, no. 5, pp. 2364–2373, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - A. H. Bhrawy and A. S. Alofi, “The operational matrix of fractional integration for shifted Chebyshev polynomials,”
*Applied Mathematics Letters*, vol. 26, no. 1, pp. 25–31, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Baleanu, A. H. Bhrawy, and T. M. Taha, “Two efficient generalized Laguerre spectral algorithms for fractional initial value problems,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 546502, 10 pages, 2013. View at Publisher · View at Google Scholar - D. Baleanu, A. H. Bhrawy, and T. M. Taha, “A modified generalized Laguerre spectral method for fractional differential equations on the half line,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 413529, 12 pages, 2013. View at Publisher · View at Google Scholar - A. H. Bhrawy and T. M. Taha, “An operational matrix of fractional integration of the Laguerre polynomials and its application on a semi-infinite interval,”
*Mathematical Sciences*, vol. 6, article 41, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - A. H. Bhrawy, M. M. Alghamdi, and T. M. Taha, “A new modified generalized Laguerre operational matrix of fractional integration for solving fractional differential equations on the half line,”
*Advances in Difference Equations*, vol. 2012, article 179, 2012. View at Publisher · View at Google Scholar - D. Rostamy, M. Alipour, H. Jafari, and D. Baleanu, “Solving multi-term orders fractional differential equations by operational matrices of BPs with convergence analysis,”
*Romanian Reports in Physics*, vol. 65, pp. 334–349, 2013. View at Google Scholar - S. Yüzbaşı, “Numerical solutions of fractional Riccati type differential equations by means of the Bernstein polynomials,”
*Applied Mathematics and Computation*, vol. 219, no. 11, pp. 6328–6343, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Baleanu, M. Alipour, and H. Jafari, “The Bernstein operational matrices for solving the fractional quadratic Riccati differential equations with the Riemann-Liouville derivative,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 461970, 7 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “A new Jacobi operational matrix: an application for solving fractional differential equations,”
*Applied Mathematical Modelling*, vol. 36, no. 10, pp. 4931–4943, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Kazem, “An integral operational matrix based on Jacobi polynomials for solving fractional-order differential equations,”
*Applied Mathematical Modelling*, vol. 37, no. 3, pp. 1126–1136, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - E. H. Doha, A. H. Bhrawy, D. Baleanu, and S. S. Ezz-Eldien, “On shifted Jacobi spectral approximations for solving fractional differential equations,”
*Applied Mathematics and Computation*, vol. 219, no. 15, pp. 8042–8056, 2013. View at Publisher · View at Google Scholar - R. P. Agarwal, V. Lakshmikantham, and J. J. Nieto, “On the concept of solution for fractional differential equations with uncertainty,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 72, no. 6, pp. 2859–2862, 2010. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - T. Allahviranloo, S. Salahshour, and S. Abbasbandy, “Explicit solutions of fractional differential equations with uncertainty,”
*Soft Computing*, vol. 16, no. 2, pp. 297–302, 2012. View at Publisher · View at Google Scholar · View at Scopus - T. Allahviranloo, Z. Gouyandeh, and A. Armand, “Fuzzy fractional differential equations under generalized fuzzy Caputo derivative,”
*Journal of Intelligent and Fuzzy Systems*. In press. - S. Salahshour, T. Allahviranloo, S. Abbasbandy, and D. Baleanu, “Existence and uniqueness results for fractional differential equations with uncertainty,”
*Advances in Difference Equations*, vol. 2012, article 112, 2012. View at Google Scholar - S. Salahshour, T. Allahviranloo, and S. Abbasbandy, “Solving fuzzy fractional differential equations by fuzzy Laplace transforms,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 17, no. 3, pp. 1372–1381, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - A. Ahmadian, M. Suleiman, S. Salahshour, and D. Baleanu, “A Jacobi operational matrix for solving a fuzzy linear fractional differential equation,”
*Advances in Difference Equations*, vol. 2013, article 104, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - A. Ahmadian, M. Suleiman, and S. Salahshour, “An operational matrix based on Legendre polynomials for solving fuzzy fractional-order differential equations,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 505903, 29 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - M. R. Balooch Shahriyar, F. Ismail, S. Aghabeigi, A. Ahmadian, and S. Salahshour, “An eigenvalue-eigenvector method for solving a system of fractional differential equations with uncertainty,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 579761, 11 pages, 2013. View at Publisher · View at Google Scholar - F. Ghaemi, R. Yunus, A. Ahmadian, S. Salahshour, M. Suleiman, and Sh. Faridah Saleh, “Application of fuzzy fractional kinetic equations to modelling of the acid hydrolysis reaction,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 610314, 19 pages, 2013. View at Publisher · View at Google Scholar - M. Mazandarani and A. Vahidian Kamyad, “Modified fractional Euler method for solving fuzzy fractional initial value problem,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 18, no. 1, pp. 12–21, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Dubois and H. Prade, “Towards fuzzy differential calculus—part 3: differentiation,”
*Fuzzy Sets and Systems*, vol. 8, no. 3, pp. 225–233, 1982. View at Google Scholar · View at Scopus - R. Goetschel Jr. and W. Voxman, “Elementary fuzzy calculus,”
*Fuzzy Sets and Systems*, vol. 18, no. 1, pp. 31–43, 1986. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - G. A. Anastassiou,
*Fuzzy Mathematics: Approximation Theory*, vol. 251 of*Studies in Fuzziness and Soft Computing*, Springer, Berlin, Germany, 2010. View at Publisher · View at Google Scholar · View at MathSciNet - G. A. Anastassiou and S. G. Gal, “On a fuzzy trigonometric approximation theorem of Weierstrass-type,”
*Journal of Fuzzy Mathematics*, vol. 9, no. 3, pp. 701–708, 2001. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - O. Kaleva, “Fuzzy differential equations,”
*Fuzzy Sets and Systems*, vol. 24, no. 3, pp. 301–317, 1987. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - Y. Chalco-Cano and H. Román-Flores, “On new solutions of fuzzy differential equations,”
*Chaos, Solitons & Fractals*, vol. 38, no. 1, pp. 112–119, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - T. Allahviranloo and M. Afshar Kermani, “Solution of a fuzzy system of linear equation,”
*Applied Mathematics and Computation*, vol. 175, no. 1, pp. 519–531, 2006. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - S. Abbasbandy and R. Ezzati, “Newton's method for solving a system of fuzzy nonlinear equations,”
*Applied Mathematics and Computation*, vol. 175, no. 2, pp. 1189–1199, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - Z. M. Odibat and N. T. Shawagfeh, “Generalized Taylor's formula,”
*Applied Mathematics and Computation*, vol. 186, no. 1, pp. 286–293, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - Y. Luke,
*The Special Functions and Their Approximations*, vol. 2, Academic Press, New York, NY, USA, 1969. - G. Szegö,
*Orthogonal Polynomials*, vol. 23 of*Colloquium Publications*, American Mathematical Society, Providence, RI, USA, 1985. - A. A. Kilbas and J. J. Trujillo, “Differential equations of fractional order: methods, results and problems. I,”
*Applicable Analysis*, vol. 78, no. 1-2, pp. 153–192, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. S. Miller and B. Ross,
*An Introduction to the Fractional Calculus and Fractional Differential Equations*, A Wiley-Interscience Publication, John Wiley & Sons, New York, NY, USA, 1993. View at MathSciNet