Skip to content
BY 4.0 license Open Access Published by De Gruyter Open Access September 26, 2022

Generalized invexity and duality in multiobjective variational problems involving non-singular fractional derivative

  • Ved Prakash Dubey , Devendra Kumar EMAIL logo , Hashim M. Alshehri , Jagdev Singh and Dumitru Baleanu
From the journal Open Physics

Abstract

In this article, we extend the generalized invexity and duality results for multiobjective variational problems with fractional derivative pertaining to an exponential kernel by using the concept of weak minima. Multiobjective variational problems find their applications in economic planning, flight control design, industrial process control, control of space structures, control of production and inventory, advertising investment, impulsive control problems, mechanics, and several other engineering and scientific problems. The proposed work considers the newly derived Caputo–Fabrizio (CF) fractional derivative operator. It is actually a convolution of the exponential function and the first-order derivative. The significant characteristic of this fractional derivative operator is that it provides a non-singular exponential kernel, which describes the dynamics of a system in a better way. Moreover, the proposed work also presents various weak, strong, and converse duality theorems under the diverse generalized invexity conditions in view of the CF fractional derivative operator.

1 Introduction

The present scenario indicates that the fractional differential equations (FDEs) and fractional variational problems (FVPs) are being used to delineate the physical models and engineering processes in a better way. The clear reason is that the standard mathematical models of integer-order derivatives incorporating models of non-linear nature do not perform efficiently in many instances according to desired results. Recently, the field of fractional calculus has portrayed a significant part in various areas of knowledge such as chemistry [1], biology [2,3], mechanics [4,5], and finance [6]. The application area of fractional modelling and fractional operators encompasses anomalous diffusion [7], physics [8], heat conduction [9], geophysics [10], epidemiology [11], fractals and fractional derivative [12], computational fractional derivative equations [13], fractional predator–prey system [14], and porous media [15]. The models related to these fields utilize fractional derivative operators frequently.

There are various types of fractional derivative operators in the literature of fractional calculus founded by so many famous mathematicians. But the most popular definitions of them are Riemann–Liouville (RL) fractional derivative and Riesz fractional derivative described in the studies of Samko et al. [16] and Podlubny [17], Caputo fractional derivative in refs. [18,19], Weyl fractional derivative [20], Hadamard fractional derivative [16,17], Jumarie’s fractional derivative propounded in works of Jumarie [21,22], Atangana–Baleanu derivative proposed in ref. [23], and Liouville–Caputo derivative described in ref. [24]. The interesting fact is that these definitions of fractional derivative operators have their own significance and their uses vary according to the structure and behaviour of particular models along with initial conditions. A wide literature is available on different perceptions of fractional derivatives. But the most celebrated fractional calculi are the Caputo fractional derivative and the RL derivative. The Caputo fractional derivative handles initial value problems efficiently in comparison to the RL derivative. The newly introduced Caputo–Fabrizio (CF) fractional derivative operator propounded by Caputo and Fabrizio in ref. [25] is actually a convolution of an exponential function and the first-order derivative. In this definition, the derivative of a constant is equal to zero like the usual Liouville–Caputo definition but it also provides the non-singular kernel which was not a characteristic of the Liouville–Caputo fractional derivative. The main purpose of the CF definition was to introduce a new fractional derivative with an exponential kernel to describe even better the dynamics of systems with memory effect.

Recently, some authors presented a new analysis on fractional modelling of real-world problems and application of fractional order Lagrangian approach towards study of problems arising in physical sciences and engineering. Some recent works related to these fields are necessary to be cited here. Jajarmi et al. [26] suggested a general fractional formulation for immunogenic tumour dynamics. Baleanu et al. [27] presented a new study on the general fractional model of COVID-19 with isolation and quarantine effects. Erturk et al. [28] utilized a new fractional-order Lagrangian to describe the dynamics of a beam on nanowire. Jajarmi et al. [29] implemented a new fractional Lagrangian approach to study the case of capacitor microphone. Dubey et al. [30] solved the fractional model of Phytoplankton–Toxic Phytoplankton–Zooplankton system with convergence analysis. Moreover, a fractional model of atmospheric dynamics of carbon dioxide gas [31] and a fractional-order hepatitis E virus model [32] were also recently investigated with efficient computational methods.

Multiobjective variational problems proficiently handle the problems of science, engineering, logistics, and economics where optimal decisions have to be decided between two or more clashing objectives. To derive the optimality conditions it is necessary to study the behaviour of functions and their derivatives at that point. In the theory of mathematical optimization, the duality principle indicates two perspectives of optimization problems: the primal problem and the dual problem. If the primal is a minimization problem, then the dual is a maximization problem, and if the primal is a maximization problem, then the dual is a minimization problem. The concept of duality considers a problem with less number of variables and constraints and so it is much advantageous regarding computational procedure. Duality results play a major role in construction of numerical algorithms for solving some specific types of optimization problems. The duality theory is applied mainly in economics, management, physics, etc. On the other hand, calculus of variations significantly deals with the solution of several problems arising in theory of variations, optimization of orbits, dynamics of rigid bodies, etc. It is closely related to optimization of functional and is expressed in terms of definite integrals pertaining to functions and their derivatives. In the past few years, a number of contributions have been made towards the duality results for multiobjective variational problems. For the first time, Hanson [33] established and developed the linkage between classical calculus of variation and mathematical programming. After that, Mond and Hanson [34] derived optimality and duality results for scalar valued variational problems in view of convexity assumptions. Chandra et al. [35] studied optimality and duality for a class of non-differentiable variational problems. In this sequence, Bector and Husain [36] investigated duality for multiobjective variational problems. Nahak and Nanda [37] and Chen [38] constructed duality results for multiobjective variational problems with invexity. Some years later, Bhatia and Mehra [39] extended further the results of Mond et al. [40] and explored the optimality conditions and duality results for multiobjective variational problems with generalized B-invexity.

The concept of invexity is of great significance in variational problems and mathematical programming. Hanson [41] introduced the notion of invexity to mathematical programming. Mishra and Mukherjee [42] presented duality results for multiobjective FVPs. Furthermore, Mond and Husain [43] also investigated sufficient optimality criteria and duality for variational problems with generalized invexity. It is clearly observed that the duality results derived for variational problems presented in refs. [40,42,43] that hold for convex functions are also well-fitted for the wide range of invex functions. Weir and Mond [44] considered the concept of weak minima to derive the duality results for multiobjective programming problems. Different scalar duality results have also been extended for multiobjective programming problems by Weir and Mond [44]. Mukherjee and Mishra [45] have considered the concept of weak minima in the continuous case and have delivered a complete generalization of the results of Weir and Mond [44] to multiobjective variational problems. Moreover, they also relaxed the generalized convexity conditions to generalized invexity conditions.

Recently, Kumar [46] extended the invexity for continuous functions to invexity of order m. They further generalized the invexity of order m to ρ-pseudoinvexity type-I of order m, ρ-pseudoinvexity type-II of order m, as well as ρ-quasi-invexity type-I and type-II of order m. In 2016, Kumar et al. [47] also analysed the multiobjective FVP under F-Kuhn–Tucker (KT) pseudoinvexity conditions. Hachimi and Aghezzaf [48] established the mixed duality results and the sufficient optimality conditions concerning multiobjective variational problems under generalized (F,α,ρ,d)-type I functions which assimilate the several concepts of generalized type-I functions successfully. Later on, Mishra et al. [49] extended the generalized type-I invexity and duality for non-differentiable multiobjective variational problems. In 2014, Wolfe-type and Mond–Weir-type duality results were formulated for multiobjective variational control models under (ϕ,ρ)-invexity conditions by Antczak [50]. More recently, Upadhyay et al. [51] presented optimality conditions and duality for multiobjective semi-infinite programming problems on Hadamard manifolds utilizing generalized geodesic convexity. Moreover, Upadhyay et al. [52] also investigated Minty’s variational principle for non-smooth multiobjective optimization problems on Hadamard manifolds. Guo et al. [53] showed applications of symmetric gH-derivative to dual interval-valued optimization problems in a very efficient way. Furthermore, optimality conditions and duality for a class of generalized convex interval-valued optimization problems are recently investigated in works of Guo et al. [54].

The main purpose of this study is to derive the weak and strong duality results for multiobjective variational problems pertaining to a CF fractional derivative operator with exponential kernel. The CF fractional derivative possesses the non-singular kernel and so is better than the Caputo and RL fractional derivative operators. The proposed work presents the derivation of duality as well as strict converse duality theorems for variational problems with CF fractional derivative by employing some propositions and theorems of fractional calculus. In this article, we propounded first the optimality conditions for the variational problem. Furthermore, we present Theorem 1 which proves that a minimizer of the variational problem is a solution of the fractional Euler–Lagrange equation containing the CF fractional derivative. Now, we derived the formula for integration by parts for the CF fractional derivative in Proposition 1. The extended invexity definitions in view of CF fractional derivative operator along with Proposition 1 and Theorem 1 have been the key motivation behind the study of the variational problems with fractional calculus approach. Theorem 2 proves the fact that if a function is convex, the solution of the fractional Euler–Lagrange equation containing the CF fractional derivative will be a minimizer of the variational problem. Theorems 3 and 4 present the results for primal variational problems having a weak minimum. Furthermore, Theorems 5–10 are concerned with weak and strong duality results depending on the CF fractional derivative. Finally, Theorems 11 and 12 provide the strict converse duality results in view of the CF fractional derivative.

In the present work, the concept of weak minima has been considered and the generalizations of weak, strong, and strict converse duality results of Mukherjee and Mishra [45] have been extended to multiobjective variational problems pertaining to the CF fractional derivative operator. The remaining part of the article is organized as follows: In Section 2, we present the elemental definitions, formulae, and theorems regarding invexity and fractional derivative operators. Section 3 derives weak and strong duality results. Section 4 presents a strict converse duality result. Finally, Section 5 records the epilogue for the proposed work.

2 Basic definitions, theorems, and symbols

We follow these definitions and symbols in the present article.

Definition 1

[55]: The left and the right RL fractional derivatives of order α are defined by

