1 Introduction

The inequality of Ostrowski gives us an estimate for the deviation of the values of a function from its mean value. More precisely, let \(f:[a,b]\rightarrow \mathbb{R}\) be a differentiable function with bounded derivative. Then the following integral inequality:

$$ \biggl\vert f(x)-\frac{1}{b-a} \int _{a}^{b}f(t)\,dt \biggr\vert \leq \biggl[ \frac{1}{4}+ \frac{(x-\frac{a+b}{2})^{2}}{(b-a)^{2}} \biggr] (b-a) \bigl\Vert f^{ \prime } \bigr\Vert _{\infty } $$
(1.1)

is valid for every \(x\in {}[ a,b]\), which is proved by Ostrowski in 1938. In addition to this, the constant \(1/4\) is the best possible.

Simpson’s inequality plays significant role in many areas of mathematics. To be more precise, the classical Simpson’s inequality is expressed as follows for four times continuously differentiable functions.

Theorem 1

([12])

Suppose for a mapping \(f: [ a,b ] \rightarrow \mathbb{R}\) which is four times continuously differentiable on \(( a,b ) \), and suppose also \(\Vert f^{ ( 4 ) } \Vert _{\infty }=\sup_{x\in ( a,b ) } \vert f^{ ( 4 ) }(x) \vert <\infty \). Then one has the inequality

$$ \biggl\vert \frac{1}{3} \biggl[ \frac{f(a)+f(b)}{2}+2f \biggl( \frac{a+b}{2} \biggr) \biggr] -\frac{1}{b-a} \int _{a}^{b}f(x)\,dx \biggr\vert \leq \frac{1}{2880} \bigl\Vert f^{ ( 4 ) } \bigr\Vert _{\infty } ( b-a ) ^{4}. $$

Over the years, many variations of Ostrowski and Simpson type inequalities have been studied for various function classes, such as convex functions, bounded functions, functions of bounded variation, and so on. Specifically, since convexity theory is an effective and powerful way to solve a large number of problems from different branches of pure and applied mathematics, many papers have been dedicated to the Simpson inequality for convex functions. For instance, Alomari et al. proved some Ostrowski type inequalities for s-convex functions in [4]. Moreover, Ostrowski type inequalities were studied for different kinds of convexities such as m-convex functions [20], \((\alpha ,m)\)-convex functions [29], and h-convex functions [26, 44]. In [41], Set first obtained the Riemann–Liouville fractional version of the Ostrowski inequality for s-convex functions. In addition to this, many researchers focused on establishing Ostrowski type inequalities for certain fractional integral operators, such as k-Riemann–Liouville fractional integrals [16], local fractional integrals [34], Raina fractional integrals [2], generalized k-g-fractional integrals [10] and ψ-Hilfer fractional integrals [5]. On the other hand, several Ostrowski inequalities for co-ordinated convex mapping in involving double Riemann integrals and double Riemann–Liouville fractional integrals are introduced in [24] and [23], respectively.

Some authors established several Simpson type inequalities for differentiable and twice differentiable convex functions [3, 3638]. Ozdemir et al. proved Simpson type inequalities for co-ordinated convex functions in [28]. In [8, 19, 42], authors obtained some new Simpson inequalities for Riemann–Liouville fractional integrals. Furthermore, a number of papers were devoted to Simpson inequalities for other fractional integrals or for functions belong to other convex classes such as p-convex function [1], h-convex function [27], preinvex functions [31], \((m,h_{1},h_{2})\)-preinvexity [31], and generalized harmonic convex functions [43].

The inequalities, introduced by C. Hermite and J. Hadamard for convex functions, are of considerable significance in the literature. These inequalities state that, if \(\digamma :I\rightarrow \mathbb{R}\) is a convex function on the interval I of real numbers and \(a,b\in I\) with \(a< b\), then

$$ \digamma \biggl( \frac{a+b}{2} \biggr) \leq \frac{1}{b-a} \int _{a}^{b}\digamma (x)\,dx\leq \frac{\digamma ( a ) +\digamma ( b ) }{2}. $$
(1.2)

If Ϝ is concave, then both inequalities in (1.2) hold in the reverse direction.

Over the years, considerable number of studies have been focused on obtaining trapezoidal and midpoint type inequalities which give bounds for the right-hand side and left-hand side of the inequality (1.2), respectively. For example, Dragomir and Agarwal first obtained trapezoidal inequalities for convex functions in [11] and Kirmacı first established midpoint inequalities for convex functions in [21]. In [40], Sarikaya et al. generalized the inequalities (1.2) for fractional integrals and the authors also proved some corresponding trapezoidal type inequalities. Iqbal et al. presented some fractional midpoint type inequalities for convex functions in [18]. On the other hand, Dragomir proved Hermite–Hadamard inequalities for co-ordinated convex mappings in [9]. The midpoint and trapezoidal type inequalities for co-ordinated convex functions were established in [22] and [39], respectively. Moreover, Sarikaya obtained fractional Hermite–Hadamard inequalities and fractional trapezoidal inequalities for functions with two variables in [32]. Tunç et al. presented some fractional midpoint type inequalities for co-ordinated convex functions in [46]. In [35], Sarikaya and Ertuğral first introduced new fractional integrals which are called generalized fractional integrals. In addition, they proved Hermite–Hadamard inequalities and several trapezoidal and midpoint type inequalities for generalized fractional integrals. Moreover, Budak et al. proved Midpoint type inequalities and extensions of Hermite–Hadamard inequalities in the papers [6] and [7], respectively. In [13], Ertuğral and Sarikaya presented some Simpson type inequalities for these fractional integral operators. For some of other papers on inequalities for generalized fractional integrals, we refer to [17, 48]. On the other hand, Turkay et al. described the generalized fractional integrals for functions with two variables. These authors presented Hermite–Hadamard and trapezoidal type inequalities for this kind of fractional integrals in [47]. For the other similar inequalities, we refer to [14, 15, 25, 30, 33, 45].

The aims of this paper is to establish some generalized inequalities for co-ordinated convex functions involving generalized fractional integrals. The general structure of the paper consists of five sections including an introduction. The remaining part of the paper proceeds as follows: In Sect. 2, we give the definitions of generalized fractional integrals and relations between generalized fractional integrals and other type of fractional integrals. In Sect. 3, an identity involving some parameters are proved for partially differentiable functions. Then we establish several generalized inequalities for mappings whose partially derivatives in absolute value are co-ordinated convex. With the special choice of the given parameters, we show that our results reproduce the results proved in the earlier work and we also give some new trapezoidal and Simpson type inequalities in Sect. 4. At the end of the paper, some conclusions and further directions of research are discussed in Sect. 5.

2 Generalized fractional Integrals

In this section, we summarize the generalized fractional integrals defined by Sarikaya and Ertuğral in [35].

Let us define a function \(\varphi :[0,\infty )\rightarrow {}[ 0,\infty )\) satisfying the following condition:

$$ \int _{0}^{1}\frac{\varphi ( t ) }{t}\,dt< \infty . $$

We consider the following left-sided and right-sided generalized fractional integral operators:

$$ _{a+}I_{\varphi }f(x)= \int _{a}^{x} \frac{\varphi ( x-t ) }{x-t}f(t)\,dt, \quad x>a, $$
(2.1)

and

$$ _{b-}I_{\varphi }f(x)= \int _{x}^{b} \frac{\varphi ( t-x ) }{t-x}f(t)\,dt, \quad x< b, $$
(2.2)

respectively.

Some forms of fractional integrals, namely Riemann–Liouville fractional integrals, k-Riemann–Liouville fractional integrals, Katugampola fractional integrals, conformable fractional integrals, and Hadamard fractional integrals, are the most significant features of generalized fractional integrals. These important special cases of the integral operators (2.1) and (2.2) are mentioned below:

Remark 1

If we choose \(\varphi ( t ) =t\), the operators (2.1) and (2.2) reduce to the Riemann integral.

Remark 2

Considering \(\varphi ( t ) = \frac{t^{\alpha }}{\Gamma ( \alpha ) } \) and \(\alpha >0\), the operators (2.1) and (2.2) reduce to the Riemann–Liouville fractional integrals \(J_{{a+}}^{\alpha }f(x) \) and \(J_{b-}^{\alpha }f(x)\), respectively. Here, Γ is Gamma function.

Remark 3

For \(\varphi ( t ) = \frac{1}{k\Gamma _{k} ( \alpha ) }t^{\frac{\alpha }{k}}\) and \(\alpha ,k>0\), the operators (2.1) and (2.2) reduce to the k-Riemann–Liouville fractional integrals \(J_{{a+,k}}^{\alpha }f(x)\) and \(J_{b-,k}^{\alpha }f(x)\), respectively. Here, \(\Gamma _{k}\) is the k-Gamma function.

Generalized double fractional integrals are given by Turkay et al. in [47], as follows.

Definition 1

The generalized double fractional integrals \(_{a+,c+}I_{\varphi ,\psi }\), \(_{a+,d-}I_{\varphi ,\psi }\), \(_{b-,c+}I_{\varphi ,\psi }\), \(_{b-,d-}I_{\varphi ,\psi }\) are defined by

$$\begin{aligned}& _{a+,c+}I_{\varphi ,\psi }f ( x,y ) = \int _{a}^{x} \int _{c}^{y} \frac{\varphi ( x-t ) }{x-t} \frac{\psi ( y-s ) }{y-s}f ( t,s )\,ds\,dt,\quad x>a, y>c, \end{aligned}$$
(2.3)
$$\begin{aligned}& _{a+,d-}I_{\varphi ,\psi }f ( x,y ) = \int _{a}^{x} \int _{y}^{d} \frac{\varphi ( x-t ) }{x-t} \frac{\psi ( s-y ) }{s-y}f ( t,s )\,ds\,dt,\quad x>a, y< d, \end{aligned}$$
(2.4)
$$\begin{aligned}& _{b-,c+}I_{\varphi ,\psi }f ( x,y ) = \int _{x}^{b} \int _{c}^{y} \frac{\varphi ( t-x ) }{t-x} \frac{\psi ( y-s ) }{y-s}f ( t,s )\,ds\,dt,\quad x< b, y>c, \end{aligned}$$
(2.5)

and

$$ _{b-,d-}I_{\varphi ,\psi }f ( x,y ) = \int _{x}^{b} \int _{y}^{d} \frac{\varphi ( t-x ) }{t-x} \frac{\psi ( s-y ) }{s-y}f ( t,s )\,ds\,dt, \quad x< b, y< d. $$
(2.6)

Here, \(f \in L_{1}([a,b]\times {}[ c,d])\) and the functions \(\varphi ,\psi : [ 0,\infty ) \rightarrow [ 0, \infty ) \) satisfy, respectively, \(\int _{0}^{1}\frac{\varphi ( t ) }{t}\,dt<\infty \) and \(\int _{0}^{1}\frac{\psi ( s ) }{s}\,ds<\infty \).

By using Definition 1, well-known fractional integrals can be obtained by some special choices. For example;

  1. (1)

    If we choose \(\varphi ( t ) =t\) and \(\psi ( s ) =s \), the operators (2.3), (2.4), (2.5) and (2.6) reduce to the double Riemann integral.

  2. (2)

    Considering \(\varphi ( t ) =\frac{t^{\alpha }}{\Gamma (\alpha )}\), \(\psi ( s ) =\frac{s^{\beta }}{\Gamma (\beta )}\), then, for \(\alpha ,\beta >0\), the operators (2.3), (2.4), (2.5) and (2.6) reduce to the Riemann–Liouville fractional integrals \(J_{a+,c+}^{\alpha ,\beta }f ( x,y ) \), \(J_{a+,d-}^{\alpha ,\beta }f ( x,y ) \) \(J_{b-,c+}^{\alpha ,\beta }f ( x,y ) \) and \(J_{b-,d-}^{\alpha ,\beta }f ( x,y ) \), respectively.

  3. (3)

    For \(\varphi ( t ) = \frac{t^{\frac{\alpha }{k}}}{k\Gamma _{k}(\alpha )}\) and \(\psi ( s ) = \frac{s^{\frac{\beta }{k}}}{k\Gamma _{k}(\beta )}\), for \(\alpha ,\beta ,k>0\), the operators (2.3), (2.4), (2.5) and (2.6) reduce to the k-Riemann–Liouville fractional integrals \(J_{a+,c+}^{\alpha ,\beta ,k}f ( x,y ) \), \(J_{a+,d-}^{\alpha ,\beta ,k}f ( x,y ) \), \(J_{b-,c+}^{\alpha ,\beta ,k}f ( x,y ) \) and \(J_{b-,d-}^{\alpha ,\beta ,k}f ( x,y ) \), respectively.

For our work to be done in two-dimensional space, we will need the following descriptions.

3 An identity for Generalized double fractional integrals

Throughout this study, we assume that \(\lambda _{i},\mu _{i}\geq 0\), \(i=1,2\) and for brevity, we define

$$ \Lambda _{1}(x,t)= \int _{0}^{t} \frac{\varphi ( ( b-x ) u ) }{u}\,du,\qquad \Delta _{1}(x,t)= \int _{0}^{t} \frac{\varphi ( ( x-a ) u ) }{u}\,du, $$
(3.1)

and

$$ \Lambda _{2}(y,s)= \int _{0}^{s} \frac{\psi ( ( d-y ) u ) }{u}\,du, \qquad \Delta _{2}(y,s)= \int _{0}^{s} \frac{\psi ( ( y-c ) u ) }{u}\,du. $$
(3.2)

Lemma 1

Let \(f:\Delta :=[a,b]\times {}[ c,d]\rightarrow \mathbb{R} \) be an absolutely continuous function on Δ such that the partial derivative of order \(\frac{\partial ^{2}f ( t,s ) }{\partial t\partial s}\) exist for all \(( t,s ) \in \Delta \). Then the following equality for generalized fractional integrals holds:

$$\begin{aligned}& \Omega (a,b,x;c,d,y) \\& \quad = ( b-x ) ( d-y ) \int _{0}^{1} \int _{0}^{1} \bigl( \lambda _{1}\Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr) \bigl( \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr) \\& \qquad {} \times \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s )\,d \bigr)\,ds\,dt \\& \qquad {} - ( b-x ) ( y-c ) \int _{0}^{1} \int _{0}^{1} \bigl( \lambda _{1}\Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr) \bigl( \mu _{2} \Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr) \\& \qquad {} \times \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s ) c \bigr)\,ds\,dt \\& \qquad {} - ( x-a ) ( d-y ) \int _{0}^{1} \int _{0}^{1} \bigl( \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr) \bigl( \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr) \\& \qquad {}\times \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s )\,d \bigr)\,ds\,dt \\& \qquad {} + ( x-a ) ( y-c ) \int _{0}^{1} \int _{0}^{1} \bigl( \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr) \bigl( \mu _{2} \Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr) \\& \qquad {}\times \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s ) c \bigr)\,ds\,dt, \end{aligned}$$

where

$$\begin{aligned}& \Omega (a,b,x;c,d,y) \\& \quad = \Lambda _{1}(x,1)\Lambda _{2}(y,1) \bigl[ ( 1-\lambda _{1} ) ( 1-\lambda _{2} ) f ( x,y ) + ( 1- \lambda _{1} ) \lambda _{2}f ( x,d ) \\& \qquad {} +\lambda _{1} ( 1- \lambda _{2} ) f ( b,y ) +\lambda _{1} \lambda _{2}f ( b,d ) \bigr] \\& \qquad {}+\Lambda _{1}(x,1)\Delta _{2}(y,1) \bigl[ ( 1-\lambda _{1} ) ( 1-\mu _{2} ) f ( x,y ) + ( 1- \lambda _{1} ) \mu _{2}f ( x,c ) \\& \qquad {}+\lambda _{1} ( 1-\mu _{2} ) f ( b,y ) +\lambda _{1}\mu _{2}f ( b,c ) \bigr] \\& \qquad {} +\Delta _{1}(x,1)\Lambda _{2}(y,1) \bigl[ ( 1-\mu _{1} ) ( 1-\lambda _{2} ) f ( x,y ) + ( 1-\mu _{1} ) \lambda _{2}f ( x,d ) \\& \qquad {}+\mu _{1} ( 1-\lambda _{2} ) f ( a,y ) +\mu _{1}\lambda _{2}f ( a,d ) \bigr] \\& \qquad {} +\Delta _{1}(x,1)\Delta _{2}(y,1) \bigl[ ( 1-\mu _{1} ) ( 1-\mu _{2} ) f ( x,y ) + ( 1-\mu _{1} ) \mu _{2}f ( x,c ) \\& \qquad {} +\mu _{1} ( 1-\mu _{2} ) f ( a,y ) +\mu _{1}\mu _{2}f ( a,c ) \bigr] \\& \qquad {} - \bigl( ( 1-\lambda _{2} ) \Lambda _{2}(y,1)+ ( 1- \mu _{2} ) \Delta _{2}(y,1) \bigr) \bigl[ _{x+}I_{\varphi }f ( b,y ) + _{x-}I_{\varphi }f ( a,y ) \bigr] \\& \qquad {} - \bigl( ( 1-\lambda _{1} ) \Lambda _{1}(x,1)+ ( 1- \mu _{1} ) \Delta _{1}(x,1) \bigr) \bigl[ _{y+}I_{\psi }f ( x,d ) + _{y-}I_{\psi }f ( x,c ) \bigr] \\& \qquad {} -\lambda _{2}\Lambda _{2}(y,1) \bigl[ _{x+}I_{\varphi }f ( b,d ) +_{x-}I_{\varphi }f ( a,d ) \bigr] -\mu _{2} \Delta _{2}(y,1) \bigl[ _{x+}I_{\varphi }f ( b,c ) +_{x-}I_{ \varphi }f ( a,c ) \bigr] \\& \qquad {} -\lambda _{1}\Lambda _{1}(x,1) \bigl[ _{y+}I_{\psi }f ( b,d ) + _{y-}I_{\psi }f ( b,c ) \bigr] -\mu _{1} \Delta _{1}(x,1) \bigl[ _{y+}I_{\psi }f ( a,d ) +_{y-}I_{\psi }f ( a,c ) \bigr] \\& \qquad {} +_{x+,y+}I_{\varphi ,\psi }f ( b,d ) +_{x+,y-}I_{ \varphi ,\psi }f ( b,c ) +_{x-,y+}I_{\varphi ,\psi }f ( a,d ) +_{x-,y-}I_{\varphi ,\psi }f ( a,c ) . \end{aligned}$$

Proof

By using integration by parts, we have

$$\begin{aligned}& H_{1} = \int _{0}^{1} \int _{0}^{1} \bigl( \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr) \bigl( \lambda _{2} \Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr) \frac{\partial ^{2}}{\partial t\partial s} \\& \hphantom{ H_{1}=}{}\times f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s )\,d \bigr)\,ds\,dt \\& \hphantom{H_{1}} = \frac{\Lambda _{1}(x,1)\Lambda _{2}(y,1)}{(b-x)(d-y)} \bigl[ ( 1-\lambda _{1} ) ( 1-\lambda _{2} ) f(x,y)+ ( 1-\lambda _{1} ) \lambda _{2}f(x,d) \\& \hphantom{ H_{1}=}{}+\lambda _{1} ( 1- \lambda _{2} ) f(b,y)+ \lambda _{1}\lambda _{2}f(b,d) \bigr] \\& \hphantom{H_{1} =} {}- \frac{ ( 1-\lambda _{2} ) \Lambda _{2}(y,1)}{(b-x)(d-y)} _{x+}I_{\varphi }f ( b,y ) - \frac{\lambda _{2}\Lambda _{2}(y,1)}{(b-x)(d-y)} _{x+}I_{\varphi }f ( b,d ) \end{aligned}$$
(3.3)
$$\begin{aligned}& \hphantom{H_{1} =}{}- \frac{ ( 1-\lambda _{1} ) \Lambda _{1}(x,1)}{(b-x)(d-y)} _{y+}I_{\psi }f ( x,d ) - \frac{\lambda _{1}\Lambda _{1}(x,1)}{(b-x)(d-y)} _{y+}I_{\psi }f ( b,d ) \\& \hphantom{ H_{1}=}{}+\frac{1}{(b-x)(d-y)} _{x+,y+}I_{\varphi ,\psi }f ( b,d ) , \\& H_{2}= \int _{0}^{1} \int _{0}^{1} \bigl( \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr) \bigl( \mu _{2} \Delta _{2}(y,1)- \Delta _{2}(y,s) \bigr) \frac{\partial ^{2}}{\partial t\partial s} \\& \hphantom{ H_{2}=}{}\times f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s ) c \bigr)\,ds\,dt \\& \hphantom{H_{2}}= -\frac{\Lambda _{1}(x,1)\Lambda _{2}(y,1)}{(b-x)(y-c)} \bigl[ ( 1-\lambda _{1} ) ( 1-\mu _{2} ) f(x,y)+ ( 1-\lambda _{1} ) \mu _{2}f(x,c) \\& \hphantom{ H_{2}=}{}+ \lambda _{1} ( 1- \mu _{2} ) f(b,y)+\lambda _{1}\mu _{2}f(b,c) \bigr] \\& \hphantom{H_{2}=} {}+\frac{ ( 1-\mu _{2} ) \Delta _{2}(y,1)}{(b-x)(y-c)} _{x+}I_{\varphi }f ( b,y ) + \frac{\mu _{2}\Delta _{2}(y,1)}{(b-x)(y-c)} _{x+}I_{\varphi }f ( b,c ) \end{aligned}$$
(3.4)
$$\begin{aligned}& \hphantom{H_{2}=} {}+ \frac{ ( 1-\lambda _{1} ) \Lambda _{1}(x,1)}{(b-x)(y-c)} _{y-}I_{\psi }f ( x,c ) + \frac{\lambda _{1}\Lambda _{1}(x,1)}{(b-x)(y-c)} _{y-}I_{\psi }f ( b,c ) \\& \hphantom{ H_{2}=}{} -\frac{1}{(b-x)(y-c)} _{x+,y-}I_{\varphi ,\psi }f ( b,c ) , \\& H_{3}= \int _{0}^{1} \int _{0}^{1} \bigl( \mu _{1} \Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr) \bigl( \lambda _{2} \Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr) \frac{\partial ^{2}}{\partial t\partial s} \\& \hphantom{ H_{3}=}{}\times f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s )\,d \bigr)\,ds\,dt \\& \hphantom{H_{3}}= -\frac{\Delta _{1}(x,1)\Lambda _{2}(y,1)}{(x-a)(d-y)} \bigl[ ( 1-\mu _{1} ) ( 1-\lambda _{2} ) f(x,y)+ ( 1-\mu _{1} ) \lambda _{2}f(x,d) \\& \hphantom{ H_{3}=}{}+ \mu _{1} ( 1- \lambda _{2} ) f(a,y)+\mu _{1}\lambda _{2}f(a,d) \bigr] \\& \hphantom{ H_{3}=} {}+ \frac{ ( 1-\lambda _{2} ) \Lambda _{2}(y,1)}{(x-a)(d-y)} _{x-}I_{\varphi }f ( a,y ) + \frac{\lambda _{2}\Lambda _{2}(y,1)}{(x-a)(d-y)} _{x-}I_{\varphi }f ( a,d ) \\& \hphantom{ H_{3}=} {}+\frac{ ( 1-\mu _{1} ) \Delta _{1}(x,1)}{(x-a)(d-y)} _{y+}I_{\psi }f ( x,d ) + \frac{\mu _{1}\Delta _{1}(x,1)}{(x-a)(d-y)} _{y+}I_{\psi }f ( a,d ) \\& \hphantom{ H_{3}=}{} - \frac{1}{(x-a)(d-y)} _{x-,y+}I_{ \varphi ,\psi }f ( a,d ) , \end{aligned}$$
(3.5)