D ξ α a y ( ξ ) = 1 Γ ( 1 α ) d d ξ a ξ ( ξ τ ) α y ( τ ) d τ ,

D b α ξ y ( ξ ) = 1 Γ ( 1 α ) d d ξ ξ b ( τ ξ ) α y ( τ ) d τ , α ( 0 , 1 ) .

Definition 2

[56]: The Caputo fractional derivative of y ( ξ ) : [ a , b ] of order α ∈ (0, 1) is stated as:

D a + α C y ( ξ ) = 1 Γ ( 1 α ) d d ξ a ξ 1 ( ξ τ ) α [ y ( τ ) y ( a ) ] d τ .

If yC 1, then

D a + α C y ( ξ ) = 1 Γ ( 1 α ) a ξ 1 ( ξ τ ) α y ( τ ) d τ .

As α → 1, D a + α C y ( ξ ) approaches to y′(ξ).

Definition 3

[25]: The new CF fractional derivative operator is described as follows:

D a + α CF y ( ξ ) = K ( α ) ( 1 α ) a ξ exp α ( ξ τ ) 1 α y ( τ ) d τ , α ( 0 , 1 ) ,

where K(α) signifies the normalization function with the property K(0) = K(1) = 1. Clearly, D a + α CF y ( ξ ) = 0 if y(ξ) is a constant function, i.e. the CF derivative of a constant function vanishes to zero same as the Caputo derivative but kernel of CF derivative does not have singularity for ξ = τ like the Caputo fractional derivative. It is remarkable that the CF fractional derivative has an exponential kernel.

Remark 1

Here we consider the value of K(α) as ( 1 α ) + α Γ ( α ) .

Remark 2

As α → 1, D a + α CF y ( ξ ) approaches to y′(ξ) and as α → 0, D a + α CF y ( ξ ) approaches to y(ξ) − y(a).

Definition 4

Abdeljawad and Baleanu [57] have defined the right CF fractional derivative as

D b α CF y ( ξ ) = K ( α ) ( 1 α ) ξ b exp α ( τ ξ ) 1 α y ( τ ) d τ , α ( 0 , 1 ) .

Definition 5

The first-order Sobolev space defined in the interval (a, b) is stated as H 1(a, b) = {xL 2(a, b) |x′ ∈ L 2(a, b)}, where x′ denotes the weak derivative of x.

Definition 6

[25]: Let yH 1(a, b), b > a, 0 < α < 1, then the CF fractional derivative is stated as in Definition (3), where K(α) specifies the normalization function with characteristic K(0) = K(1) = 1. If the function yH 1(a, b), then the derivative is formulated as follows:

D a + α CF y ( ξ ) = α K ( α ) ( 1 α ) a ξ exp α ( ξ τ ) 1 α [ y ( ξ ) y ( τ ) ] d τ .

Here, the CF fractional derivative has an exponential kernel.

Definition 7

[57]: Let x be a function in such a way that xH 1(a, b) a < b. The left Riemann fractional derivative of order α in the CF sense is given by

D a + α CFR y ( ξ ) = K ( α ) 1 α d d ξ a ξ exp α 1 α ( ξ τ ) y ( τ ) d τ ,

where a ξ , α ( 0 < α < 1 ) is a real number and K ( α ) is a normalization function depending on α with K ( 0 ) = K ( 1 ) = 1 .

Similarly, the right Riemann fractional derivative of order α in the CF sense can be written as follows:

D b α CFR y ( ξ ) = K ( α ) 1 α d d ξ ξ b exp α 1 α ( τ ξ ) y ( τ ) d τ ,

where ξ b .

Remark 3

When α 0 , lim α 0 D a + α CFR y ( ξ ) = d d ξ a ξ y ( τ ) d τ = y ( ξ ) .

Proposition 1

Let α ( 0 , 1 ) and y , z : [ a , b ] be two continuous functions of class C 1 [ a , b ] . Then the following formula for integration by parts holds:

a b y ( ξ ) D a + α CF z ( ξ ) d ξ = [ z ( ξ ) I b 1 α y ( ξ ) ] ξ = a ξ = b + a b z ( ξ ) D b α CFR y ( ξ ) d ξ .

Proof

We define the left and right auxiliary fractional integrals as

(1) I a + 1 α z ( ξ ) = K ( α ) ( 1 α ) a ξ exp α 1 α ( ξ τ ) z ( τ ) d τ ,

(2) I b 1 α z ( ξ ) = K ( α ) ( 1 α ) ξ b exp α 1 α ( τ ξ ) z ( τ ) d τ .

Now in view of Definition (3) and Eq. (1), it is concluded that

(3) D a + α CF z ( ξ ) = I a + 1 α d d ξ z ( ξ ) .

In the next step, we evaluate the integral a b y ( ξ ) D a + α CF z ( ξ ) d ξ as follows.

Using Eq. (3) along with further utilization of Theorem 1 of ref. [57] and integration by parts for classical derivatives, we obtain

a b y ( ξ ) D a + α CF z ( ξ ) d ξ = a b y ( ξ ) I a + 1 α d d ξ z ( ξ ) d ξ = a b d d ξ z ( ξ ) I b 1 α y ( ξ ) d ξ = a b I b 1 α y ( ξ ) d d ξ z ( ξ ) d ξ = [ z ( ξ ) I b 1 α y ( ξ ) ] ξ = a ξ = b a b z ( ξ ) × K ( α ) ( 1 α ) d d ξ ξ b exp α 1 α ( ξ τ ) y ( τ ) d τ d ξ .

Now in view of Definition 7, we obtain

a b y ( ξ ) D a + α CF z ( ξ ) d ξ = [ z ( ξ ) I b 1 α y ( ξ ) ] ξ = a ξ = b + a b z ( ξ ) D b α CFR y ( ξ ) d ξ .

Definition 8

(Optimality conditions for variational problems):

The following variational problem with the CF fractional derivative is considered here for given y C 1 ( a , b ) ,

(4) min V ( y ) = a b Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ ,

with y ( a ) = y a and y ( b ) = y b , where y ( a ) , y ( b ) . The assumptions are as follows:

  1. Q : [ a , b ] × 2 is continuously differentiable w.r.t. the second and third arguments.

  2. Given any x , the map ξ D b α CFR ( 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) ) = 0 is continuous.

Here, we denote i g ( y 1 , y 2 , y n ) = g y i ( y 1 , y 2 , y n ) for a function g : T n .

Theorem 1

Let y be a minimizer of the variational V(y) defined on E = { y C 1 ( a , b ) : y ( a ) = y a , y ( b ) = y b } , where y a , y b are fixed. Then y is a solution of the following fractional Euler–Lagrange equation:

(5) 2 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) + D b α CFR ( 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) ) = 0 , ξ [ a , b ] .

Proof

Let y be a solution for the functional V ( y ) . Assume y + δ ω be a variation of y with δ 1 , and ω : [ a , b ] R be a function of class C 1 [ a , b ] in such a way that the conditions ω ( a ) = ω ( b ) = 0 hold. Let ϑ ( δ ) = V ( y + δ ω ) . Since y satisfies Eq. (4) as a solution, the first variation of V must vanish, and hence ϑ ( 0 ) = 0 . Now, computing ϑ ( δ ) δ = 0 , equating to zero, and further employing Proposition 1, we have

a b 2 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) ω ( ξ ) d ξ + a b 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) D a + α CF ω ( ξ ) d ξ = a b 2 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) ω ( ξ ) d ξ + a b ω ( ξ ) D b α CFR ( 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) ) d ξ + [ ω ( ξ ) I b 1 α ( 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) ) ] ξ = a ξ = b = a b [ 2 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) + D b α CFR ( 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) ) ] ω ( ξ ) d ξ , ( ω ( a ) = ω ( b ) = 0 ) = 0 .

Now utilizing the boundary conditions ω ( a ) = ω ( b ) = 0 along with the assumption of arbitrariness of ω , we obtain the desired equation as:

2 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) + D b α CFR ( 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) ) = 0 ξ [ a , b ] .

Remark 4

Eq. (5) is called the Euler–Lagrange equation associated with the variational V ( y ) and the solutions for this equation are termed as extremals.

Remark 5

It is notable that Eq. (5) provides necessary criterion only. Now to obtain sufficient criterion, the concept of convex function is necessary to recall.

Definition 9

A function Q ( ξ , , ) is said to be convex in T R 3 if Q possesses continuous derivatives in respect of the second and third arguments and also satisfies the following inequality:

Q ( ξ , + 1 , + 1 ) Q ( ξ , , ) 2 Q ( ξ , , ) 1 + 3 Q ( ξ , , ) 1 , ( ξ , , ) , ( ξ , + 1 , + 1 ) T .

Theorem 2

If the function Q as described in Eq. ( 4 ) is convex in [ a , b ] × 2 , then each solution of the fractional Euler–Lagrange Eq. (5) minimizes V in E.

Proof

Let y be a solution for the fractional Euler–Lagrange Eq. (5). Assume y + δ ω to be a variation of y with δ 1 , and ω : [ a , b ] is a function that belongs to C 1 [ a , b ] such that the boundary conditions ω ( a ) = ω ( b ) = 0 hold. Now, we compute V ( y + δ ω ) V ( y ) in view of Definition 9 and Proposition 1 as follows:

(6) V ( y + δ ω ) V ( y ) = a b [ Q ( ξ , y ( ξ ) + δ ω ( ξ ) , D a + α CF y ( ξ ) + δ D a + α CF ω ( ξ ) ) Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) ] d ξ a b [ 2 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) δ ω ( ξ ) + 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) δ D a + α CF ω ( ξ ) ] d ξ = a b 2 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) δ ω ( ξ ) d ξ + a b 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) δ D a + α CF ω ( ξ ) d ξ = a b 2 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) δ ω ( ξ ) d ξ + δ [ ω ( ξ ) I b 1 α 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) ] ξ = a ξ = b + a b δ ω ( ξ ) D b α CFR 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ = a b 2 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) δ ω ( ξ ) d ξ + a b δ ω ( ξ ) D b α CFR 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ ( ω ( a ) = ω ( b ) = 0 ) , = a b [ 2 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) + D b α CFR 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) ] δ ω ( ξ ) d ξ .

Since y is a solution of the fractional Euler–Lagrange Eq. (5), we have

(7) 2 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) + D b α CFR 3 Q ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) = 0 .

Hence in view of Eq. (7), Eq. (6) provides the inequality as follows:

(8) V ( y + δ ω ) V ( y ) 0 .

Consequently, V ( y + δ ω ) V ( y ) , which implies that y is a local minimizer of V. □

Definition 10

Invexity definitions

Let Ω = [ a , b ] be a real interval. Let g : Ω × n × n be a continuously differentiable function. Consider g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) , where y : Ω n is a function of class C 1 [ a , b ] and D a + α CF y represents the CF derivative of order 0 < α < 1 of a function y. We denote the partial derivatives of g by

g ξ = g ξ , g y = g y 1 , g y 2 , g y 3 , , g y n ,

g D a + α CF y = g ( D a + α CF y 1 ) , g ( D a + α CF y 2 ) , g ( D a + α CF y 3 ) , , g ( D a + α CF y n ) .

Let Y be the space of piecewise smooth functions y : Ω n along with the norm y = y + D y , where the differential operator D is described as follows:

v = D y y ( ξ ) = y 0 + a ξ v ( s ) d s ,

where y 0 signifies the boundary value.

Let G : Y defined by G ( y ) = a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ be Fréchet differentiable. For notational convenience g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) will be written as g ( ξ , y , D a + α CF y ) . Here, it is assumed that g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) is convex in 3 if 2 g and 3 g exist and are continuous, and the condition

g ( ξ , y + y 1 , D a + α CF y + D a + α CF y 2 ) g ( ξ , y , D a + α CF y ) g y ( ξ , y , D a + α CF y ) y 1 + g D a + α CF y ( ξ , y , D a + α CF y ) D a + α CF y 2 ,

holds for every ( ξ , y , D a + α CF y ) , ( ξ , y + y 1 , D a + α CF y + D a + α CF y 2 ) 3 . Here, 2 g and 3 g denote the partial derivatives of g with respect to y and D a + α CF y , respectively.

Let y ¯ be a solution of the variational functional G ( y ¯ ) = a b g ( ξ , y ¯ ( ξ ) , D a + α CF y ¯ ( ξ ) ) d ξ y ¯ Y and ξ [ a , b ] . Let y = y ¯ + η ( ξ , y , y ¯ ) and η C 1 [ a , b ] with η ( ξ , y , y ¯ ) ξ = a = η ( ξ , y , y ¯ ) ξ = b . Clearly also η ( ξ , y , y ) = 0 .

Now utilizing the linearity property of the CF derivative operator and further the convexity assumption of g ( ξ , y ¯ ( ξ ) , D a + α CF y ¯ ( ξ ) ) , we have

G ( y ) G ( y ¯ ) = G ( y ¯ + η ) G ( y ¯ ) = a b [ g ( ξ , y ¯ ( ξ ) + η ( ξ , y , y ¯ ) , D a + α CF y ¯ ( ξ ) + D a + α CF η ( ξ , y , y ¯ ) ) g ( ξ , y ¯ ( ξ ) , D a + α CF y ¯ ( ξ ) ) ] d ξ a b [ g y ¯ ( ξ , y ¯ ( ξ ) , D a + α CF y ¯ ( ξ ) ) η ( ξ , y , y ¯ ) + g D a + α CF y ¯ ( ξ , y ¯ ( ξ ) , D a + α CF y ¯ ( ξ ) ) D a + α CF η ( ξ , y , y ¯ ) ] d ξ = a b [ η ( ξ , y , y ¯ ) g y ¯ ( ξ , y ¯ ( ξ ) , D a + α CF y ¯ ( ξ ) ) + ( D a + α CF η ( ξ , y , y ¯ ) ) g D a + α CF y ¯ ( ξ , y ¯ ( ξ ) , D a + α CF y ¯ ( ξ ) ) ] d ξ .

Clearly,

G ( y ) G ( y ¯ ) a b { η ( ξ , y , y ¯ ) g y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) + ( D a + α CF η ( ξ , y , y ¯ ) ) g D a + α CF y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) } d ξ , y , y ¯ Y .

Clearly as α 1 , the above obtained inequality reduces to

a b g ξ , y ( ξ ) , d d ξ y ( ξ ) d ξ a b g ξ , y ¯ ( ξ ) , d d ξ y ¯ ( ξ ) d ξ a b η ( ξ , y , y ¯ ) g y ¯ ξ , y ¯ , d d ξ y ¯ + d d ξ η ( ξ , y , y ¯ ) g d d ξ y ¯ ξ , y ¯ , d d ξ y ¯ d ξ ,

which is the definition of invexity in the continuous case extended by Mond et al. [40]. It is notable that if the function g is independent of ξ, the above given definition of invexity transforms to the inequality g ( y ) g ( y ¯ ) η ( y , y ¯ ) g y ¯ ( y ¯ ) , which is the fundamental interpretation of invexity prescribed by Hanson [41].

Example 1

The proposed inequality which is derived earlier is given as follows:

(9) G ( y ) G ( y ¯ ) a b { η ( ξ , y , y ¯ ) g y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) + ( D a + α CF η ( ξ , y , y ¯ ) ) g D a + α CF y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) } d ξ .

Let y ¯ = ξ , y = 2 ξ , η ( ξ , y , y ¯ ) = y y ¯ , g ( ξ , y ¯ , D a + α CF y ¯ ) = ξ + y ¯ + D a + α CF y ¯ , and g ( ξ , y , D a + α CF y ) = ξ + y + D a + α CF y .

Then

(10) η ( ξ , y , y ¯ ) = y y ¯ = ξ , D a + α CF η ( ξ , y , y ¯ ) = D a + α CF ξ ,

(11) g ( ξ , y ¯ , D a + α CF y ¯ ) = ξ + y ¯ + D a + α CF y ¯ = 2 ξ + D a + α CF ξ .

Now utilizing the formula of CF derivative, we obtain

(12) D a + α CF ξ = K ( α ) ( 1 α ) a ξ exp α 1 α ( ξ τ ) d τ = K ( α ) α 1 exp α 1 α ( ξ a ) ,

where K ( α ) = ( 1 α ) + α Γ ( α ) signifies the normalization function with the property K ( 0 ) = K ( 1 ) = 1 .

Thus,

(13) D a + α CF η ( ξ , y , y ¯ ) = K ( α ) α 1 exp α 1 α ( ξ a ) ,

and

(14) g ( ξ , y ¯ , D a + α CF y ¯ ) = 2 ξ + K ( α ) α 1 exp α 1 α ( ξ a ) .

Now

(15) g y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) = 2 + K ( α ) ( 1 α ) exp α 1 α ( ξ a ) , g D a + α CF y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) = 1 .

Now, we evaluate the term G ( y ) G ( y ¯ ) of the aforementioned proposed inequality (9) in view of Eq. (14) as follows:

(16) G ( y ) G ( y ¯ ) = a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ a b g ( ξ , y ¯ ( ξ ) , D a + α CF y ¯ ( ξ ) ) d ξ = 1 2 ( b 2 a 2 ) + K ( α ) α ( b a ) + K ( α ) ( 1 α ) α 2 exp α ( b a ) 1 α K ( α ) ( 1 α ) α 2 .

In this sequence, we also evaluate the integral term of the proposed inequality (9) in view of Eqs. (10) and (15) as follows:

(17) a b { η ( ξ , y , y ¯ ) g y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) + ( D a + α CF η ( ξ , y , y ¯ ) ) g D a + α CF y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) } d ξ = ( b 2 a 2 ) b K ( α ) α exp α ( b a ) 1 α + 2 K ( α ) ( 1 α ) α 2 exp α ( b a ) 1 α 2 K ( α ) ( 1 α ) α 2 + b K ( α ) α .

Case I: For a = 0 , b = 1 , α = 0.5 , K ( α ) = ( 1 α ) + α Γ ( α ) , the proposed inequality is satisfied.

Case II: For a = 0 , b = 1 , α = 0.2 , K ( α ) = ( 1 α ) + α Γ ( α ) , the proposed inequality is also satisfied.

Consequently, it is concluded that the proposed inequality (9) with CF fractional derivative holds well for the function g ( ξ , y ¯ , D a + α CF y ¯ ) = 2 ξ + 1 α α + 1 Γ ( α ) 1 exp α 1 α ξ , where 0 < α < 1 .

Now, we extend the definitions of invex, pseudoinvex (PIX), strictly pseudoinvex (SPIX), and quasi-invex (QIX) as described in ref. [58] with the CF fractional derivative of order 0 < α < 1 in the following way:

Definition 11

Invex

The functional G is stated as invex with respect to η if there exists a differentiable vector function η ( ξ , y , y ¯ ) C 1 [ a , b ] with η ( ξ , y , y ) = 0 such that y , y ¯ Y ,

G ( y ) G ( y ¯ ) a b { η ( ξ , y , y ¯ ) g y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) + ( D a + α CF η ( ξ , y , y ¯ ) ) g D a + α CF y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) } d ξ .

Definition 12

PIX

The functional G is stated as PIX w.r.t. η if a differentiable vector function η ( ξ , y , y ¯ ) C 1 [ a , b ] with η ( ξ , y , y ) = 0 such that y , y ¯ Y ,

a b { η ( ξ , y , y ¯ ) g y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) + ( D a + α CF η ( ξ , y , y ¯ ) ) g D a + α CF y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) } d ξ 0 G ( y ) G ( y ¯ ) ,

or equivalently,

G ( y ) < G ( y ¯ ) a b { η ( ξ , y , y ¯ ) g y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) + ( D a + α CF η ( ξ , y , y ¯ ) ) g D a + α CF y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) } d ξ < 0 .

Definition 13

SPIX

The functional G is stated as SPIX w.r.t. η if a differentiable vector function η ( ξ , y , y ¯ ) C 1 [ a , b ] with η ( ξ , y , y ) = 0 such that y , y ¯ Y ,