and

$$\begin{aligned} H_{4} =& \int _{0}^{1} \int _{0}^{1} \bigl( \mu _{1} \Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr) \bigl( \mu _{2} \Delta _{2}(y,1)- \Delta _{2}(y,s) \bigr) \frac{\partial ^{2}}{\partial t\partial s} \\ &{}\times f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s ) c \bigr)\,ds\,dt \\ =&\frac{\Delta _{1}(x,1)\Delta _{2}(y,1)}{(x-a)(y-c)} \bigl[ ( 1- \mu _{1} ) ( 1-\mu _{2} ) f(x,y)+ ( 1-\mu _{1} ) \mu _{2}f(x,c) \\ &{}+\mu _{1} ( 1- \mu _{2} ) f(a,y)+ \mu _{1}\mu _{2}f(a,c) \bigr] \\ &{}-\frac{ ( 1-\mu _{2} ) \Delta _{2}(y,1)}{(x-a)(y-c)} _{x-}I_{\varphi }f ( a,y ) - \frac{\mu _{2}\Delta _{2}(y,1)}{(x-a)(y-c)} _{x-}I_{\varphi }f ( a,c ) \\ &{}-\frac{ ( 1-\mu _{1} ) \Delta _{1}(x,1)}{(x-a)(y-c)} _{y-}I_{\psi }f ( x,c ) - \frac{\mu _{1}\Delta _{1}(x,1)}{(x-a)(y-c)} _{y-}I_{\psi }f ( a,c ) \\ &{} + \frac{1}{(x-a)(y-c)} _{x-,y-}I_{ \varphi ,\psi }f ( a,c ) . \end{aligned}$$
(3.6)

By using Eqs. (3.3)–(3.6), we have

$$\begin{aligned}& ( b-x ) ( d-y ) H_{1}- ( b-x ) ( y-c ) H_{2}- ( x-a ) ( d-y ) H_{3}+ ( x-a ) ( y-c ) H_{4} \\& \quad =\Omega (a,b,x;c,d,y), \end{aligned}$$

which completes the proof of Lemma 1. □

Corollary 1

In Lemma 1, if we choose \(\varphi (t)=t\) and \(\psi (s)=s\) for all \(( t,s ) \in \triangle \), then we obtain the equality for the Riemann integral

$$\begin{aligned}& \aleph (a,b,x;c,d,y) \\& \quad = ( b-x ) ^{2} ( d-y ) ^{2} \int _{0}^{1} \int _{0}^{1} ( \lambda _{1}-t ) ( \lambda _{2}-s ) \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s )\,d \bigr)\,ds\,dt \\& \qquad {}- ( b-x ) ^{2} ( y-c ) ^{2} \int _{0}^{1} \int _{0}^{1} ( \lambda _{1}-t ) ( \mu _{2}-s ) \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s ) c \bigr) \,ds\,dt \\& \qquad {}- ( x-a ) ^{2} ( d-y ) ^{2} \int _{0}^{1} \int _{0}^{1} ( \mu _{1}-t ) ( \lambda _{2}-s ) \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s )\,d \bigr) \,ds\,dt \\& \qquad {}+ ( x-a ) ^{2} ( y-c ) ^{2} \int _{0}^{1} \int _{0}^{1} ( \mu _{1}-t ) ( \mu _{2}-s ) \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s ) c \bigr) \,ds\,dt. \end{aligned}$$

Here,

$$\begin{aligned}& \aleph (a,b,x;c,d,y) \\& \quad = ( b-x ) ( d-y ) \bigl[ ( 1-\lambda _{1} ) ( 1-\lambda _{2} ) f ( x,y ) + ( 1- \lambda _{1} ) \lambda _{2}f ( x,d ) \\& \qquad {} +\lambda _{1} ( 1-\lambda _{2} ) f ( b,y ) +\lambda _{1} \lambda _{2}f ( b,d ) \bigr] \\& \qquad {}+ ( b-x ) ( y-c ) \bigl[ ( 1-\lambda _{1} ) ( 1-\mu _{2} ) f ( x,y ) + ( 1- \lambda _{1} ) \mu _{2}f ( x,c ) \\& \qquad {} +\lambda _{1} ( 1-\mu _{2} ) f ( b,y ) +\lambda _{1}\mu _{2}f ( b,c ) \bigr] \\& \qquad {}+ ( x-a ) ( d-y ) \bigl[ ( 1-\mu _{1} ) ( 1-\lambda _{2} ) f ( x,y ) + ( 1- \mu _{1} ) \lambda _{2}f ( x,d ) \\& \qquad {} +\mu _{1} ( 1- \lambda _{2} ) f ( a,y ) +\mu _{1}\lambda _{2}f ( a,d ) \bigr] \\& \qquad {}+ ( x-a ) ( y-c ) \bigl[ ( 1-\mu _{1} ) ( 1-\mu _{2} ) f ( x,y ) + ( 1- \mu _{1} ) \mu _{2}f ( x,c ) \\& \qquad {}+\mu _{1} ( 1-\mu _{2} ) f ( a,y ) +\mu _{1}\mu _{2}f ( a,c ) \bigr] \\& \qquad {}- \bigl( ( 1-\lambda _{2} ) ( d-y ) + ( 1-\mu _{2} ) ( y-c ) \bigr) \int _{a}^{b}f ( t,y )\,dt- \bigl( ( 1-\lambda _{1} ) ( b-x ) \\& \qquad {} + ( 1-\mu _{1} ) ( x-a ) \bigr) \int _{c}^{d}f ( x,s )\,ds \\& \qquad {}-\lambda _{2} ( d-y ) \int _{a}^{b}f ( t,d )\,dt-\mu _{2} ( y-c ) \int _{a}^{b}f ( t,c )\,dt \\& \qquad {}-\lambda _{1} ( b-x ) \int _{c}^{d}f ( b,s )\,dt-\mu _{1} ( x-a ) \int _{c}^{d}f ( a,s )\,ds+ \int _{a}^{b} \int _{c}^{d}f ( t,s )\,ds\,dt. \end{aligned}$$

Corollary 2

In Lemma 1, let us consider \(\varphi (t)=\frac{t^{\alpha }}{\Gamma (\alpha )}\) and \(\psi (s)=\frac{s^{\beta }}{\Gamma (\beta )}\) for all \(( t,s ) \in \triangle \). Then we obtain the equality for Riemann–Liouville fractional integrals

$$\begin{aligned}& \Phi (a,b,x;c,d,y) \\& \quad = ( b-x ) ^{\alpha +1} ( d-y ) ^{\beta +1} \int _{0}^{1} \int _{0}^{1} \bigl( \lambda _{1}-t^{ \alpha } \bigr) \bigl( \lambda _{2}-s^{\beta } \bigr) \frac{\partial ^{2}}{\partial t\partial s} \\& \qquad {}\times f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s )\,d \bigr)\,ds\,dt \\& \qquad {}- ( b-x ) ^{\alpha +1} ( y-c ) ^{\beta +1} \int _{0}^{1} \int _{0}^{1} \bigl( \lambda _{1}-t^{ \alpha } \bigr) \bigl( \mu _{2}-s^{\beta } \bigr) \frac{\partial ^{2}}{\partial t\partial s} \\& \qquad {}\times f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s ) c \bigr)\,ds\,dt \\& \qquad {}- ( x-a ) ^{\alpha +1} ( d-y ) ^{\beta +1} \int _{0}^{1} \int _{0}^{1} \bigl( \mu _{1}-t^{\alpha } \bigr) \bigl( \lambda _{2}-s^{\beta } \bigr) \frac{\partial ^{2}}{\partial t\partial s} \\& \qquad {}\times f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s )\,d \bigr)\,ds\,dt \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( y-c ) ^{\beta +1} \int _{0}^{1} \int _{0}^{1} \bigl( \mu _{1}-t^{\alpha } \bigr) \bigl( \mu _{2}-s^{\beta } \bigr) \frac{\partial ^{2}}{\partial t\partial s}\\& \qquad {}\times f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s ) c \bigr)\,ds\,dt, \end{aligned}$$

where

$$\begin{aligned}& \Phi (a,b,x;c,d,y) \\& \quad = ( b-x ) ^{\alpha } ( d-y ) ^{\beta } \bigl[ ( 1-\lambda _{1} ) ( 1-\lambda _{2} ) f ( x,y ) + ( 1-\lambda _{1} ) \lambda _{2}f ( x,d ) \\& \qquad {} +\lambda _{1} ( 1- \lambda _{2} ) f ( b,y ) +\lambda _{1}\lambda _{2}f ( b,d ) \bigr] \\& \qquad {}+ ( b-x ) ^{\alpha } ( y-c ) ^{\beta } \bigl[ ( 1-\lambda _{1} ) ( 1-\mu _{2} ) f ( x,y ) + ( 1-\lambda _{1} ) \mu _{2}f ( x,c ) \\& \qquad {}+ \lambda _{1} ( 1-\mu _{2} ) f ( b,y ) +\lambda _{1} \mu _{2}f ( b,c ) \bigr] \\& \qquad {}+ ( x-a ) ^{\alpha } ( d-y ) ^{\beta } \bigl[ ( 1-\mu _{1} ) ( 1-\lambda _{2} ) f ( x,y ) + ( 1-\mu _{1} ) \lambda _{2}f ( x,d ) \\& \qquad {} + \mu _{1} ( 1-\lambda _{2} ) f ( a,y ) +\mu _{1} \lambda _{2}f ( a,d ) \bigr] \\& \qquad {}+ ( x-a ) ^{\alpha } ( y-c ) ^{\beta } \bigl[ ( 1-\mu _{1} ) ( 1-\mu _{2} ) f ( x,y ) + ( 1-\mu _{1} ) \mu _{2}f ( x,c ) \\& \qquad {}+ \mu _{1} ( 1-\mu _{2} ) f ( a,y ) +\mu _{1}\mu _{2}f ( a,c ) \bigr] \\& \qquad {}-\Gamma (\alpha +1) \bigl( ( 1-\lambda _{2} ) ( d-y ) ^{\beta }+ ( 1-\mu _{2} ) ( y-c ) ^{ \beta } \bigr) \bigl[ J_{x+}^{\alpha }f ( b,y ) +J_{x-}^{\alpha }f ( a,y ) \bigr] \\& \qquad {}-\Gamma ( \beta +1 ) \bigl( ( 1-\lambda _{1} ) ( b-x ) ^{\alpha }+ ( 1-\mu _{1} ) ( x-a ) ^{\alpha } \bigr) \bigl[ J_{y+}^{\beta }f ( x,d ) +J_{y-}^{\beta }f ( x,c ) \bigr] \\& \qquad {}-\lambda _{2}\Gamma (\alpha +1) ( d-y ) ^{\beta } \bigl[ J_{x+}^{ \alpha }f ( b,d ) +J_{x-}^{\alpha }f ( a,d ) \bigr] \\& \qquad {}-\mu _{2}\Gamma ( \alpha +1 ) ( y-c ) ^{\beta } \bigl[ J_{x+}^{\alpha }f ( b,c ) +J_{x-}^{\alpha }f ( a,c ) \bigr] \\& \qquad {}-\lambda _{1}\Gamma (\beta +1) ( b-x ) ^{\alpha } \bigl[ J_{y+}^{ \beta }f ( b,d ) +J_{y-}^{\beta }f ( b,c ) \bigr] \\& \qquad {}-\mu _{1}\Gamma ( \beta +1 ) ( x-a ) ^{ \alpha } \bigl[ J_{y+}^{\beta }f ( a,d ) +J_{y-}^{\beta }f ( a,c ) \bigr] \\& \qquad {}+\Gamma (\alpha +1)\Gamma ( \beta +1 ) \bigl[ J_{x+,y+}^{ \alpha ,\beta }f ( b,d ) +J_{x+,y-}^{\alpha ,\beta }f ( b,c ) +J_{x-,y+}^{\alpha ,\beta }f ( a,d ) +J_{x-,y-}^{ \alpha ,\beta }f ( a,c ) \bigr] . \end{aligned}$$