a b { η ( ξ , y , y ¯ ) g y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) + ( D a + α CF η ( ξ , y , y ¯ ) ) g D a + α CF y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) } d ξ 0 G ( y ) > G ( y ¯ ) ,

or equivalently,

G ( y ) G ( y ¯ ) a b { η ( ξ , y , y ¯ ) g y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) + ( D a + α CF η ( ξ , y , y ¯ ) ) g D a + α CF y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) } d ξ < 0 .

Definition 14

QIX

The functional G is stated as QIX w.r.t. η if a differentiable vector function η ( ξ , y , y ¯ ) C 1 [ a , b ] with η ( ξ , y , y ) = 0 such that y , y ¯ Y ,

a b { η ( ξ , y , y ¯ ) g y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) + ( D a + α CF η ( ξ , y , y ¯ ) ) g D a + α CF y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) } d ξ > 0 G ( y ) > G ( y ¯ ) ,

or equivalently,

G ( y ) G ( y ¯ ) a b { η ( ξ , y , y ¯ ) g y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) + ( D a + α CF η ( ξ , y , y ¯ ) ) g D a + α CF y ¯ ( ξ , y ¯ , D a + α CF y ¯ ) } d ξ 0 .

In the aforementioned definitions, D a + α CF η ( ξ , y , y ¯ ) is the vector whose ith component is ( d α / d ξ α ) η i ( ξ , y , y ¯ ) . Let g ( ξ , y , D a + α CF y ( ξ ) ) be a real scalar function and h ( ξ , y , D a + α CF y ( ξ ) ) be an m-dimensional function with continuous derivatives up to the second order with respect to each of its arguments. Here, y is an n-dimensional function of ξ and D a + α CF y ( ξ ) denotes the CF fractional derivative of order α with respect to ξ where 0 < α < 1 .

Now we deal with the multiobjective variational primal problem, as discussed in the work of Mukherjee and Mishra [45], with the CF fractional derivative operator in the following way:

(P) Minimize

a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ ,

subject to y ( a ) = y 0 , y ( b ) = y 1 ,

h ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) 0 ,

where g : [ a , b ] × n × n p and h : [ a , b ] × n × n m , 0 < α < 1 .

For the primal problem (P), a point y 0 is referred to as a weak minimum if there exists no other feasible point y for which the following inequality will hold

(18) a b g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ > a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ .

Now, we frame the continuous versions of Theorems 2.1 and 2.2, as described in the work of Weir and Mond [44], involving fractional derivative operators with exponential kernel in the following way:

Theorem 3

Let y = y 0 be a weak minimum for the primal problem (P). Then λ p , z m such that

(19) λ T g y ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) + z ( ξ ) T h y ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) = D b α CFR [ λ T g D a + α CF y ( ξ ) ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) + z ( ξ ) T h D a + α CF y ( ξ ) ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) ] ,

(20) z ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) = 0 ,

(21) ( λ , y ) 0 .

Proof

The proof is easily established through Theorem 1.□

Theorem 4

Let the primal problem (P) have a weak minimum at a point y 0 , which satisfies the KT constraint qualification. Then λ p , z m such that

(22) λ T g y ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) + z ( ξ ) T h y ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) = D b α CFR [ λ T g D a + α CF y ( ξ ) ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) + z ( ξ ) T h D a + α CF y ( ξ ) ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) ] ,

(23) z ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) = 0 ,

(24) z ( ξ ) 0 ,

(25) λ ( ξ ) 0 , λ T e = 1 ,

where e = ( 1 , .... . , 1 ) p .

Proof

The proof is easily established through Theorem 1.□

3 Duality

In relation to the primal problem (P), the dual problem (D) as discussed in ref. [45] is considered with fractional derivative operators pertaining to exponential kernel in the following way:

(D) Maximize

a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ ,

subject to

(26) y ( a ) = y 0 , y ( b ) = y 1 ,

(27) g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e = D b α CFR [ g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e ] ,

(28) z 0 .

λ Λ ,

where

(29) Λ = { λ p : λ 0 , λ T e = 1 } .

In upcoming steps, we discuss the duality theorems, as discussed in the work of Mukherjee and Mishra [45], with the CF fractional derivative operator of order 0 < α < 1 as follows:

Theorem 5

(Weak duality): If, for all feasible ( y , v , z , λ ) ,

  1. a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ is PIX or

  2. a b { λ T g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ is PIX, then

a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ .

Proof

(a) Let y be feasible for (P) and ( v , z , λ ) feasible for (D).

From Eq. (27), we have

(30) a b η ( ξ , y , v ) { g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ = a b η ( ξ , y , v ) D b α CFR { g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ .

Suppose contrary to the result given in statement of Theorem 5, i.e.

(31) a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ < a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ a b { g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) + z ( ξ ) T h ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) e } d ξ

(31) < a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ .

Now, use of pseudoinvexity of a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ along with the above obtained inequality (31) provides

a b { η ( ξ , y , v ) [ g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e ] + D a + α CF η ( ξ , y , v ) [ g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e ] } d ξ < 0 ,

or

(32) a b η ( ξ , y , v ) [ g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e ] d ξ + a b [ g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e ] D a + α CF η ( ξ , y , v ) d ξ < 0 .

Now using Proposition 1 in the above obtained inequality (32), we have

a b η ( ξ , y , v ) [ g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e ] d ξ + [ η ( ξ , y , v ) I b 1 α { g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } ] ξ = a ξ = b + a b η ( ξ , y , v ) D b α CFR { g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ < 0 .

η ( ξ , y , y ) = 0 ,

we obtain

(33) a b η ( ξ , y , v ) [ g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e ] d ξ + a b η ( ξ , y , v ) D b α CFR { g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ < 0 .

Now in view of Eq. (30) and inequality (33), we have

a b η ( ξ , y , v ) [ g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e ] d ξ + a b η ( ξ , y , v ) D b α CFR { g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ) e } d ξ = 0 < 0 ,

which is a contradiction.

Hence,

a b { g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) + z ( ξ ) T h ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) e } d ξ a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ .

Thus, the supposition (31), which is contrary to the result given in the statement of Theorem 5, is wrong and consequently

a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ .

(b) Let y be feasible for (P) and ( v , z , λ ) feasible for (D).

Multiplying Eq. (30) by λ T and using λ T e = 1 , we obtain

(34) a b η ( ξ , y , v ) { λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ = a b η ( ξ , y , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ .

Suppose contrary to the result given in the statement of Theorem 5, i.e.

(35) a b { g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) + z ( ξ ) T h ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) e } d ξ < a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ .

Multiplying the aforementioned obtained inequality (35) by λ T and using λ T e = 1 , we obtain

(36) a b { λ T g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) + z ( ξ ) T h ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) } d ξ < a b { λ T g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ .

Now, the pseudoinvexity of a b { λ T g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ along with the above obtained inequality (36) provides

a b { η ( ξ , y , v ) [ λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] + D a + α CF η ( ξ , y , v ) [ λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] } d ξ < 0 ,

or

(37) a b η ( ξ , y , v ) [ λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] d ξ + a b [ λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] D a + α CF η ( ξ , y , v ) d ξ < 0 .

Now by using Proposition 1 in the above obtained inequality (37), we have

(38) a b η ( ξ , y , v ) [ λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] d ξ + [ η ( ξ , y , v ) I b 1 α { λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } ] ξ = a ξ = b + a b η ( ξ , y , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ < 0 .

η ( ξ , y , y ) = 0 , the above obtained inequality (38) reduces to

(39) a b η ( ξ , y , v ) [ λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] d ξ + a b η ( ξ , y , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ < 0 .

Now in view of Eq. (34) and inequality (39), we have

a b η ( ξ , y , v ) [ λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] d ξ + a b η ( ξ , y , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } d ξ = 0 < 0 ,

which is a contradiction. Thus, the supposition (35) is false. Consequently,

a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ .

Theorem 6

(Strong duality): Let the primal (P) have a weak minimum at y 0 , which satisfies the KT constraint qualification. Then there exists ( z , λ ) in such a way that ( y 0 , z , λ ) is feasible for (D) and the objective values of (P) and (D) are equal. If also,

  1. a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ is PIX or

  2. a b { λ T g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ is PIX,

then ( y 0 , z , λ ) is a weak maximum for (D).

Proof

Since the primal (P) has a weak minimum at y 0 , which satisfies the KT constraint qualification, by Theorem 4, z 0 , λ 0 , λ T e = 1 such that

λ T g y ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) + z ( ξ ) T h y ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) = D b α CFR [ λ T g D a + α CF y ( ξ ) ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) + z ( ξ ) T h D a + α CF y ( ξ ) ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) ] ,

and z ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) = 0 . Hence, ( y 0 , z , λ ) is feasible for dual (D) and the objective values of primal (P) and dual (D) are equal.

If ( y 0 , z , λ ) is not a weak maximum for (D), then a feasible ( v , z , λ ) for dual (D) occurs so that

(40) a b { g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) + z ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) e } d ξ < a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ .

(a) The pseudoinvexity of a b { g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e } d ξ together with the above obtained inequality (40) gives

(41) a b { η ( ξ , y 0 , v ) [ g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e ] + D a + α CF η ( ξ , y 0 , v ) [ g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) e ] } d ξ < 0 .

Since λ 0 thus multiplying inequality (41) by λ and using λ ( ξ ) T e = 1 , we obtain

(42) a b { η ( ξ , y 0 , v ) [ λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] + D a + α CF η ( ξ , y 0 , v ) [ λ ( ξ ) T g D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) ] } d ξ < 0 .

Now by using Proposition 1 in inequality (42), we have

(43) a b η ( ξ , y 0 , v ) [ λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] d ξ + [ η ( ξ , y 0 , v ) I b 1 α { λ ( ξ ) T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } ] ξ = a ξ = b + a b η ( ξ , y 0 , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } d ξ < 0 .

η ( ξ , y , y ) = 0 , inequality (43) reduces to