4 New inequalities for Generalized Fractional Integrals

Theorem 2

Assume that the assumptions of Lemma 1hold. Assume also that the mapping \(\vert \frac{\partial ^{2}f}{\partial t\partial s} \vert \) is co-ordinated convex on Δ. Then we obtain the following inequality for generalized fractional integrals:

$$\begin{aligned}& \bigl\vert \Omega (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ( d-y ) \biggl[ A_{1}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert +A_{1}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert \\& \qquad {}+ A_{2}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert +A_{2}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert \biggr] \\& \qquad {}+ ( b-x ) ( y-c ) \biggl[ A_{1}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert +A_{1}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert \\& \qquad {}+ A_{2}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert +A_{2}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert \biggr] \\& \qquad {}+ ( x-a ) ( d-y ) \biggl[ A_{4}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert +A_{4}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert \\& \qquad {}+ A_{3}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert +A_{3}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert \biggr] \\& \qquad {}+ ( x-a ) ( y-c ) \biggl[ A_{4}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert +A_{4}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert \\& \qquad {}+ A_{3}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert +A_{3}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert \biggr]. \end{aligned}$$

Here,

$$ \textstyle\begin{cases} A_{1}=\int _{0}^{1}t \vert \lambda _{1}\Lambda _{1}(x,1)- \Lambda _{1}(x,t) \vert \,dt, \\ A_{2}=\int _{0}^{1} ( 1-t ) \vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \vert \,dt, \\ A_{3}=\int _{0}^{1} ( 1-t ) \vert \mu _{1} \Delta _{1}(x,1)-\Delta _{1}(x,t) \vert \,dt, \\ A_{4}=\int _{0}^{1}t \vert \mu _{1}\Delta _{1}(x,1)- \Delta _{1}(x,t) \vert \,dt,\end{cases} $$
(4.1)

and

$$ \textstyle\begin{cases} B_{1}=\int _{0}^{1}s \vert \lambda _{2}\Lambda _{2}(y,1)- \Lambda _{2}(y,s) \vert \,ds, \\ B_{2}=\int _{0}^{1} ( 1-s ) \vert \lambda _{2} \Lambda _{2}(y,1)-\Lambda _{2}(y,s) \vert \,ds, \\ B_{3}=\int _{0}^{1} ( 1-s ) \vert \mu _{2} \Delta _{2}(y,1)-\Delta _{2}(y,s) \vert \,ds, \\ B_{4}=\int _{0}^{1}s \vert \mu _{2}\Delta _{2}(y,1)- \Delta _{2}(y,s) \vert \,ds.\end{cases} $$
(4.2)

Proof

By taking the modulus in Lemma 1, we have

$$\begin{aligned}& \bigl\vert \Omega (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ( d-y ) \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)- \Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)- \Lambda _{2}(y,s) \bigr\vert \\& \qquad {} \times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s )\,d \bigr) \biggr\vert \,ds\,dt \\& \qquad {}+ ( b-x ) ( y-c ) \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \mu _{2}\Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert \\& \qquad {} \times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s ) c \bigr) \biggr\vert \,ds\,dt \\& \qquad {}+ ( x-a ) ( d-y ) \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \\& \qquad {} \times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s )\,d \bigr) \biggr\vert \,ds\,dt \\& \qquad {}+ ( x-a ) ( y-c ) \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \mu _{2}\Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert \\& \qquad {} \times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s ) c \bigr) \biggr\vert \,ds\,dt. \end{aligned}$$
(4.3)

Since \(\vert \frac{\partial ^{2}f}{\partial t\partial s} \vert \) is a co-ordinated convex, we get

$$\begin{aligned}& \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2} \Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \\& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s )\,d \bigr) \biggr\vert \,ds\,dt \\& \quad \leq \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2} \Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \\& \qquad {} \times \biggl( ts \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert +t ( 1-s ) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert \\& \qquad {}+ ( 1-t ) s \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert + ( 1-t ) ( 1-s ) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert \biggr)\,ds \,dt \\& \quad = A_{1}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert +A_{1}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert +A_{2}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert +A_{2}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert . \end{aligned}$$
(4.4)

Similarly, we have

$$\begin{aligned}& \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \mu _{2} \Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert \end{aligned}$$
(4.5)
$$\begin{aligned}& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s ) c \bigr) \biggr\vert \,ds\,dt \\& \quad \leq A_{1}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert +A_{1}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert +A_{2}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert +A_{2}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert , \\& \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1} \Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2} \Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \\& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s )\,d \bigr) \biggr\vert \,ds\,dt \\& \quad \leq A_{4}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert +A_{4}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert +A_{3}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert +A_{3}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert , \end{aligned}$$
(4.6)

and

$$\begin{aligned}& \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1} \Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \mu _{2} \Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert \\& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s ) c \bigr) \biggr\vert \,ds\,dt \\& \quad \leq A_{4}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert +A_{4}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert +A_{3}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert +A_{3}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert . \end{aligned}$$
(4.7)

If we substitute the inequalities (4.4)–(4.7) in (4.3), we obtain the desired result. This ends the proof of Theorem 2. □

Corollary 3

In Theorem 2, if we assign \(\varphi (t)=t\) and \(\psi (s)=s\) for all \(( t,s ) \in \triangle \), then we obtain the following inequality for Riemann integrals:

$$\begin{aligned}& \bigl\vert \aleph (a,b,x;c,d,y) \bigr\vert \\& \quad \leq \frac{ ( b-x ) ^{2} ( d-y ) ^{2}}{36} \biggl[ \bigl( 2\lambda _{1}^{3}-3 \lambda _{1}+2 \bigr) \bigl( 2 \lambda _{2}^{3}-3 \lambda _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {}+ \bigl( 2\lambda _{1}^{3}-3\lambda _{1}+2 \bigr) \bigl( -2 \lambda _{2}^{3}+6\lambda _{2}^{2}-3\lambda _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert \\& \qquad {}+ \bigl( -2\lambda _{1}^{3}+6\lambda _{1}^{2}-3 \lambda _{1}+1 \bigr) \bigl( 2\lambda _{2}^{3}-3 \lambda _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert \\& \qquad {} + \bigl( -2\lambda _{1}^{3}+6\lambda _{1}^{2}-3\lambda _{1}+1 \bigr) \bigl( -2\lambda _{2}^{3}+6\lambda _{2}^{2}-3\lambda _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert \biggr] \\& \qquad {}+\frac{ ( b-x ) ^{2} ( y-c ) ^{2}}{36} \biggl[ \bigl( 2\lambda _{1}^{3}-3 \lambda _{1}+2 \bigr) \bigl( 2\mu _{2}^{3}-3 \mu _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {}+ \bigl( 2\lambda _{1}^{3}-3\lambda _{1}+2 \bigr) \bigl( -2\mu _{2}^{3}+6 \mu _{2}^{2}-3 \mu _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert \\& \qquad {}+ \bigl( -2\lambda _{1}^{3}+6\lambda _{1}^{2}-3 \lambda _{1}+1 \bigr) \bigl( 2\mu _{2}^{3}-3\mu _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert \\& \qquad {} + \bigl( -2\lambda _{1}^{3}+6\lambda _{1}^{2}-3\lambda _{1}+1 \bigr) \bigl( -2\mu _{2}^{3}+6\mu _{2}^{2}-3\mu _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert \biggr] \\& \qquad {}+\frac{ ( x-a ) ^{2} ( d-y ) ^{2}}{36} \biggl[ \bigl( 2\mu _{1}^{3}-3\mu _{1}+2 \bigr) \bigl( 2\lambda _{2}^{3}-3 \lambda _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {}+ \bigl( 2\mu _{1}^{3}-3\mu _{1}+2 \bigr) \bigl( -2\lambda _{2}^{3}+6 \lambda _{2}^{2}-3 \lambda _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert \\& \qquad {}+ \bigl( -2\mu _{1}^{3}+6\mu _{1}^{2}-3 \mu _{1}+1 \bigr) \bigl( 2 \lambda _{2}^{3}-3 \lambda _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert \\& \qquad {} + \bigl( -2\mu _{1}^{3}+6\mu _{1}^{2}-3 \mu _{1}+1 \bigr) \bigl( -2\lambda _{2}^{3}+6 \lambda _{2}^{2}-3\lambda _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert \biggr] \\& \qquad {}+\frac{ ( x-a ) ^{2} ( y-c ) ^{2}}{36} \biggl[ \bigl( 2\mu _{1}^{3}-3\mu _{1}+2 \bigr) \bigl( 2\mu _{2}^{3}-3\mu _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {}+ \bigl( 2\mu _{1}^{3}-3\mu _{1}+2 \bigr) \bigl( -2\mu _{2}^{3}+6 \mu _{2}^{2}-3\mu _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert \\& \qquad {}+ \bigl( -2\mu _{1}^{3}+6\mu _{1}^{2}-3 \mu _{1}+1 \bigr) \bigl( 2 \mu _{2}^{3}-3\mu _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert \\& \qquad {} + \bigl( -2\mu _{1}^{3}+6\mu _{1}^{2}-3 \mu _{1}+1 \bigr) \bigl( -2\mu _{2}^{3}+6\mu _{2}^{2}-3\mu _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert \biggr] , \end{aligned}$$

where \(\aleph (a,b,x;c,d,y)\) is defined as in Corollary 1.

Corollary 4

In Theorem 2, let us now consider \(\varphi (t)=\frac{t^{\alpha }}{\Gamma (\alpha )}\) and \(\psi (s)=\frac{s^{\beta }}{\Gamma (\beta )}\) for all \(( t,s ) \in \triangle \). Then we have the following inequality for Riemann–Liouville fractional integrals:

$$\begin{aligned}& \bigl\vert \Phi (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl[ C_{1}( \alpha ,\lambda _{1})D_{1}(\beta ,\lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {}+C_{1}(\alpha , \lambda _{1})D_{2}(\beta ,\lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert \\& \qquad {}+ C_{2}(\alpha ,\lambda _{1})D_{1}(\beta ,\lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert +C_{2}(\alpha ,\lambda _{1})D_{2}(\beta , \lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert \biggr] \\& \qquad {}+ ( b-x ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl[ C_{1}( \alpha ,\lambda _{1})D_{1}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {}+C_{1}(\alpha ,\lambda _{1})D_{2}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert \\& \qquad {}+ C_{2}(\alpha ,\lambda _{1})D_{1}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert +C_{2}(\alpha ,\lambda _{1})D_{2}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert \biggr] \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl[ C_{1}( \alpha ,\mu _{1})D_{1}(\beta ,\lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {} +C_{1}(\alpha ,\mu _{1})D_{2}(\beta ,\lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert \\& \qquad {}+ C_{2}(\alpha ,\mu _{1})D_{1}(\beta , \lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert +C_{2}(\alpha ,\mu _{1})D_{2}(\beta ,\lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert \biggr] \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl[ C_{1}( \alpha ,\mu _{1})D_{1}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {} +C_{1}( \alpha ,\mu _{1})D_{2}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert \\& \qquad {}+ C_{2}(\alpha ,\mu _{1})D_{1}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert +C_{2}(\alpha ,\mu _{1})D_{2}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert \biggr]. \end{aligned}$$

Here, \(\Phi (a,b,x;c,d,y)\) is defined as in Corollary 2and

$$ \textstyle\begin{cases} C_{1}(\alpha ,\varsigma )=\frac{\alpha }{\alpha +2}\varsigma ^{1+ \frac{2}{\alpha }}+\frac{1}{\alpha +2}-\frac{\varsigma }{2}, \\ C_{2}(\alpha ,\varsigma )=\frac{2\alpha }{\alpha +1}\varsigma ^{1+ \frac{1}{\alpha }}-\frac{\alpha }{\alpha +2}\varsigma ^{1+\frac{2}{\alpha }}+ \frac{1}{ ( \alpha +1 ) ( \alpha +2 ) }- \frac{\varsigma }{2}, \\ D_{1}(\beta ,\varsigma )=\frac{\beta }{\beta +2}\varsigma ^{1+ \frac{2}{\beta }}+\frac{1}{\beta +2}-\frac{\varsigma }{2}, \\ D_{2}(\beta ,\varsigma )=\frac{2\beta }{\beta +1}\varsigma ^{1+ \frac{1}{\alpha }}-\frac{\beta }{\beta +2}\varsigma ^{1+\frac{2}{\beta }}+ \frac{1}{ ( \beta +1 ) ( \beta +2 ) }- \frac{\varsigma }{2}.\end{cases} $$
(4.8)

Theorem 3

Suppose that the assumptions of Lemma 1hold. Suppose also that the mapping \(\vert \frac{\partial ^{2}f}{\partial t\partial s} \vert ^{q}\), \(q>1\) is co-ordinated convex on Δ. Then we get the following inequality for generalized fractional integrals:

$$\begin{aligned}& \bigl\vert \Omega (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ( d-y ) \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)- \Lambda _{1}(x,t) \bigr\vert ^{p} \bigl\vert \lambda _{2}\Lambda _{2}(y,1)- \Lambda _{2}(y,s) \bigr\vert ^{p}\,ds\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {} \times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ( y-c ) \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)- \Lambda _{1}(x,t) \bigr\vert ^{p} \bigl\vert \mu _{2}\Delta _{2}(y,1)- \Delta _{2}(y,s) \bigr\vert ^{p}\,ds\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {} \times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ( d-y ) \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert ^{p} \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert ^{p}\,ds\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {} \times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ( y-c ) \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert ^{p} \bigl\vert \mu _{2}\Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert ^{p}\,ds\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}}, \end{aligned}$$

where \(\frac{1}{p}+\frac{1}{q}=1\) and \(\Omega (a,b,x;c,d,y)\) are defined as in Lemma 1.

Proof

With the help of the Hölder inequality and co-ordinated convexity of \(\vert \frac{\partial ^{2}f}{\partial t\partial s} \vert ^{q}\), we have

$$\begin{aligned}& \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2} \Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \\& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s ) \,d \bigr) \biggr\vert \,ds\,dt \\& \quad \leq \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert ^{p} \bigl\vert \lambda _{2}\Lambda _{2}(y,1)- \Lambda _{2}(y,s) \bigr\vert ^{p}\,ds\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {}\times \biggl( \int _{0}^{1} \int _{0}^{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s )\,d \bigr) \biggr\vert ^{q}\,ds\,dt \biggr) ^{ \frac{1}{q}} \\& \quad \leq \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert ^{p} \bigl\vert \lambda _{2}\Lambda _{2}(y,1)- \Lambda _{2}(y,s) \bigr\vert ^{p}\,ds\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {} \times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}}. \end{aligned}$$
(4.9)

Similarly, we get

$$\begin{aligned}& \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \mu _{2} \Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert \end{aligned}$$
(4.10)
$$\begin{aligned}& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s ) c \bigr) \biggr\vert \,ds\,dt \\& \quad \leq \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert ^{p} \bigl\vert \mu _{2}\Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert ^{p}\,ds\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}}, \\& \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1} \Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2} \Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \\& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s )\,d \bigr) \biggr\vert \,ds\,dt \\& \quad \leq \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert ^{p} \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert ^{p}\,ds\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {} \times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}}, \end{aligned}$$
(4.11)

and

$$\begin{aligned}& \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1} \Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \mu _{2} \Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert \\& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s ) c \bigr) \biggr\vert \,ds\,dt \\& \quad \leq \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert ^{p} \bigl\vert \mu _{2}\Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert ^{p}\,ds\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {} \times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}}. \end{aligned}$$
(4.12)

By substituting the inequalities (4.9)–(4.12) in (4.3), we establish required result. □

Corollary 5

In Theorem 3, if we select \(\varphi (t)=t\) and \(\psi (s)=s\) for all \(( t,s ) \in \triangle \), then we obtain the following inequality for the Riemann integral:

$$\begin{aligned}& \bigl\vert \aleph (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ^{2} ( d-y ) ^{2} \biggl( \frac{ ( 1-\lambda _{1} ) ^{p+1}+\lambda _{1}^{p+1}}{p+1} \biggr) ^{ \frac{1}{p}} \biggl( \frac{ ( 1-\lambda _{2} ) ^{p+1}+\lambda _{2}^{p+1}}{p+1} \biggr) ^{\frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ^{2} ( y-c ) ^{2} \biggl( \frac{ ( 1-\lambda _{1} ) ^{p+1}+\lambda _{1}^{p+1}}{p+1} \biggr) ^{\frac{1}{p}} \biggl( \frac{ ( 1-\mu _{2} ) ^{p+1}+\mu _{2}^{p+1}}{p+1} \biggr) ^{ \frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( d-y ) ^{2} \biggl( \frac{ ( 1-\mu _{1} ) ^{p+1}+\mu _{1}^{p+1}}{p+1} \biggr) ^{ \frac{1}{p}} \biggl( \frac{ ( 1-\lambda _{2} ) ^{p+1}+\lambda _{2}^{p+1}}{p+1} \biggr) ^{\frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( y-c ) ^{2} \biggl( \frac{ ( 1-\mu _{1} ) ^{p+1}+\mu _{1}^{p+1}}{p+1} \biggr) ^{ \frac{1}{p}} \biggl( \frac{ ( 1-\mu _{2} ) ^{p+1}+\mu _{2}^{p+1}}{p+1} \biggr) ^{ \frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}}. \end{aligned}$$

Here, \(\aleph (a,b,x;c,d,y)\) is defined as in Corollary 1.

Corollary 6

In Theorem 3, Let us consider \(\varphi (t)=\frac{t^{\alpha }}{\Gamma (\alpha )}\) and \(\psi (s)=\frac{s^{\beta }}{\Gamma (\beta )}\) for all \(( t,s ) \in \triangle \). Then we have the inequality for Riemann–Liouville fractional integrals

$$\begin{aligned}& \bigl\vert \Phi (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl( \int _{0}^{1} \bigl\vert \lambda _{1}-t^{\alpha } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \bigl\vert \lambda _{2}-s^{\beta } \bigr\vert ^{p}\,dt \biggr) ^{ \frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl( \int _{0}^{1} \bigl\vert \lambda _{1}-t^{\alpha } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \bigl\vert \mu _{2}-s^{\beta } \bigr\vert ^{p}\,dt \biggr) ^{ \frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl( \int _{0}^{1} \bigl\vert \mu _{1}-t^{\alpha } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \bigl\vert \lambda _{2}-s^{\beta } \bigr\vert ^{p}\,dt \biggr) ^{ \frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl( \int _{0}^{1} \bigl\vert \mu _{1}-t^{\alpha } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \bigl\vert \mu _{2}-s^{\beta } \bigr\vert ^{p}\,dt \biggr) ^{ \frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}}, \end{aligned}$$

where \(\Phi (a,b,x;c,d,y)\) is defined as in Corollary 2.

Theorem 4

Assume that the assumptions of Lemma 1hold. If the mapping \(\vert \frac{\partial ^{2}f}{\partial t\partial s} \vert ^{q}\), \(q\geq 1\) is co-ordinated convex on Δ, then we get the following inequality for generalized fractional integrals:

$$\begin{aligned}& \bigl\vert \Omega (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ( d-y ) \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)- \Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)- \Lambda _{2}(y,s) \bigr\vert \,ds\,dt \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl( A_{1}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+A_{1}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q}+A_{2}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {}+A_{2}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert ^{q} \biggr) ^{ \frac{1}{q}} \\& \qquad {}+ ( b-x ) ( y-c ) \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)- \Lambda _{1}(x,t) \bigr\vert \bigl\vert \mu _{2}\Delta _{2}(y,1)- \Delta _{2}(y,s) \bigr\vert \,ds\,dt \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl( A_{1}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+A_{1}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q}+A_{2}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {}+A_{2}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert ^{q} \biggr) ^{ \frac{1}{q}} \\& \qquad {}+ ( x-a ) ( d-y ) \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \,ds\,dt \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl( A_{4}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+A_{4}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q}+A_{3}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q} \\& \qquad {}+A_{3}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert ^{q} \biggr) ^{ \frac{1}{q}} \\& \qquad {}+ ( x-a ) ( y-c ) \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \mu _{2}\Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert \,ds\,dt \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl( A_{4}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+A_{4}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q}+A_{3}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q} \\& \qquad {}+A_{3}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert ^{q} \biggr) ^{ \frac{1}{q}}. \end{aligned}$$

Here, \(\Omega (a,b,x;c,d,y)\) is defined as in Lemma 1, \(A_{i}\), \(i=1,2,3,4\) are defined as in (4.1) and \(B_{i}\), \(i=1,2,3,4\) are defined as in (4.2).

Proof

Power mean inequality and co-ordinated convexity of \(\vert \frac{\partial ^{2}f}{\partial t\partial s} \vert ^{q}\) yield

$$\begin{aligned}& \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2} \Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \\& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s ) \,d \bigr) \biggr\vert \,ds\,dt \\& \quad \leq \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \,ds\,dt \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \\& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s )\,d \bigr) \biggr\vert ^{q}\,ds\,dt \biggr) ^{\frac{1}{q}} \\& \quad \leq \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \,ds\,dt \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl[ \int _{0}^{1} \int _{0}^{1} \bigg( \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert ts \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {}+ \bigl\vert \lambda _{1}\Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert t ( 1-s ) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {}+ \bigl\vert \lambda _{1}\Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert ( 1-t ) s \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {}+ \bigl\vert \lambda _{1}\Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1) \\& \qquad {}-\Lambda _{2}(y,s) \bigr\vert ( 1-t ) ( 1-s ) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert ^{q} \bigg)\,ds\,dt \biggr] ^{\frac{1}{q}} \\& \quad = \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \,ds\,dt \biggr) ^{1-\frac{1}{q}} \\& \qquad {} \times \biggl( A_{1}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+A_{1}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {}+A_{2}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q}+A_{2}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert ^{q} \biggr) ^{ \frac{1}{q}}. \end{aligned}$$
(4.13)

Similarly, we obtain

$$\begin{aligned}& \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \mu _{2} \Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert \end{aligned}$$
(4.14)
$$\begin{aligned}& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) b,sy+ ( 1-s ) c \bigr) \biggr\vert \,ds\,dt \\& \quad \leq \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \lambda _{1} \Lambda _{1}(x,1)-\Lambda _{1}(x,t) \bigr\vert \bigl\vert \mu _{2}\Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert \,ds\,dt \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl( A_{1}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+A_{1}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q}+A_{2}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {}+A_{2}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert ^{q} \biggr) ^{ \frac{1}{q}}, \\& \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1} \Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2} \Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \\& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s )\,d \bigr) \biggr\vert \,ds\,dt \\& \quad \leq \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \lambda _{2}\Lambda _{2}(y,1)-\Lambda _{2}(y,s) \bigr\vert \,ds\,dt \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl( A_{4}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+A_{4}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q}+A_{3}B_{1} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q} \\& \qquad {}+A_{3}B_{2} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert ^{q} \biggr) ^{ \frac{1}{q}}, \end{aligned}$$
(4.15)

and