(44) a b η ( ξ , y 0 , v ) [ λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] d ξ + a b η ( ξ , y 0 , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } d ξ < 0 .

Since ( v , z , λ ) is feasible for (D) thus from Eq. (34), we have

(45) a b η ( ξ , y 0 , v ) { λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ

= a b η ( ξ , y 0 , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } d ξ .

Using Eq. (45) in inequality (44), we obtain

a b η ( ξ , y 0 , v ) [ λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] d ξ + a b η ( ξ , y 0 , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } d ξ = 0 < 0 .

This is clearly a contradiction. Hence, the supposition of the existence of feasibility of ( v , z , λ ) for (D) is false and consequently ( y 0 , z , λ ) is a weak maximum for (D). □

(b) The proof of part (b) is very similar to that of part (a).

We now consider the following dual problem (D1) in relation to the primal problem (P).

(D1) Maximize

(46) a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ subject to y ( a ) = y 0 , y ( b ) = y 1 ,

(47) λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e = D b α CFR [ λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] ,

(48) z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) 0 ,

(49) z 0 .

λ Λ , where

(50) Λ = { λ p : λ 0 , λ T e = 1 } .

The next step is to establish the weak and strong duality theorems, as discussed in the work of Mukherjee and Mishra [45], with the CF fractional derivative operators.□

Theorem 7

(Weak duality): If for all feasible ( y , v , z , λ ) ,

  1. a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is PIX and a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is QIX or

  2. a b λ T g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is PIX and a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is QIX or

  3. a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is QIX and a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is SPIX or

  4. a b λ T g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is QIX and a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is SPIX.

Then

a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ .

Proof

(a) Let y be feasible for primal (P) and ( v , z , λ ) feasible for dual (D1).

Suppose contrary to the result, i.e.

(51) a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ < a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ .

Thus, in view of the pseudoinvexity of a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ , we have

a b { η ( ξ , y , v ) g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , y , v ) ) g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ < 0 .

λ 0 , we have

(52) a b { η ( ξ , y , v ) λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , y , v ) ) λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ < 0 .

From the constraint of the primal problem (P), we have

h ( ξ , y ( ξ ) , D a + α C F y ( ξ ) ) 0 .

z ( ξ ) 0 therefore

(53) z ( ξ ) T h ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) 0 .

Now from the constraint of the dual problem (D1), we have

(54) z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) 0 .

Combining both inequalities (53) and (54), we obtain

z ( ξ ) T h ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) 0 z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ,

or

z ( ξ ) T h ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ,

or

(55) a b z ( ξ ) T h ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ .

Now, inequality (55) together with the quasi-invexity of a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ implies that

(56) a b { η ( ξ , y , v ) z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , y , v ) ) z ( t ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } d ξ 0 .

Now adding inequalities (52) and (56), we obtain

(57) a b { η ( ξ , y , v ) [ λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] } d ξ + a b { ( D a + α CF η ( ξ , y , v ) ) [ λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) ] } d ξ < 0 .

Now by using Proposition 1 in inequality (57), we have

(58) a b { η ( ξ , y , v ) [ λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] } d ξ + [ η ( ξ , y , v ) I b 1 α { λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } ] ξ = a ξ = b + a b η ( ξ , y , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ < 0 .

Using Eq. (47) and η ( ξ , y , y ) = 0 , inequality (58) reduces to

a b { η ( ξ , y , v ) [ λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] } d ξ + a b η ( ξ , y , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } d ξ = 0 < 0 ,

which is a contradiction. Thus,

a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ .

(b) The proof of part (b) is very similar to the proof of part (a).

(c) Let y be feasible for (P) and ( v , z , λ ) feasible for (D1).

Suppose contrary to the result, i.e.

(59) a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ < a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ .

Then in view of the quasi-invexity of a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ , we have

a b { η ( ξ , y , v ) g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , y , v ) ) g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ 0 .

λ 0 , we have

(60) a b { η ( ξ , y , v ) λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , y , v ) ) λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ 0 .

(61) a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ .

Thus, the strict pseudoinvexity of a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ along with inequality (61) provides

(62) a b { η ( ξ , y , v ) z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , v , v ) ) z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } d ξ < 0 .

Now adding inequalities (60) and (62), we obtain

(63) a b { η ( ξ , y , v ) [ λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] } d ξ + a b { ( D a + α CF η ( ξ , y , v ) ) [ λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) ] } d ξ < 0 .

Now using Proposition 1 in inequality (63) along with η ( ξ , y , y ) = 0 , and from Eq. (47), we obtain

a b { η ( ξ , y , v ) [ λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] } d ξ + a b η ( ξ , y , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } d ξ = 0 < 0 ,

which is a contradiction. Consequently,

a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ .

  1. The proof of part (d) is very similar to the proof of part (c).□

Theorem 8

(Strong duality): Let the primal (P) have a weak minimum at the point y 0 , which satisfies the KT constraint qualification. Then ( z , λ ) so that ( y 0 , z , λ ) is feasible for (D1) and the objective values of primal (P) and dual (D1) are equal. If also,

  1. a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is PIX and a b z ( t ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is QIX or

  2. a b λ T g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is PIX and a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is QIX or

  3. a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is QIX and a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is SPIX or

  4. a b λ T g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is QIX and a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is SPIX,

then ( y 0 , z , λ ) is a weak maximum for (D1).

Proof

Since the primal (P) have a weak minimum at y 0 for which the KT constraint qualification is satisfied, then by Theorem 4, there exists z 0 , λ 0 , λ T e = 1 such that

λ T g y ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) + z ( ξ ) T h y ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) = D b α CFR [ λ T g D a + α CF y ( ξ ) ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) + z ( ξ ) T h D a + α CF y ( ξ ) ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) ] ,

and z ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) = 0 . Therefore, ( y 0 , z , λ ) is feasible for (D1) and the objective values of primal (P) and dual (D1) are equal.

If ( y 0 , z , λ ) is not a weak maximum for (D1), then a feasible solution ( v , z , λ ) for dual (D1) occurs such that

(64) a b g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ < a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ .

(a) The pseudoinvexity of a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ together with inequality (64) gives

(65) a b { η ( ξ , y 0 , v ) g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , y 0 , v ) ) g D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } d ξ < 0 .

Since λ 0 thus multiplying the above obtained inequality (65) by λ , and using λ ( ξ ) T e = 1 , we obtain

(66) a b { η ( ξ , y 0 , v ) λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , y 0 , v ) ) λ ( ξ ) T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ < 0 .

From the constraint of the primal problem (P), we have

h ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) 0 .

Since y 0 is feasible for primal problem (P) thus it satisfies the aforementioned constraint and so h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) 0 .

z ( t ) 0

thus

(67) z ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) 0 .

Now from the constraint of the dual problem (D1), we have

z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) 0 .

Since ( v , z , λ ) is feasible for (D1) thus

(68) z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) 0 .

Combining both inequalities (67) and (68), we obtain

z ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) 0 z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ,

or

z ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ,

or

(69) a b z ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ .

The quasi-invexity of a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ together with the above obtained inequality (69) provides

(70) a b { η ( ξ , y 0 , v ) z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , y 0 , v ) ) z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ 0 .

Now adding inequalities (66) and (70), we obtain

(71) a b { η ( ξ , y 0 , v ) [ λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] + ( D a + α CF η ( ξ , y 0 , v ) ) [ λ T g D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) ] } d ξ < 0 .

Now using Proposition 1 in the above obtained inequality (71) and η ( ξ , y 0 , y 0 ) = 0 , we obtain

(72) a b η ( ξ , y 0 , v ) [ λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] d ξ + a b ( η ( ξ , y 0 , v ) ) D b α CFR [ λ T g D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) ] d ξ < 0 .

Since ( v , z , λ ) is feasible for (D1) thus from Eq. (47), we have

(73) a b η ( ξ , y 0 , v ) { λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ = a b η ( ξ , y 0 , v ) D b α CFR { λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) } d ξ .

Using Eq. (73) in inequality (72), we obtain

a b η ( ξ , y 0 , v ) { λ T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ + a b ( η ( ξ , y 0 , v ) ) D b α CFR [ λ T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) ] d ξ = 0 < 0 ,

which is a discrepancy and consequently contradicts the feasibility of ( v , z , λ ) . □

(b) The proof of part (b) is very similar to the proof of part (a).

(c) The quasi-invexity of a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ together with inequality (64) gives

(74) a b { η ( ξ , y 0 , v ) g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , y 0 , v ) ) g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ 0 .

Since λ 0 thus multiplying inequality (74) by λ and using λ ( ξ ) T e = 1 , we obtain

(75) a b { η ( ξ , y 0 , v ) λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , y 0 , v ) ) λ ( ξ ) T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ 0 .

Since

(76) a b z ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ .

( a b z ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) 0 and a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) 0 ).

Thus, in view of strictly pseudoinvexity of a b z ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ , we have

(77) a b { η ( ξ , y 0 , v ) z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + ( D a + α CF η ( ξ , y 0 , v ) ) z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) } d ξ < 0 .

Now combining inequalities (75) and (77), we obtain

(78) a b { η ( ξ , y 0 , v ) [ λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] + ( D a + α CF η ( ξ , y 0 , v ) ) [ λ T g D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) ] } d ξ < 0 .

Now using Proposition 1 in inequality (78) and utilizing η ( ξ , y 0 , y 0 ) = 0 , we obtain

(79) a b η ( ξ , y 0 , v ) [ λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] d ξ + a b ( η ( ξ , y 0 , v ) ) D b α CFR [ λ T g D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) ] d ξ < 0 .

Now in view of Eq. (73), inequality (79) can be expressed as

a b η ( ξ , y 0 , v ) [ λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] d ξ + a b ( η ( ξ , y 0 , v ) ) D b α CFR [ λ T g D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) + z T h D a + α CF v ( ξ ) ( ξ , v , D a + α CF v ) ] d ξ = 0 < 0 ,

which clearly contradicts the feasibility of ( v , z , λ ) . □

(d) The proof of part (d) is very much identical to the proof of part (c).

Now we state the continuous form of a general primal (CP) and dual (CD) for the multiobjective variational optimization problem, as discussed in the work of Mukherjee and Mishra [45], involving fractional derivative operator with exponential kernel. Consider the problem given as follows:

(CP) Minimize

a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ ,

subject to y ( a ) = y 0 , y ( b ) = y 1

h ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) 0 ,

ϑ ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) = 0 ,

where g : I × n × n k , h : I × n × n m , and ϑ : I × n × n l are all differentiable.

Let K = { 1 , 2 , , k } , U = { 1 , 2 , , u } , I β K , β = 0 , 1 , , κ with I β I γ = φ , β γ , and β = 0 κ I β = K and J β U , β = 0 , 1 , , κ with J β J γ = φ , β γ , and β = 0 κ J β = U . It is notable that any particular I β or J β may be empty. Hence, if K has κ 1 disjoint subsets and U has κ 2 disjoint subsets, κ = max [ κ 1 , κ 2 ] . So, if κ 1 > κ 2 , then J β , β > κ 2 is empty.

In connection to (CP), we investigate the problem:

(CD) Maximize

a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + i I 0 z i ( ξ ) T h i ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e + j J 0 w j ( ξ ) T ϑ j ( ξ , v , D a + α CF v ) e d ξ ,

subject to

y ( a ) = y 0 , y ( b ) = y 1

λ ( ξ ) T g v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + w ( ξ ) T ϑ v ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) = D b α CFR [ λ ( ξ ) T g D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + z ( ξ ) T h D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + w ( ξ ) T ϑ D a + α CF v ( ξ ) ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) ] ,

i I β z i ( ξ ) T h i ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + j J β w j ( ξ ) T ϑ j ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) 0 ,

β = 0 , 1 , , κ , y 0 , λ Λ , where Λ = { λ p : λ 0 , λ T e = 1 } .

The weak and strong duality theorems, as described in ref. [45], are articulated here with the CF fractional derivative operators without proof because their proof may be delivered in a very identical manner to that of Theorems (5)–(8).□

Theorem 9

(Weak duality): If, for all feasible ( y , v , z , w , λ ) ,

  1. a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + i I 0 z i ( ξ ) T h i ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e + j J 0 w j ( ξ ) T ϑ j ( ξ , v , D a + α CF v ) e d ξ is PIX and a b i I α z i ( ξ ) T h i ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + j J α w j ( ξ ) T ϑ j ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ , β = 1 , 2 , , κ , is QIX; or

  2. a b λ ( t ) T g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + i I 0 z i ( ξ ) T h i ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + j J 0 w j ( ξ ) T ϑ j ( ξ , v , D a + α CF v ) d ξ is PIX and a b i I α z i ( ξ ) T h i ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + j J α w j ( ξ ) T ϑ j ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ , β = 1 , 2 , , κ , is QIX; or

  3. I 0 K and J 0 U ,

    a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + i I 0 z i ( ξ ) T h i ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e + j J 0 w j ( ξ ) T ϑ j ( ξ , v , D a + α CF v ) e d ξ ,

    is QIX and a b i I α z i ( ξ ) T h i ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + j J α w j ( ξ ) T ϑ j ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ , α = 1 , 2 , , v , is SPIX; or

  4. I 0 K and J 0 U ,

a b λ ( t ) T g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + i I 0 z i ( ξ ) T h i ( ξ , v , D a + α CF v ) + j J 0 w j ( ξ ) T ϑ j ( ξ , v , D a + α CF v ) d ξ ,

is QIX and a b i I α z i ( ξ ) T h i ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + j J α w j ( ξ ) T ϑ j ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ , α = 1 , 2 , , v , is SPIX,

then

a b g ( ξ , y ( ξ ) , D a + α CF y ( ξ ) ) d ξ a b g ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) + i I 0 z i ( ξ ) T h i ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e + j J 0 w j ( ξ ) T ϑ j ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) e .

Theorem 10

(Strong duality): Let (CP) have a weak minimum at y 0 , which satisfies the KT constraint qualification. Then ( z , w , λ ) , so that ( y 0 , z , w , λ ) is feasible for dual (CD) and the objective values of primal (CP) and dual (CD) are equal. If one of the presumptions (a), (b), (c), or (d) of Theorem 9 is fulfilled, then ( y 0 , z , w , λ ) is a weak maximum for (CD).

Now, we establish the strict converse duality theorem as stated in ref. [45] with CF fractional derivative operator in the forthcoming section.

4 Strict converse duality

Theorem 11

Let y 0 be a weak minimum for (P) and ( v 0 , z 0 , λ 0 ) be a weak maximum for (D1) such that a b λ 0 T ( ξ ) g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ a b λ 0 T g ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) d ξ .

Assume that

  1. a b λ 0 T ( ξ ) g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ is SPIX at v 0 and a b z 0 ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is QIX at v 0 ; or

  2. a b λ 0 T ( ξ ) g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ is QIX at v 0 and a b z 0 ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ is SPIX at v 0;

then y 0 = v 0 , i.e. v 0 is a weak minimum for (P).

Proof

It is assumed that y 0 v 0 . Since y 0 and ( v 0 , z 0 , λ 0 ) are feasible for (P) and (D1), respectively. Thus, h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) 0 and z 0 ( ξ ) T h ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) 0 .

Now, the aforementioned inequalities imply that a b z 0 ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ 0 , and a b z 0 ( ξ ) T h ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) d ξ 0 . Finally, they can be written in the combined form as

(80) a b z 0 ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ a b z 0 ( ξ ) T h ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) d ξ .

(a) Now, inequality (80) in view of quasi-invexity of a b z 0 ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ at v 0 provides the following inequality:

(81) a b { η ( ξ , y 0 , v 0 ) z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + ( D a + α CF η ( ξ , y 0 , v 0 ) ) z 0 ( ξ ) T h D a + α CF v 0 ( ξ ) ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) } d ξ 0 .

Since it is given that a b λ 0 T ( ξ ) g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ a b λ 0 T g ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) d ξ , thus in view of strict pseudoinvexity of a b λ 0 ( ξ ) T g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ , we have

(82) a b { η ( ξ , y 0 , v 0 ) λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + ( D a + α CF η ( ξ , y 0 , v 0 ) ) λ 0 ( ξ ) T g D a + α CF v 0 ( ξ ) ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) } d ξ < 0 .

Now combining inequalities (81) and (82), we obtain

(83) a b { η ( ξ , y 0 , v 0 ) [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) ] + ( D a + α CF η ( ξ , y 0 , v 0 ) ) [ λ 0 ( ξ ) T g D a + α C F v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) + z 0 ( ξ ) T h D a + α C F v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) ] } d ξ 0 .

Now using Proposition 1 in inequality (83) and η ( ξ , y 0 , y 0 ) = 0 , we obtain

(84) a b { η ( ξ , y 0 , y 0 ) [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) ] d ξ + a b η ( ξ , y 0 , v 0 ) D b α CFR [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) ] d ξ 0 .

Since ( v 0 , z 0 , λ 0 ) is feasible for (D1) thus from Eq. (47), we have

(85) a b η ( ξ , y 0 , v 0 ) { λ 0 T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) } d ξ = a b η ( ξ , y 0 , v 0 ) D b α CFR { λ 0 T g D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) + z 0 ( ξ ) T h D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) } d ξ .

In view of Eq. (85), inequality (84) can be written as:

a b η ( ξ , y 0 , v 0 ) [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) ] d ξ + a b η ( ξ , y 0 , v 0 ) D b α CFR [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 , D a + α CF v 0 ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 , D a + α CF v 0 ) ] d ξ = 0 0 ,

which contradicts the assumption y 0 v 0 . Consequently, y 0 = v 0 and since y 0 is a weak minimum for (P), thus v 0 is a weak minimum for (P). □

(b) Since a b z 0 ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ a b z 0 ( ξ ) T h ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) d ξ thus in view of the strictly pseudoinvexity of a b z 0 ( ξ ) T h ( ξ , v ( ξ ) , D a + α CF v ( ξ ) ) d ξ at v 0, we have

(86) a b { η ( ξ , y 0 , v 0 ) z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + ( D a + α CF η ( ξ , y 0 , v 0 ) ) z 0 ( ξ ) T h D a + α CF v 0 ( ξ ) ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) } d ξ < 0 .

Since it is given that a b λ 0 T ( ξ ) g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ a b λ 0 T g ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) d ξ , thus in view of the quasi-invexity of a b λ 0 ( ξ ) T g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ , we have

(87) a b { η ( ξ , y 0 , v 0 ) λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + ( D a + α CF η ( ξ , y 0 , v 0 ) ) λ 0 ( ξ ) T g D a + α CF v 0 ( ξ ) ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) } d ξ 0 .

Now combining inequalities (86) and (87), we obtain

(88) a b { η ( ξ , y 0 , v 0 ) [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) ] + ( D a + α CF η ( ξ , y 0 , v 0 ) ) [ λ 0 ( ξ ) T g D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) + z 0 ( ξ ) T h D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) ] } d ξ 0 .

Now using Proposition 1 in inequality (88) along with η ( ξ , y 0 , y 0 ) = 0 , we obtain

(89) a b η ( ξ , y 0 , v 0 ) [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) ] d ξ + a b η ( ξ , y 0 , v 0 ) D b α CFR [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) ] d ξ 0 .

In view of Eq. (85), inequality (89) can be written as

a b η ( ξ , y 0 , v 0 ) [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) ] d ξ + a b η ( ξ , y 0 , v 0 ) D b α CFR [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 , D a + α CF v 0 ) ] d ξ = 0 0 ,

which contradicts the assumption y 0 v 0 . Consequently, y 0 = v 0 and v 0 is a weak minimum for (P).□

In the next step, we state the theorem as stated in ref. [45] with the CF fractional derivative operator in the following way:

Theorem 12

Let y 0 be a weak minimum for (CP) and ( v 0 , z 0 , w 0 , λ 0 ) be a weak maximum for (CD) such that