$$\begin{aligned}& \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1} \Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \mu _{2} \Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert \\& \qquad {}\times \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f \bigl( tx+ ( 1-t ) a,sy+ ( 1-s ) c \bigr) \biggr\vert \,ds\,dt \\& \leq \biggl( \int _{0}^{1} \int _{0}^{1} \bigl\vert \mu _{1}\Delta _{1}(x,1)-\Delta _{1}(x,t) \bigr\vert \bigl\vert \mu _{2}\Delta _{2}(y,1)-\Delta _{2}(y,s) \bigr\vert \,ds\,dt \biggr) ^{1- \frac{1}{q}} \\& {}\times \biggl( A_{4}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+A_{4}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q} \\& \qquad {}+A_{3}B_{4} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q}+A_{3}B_{3} \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert ^{q} \biggr) ^{ \frac{1}{q}}. \end{aligned}$$
(4.16)

If we substitute the inequalities (4.13)–(4.16) in (4.3), then we establish the desired result. □

Corollary 7

In Theorem 4, if we take \(\varphi (t)=t\) and \(\psi (s)=s\) for all \(( t,s ) \in \triangle \), then we have the following inequality for the Riemann integral:

$$\begin{aligned}& \bigl\vert \aleph (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ^{2} ( d-y ) ^{2} \frac{ ( \lambda _{1}^{2}-\lambda _{1}+\frac{1}{2} ) ^{1-\frac{1}{q}} ( \lambda _{2}^{2}-\lambda _{2}+\frac{1}{2} ) ^{1-\frac{1}{q}}}{6^{\frac{2}{q}}} \\& {}\times \biggl[ \bigl( 2\lambda _{1}^{3}-3\lambda _{1}+2 \bigr) \bigl( 2\lambda _{2}^{3}-3\lambda _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {}+ \bigl( 2\lambda _{1}^{3}-3\lambda _{1}+2 \bigr) \bigl( -2 \lambda _{2}^{3}+6\lambda _{2}^{2}-3\lambda _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {}+ \bigl( -2\lambda _{1}^{3}+6\lambda _{1}^{2}-3 \lambda _{1}+1 \bigr) \bigl( 2\lambda _{2}^{3}-3 \lambda _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {} + \bigl( -2\lambda _{1}^{3}+6\lambda _{1}^{2}-3\lambda _{1}+1 \bigr) \bigl( -2\lambda _{2}^{3}+6\lambda _{2}^{2}-3\lambda _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ^{2} ( y-c ) ^{2} \frac{ ( \lambda _{1}^{2}-\lambda _{1}+\frac{1}{2} ) ^{1-\frac{1}{q}} ( \mu _{2}^{2}-\mu _{2}+\frac{1}{2} ) ^{1-\frac{1}{q}}}{6^{\frac{2}{q}}} \\& \qquad {}\times \biggl[ \bigl( 2\lambda _{1}^{3}-3\lambda _{1}+2 \bigr) \bigl( 2\mu _{2}^{3}-3\mu _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {}+ \bigl( 2\lambda _{1}^{3}-3\lambda _{1}+2 \bigr) \bigl( -2\mu _{2}^{3}+6 \mu _{2}^{2}-3 \mu _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q} \\& \qquad {}+ \bigl( -2\lambda _{1}^{3}+6\lambda _{1}^{2}-3 \lambda _{1}+1 \bigr) \bigl( 2\mu _{2}^{3}-3\mu _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {} + \bigl( -2\lambda _{1}^{3}+6\lambda _{1}^{2}-3\lambda _{1}+1 \bigr) \bigl( -2\mu _{2}^{3}+6\mu _{2}^{2}-3\mu _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( d-y ) ^{2} \frac{ ( \mu _{1}^{2}-\mu _{1}+\frac{1}{2} ) ^{1-\frac{1}{q}} ( \lambda _{2}^{2}-\lambda _{2}+\frac{1}{2} ) ^{1-\frac{1}{q}}}{6^{\frac{2}{q}}} \\& \qquad {}\times \biggl[ \bigl( 2\mu _{1}^{3}-3\mu _{1}+2 \bigr) \bigl( 2 \lambda _{2}^{3}-3\lambda _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {}+ \bigl( 2\mu _{1}^{3}-3\mu _{1}+2 \bigr) \bigl( -2\lambda _{2}^{3}+6 \lambda _{2}^{2}-3 \lambda _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {}+ \bigl( -2\mu _{1}^{3}+6\mu _{1}^{2}-3 \mu _{1}+1 \bigr) \bigl( 2 \lambda _{2}^{3}-3 \lambda _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q} \\& \qquad {} + \bigl( -2\mu _{1}^{3}+6\mu _{1}^{2}-3 \mu _{1}+1 \bigr) \bigl( -2\lambda _{2}^{3}+6 \lambda _{2}^{2}-3\lambda _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( y-c ) ^{2} \frac{ ( \mu _{1}^{2}-\mu _{1}+\frac{1}{2} ) ^{1-\frac{1}{q}} ( \mu _{2}^{2}-\mu _{2}+\frac{1}{2} ) ^{1-\frac{1}{q}}}{6^{\frac{2}{q}}} \\& \qquad {}\times \biggl[ \bigl( 2\mu _{1}^{3}-3\mu _{1}+2 \bigr) \bigl( 2 \mu _{2}^{3}-3\mu _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {}+ \bigl( 2\mu _{1}^{3}-3\mu _{1}+2 \bigr) \bigl( -2\mu _{2}^{3}+6 \mu _{2}^{2}-3\mu _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q} \\& \qquad {}+ \bigl( -2\mu _{1}^{3}+6\mu _{1}^{2}-3 \mu _{1}+1 \bigr) \bigl( 2 \mu _{2}^{3}-3\mu _{2}+2 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q} \\& \qquad {} + \bigl( -2\mu _{1}^{3}+6\mu _{1}^{2}-3 \mu _{1}+1 \bigr) \bigl( -2\mu _{2}^{3}+6\mu _{2}^{2}-3\mu _{2}+1 \bigr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}}, \end{aligned}$$

where \(\aleph (a,b,x;c,d,y)\) is defined as in Corollary 1.

Corollary 8

In Theorem 4, let us consider \(\varphi (t)=\frac{t^{\alpha }}{\Gamma (\alpha )}\) and \(\psi (s)=\frac{s^{\beta }}{\Gamma (\beta )}\) for all \(( t,s ) \in \triangle \). Then we obtain the following inequality for Riemann–Liouville fractional integrals:

$$\begin{aligned}& \bigl\vert \Phi (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ^{\alpha +1} ( d-y ) ^{\beta +1} \bigl( C_{1}( \alpha ,\lambda _{1})+C_{2}(\alpha ,\lambda _{1}) \bigr) ^{1-\frac{1}{q}} \bigl( D_{1}(\beta ,\lambda _{2})+D_{2}( \beta ,\lambda _{2}) \bigr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl[ C_{1}(\alpha ,\lambda _{1})D_{1}( \beta ,\lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+C_{1}(\alpha ,\lambda _{1})D_{2}( \beta , \lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {}+ C_{2}(\alpha ,\lambda _{1})D_{1}(\beta ,\lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q}+C_{2}(\alpha ,\lambda _{1})D_{2}( \beta , \lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ^{\alpha +1} ( y-c ) ^{\beta +1} \bigl( C_{1}( \alpha ,\lambda _{1})+C_{2}(\alpha ,\lambda _{1}) \bigr) ^{1-\frac{1}{q}} \bigl( D_{1}(\beta ,\mu _{2})+D_{2}( \beta ,\mu _{2}) \bigr) ^{1- \frac{1}{q}} \\& \qquad {}\times \biggl[ C_{1}(\alpha ,\lambda _{1})D_{1}( \beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+C_{1}(\alpha ,\lambda _{1})D_{2}( \beta , \mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q} \\& \qquad {}+ C_{2}(\alpha ,\lambda _{1})D_{1}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q}+C_{2}(\alpha ,\lambda _{1})D_{2}( \beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( d-y ) ^{\beta +1} \bigl( C_{1}( \alpha ,\mu _{1})+C_{2}(\alpha ,\mu _{1}) \bigr) ^{1- \frac{1}{q}} \bigl( D_{1}(\beta ,\lambda _{2})+D_{2}( \beta ,\lambda _{2}) \bigr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl[ C_{1}(\alpha ,\mu _{1})D_{1}( \beta ,\lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+C_{1}(\alpha ,\mu _{1})D_{2}( \beta , \lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {}+ C_{2}(\alpha ,\mu _{1})D_{1}(\beta , \lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q}+C_{2}(\alpha ,\mu _{1})D_{2}(\beta ,\lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( y-c ) ^{\beta +1} \bigl( C_{1}( \alpha ,\mu _{1})+C_{2}(\alpha ,\mu _{1}) \bigr) ^{1- \frac{1}{q}} \bigl( D_{1}(\beta ,\mu _{2})+D_{2}( \beta ,\mu _{2}) \bigr) ^{1- \frac{1}{q}} \\& \qquad {}\times \biggl[ C_{1}(\alpha ,\mu _{1})D_{1}( \beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+C_{1}(\alpha ,\mu _{1})D_{2}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q} \\& \qquad {}+ C_{2}(\alpha ,\mu _{1})D_{1}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q}+C_{2}(\alpha ,\mu _{1})D_{2}(\beta ,\mu _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}}. \end{aligned}$$

Here, \(\Phi (a,b,x;c,d,y)\) is defined as in Corollary 2and \(C_{1} \), \(C_{2}\) are defined as in (4.8).

5 Special Cases

In this section, some special cases of our results are presented and we show that our results reduce to inequalities obtained in earlier work.

Remark 4

In Corollary 3:

  1. (1)

    If we assign \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=0\) and if \(\vert \frac{\partial ^{2}}{\partial t\partial s}f ( t,s ) \vert \leq M\) for all \(( t,s ) \in \Delta \), then Corollary 3 reduces to [24, Theorem 3].

  2. (2)

    Let us note that \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=1\), \(x=\frac{a+b}{2}\) and \(y=\frac{c+d}{2}\). Then Corollary 3 reduces to [39, Theorem 2].

  3. (3)

    For \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=\frac{1}{3}\), \(x=\frac{a+b}{2}\) and \(y=\frac{c+d}{2}\), Corollary 3 is equal to [28, Theorem 3].

Remark 5

In Corollary 4;

  1. (1)

    Assume \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=0\) and assume also \(\vert \frac{\partial ^{2}}{\partial t\partial s}f ( t,s ) \vert \leq M\) for all \(( t,s ) \in \Delta \). Then Corollary 4 reduces to [23, Theorem 3].

  2. (2)

    If we put \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=1\), then we have the following trapezoidal type inequalities for Riemann–Liouville fractional integrals:

    $$\begin{aligned}& \bigl\vert \Im (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl[ \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {} +\frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert \\& \qquad {}+ \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert \\& \qquad {} + \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert \biggr] \\& \qquad {}+ ( b-x ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl[ \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {} +\frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert \\& \qquad {}+ \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert \\& \qquad {}+ \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert \biggr] \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl[ \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {} +\frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert \\& \qquad {}+ \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert \\& \qquad {} + \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert \biggr] \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl[ \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {} +\frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert \\& \qquad {}+ \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert \\& \qquad {}+\frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert \biggr] , \end{aligned}$$

    where

    $$\begin{aligned}& \Im (a,b,x;c,d,y) \\& \quad = ( b-x ) ^{\alpha } ( d-y ) ^{\beta } \bigl[ f ( b,d ) \bigr] + ( b-x ) ^{\alpha } ( y-c ) ^{\beta } \bigl[ f ( b,c ) \bigr] \\& \qquad {}+ ( x-a ) ^{\alpha } ( d-y ) ^{\beta } \bigl[ f ( a,d ) \bigr] + ( x-a ) ^{\alpha } ( y-c ) ^{\beta } \bigl[ f ( a,c ) \bigr] \\& \qquad {}-\Gamma (\alpha +1) ( d-y ) ^{\beta } \bigl[ J_{x+}^{ \alpha }f ( b,d ) +J_{x-}^{\alpha }f ( a,d ) \bigr] \\& \qquad {} -\Gamma ( \alpha +1 ) ( y-c ) ^{ \beta } \bigl[ J_{x+}^{\alpha }f ( b,c ) +J_{x-}^{\alpha }f ( a,c ) \bigr] \\& \qquad {}-\Gamma (\beta +1) ( b-x ) ^{\alpha } \bigl[ J_{y+}^{ \beta }f ( b,d ) +J_{y-}^{\beta }f ( b,c ) \bigr] \\& \qquad {}-\Gamma ( \beta +1 ) ( x-a ) ^{ \alpha } \bigl[ J_{y+}^{\beta }f ( a,d ) +J_{y-}^{\beta }f ( a,c ) \bigr] \\& \qquad {}+\Gamma (\alpha +1)\Gamma ( \beta +1 ) \bigl[ J_{x+,y+}^{ \alpha ,\beta }f ( b,d ) +J_{x+,y-}^{\alpha ,\beta }f ( b,c ) +J_{x-,y+}^{\alpha ,\beta }f ( a,d ) +J_{x-,y-}^{ \alpha ,\beta }f ( a,c ) \bigr] . \end{aligned}$$
  3. (3)

    Let us consider \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=\frac{1}{3}\). Then we obtain the following Simpson type inequalities for Riemann–Liouville fractional integrals:

$$\begin{aligned}& \bigl\vert \chi (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl[ C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta , \frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {}+C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert \\& \qquad {}+ C_{2} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert +C_{2}\biggl(\alpha , \frac{1}{3}\biggr)D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert \biggr] \\& \qquad {}+ ( b-x ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl[ C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta , \frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {} +C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert \\& \qquad {}+ C_{2} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert +C_{2} \biggl( \alpha , \frac{1}{3} \biggr) D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert \biggr] \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl[ C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta , \frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {} +C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert \\& \qquad {}+ C_{2} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert +C_{2} \biggl( \alpha , \frac{1}{3} \biggr) D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert \biggr] \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl[ C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta , \frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert \\& \qquad {} +C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert \\& \qquad {}+ C_{2} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert +C_{2} \biggl( \alpha , \frac{1}{3} \biggr) D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert \biggr]. \end{aligned}$$

Here,

$$\begin{aligned}& \chi (a,b,x;c,d,y) \\& \quad = \frac{ ( b-x ) ^{\alpha } ( d-y ) ^{\beta }}{9} \bigl[ 4f ( x,y ) +2f ( x,d ) +2f ( b,y ) +f ( b,d ) \bigr] \\& \qquad {}+ \frac{ ( b-x ) ^{\alpha } ( y-c ) ^{\beta }}{9} \bigl[ 4f ( x,y ) +2f ( x,c ) +2f ( b,y ) +f ( b,c ) \bigr] \\& \qquad {}+ \frac{ ( x-a ) ^{\alpha } ( d-y ) ^{\beta }}{9} \bigl[ 4f ( x,y ) +2f ( x,d ) +2f ( a,y ) +f ( a,d ) \bigr] \\& \qquad {}+ \frac{ ( x-a ) ^{\alpha } ( y-c ) ^{\beta }}{9} \bigl[ 4f ( x,y ) +2f ( x,c ) +2f ( a,y ) +f ( a,c ) \bigr] \\& \qquad {}-\frac{2}{3}\Gamma (\alpha +1) \bigl( ( d-y ) ^{\beta }+ ( y-c ) ^{\beta } \bigr) \bigl[ J_{x+}^{\alpha }f ( b,y ) +J_{x-}^{\alpha }f ( a,y ) \bigr] \\& \qquad {}-\frac{2}{3}\Gamma ( \beta +1 ) \bigl( ( b-x ) ^{\alpha }+ ( x-a ) ^{\alpha } \bigr) \bigl[ J_{y+}^{ \beta }f ( x,d ) +J_{y-}^{\beta }f ( x,c ) \bigr] \\& \qquad {}-\Gamma (\alpha +1)\frac{ ( d-y ) ^{\beta }}{3} \bigl[ J_{x+}^{ \alpha }f ( b,d ) +J_{x-}^{\alpha }f ( a,d ) \bigr] \\& \qquad {} -\Gamma ( \alpha +1 ) \frac{ ( y-c ) ^{\beta }}{3} \bigl[ J_{x+}^{\alpha }f ( b,c ) +J_{x-}^{\alpha }f ( a,c ) \bigr] \\& \qquad {}-\Gamma (\beta +1)\frac{ ( b-x ) ^{\alpha }}{3} \bigl[ J_{y+}^{ \beta }f ( b,d ) +J_{y-}^{\beta }f ( b,c ) \bigr] \\& \qquad {}-\Gamma ( \beta +1 ) \frac{ ( x-a ) ^{\alpha }}{3} \bigl[ J_{y+}^{\beta }f ( a,d ) +J_{y-}^{\beta }f ( a,c ) \bigr] \\& \qquad {}+\Gamma (\alpha +1)\Gamma ( \beta +1 ) \bigl[ J_{x+,y+}^{ \alpha ,\beta }f ( b,d ) +J_{x+,y-}^{\alpha ,\beta }f ( b,c ) +J_{x-,y+}^{\alpha ,\beta }f ( a,d ) +J_{x-,y-}^{ \alpha ,\beta }f ( a,c ) \bigr] . \end{aligned}$$

Remark 6

In Corollary 5;

  1. (1)

    Suppose \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=0\) and suppose also \(\vert \frac{\partial ^{2}}{\partial t\partial s}f ( t,s ) \vert \leq M\) for all \(( t,s ) \in \Delta \). Then Corollary 5 equals [24, Theorem 4].

  2. (2)

    If we take \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=1\), then we get the following trapezoidal type inequalities for the Riemann integral:

    $$\begin{aligned}& \biggl\vert ( b-x ) ( d-y ) f ( b,d ) + ( b-x ) ( y-c ) f ( b,c ) + ( x-a ) ( d-y ) f ( a,d ) \\& \qquad {} + ( x-a ) ( y-c ) f ( a,c ) \\& \qquad {}- ( d-y ) \int _{a}^{b}f ( t,d )\,dt- ( y-c ) \int _{a}^{b}f ( t,c )\,dt \\& \qquad {} - ( b-x ) \int _{c}^{d}f ( b,s )\,dt- ( x-a ) \int _{c}^{d}f ( a,s )\,ds+ \int _{a}^{b} \int _{c}^{d}f ( t,s )\,ds\,dt \biggr\vert \\& \quad \leq ( b-x ) ^{2} ( d-y ) ^{2} \biggl( \frac{1}{p+1} \biggr) ^{\frac{2}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ^{2} ( y-c ) ^{2} \biggl( \frac{1}{p+1} \biggr) ^{\frac{2}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( d-y ) ^{2} \biggl( \frac{1}{p+1} \biggr) ^{\frac{2}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( y-c ) ^{2} \biggl( \frac{1}{p+1} \biggr) ^{\frac{2}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}}. \end{aligned}$$
  3. (3)

    For \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=\frac{1}{3}\), we have the following Simpson type inequalities for the Riemann integral:

    $$\begin{aligned}& \biggl\vert \frac{ ( b-x ) ( d-y ) }{9} \bigl[ 4f ( x,y ) +2_{2}f ( x,d ) +2f ( b,y ) +f ( b,d ) \bigr] \\& \qquad {}+\frac{ ( b-x ) ( y-c ) }{9} \bigl[ 4f ( x,y ) +2f ( x,c ) +2f ( b,y ) +f ( b,c ) \bigr] \\& \qquad {}+\frac{ ( x-a ) ( d-y ) }{9} \bigl[ 4f ( x,y ) +2f ( x,d ) +2f ( a,y ) +f ( a,d ) \bigr] \\& \qquad {}+\frac{ ( x-a ) ( y-c ) }{9} \bigl[ 4f ( x,y ) +2f ( x,c ) +2f ( a,y ) +f ( a,c ) \bigr] \\& \qquad {}-\frac{2}{3} \bigl( ( d-y ) + ( y-c ) \bigr) \int _{a}^{b}f ( t,y )\,dt-\frac{2}{3} \bigl( ( b-x ) + ( x-a ) \bigr) \int _{c}^{d}f ( x,s )\,ds \\& \qquad {}-\frac{ ( d-y ) }{3} \int _{a}^{b}f ( t,d )\,dt-\frac{ ( y-c ) }{3} \int _{a}^{b}f ( t,c )\,dt \\& \qquad {} -\frac{ ( b-x ) }{3} \int _{c}^{d}f ( b,s )\,dt-\frac{ ( x-a ) }{3} \int _{c}^{d}f ( a,s )\,ds+ \int _{a}^{b} \int _{c}^{d}f ( t,s )\,ds\,dt \biggr\vert \\& \quad \leq ( b-x ) ^{2} ( d-y ) ^{2} \biggl( \frac{ ( \frac{2}{3} ) ^{p+1}+ ( \frac{1}{3} ) ^{p+1}}{p+1} \biggr) ^{ \frac{2}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ^{2} ( y-c ) ^{2} \biggl( \frac{ ( \frac{2}{3} ) ^{p+1}+ ( \frac{1}{3} ) ^{p+1}}{p+1} \biggr) ^{ \frac{2}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( d-y ) ^{2} \biggl( \frac{ ( \frac{2}{3} ) ^{p+1}+ ( \frac{1}{3} ) ^{p+1}}{p+1} \biggr) ^{ \frac{2}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( y-c ) ^{2} \biggl( \frac{ ( \frac{2}{3} ) ^{p+1}+ ( \frac{1}{3} ) ^{p+1}}{p+1} \biggr) ^{ \frac{2}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}}. \end{aligned}$$

Remark 7

In Corollary 6;

  1. (1)

    If we assign \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=0\) and if \(\vert \frac{\partial ^{2}}{\partial t\partial s}f ( t,s ) \vert \leq M\) for all \(( t,s ) \in \Delta \), then Corollary 6 reduces to [23, Theorem 4].

  2. (2)

    For \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=1\), the following trapezoidal type inequality for Riemann–Liouville fractional integrals holds:

    $$\begin{aligned}& \bigl\vert \Im (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl( \int _{0}^{1} \bigl\vert 1-t^{\alpha } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \bigl\vert 1-s^{\beta } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl( \int _{0}^{1} \bigl\vert 1-t^{\alpha } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \bigl\vert 1-s^{\beta } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl( \int _{0}^{1} \bigl\vert 1-t^{\alpha } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \bigl\vert 1-s^{\beta } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl( \int _{0}^{1} \bigl\vert 1-t^{\alpha } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \bigl\vert 1-s^{\beta } \bigr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \\& \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}}, \end{aligned}$$

    where \(\Im (a,b,x;c,d,y)\) is defined as in Remark 5.

  3. (3)

    Considering \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=\frac{1}{3}\), we obtain the following Simpson type inequality for Riemann–Liouville fractional integrals:

    $$\begin{aligned}& \bigl\vert \chi (a,b,x;c,d,y) \bigr\vert \\ & \quad \leq ( b-x ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl( \int _{0}^{1} \biggl\vert \frac{1}{3}-t^{\alpha } \biggr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \biggl\vert \frac{1}{3}-s^{\beta } \biggr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \\ & \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\ & \qquad {}+ ( b-x ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl( \int _{0}^{1} \biggl\vert \frac{1}{3}-t^{\alpha } \biggr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \biggl\vert \frac{1}{3}-s^{\beta } \biggr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \\ & \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\ & \qquad {}+ ( x-a ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl( \int _{0}^{1} \biggl\vert \frac{1}{3}-t^{\alpha } \biggr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \biggl\vert \frac{1}{3}-s^{\beta } \biggr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \\ & \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}} \\ & \qquad {}+ ( x-a ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl( \int _{0}^{1} \biggl\vert \frac{1}{3}-t^{\alpha } \biggr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \biggl( \int _{0}^{1} \biggl\vert \frac{1}{3}-s^{\beta } \biggr\vert ^{p}\,dt \biggr) ^{\frac{1}{p}} \\ & \qquad {}\times \biggl( \frac{ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \vert ^{q}+ \vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \vert ^{q}}{4} \biggr) ^{\frac{1}{q}}. \end{aligned}$$

    Here, \(\chi (a,b,x;c,d,y)\) is defined as in Remark 5.

Remark 8

In Corollary 7;

  1. (1)

    Let us consider \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=0\) and \(\vert \frac{\partial ^{2}}{\partial t\partial s}f ( t,s ) \vert \leq M\) for all \(( t,s ) \in \Delta \). Then Corollary 7 equals [24, Theorem 5].

  2. (2)

    For \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=1\), we have the following trapezoidal type inequalities for the Riemann integral:

    $$\begin{aligned}& \biggl\vert ( b-x ) ( d-y ) f ( b,d ) + ( b-x ) ( y-c ) f ( b,c ) + ( x-a ) ( d-y ) f ( a,d ) \\& \qquad {} + ( x-a ) ( y-c ) f ( a,c ) \\& \qquad {}- ( d-y ) \int _{a}^{b}f ( t,d )\,dt- ( y-c ) \int _{a}^{b}f ( t,c )\,dt \\& \qquad {} - ( b-x ) \int _{c}^{d}f ( b,s )\,dt- ( x-a ) \int _{c}^{d}f ( a,s )\,ds+ \int _{a}^{b} \int _{c}^{d}f ( t,s )\,ds\,dt \biggr\vert \\& \quad \leq ( b-x ) ^{2} ( d-y ) ^{2} \frac{1}{2^{2}\cdot 3^{\frac{2}{q}}} \biggl[ \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+2 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {} +2 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q}+4 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ^{2} ( y-c ) ^{2} \frac{1}{2^{2}\cdot 3^{\frac{2}{q}}} \biggl[ \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+2 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q} \\& \qquad {} +2 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q}+4 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( d-y ) ^{2} \frac{1}{2^{2}\cdot 3^{\frac{2}{q}}} \biggl[ \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+2 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {} +2 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q}+4 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( y-c ) ^{2} \frac{1}{2^{2}\cdot 3^{\frac{2}{q}}} \biggl[ \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+2 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q} \\& \qquad {} +2 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q}+4 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}}. \end{aligned}$$
  3. (3)

    If we elect \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=\frac{1}{3}\), then we have he following Simpson type inequalities for the Riemann integral:

    $$\begin{aligned}& \biggl\vert \frac{ ( b-x ) ( d-y ) }{9} \bigl[ 4f ( x,y ) +2_{2}f ( x,d ) +2f ( b,y ) +f ( b,d ) \bigr] \\& \qquad {}+\frac{ ( b-x ) ( y-c ) }{9} \bigl[ 4f ( x,y ) +2f ( x,c ) +2f ( b,y ) +f ( b,c ) \bigr] \\& \qquad {}+\frac{ ( x-a ) ( d-y ) }{9} \bigl[ 4f ( x,y ) +2f ( x,d ) +2f ( a,y ) +f ( a,d ) \bigr] \\& \qquad {}+\frac{ ( x-a ) ( y-c ) }{9} \bigl[ 4f ( x,y ) +2f ( x,c ) +2f ( a,y ) +f ( a,c ) \bigr] \\& \qquad {}-\frac{2}{3} \bigl( ( d-y ) + ( y-c ) \bigr) \int _{a}^{b}f ( t,y )\,dt-\frac{2}{3} \bigl( ( b-x ) + ( x-a ) \bigr) \int _{c}^{d}f ( x,s )\,ds \\& \qquad {}-\frac{ ( d-y ) }{3} \int _{a}^{b}f ( t,d )\,dt-\frac{ ( y-c ) }{3} \int _{a}^{b}f ( t,c )\,dt \\& \qquad {} -\frac{ ( b-x ) }{3} \int _{c}^{d}f ( b,s )\,dt-\frac{ ( x-a ) }{3} \int _{c}^{d}f ( a,s )\,ds+ \int _{a}^{b} \int _{c}^{d}f ( t,s )\,ds \,dt \biggr\vert \\& \quad \leq ( b-x ) ^{2} ( d-y ) ^{2} \frac{5^{2}}{18^{2} ( 45 ) ^{\frac{2}{q}}} \\& \qquad {}\times \biggl[ 29\cdot 29 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+29\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q}+29\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {} +16\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ^{2} ( y-c ) ^{2} \frac{5^{2}}{18^{2} ( 45 ) ^{\frac{2}{q}}} \\& \qquad {}\times \biggl[ 29\cdot 29 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+29\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q}+29\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {} +16\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( d-y ) ^{2} \frac{5^{2}}{18^{2} ( 45 ) ^{\frac{2}{q}}} \\& \qquad {}\times \biggl[ 29\cdot 29 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+29\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q}+29\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q} \\& \qquad {} +16\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{2} ( y-c ) ^{2} \frac{5^{2}}{18^{2} ( 45 ) ^{\frac{2}{q}}} \\& \qquad {}\times \biggl[ 29\cdot 29 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q}+29\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q}+29\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q} \\& \qquad {} +16\cdot 16 \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}}. \end{aligned}$$

Remark 9

In Corollary 8;

  1. (1)

    Assume \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=0\) and assume also \(\vert \frac{\partial ^{2}}{\partial t\partial s}f ( t,s ) \vert \leq M\) for all \(( t,s ) \in \Delta \). Then Corollary 8 reduces to [23, Theorem 5].

  2. (2)

    If \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=1\) is chosen, then we get the following trapezoidal type inequality for Riemann–Liouville fractional integrals:

    $$\begin{aligned}& \bigl\vert \Im (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl( \frac{\alpha }{\alpha +1} \biggr) ^{1-\frac{1}{q}} \biggl( \frac{\beta }{\beta +1} \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl[ \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {} + \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {}+ \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {} + \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl( \frac{\alpha }{\alpha +1} \biggr) ^{1-\frac{1}{q}} \biggl( \frac{\beta }{\beta +1} \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl[ \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {} + \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q} \\& \qquad {}+ \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {} + \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl( \frac{\alpha }{\alpha +1} \biggr) ^{1-\frac{1}{q}} \biggl( \frac{\beta }{\beta +1} \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl[ \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {} + \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {}+ \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q} \\& \qquad {} + \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl( \frac{\alpha }{\alpha +1} \biggr) ^{1-\frac{1}{q}} \biggl( \frac{\beta }{\beta +1} \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl[ \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {} + \frac{\alpha }{2 ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q} \\& \qquad {}+ \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta }{2 ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q} \\& \qquad {} + \frac{\alpha ^{2}+3\alpha }{2 ( \alpha +1 ) ( \alpha +2 ) } \frac{\beta ^{2}+3\beta }{2 ( \beta +1 ) ( \beta +2 ) } \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}}, \end{aligned}$$

    where \(\Im (a,b,x;c,d,y)\) is defined as in Remark 5.

  3. (3)

    For \(\lambda _{1}=\lambda _{2}=\mu _{1}=\mu _{2}=\frac{1}{3}\), the following Simpson type inequality for Riemann–Liouville fractional integrals holds:

    $$\begin{aligned}& \bigl\vert \chi (a,b,x;c,d,y) \bigr\vert \\& \quad \leq ( b-x ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl( C_{1} \biggl( \alpha ,\frac{1}{3} \biggr)+C_{2} \biggl( \alpha , \frac{1}{3} \biggr) \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl( D_{1} \biggl( \beta , \frac{1}{3} \biggr) +D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggr) ^{1- \frac{1}{q}} \\& \qquad {}\times \biggl[ C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {}+C_{1} \biggl( \alpha , \frac{1}{3} \biggr) D_{2} \biggl( \beta , \frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {}+ C_{2} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1}(\beta , \lambda _{2}) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {}+C_{2} \biggl( \alpha , \frac{1}{3} \biggr) D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,d ) \biggr\vert ^{q} \biggr] ^{\frac{1}{q}} \\& \qquad {}+ ( b-x ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl( C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) +C_{2} \biggl( \alpha , \frac{1}{3} \biggr) \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl( D_{1} \biggl( \beta , \frac{1}{3} \biggr) +D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggr) ^{1- \frac{1}{q}} \\& \qquad {}\times \biggl[ C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {}+C_{1} \biggl( \alpha , \frac{1}{3} \biggr) D_{2} \biggl( \beta , \frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q} \\& \qquad {}+ C_{2} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,y ) \biggr\vert ^{q} \\& \qquad {}+C_{2} \biggl( \alpha , \frac{1}{3} \biggr) D_{2} \biggl( \beta , \frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( b,c ) \biggr\vert ^{q} \biggr] ^{ \frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( d-y ) ^{\beta +1} \biggl( C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) +C_{2} \biggl( \alpha , \frac{1}{3} \biggr) \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl( D_{1} \biggl( \beta , \frac{1}{3} \biggr) +D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggr) ^{1- \frac{1}{q}} \\& \qquad {}\times \biggl[ C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {}+C_{1} \biggl( \alpha , \frac{1}{3} \biggr) D_{2} \biggl( \beta , \frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,d ) \biggr\vert ^{q} \\& \qquad {}+ C_{2} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q} \\& \qquad {}+C_{2} \biggl( \alpha , \frac{1}{3} \biggr) D_{2} \biggl( \beta , \frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,d ) \biggr\vert ^{q} \biggr] ^{ \frac{1}{q}} \\& \qquad {}+ ( x-a ) ^{\alpha +1} ( y-c ) ^{\beta +1} \biggl( C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) +C_{2} \biggl( \alpha , \frac{1}{3} \biggr) \biggr) ^{1-\frac{1}{q}} \\& \qquad {}\times \biggl( D_{1} \biggl( \beta , \frac{1}{3} \biggr) +D_{2} \biggl( \beta ,\frac{1}{3} \biggr) \biggr) ^{1- \frac{1}{q}} \\& \qquad {}\times \biggl[ C_{1} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,y ) \biggr\vert ^{q} \\& \qquad {}+C_{1} \biggl( \alpha , \frac{1}{3} \biggr) D_{2} \biggl( \beta , \frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( x,c ) \biggr\vert ^{q} \\& \qquad {}+ C_{2} \biggl( \alpha ,\frac{1}{3} \biggr) D_{1} \biggl( \beta ,\frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,y ) \biggr\vert ^{q} \\& \qquad {}+C_{2} \biggl( \alpha , \frac{1}{3} \biggr) D_{2} \biggl( \beta , \frac{1}{3} \biggr) \biggl\vert \frac{\partial ^{2}}{\partial t\partial s}f ( a,c ) \biggr\vert ^{q} \biggr] ^{ \frac{1}{q}}. \end{aligned}$$

    Here, \(\chi (a,b,x;c,d,y)\) is defined as in Remark 5.

6 Conclusion

In this paper, we used the concepts of generalized fractional integrals and proved some new generalized inequalities for partially differentiable co-ordinated convex mappings. In addition, we discussed the special cases of the main results. Furthermore, several new inequalities of trapezoidal and Simpson type for partially differentiable co-ordinated convex functions via Riemann and Riemann–Liouville fractional integrals. It is an interesting and new problem that the upcoming researchers can address similar inequalities for other type co-ordinated convex functions in their future research. By using the Jensen inequality and a coordinated concave function in our main lemma, one can obtain new inequalities for generalized fractional integrals