a b λ 0 ( ξ ) T g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ a b { λ 0 ( ξ ) T g ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + i I 0 z 0 i ( ξ ) T × h i ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + j J 0 w 0 j ( ξ ) T × ϑ j ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) d ξ .

If

  1. a b λ 0 ( ξ ) T g ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + i I 0 z 0 i ( ξ ) T × h i ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + j J 0 w 0 j ( ξ ) T × ϑ j ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) d ξ is SPIX at v 0 and each a b i I α z 0 i ( ξ ) T × h i ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + j J α w 0 j ( ξ ) T × ϑ j ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) d ξ , β = 1 , 2 , 3 , κ , is QIX at v 0 ; or

  2. a b λ 0 ( ξ ) T g ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + i I 0 z 0 i ( ξ ) T × h i ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + j J 0 w 0 j ( ξ ) T × ϑ j ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) d ξ

is QIX at v 0 and each a b i I α z 0 i ( ξ ) T × h i ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + j J α w 0 j ( ξ ) T × ϑ j ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) d ξ , β = 1 , 2 , 3 , κ , is SPIX at v 0 then y 0 = v 0 , i.e. v 0 is a weak minimum for (P).

Corollary

Suppose y 0 is a weak minimum for (P) and ( v 0 , z 0 , λ 0 ) is a weak maximum for (D) so that a b λ 0 ( ξ ) T g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ a b { λ 0 ( ξ ) T g ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) } d ξ .

If a b { λ 0 ( ξ ) T g ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) } d ξ is SPIX at v 0 , then y 0 = v 0 , i.e. v 0 is a weak minimum for (P).

Proof

It is assumed that y 0 v 0 . Since y 0 is feasible for (P) thus h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) 0 , and ( v 0 , z 0 , λ 0 ) is feasible for (D) thus from Eq. (27), we have

(90) a b η ( ξ , y 0 , v 0 ) { g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) e } d ξ = a b η ( ξ , y 0 , v 0 ) D b α CFR { g D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) + z 0 ( ξ ) T h D a + α CF v 0 ( ξ ) ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) e } d ξ .

Multiplying Eq. (90) by λ 0 ( ξ ) T and further using λ 0 ( ξ ) T e = 1 , we obtain

(91) a b η ( ξ , y 0 , v 0 ) { λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) } d ξ = a b η ( ξ , y 0 , v 0 ) D b α CFR { λ 0 ( ξ ) T g D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) + z 0 ( ξ ) T h D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) } d ξ .

Since

(92) a b λ 0 ( ξ ) T g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ a b { λ 0 ( ξ ) T g ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) } d ξ ,

and also h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) 0 , which implies a b z 0 ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) d ξ 0 , thus inequality (92) can be written as

a b { λ 0 ( ξ ) T g ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) + z 0 ( ξ ) T h ( ξ , y 0 ( ξ ) , D a + α CF y 0 ( ξ ) ) } d ξ a b { λ 0 ( ξ ) T g ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) } d ξ ,

and since a b { λ 0 ( ξ ) T g ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) } d ξ is SPIX at v 0 , so we have

(93) a b { η ( ξ , y 0 , v 0 ) [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) ] d ξ + ( D a + α CF η ( ξ , y 0 , v 0 ) ) [ λ 0 ( ξ ) T g D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) + z 0 ( ξ ) T h D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) ] } d ξ < 0 .

Now using Proposition 1 in inequality (93) and η ( ξ , y 0 , y 0 ) = 0 , we obtain

(94) a b η ( ξ , y 0 , v 0 ) [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) ] d ξ + a b η ( ξ , y 0 , v 0 ) D b α C F R [ λ 0 ( ξ ) T g D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) + z 0 ( ξ ) T h D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) ] d ξ < 0 .

Now using Eq. (91) in inequality (94), we have

(95) a b η ( ξ , y 0 , v 0 ) [ λ 0 ( ξ ) T g v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) + z 0 ( ξ ) T h v 0 ( ξ , v 0 ( ξ ) , D a + α CF v 0 ( ξ ) ) ] d ξ + a b η ( ξ , y 0 , v 0 ) D b α CFR [ λ 0 T g D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) + z 0 ( ξ ) T h D a + α CF v 0 ( ξ ) ( ξ , v 0 , D a + α CF v 0 ) ] d ξ = 0 < 0 ,

which is a contradiction. Thus, the assumption y 0 v 0 is false. Consequently, y 0 = v 0 and v 0 is a weak minimum for (P).□

5 Conclusions

The proposed work extends and derives the generalized invexity and duality results for multiobjective variational problems with the framing of a non-singular fractional derivative pertaining to the exponential kernel by utilizing the concept of weak minima. This work considers the CF fractional derivative operator possessing a non-singular exponential kernel. Moreover, several duality results of weak, strong, and converse categories have also been derived for various types of generalized invexity conditions in view of the CF fractional derivative operator. Some basic theorems and formulas for integration by parts for the fractional derivative with an exponential kernel have played a significant role in proving the weak, strong, and converse duality theorems. This article also presents the derivation of strict converse duality theorems for multiobjective variational problems with the CF fractional derivative by employing some propositions and theorems of fractional calculus. The variational problems with CF fractional derivative may be helpful in analysing the optimization problems and physical processes. Problems related to production planning, oil refinery scheduling, portfolio selection, management sciences, and economics can be modelled successfully in the form of multiobjective variational problems. The results derived in this article are important for the growth of generalized invexity and duality results for a class of multiobjective variational problems involving a non-singular fractional derivative. However, it is difficult to illustrate any practical application on the basis of derived results. But there is an ample scope to explore the optimality and duality results for multiobjective variational problems within the scope of fractional calculus. This work can be further extended to study the multiobjective variational problems and non-differentiable multiobjective variational problems involving other kinds of fractional derivatives. As a future scope of the work, FVPs can also be studied with the theorems and propositions applied in this article.


tel: +91-9460905223

Acknowledgments

The authors are grateful to Prof. S. K. Mishra, Department of Mathematics, Institute of Science, Banaras Hindu University, Varanasi-221005, India, for fruitful ideas and valuable discussions at initial stages of their work on this article.

  1. Funding information: The authors state no funding involved.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

  4. Data availability statement: Data sharing is not applicable to this article as no datasets were generated or analysed during the current study.

References

[1] Sheybak M, Tajadodi H. Numerical solutions of fractional chemical kinetics system. Nonlinear Dyn Syst Theory. 2019;19:200–8.Search in Google Scholar

[2] Khader M. On the numerical solutions to nonlinear biochemical reaction model using Picard-Padé technique. World J Model Simul. 2013;9(1):38–46.Search in Google Scholar

[3] Malik SA, Qureshi IM, Amir M, Haq I. Numerical solution to nonlinear biochemical reaction model using hybrid polynomial basis differential evolution technique. Adv Stud Biol. 2014;6(3):99–113.10.12988/asb.2014.4520Search in Google Scholar

[4] Rossikhin YA, Shiticova MV. Application of fractional calculus for dynamic problems of solid mechanics: Novel trend and recent results. Appl Mech Rev. 2010;63(1):010801. 10.1115/1.4000563.Search in Google Scholar

[5] Varieschi GU. Applications of fractional calculus to Newtonian mechanics. J Appl Math Phys. 2018;6(6):1–11. 10.4236/jamp.2018.66105.Search in Google Scholar

[6] Kumar S, Yildirim A, Khan Y, Jafari H, Sayevand K, Wei L. Analytical solution of fractional Black-Scholes European option pricing equation by using Laplace transform. J Frac Calc Appl. 2012;2(8):1–9.Search in Google Scholar

[7] Singh J, Kumar D, Sushila, Gupta S. Application of homotopy perturbation transform method to linear and nonlinear space-time fractional reaction diffusion equations. J Math Comput Sci. 2012;5(1):40–52. 10.22436/jmcs.05.01.05.Search in Google Scholar

[8] Yildirim A. Application of the homotopy perturbation method for the Fokker-Planck equation. Int J Numer Method Biomed Eng. 2010;26(9):1144–54. 10.1002/cnm.1200.Search in Google Scholar

[9] Zhao D, Singh J, Kumar D, Rathore S, Yang X-J. An efficient computational technique for local fractional heat conduction equations in fractal media. J Nonlinear Sci Appl. 2017;10:1478–86. 10.22436/jnsa.010.04.17.Search in Google Scholar

[10] Atangana A, Botha JF. Generalized groundwater flow equation using the concept of variable order derivative. Bound Value Probl. 2013;2013:1–11. 10.1186/1687-2770-2013-53.Search in Google Scholar

[11] Singh J, Kumar D, Hammouch Z, Atangana A. A fractional epidemiological model for computer viruses pertaining to a new fractional derivative. Appl Math Comput. 2018;316:504–15. 10.1016/j.amc.2017.08.048.Search in Google Scholar

[12] Zaslavsky GM. Hamiltonian chaos and fractional dynamics. Oxford: Oxford University Press; 2004.Search in Google Scholar

[13] Atangana A, Secer A. The time-fractional coupled-the Korteweg-de Vries equations. Abstr Appl Anal. 2013;2013:1–8. 10.1155/2013/947986.Search in Google Scholar

[14] Ahmed E, El-Sayed AMA, El-Saka HAA. Equilibrium points, stability and numerical solutions of fractional order predator-prey and rabies models. J Math Anal Appl. 2007;325(1):542–53. 10.1016/j.jmaa.2006.01.087.Search in Google Scholar

[15] Chaudhary A, Kumar D, Singh J. A fractional model of fluid flow through porous media with mean capillary pressure. J Assoc Arab Univ Basic Appl Sci. 2016;21:59–63. 10.1016/j.jaubas.2015.01.002.Search in Google Scholar

[16] Samko SG, Kilbas AA, Maritchev OI. Integrals and derivatives of fractional order and some of their applications. Nauka i Tekhnika, Minsk, (1987) (Russian); English transl.: Fractional integrals and derivatives: Theory and application. London: Gordon and Breach; 1993.Search in Google Scholar

[17] Podlubny I. Geometric and physical interpretation of fractional integration and fractional differentiation. Frac Calc Appl Anal. 2002;5(4):367–86.Search in Google Scholar

[18] Caputo M. Linear model of dissipation whose Q is almost frequency independent-II. Geophys J Astronom Soc. 1967;13:529–39.10.1111/j.1365-246X.1967.tb02303.xSearch in Google Scholar

[19] Atangana A, Kilicman A. Analytical solutions of the space time-fractional derivative of advection dispersion equation. Math Probl Eng. 2013;5:1–9. 10.1155/2013/853127.Search in Google Scholar

[20] Miller KS. Fractional calculus and its applications. Conference Proceedings of the International Conference held at the University of New Haven, June 1974. Part of the Book Series: Lecture Notes in Mathematics. Vol. 457; 1975. p. 80–9.Search in Google Scholar

[21] Jumarie G. On the solution of the stochastic differential equation of exponential growth driven by fractional Brownian motion. Appl Math Lett. 2005;18(7):817–26. 10.1016/j.aml.2004.09.012.Search in Google Scholar

[22] Jumarie G. Modified Riemann-Liouville derivative and fractional Taylor series of non-differentiable functions further results. Comput Math Appl. 2006;51(9–10):1367–76. 10.1016/j.camwa.2006.02.001.Search in Google Scholar

[23] Atangana A, Baleanu D. New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model. Therm Sci. 2016;20(2):763–69. 10.2298/TSCI160111018A.Search in Google Scholar

[24] Zuñiga-Aguilar CJ, Gómez-Aguilar JF, Escobar-Jiménez RF, Romero-Ugalde HM. Robust control for factional variable-order chaotic systems with non-singular kernel. Eur Phys J Plus. 2018;133:1–13. 10.1140/epjp/i2018-11853-y.Search in Google Scholar

[25] Caputo M, Fabrizio M. A new definition of fractional derivative without singular kernel. Prog Fract Differ Appl. 2015;1(2):73–85. 10.12785/pfda/010201.Search in Google Scholar

[26] Jajarmi A, Baleanu D, Vahid KZ, Mobayen S. A general fractional formulation and tracking control for immunogenic tumor dynamics. Math Methods Appl Sci. 2022;45(2):667–80. 10.1002/mma.780.4.Search in Google Scholar

[27] Baleanu D, Abadi MH, Jajarmi A, Vahid KZ, Nieto JJ. A new comparative study on the general fractional model of COVID-19 with isolation and quarantine effects. Alex Eng J. 2022;61(6):4779–91. 10.1016/j.aej.2021.10.030.Search in Google Scholar

[28] Erturk VS, Godwe E, Baleanu D, Kumar P, Asad J, Jajarmi A. Novel fractional-order Lagrangian to describe motion of beam on nanowire. Acta Physica Polonica A. 2021;140(3):265–72. 10.12693/APhysPolA.140.265.Search in Google Scholar

[29] Jajarmi A, Baleanu D, Vahid KZ, Pirouz HM, Asad JH. A new and general fractional Lagrangian approach: A capacitor microphone case study. Results Phys. 2021;31:104950. 10.1016/j.rinp.2021.104950.Search in Google Scholar

[30] Dubey VP, Singh J, Alshehri AM, Dubey S, Kumar D. Numerical investigation of fractional model of Phytoplankton–Toxic Phytoplankton–Zooplankton system with convergence analysis. Int J Biomath. 2022;15(4):2250006. 10.1142/S1793524522500061.Search in Google Scholar

[31] Dubey VP, Dubey S, Kumar D, Singh J. A computational study of fractional model of atmospheric dynamics of carbon dioxide gas. Chaos, Solitons & Fractals. 2021;142:110375. 10.1016/j.chaos.2020.110375.Search in Google Scholar

[32] Dubey VP, Kumar D, Dubey S. A modified computational scheme and convergence analysis for fractional order hepatitis E virus model. Advanced numerical methods for differential equations: applications in science and engineering. USA: CRC Press (Taylor & Francis Group); 2021. p. 279–312.10.1201/9781003097938-11Search in Google Scholar

[33] Hanson MA. Bounds for functionally convex optimal control problems. J Math Anal Appl. 1964;8(1):84–9. 10.1016/0022-247X(64)90086-1.Search in Google Scholar

[34] Mond B, Hanson MA. Duality for variational problems. J Math Anal Appl. 1967;18(2):355–64. 10.1016/0022-247X(67)90063-7.Search in Google Scholar

[35] Chandra S, Craven BD, Husain I. A class of nondifferentiable continuous programming problems. J Math Anal Appl. 1985;107(1):122–31. 10.1016/0022-247X(85)90357-9.Search in Google Scholar

[36] Bector CR, Husain I. Duality for multiobjective variational problems. J Math Anal Appl. 1992;166(1):214–29. 10.1016/0022-247X(92)90337-D.Search in Google Scholar

[37] Nahak C, Nanda S. Duality for multiobjective variational problems with invexity. Optimization. 1996;36(3):235–48.10.1080/02331939608844181Search in Google Scholar

[38] Chen XH. Duality for multiobjective variational problems with invexity. J Math Anal Appl. 1996;203(1):236–53. 10.1006/jmaa.1996.0377.Search in Google Scholar

[39] Bhatia D, Mehra A. Optimality conditions and duality for multiobjective variational problems with generalized B-invexity. J Math Anal Appl. 1999;234(2):341–60. 10.1006/jmaa.1998.6256.Search in Google Scholar

[40] Mond B, Chandra S, Husain I. Duality for variational problems with invexity. J Math Anal Appl. 1988;134(2):322–8. 10.1016/0022-247X(88)90026-1.Search in Google Scholar

[41] Hanson MA. On sufficiency of the Kuhn-Tucker conditions. J Math Anal Appl. 1981;80(2):545–50. 10.1016/0022-247X(81)90123-2.Search in Google Scholar

[42] Mishra SK, Mukherjee RN. Duality for multiobjective fractional variational problems. J Math Anal Appl. 1994;186(3):711–25. 10.1006/jmaa.1994.1328.Search in Google Scholar

[43] Mond B, Husain I. Sufficient optimality criteria and duality for variational problems with generalized invexity. J Austral Math Soc Ser B. 1989;31(1):108–21.10.1017/S0334270000006512Search in Google Scholar

[44] Weir T, Mond B. Generalized convexity and duality in multiple objective programming. Bull Austral Math Soc. 1989;39(2):287–99. 10.1017/S000497270000277X.Search in Google Scholar

[45] Mukherjee RN, Mishra SK. Generalized invexity and duality in multiple objective variational problems. J Math Anal Appl. 1995;195(2):307–22. 10.1006/jmaa.1995.1357.Search in Google Scholar

[46] Kumar PJ. Generalized invexity of higher order and its applications in variational problems. Appl Math. 2015;6(9):1638–48. 10.4236/am.2015.69146.Search in Google Scholar

[47] Kumar P, Jyoti, Sharma B. Characterization of generalized invexity in multi-objective fractional variational problem. Stat Optim Inf Comput. 2016;4:342–9. 10.19139/soic.v4i4.208. Search in Google Scholar

[48] Hachimi M, Aghezzaf B. Sufficiency and duality in multiobjective variational problems with generalized type I functions. J Glob Optim. 2006;34(2):191–218. 10.1007/s10898-005-1653-2.Search in Google Scholar

[49] Mishra SK, Wang SY, Lai KK. Symmetric duality for a class of nondifferentiable multi-objective fractional variational problems. J Math Anal Appl. 2007;333(2):1093–110. 10.1016/j.jmaa.2006.11.054.Search in Google Scholar

[50] Antczak T. Duality for multiobjective variational control problems with (ϕ,ρ)-invexity. Calcolo. 2014;51(3):393–421.10.1007/s10092-013-0092-6Search in Google Scholar

[51] Upadhyay BB, Ghosh A, Mishra P, Treanta S. Optimality conditions and duality for multiobjective semi-infinite programming problems on Hadamard manifolds using generalized geodesic convexity. RAIRO-Oper Res. 2022;56(4):2037–65. 10.1051/ro/2022098.Search in Google Scholar

[52] Upadhyay BB, Treanta S, Mishra P. On Minty variational principle for nonsmooth multiobjective optimization problems on Hadamard manifolds. Optimization. 2022. 10.1080/02331934.2022.2088369. Search in Google Scholar

[53] Guo Y, Ye G, Liu W, Zhao D, Treanta S. On symmetric gH-derivative: Applications to dual interval-valued optimization problems. Chaos Solitons & Fractals. 2022;158(2):112068. org/10.1016/j.chaos.2022.112068.Search in Google Scholar

[54] Guo Y, Ye G, Liu W, Zhao D, Treanta S. Optimality conditions and duality for a class of generalized convex interval-valued optimization problems. Mathematics. 2021;9(22):2979. 10.3390/math9222979.Search in Google Scholar

[55] Podlubny I. Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications, mathematics in science and engineering. Vol. 198. New York, London, Sydney, Tokyo and Toronto: Academic Press; 1999.Search in Google Scholar

[56] Kilbas AA, Srivastava HM, Trujillo JJ. Theory and applications of fractional differential equations. North-Holland Math Stud Book Ser. 2006;204:1–523.Search in Google Scholar

[57] Abdeljawad T, Baleanu D. On fractional derivatives with exponential kernel and their discrete versions. Rep Math Phys. 2017;80(1):11–27. 10.1016/S0034-4877(17)30059-9.Search in Google Scholar

[58] Mishra SK, Giorgi G. Invexity and optimization. Part of the book series: nonconvex optimization and its applications. Vol. 88. Berlin Heidelberg, Germany: Springer-Verlag; 2008. p. 1–266. 10.1007/978-3-540-78562-0. Search in Google Scholar

Received: 2022-06-03
Revised: 2022-07-21
Accepted: 2022-08-07
Published Online: 2022-09-26

© 2022 Ved Prakash Dubey et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 15.5.2024 from https://www.degruyter.com/document/doi/10.1515/phys-2022-0195/html
Scroll to top button