Kumwimba Seya Didier^*, Kasongo Ntambwe Dany, Panga Lutanda Grégoire, Mayuke Katshongo Jean Paul, Lunda Ngoy Jean Pierre
Département de Mathématiques et Informatique, Faculté des Sciences, Université de Lubumbashi, Lubumbashi, Democratic Republic of the Congo.
DOI: 10.4236/jamp.2024.123056   PDF    HTML   XML   87 Downloads   319 Views   Citations

Abstract

In this article, we propose by using the Hausdorff distance Simpson’s rule for the triple integral of a fuzzy-valued function and the error bound of this method, one of the variables of which is fuzzy. In addition, thin δ-fine partitions are introduced. The integration domain is a quasi-fuzzy parallelipiped. A numerical example is presented in order to show the application and the significance of the method.

Keywords

Fuzzy-Valued Function, Hausdorff Distance, Quasi-Fuzzy Parallelipiped, Henstock Integral

Share and Cite:

1. Introduction

The term fuzzy integral was introduced by Sugeno [1] . In order to evaluate a fuzzy set, some methods digital ones have recently been proposed. Wu [2] [3] , Allahviranloo [4] and Fariborzi [5] [6] have developed some numerical methods to evaluate fuzzy integrals using quadratic methods and the definition of level sets. Wu and Gong [7] proposed the Henstock integral of a fuzzy-valued function and developed this work by applying the concept of fuzzy function differentiability. Bede and Gal [8] applied the quadrature rule to evaluate the integral of a fuzzy-valued function.

In our previous paper we developed the Henstock-Kurzweil triple integral of a fuzzy function on a classic parallelepiped [9] .

In the present article, we develop this idea for a fuzzy-valued function (function with three variables of which one is fuzzy) by applying Simpson’s rule which is powerful tool for numerical integration, especially when dealing with curves and polynomial functions and introducing Henstock-Kurzweil’s triple integral. Certain domains of space can present random deformations. Thus, determining the centers of gravity of these domains will require taking into account the different variations. These variations can be done according to the three components of space. In this paper, we model these transformations by fuzzy variables and a particular case is treated. It involves taking into account a fuzzy variable among the three variables and treating the other two classically.

In section two, we present some preliminary notions on fuzzy sets as well as some fundamental theorems that we will use later.

In section three, we introduce Simpson’s rule to compute the fuzzy triple integral of Henstock on a quasi-fuzzy parallelepiped. The calculations will be based on Simpson’s rule for the fuzzy triple integral over an almost fuzzy three-dimensional domain. This rule generalizes that used for the calculation of the double fuzzy integral on a quasi-fuzzy two-dimensional domain, Didier and Zerbo [10] . The method involves the approximation of a fuzzy integral on a quasi-fuzzy parallelepiped. To calculate this integral, we use a fuzzy approximation formula based on fuzzy quadratic interpolation polynomials in 3 dimensions, one of which is fuzzy and the other two crisp.

Finally, in order to explain an application of the proposed method, in Section 4, a fuzzy triple integral of a fuzzy function which depends on three variables one of which is fuzzy is evaluated in order to show the effectiveness of the mentioned method.

2. Preliminaries

In this section, we talk about some basic definitions of fuzzy sets theory which are being used in the following.

Definition 2.1. [11] Let X be a non-empty reference set. A nonempty subset $\left\{\left(x\mathrm{,}\tilde{A}\left(x\right)\right)\mathrm{<mo>\; :\; </mo>}x\in X\right\}$ of $X\times \left[\mathrm{0,1}\right]$ such that $\tilde{A}\mathrm{<mo>\; :\; </mo>}X\to \left[\mathrm{0,1}\right]$ is a fuzzy subset (fuzzy set) of X. The function $\tilde{A}$ is itself called fuzzy set. $\tilde{A}\left(x\right)$ denotes the degree of membership of the element x in the fuzzy set $\tilde{A}$ .

We denote by X_F the collection of all fuzzy-subsets of X.

Definition 2.2. [11] [12] Let $ℝ$ be a real set. Given a fuzzy-subset $\tilde{u}\mathrm{<mo>\; :\; </mo>}ℝ\to \left[\mathrm{0,1}\right]$ satisfying the properties below:

1) $\tilde{u}$ is normal, i.e. $\exists x_0\in ℝ$ such that $\tilde{u}\left(x_0\right)=1$ ,

2) $\tilde{u}$ is a convex fuzzy set, i.e.,

$\tilde{u}\left(\lambda x+\left(1-\lambda \right)y\right)\ge \min \left\{\tilde{u}\left(x\right),\tilde{u}\left(y\right)\right\}\forall x,y\in ℝ,\lambda \in \left[0,1\right],$

3) $\tilde{u}$ is upper semi-continuous on $ℝ$ , i.e. $\forall \epsilon >0$ , $\exists \delta >0$ such that $\left|x-x_0\right|<\delta \Rightarrow \tilde{u}\left(x\right)-\tilde{u}\left(x_0\right)<\epsilon$ .

4) The set $\overline{\left\{x\in ℝ:\stackrel{˜}{u}\left(x\right)>0\right\}}$ is compact, where $\bar{B}$ denotes the closure of B.

This function $\tilde{u}$ is called a fuzzy number.

We denote by ${ℝ}_F$ the set of all fuzzy real numbers. We define ${\left[\tilde{u}\right]}^{\alpha }=\left\{x\in ℝ\mathrm{<mo>\; :\; </mo>}\tilde{u}\left(x\right)\ge \alpha \right\}$ and ${\left[\tilde{u}\right]}^0=\left\{x\in ℝ:\tilde{u}\left(x\right)>0\right\}$ , for $0<\alpha \le 1$ , as the α-cut and respectively the support of a fuzzy number such as $\tilde{u}$ . Moreover, we define ${\tilde{u}}_{\alpha }^{-}=\inf {\left[\tilde{u}\right]}^{\alpha }$ and ${\tilde{u}}_{\alpha }^{+}=\sup {\left[\tilde{u}\right]}^{\alpha }$ .

A triangular fuzzy number $\tilde{u}=\left(a\mathrm{,}b\mathrm{,}c\right)$ where, $a<b<c$ and $a\mathrm{,}b\mathrm{,}c\in ℝ$ is defined by ${\tilde{u}}_{\alpha }^{-}=a+\left(b-a\right)\alpha$ and ${\tilde{u}}_{\alpha }^{+}=c-\left(c-b\right)\alpha$ .

Definition 2.3. [11] [12] A set of level $\alpha$ of a fuzzy set satisfies the following properties $\forall \tilde{u}\mathrm{,}\tilde{v}\in {ℝ}_F^n$ and $\alpha ,\beta \in I=\left(0,1\right]$ ,

1) ${\tilde{u}}_0={ℝ}^n$ ,

2) $\alpha \le \beta \Rightarrow {\tilde{u}}_{\beta }\subset {\tilde{u}}_{\alpha }$ ,

3) ${\tilde{u}}_{\alpha }={\displaystyle {\cap }_{\beta <\alpha }}{\tilde{u}}_{\beta }$ ,

4) ${\left(\tilde{u}\vee \tilde{v}\right)}_{\alpha ={\tilde{u}}_{\alpha }\cup {\tilde{v}}_{\alpha }}$ and ${\left(\tilde{u}\wedge \tilde{v}\right)}_{\alpha ={\tilde{u}}_{\alpha }\cap {\tilde{v}}_{\alpha }}$ .

The theorem below shows that any fuzzy set can be represented by a family of its $\alpha$ cuts $\left\{{\tilde{u}}_{\alpha }\mathrm{<mo>\; :\; </mo>}\alpha \in I\right\}$ , and it can be represented by its countable sets of level $\alpha$ denoted by $\left\{{\tilde{u}}_{\alpha }\mathrm{<mo>\; :\; </mo>}\alpha \in I\cap \mathbb{Q}\right\}$ .

Theorem 2.4. [12] Let $\tilde{u}\in {ℝ}_F^n$ and let $\left\{{\tilde{u}}_{\alpha }\mathrm{<mo>\; :\; </mo>}\alpha \in I\right\}$ be a family of its sets of level $\alpha$ . For all $x\in {ℝ}^n$ , we have:

$\tilde{u}\left(x\right)=\underset{\alpha \in I}{\sup }\left[\alpha \cdot I_{{\tilde{u}}_{\alpha }}\left(x\right)\right]=\underset{\alpha \in I\cap \mathbb{Q}}{\sup }\left[\alpha \cdot I_{{\tilde{u}}_{\alpha }}\left(x\right)\right]$

which can also be written

$\tilde{u}\left(x\right)=\sup \left\{\alpha \in I:x\in {\tilde{u}}_{\alpha }\right\}=\sup \left\{\alpha \in I\cap \mathbb{Q}:x\in {\tilde{u}}_{\alpha }\right\}$

Let $\left\{M_{\alpha }\mathrm{<mo>\; :\; </mo>}\alpha \in I\cap \mathbb{Q}\right\}$ be a family of nonempty closed sets of ${ℝ}^n$ such that $M_{\alpha }\supset M_{\beta }\forall \alpha ,\beta \in I\cap \mathbb{Q}$ with $\alpha <\beta$ .

Then the function $\tilde{u}$ defined by

$\tilde{u}\left(x\right)=\sup \left\{\alpha \in I\cap \mathbb{Q}:x\in M_{\alpha }\right\}$ ,

is upper semi-continuous. Additionally, it checks

${\tilde{u}}_{\alpha }={\displaystyle \underset{\beta <\alpha ,\beta \in \mathbb{Q}\cap \left(0,1\right]}{\cap }}{M}_{\beta },\alpha \in \left(0,1\right].$

Remark 2.5. A fuzzy set $\tilde{u}$ is said to be convex if for all $\alpha \in I$ , the level set $\alpha$ denoted by ${\tilde{u}}_{\alpha }$ is a sub convex set of ${ℝ}^n$ .

We denote by ${ℝ}_{F_c}^n$ the family of all convex fuzzy sets.

Definition 2.6. [11] [12] Let $\tilde{u}\mathrm{,}\tilde{v}\in {ℝ}_F^n$ we define for all $x\in {ℝ}^n$

$\left(\tilde{u}\oplus \tilde{v}\right)\left(x\right)=\sup \left\{\alpha \in I\mathrm{<mo>\; :\; </mo>}x\in {\tilde{u}}_{\alpha }{\oplus }_{int}{\tilde{v}}_{\alpha }\right\}$

and

$\left(\lambda \odot \tilde{u}\right)\left(x\right)=\{\begin{array}{ll}\tilde{u}\left(\frac{x}{\lambda }\right)\hfill & \text{si}\lambda \ne 0\hfill \\ I_0\left(x\right)\hfill & \text{si}\lambda =0\hfill \end{array}$

where $I_0$ is the characteristic function of the singleton {0}.

For $\tilde{u}\mathrm{,}\tilde{v}\in {ℝ}_F^n$ and $\lambda \in ℝ$ , we can define the sum $\tilde{u}\oplus \tilde{v}$ and the product $\lambda \odot \tilde{u}$ by

${\left[\tilde{u}\oplus \tilde{v}\right]}^{\alpha }={\left[\tilde{u}\right]}^{\alpha }{\oplus }_{int}{\left[\tilde{v}\right]}^{\alpha }$ and ${\left[\lambda \odot \tilde{u}\right]}^{\alpha }=\lambda {\odot }_{int}{\left[\tilde{u}\right]}^{\alpha }$ $\forall \alpha \in \left[\mathrm{0,1}\right]$ .

Given $\tilde{u}\mathrm{,}\tilde{v}$ two convex fuzzy sets, another definition of addition is given by

$\left(\tilde{u}\oplus \tilde{v}\right)\left(x\right)=\underset{x=x_1+x_2}{\sup }\min \left\{\tilde{u}\left(x_1\right)\mathrm{,}\tilde{v}\left(x_2\right)\right\}\mathrm{,}x\in {ℝ}^n\mathrm{.}$

Definition 2.7 [11] The three metrics below generalize the Hausdorff metric: for all $\tilde{u}\mathrm{,}\tilde{v}\in {ℝ}_{F_c}^n$

$D_{\infty }\left(\tilde{u}\mathrm{,}\tilde{v}\right)=\underset{\alpha \in \left(\mathrm{0,1}\right]}{\sup }D\left({\tilde{u}}_{\alpha }\mathrm{,}{\tilde{v}}_{\alpha }\right)\mathrm{<mo>\; ;\; </mo>}$

$D_1\left(\tilde{u}\mathrm{,}\tilde{v}\right)={\displaystyle {\int }_0^{1}D}\left({\tilde{u}}_{\alpha }\mathrm{,}{\tilde{v}}_{\alpha }\right)\text{d}\alpha \mathrm{<mo>\; ;\; </mo>}$

and

$D_p\left(\tilde{u}\mathrm{,}\tilde{v}\right)={\left\{{\displaystyle {\int }_0^1}{\left[D\left({\tilde{u}}_{\alpha }\mathrm{,}{\tilde{v}}_{\alpha }\right)\right]}^p\text{d}\alpha \right\}}^{\frac{1}{p}}$ ; for $p>1$ .

We note by

${\Vert \tilde{u}\Vert }_F=D_{\infty }\left(\tilde{u}\mathrm{,}{I}_0\right)=\underset{\alpha >0}{\sup }{\Vert {\tilde{u}}_{\alpha }\Vert }_K$

3. Triple Simpson’s Rule for the Fuzzy Henstock-Kurzweil Triple Integrals

The concept of the Henstock integral for a fuzzy number-valued function were introduced by Wu and Gong [7] . We introduce this definition for a three-dimensional fuzzy number-valued function in which one of the three variables is fuzzy.

Let $\tilde{f}\mathrm{<mo>\; :\; </mo>}\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\to {ℝ}_F$ and ${\Delta }_m:a=x_0<x_1<x_2<\cdots <x_m=b$ , ${\Delta }_n:c=y_0<y_1<y_2<\cdots <y_n=d$ and ${\Delta }_s:\tilde{p}={\tilde{z}}_0\prec {\tilde{z}}_1\prec {\tilde{z}}_2\prec \cdots \prec {\tilde{z}}_s=\tilde{q}$ be the partitions of the intervals $\left[a\mathrm{,}b\right]\mathrm{,}\left[c\mathrm{,}d\right]$ and $\left[\tilde{p}\mathrm{,}\tilde{q}\right]$ respectively.

Consider the points ${\xi }_i\in \left[x_{i-1},x_i\right]$ , $i=1,2,\cdots ,m$ ; ${\eta }_j\in \left[y_{j-1}\mathrm{,}{y}_j\right]$ , $j=\mathrm{1,2,}\cdots \mathrm{,}n$ ; ${\tilde{\zeta }}_k\in \left[{\tilde{z}}_{k-1}\mathrm{,}{\tilde{z}}_k\right]$ , $k=\mathrm{1,2,}\cdots \mathrm{,}s$ .

Let $\pi$ be the partition defined by

$\pi =\left\{\left(\left({\xi }_i,{\eta }_j,{\tilde{\zeta }}_k\right),\left[x_{i-1},x_i\right]\times \left[y_{j-1},y_j\right]\otimes \left[{\tilde{z}}_{k-1},{\tilde{z}}_k\right]\right)\right\}$

In what follows we will define the fuzzy Henstock triple integral.

Let the fuzzy gauge function defined on the type 1 quasi-fuzzy parallelepiped al domain $\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]$ by

$\tilde{\delta }\mathrm{<mo>\; :\; </mo>}\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\to {ℝ}_F\mathrm{.}$

The partition $\pi$ is said to be $\tilde{\delta }$ -fine if

$\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\in \left[x_{i-1}\mathrm{,}{x}_i\right]\times \left[y_{j-1}\mathrm{,}{y}_j\right]\otimes \left[{\tilde{z}}_{k-1}\mathrm{,}{\tilde{z}}_k\right]\subseteq \Delta \left(\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\mathrm{,}\tilde{\delta }\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\right)\mathrm{,}$

where $\Delta \left(\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\mathrm{,}\tilde{\delta }\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\right)$ is the open quasi-fuzzy ball with center $\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)$ and positive radius $\tilde{\delta }\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)$ .

The function $\tilde{f}\mathrm{<mo>\; :\; </mo>}\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\to {ℝ}_F$ is said to be Henstock integrable if $\exists \tilde{I}\in {ℝ}_F$ such that $\forall \epsilon >0$ , there exists $\tilde{\delta }$ (called fuzzy gauge function) such as for any subdivision ${\tilde{\zeta }}_k$ -fine (cfr the partition $\pi$ defined above), we have:

$D\left(\underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right)\odot \left({\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\right)\otimes \tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\mathrm{,}\tilde{I}\right)<\epsilon \mathrm{,}$

where $\tilde{I}=\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x$ is the Fuzzy Henstock Triple Integral.

Definition 3.1. Let $\tilde{f}\mathrm{<mo>\; :\; </mo>}\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\to {ℝ}_F$ a bounded fuzzy function. So the function ${\omega }_{\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]}\left(\tilde{f}\mathrm{,.,.,.}\right)\mathrm{<mo>\; :\; </mo>}{ℝ}^{+}\cup \left\{\infty \right\}\to {ℝ}^{+}$ such that

$\begin{array}{l}{\omega }_{\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]}\left(\tilde{f}\mathrm{,}{\delta }_1\mathrm{,}{\delta }_2\mathrm{,}{\delta }_3\right)\\ =\sup \{D\left(\tilde{f}\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)\mathrm{,}\tilde{f}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\right)\mathrm{<mo>\; ;\; </mo>}\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)\mathrm{,}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\in \left[a\mathrm{,}b\right]\\ \times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\mathrm{<mo>\; :\; </mo>}\left|x_1-x_2\right|\le {\delta }_1\mathrm{,}\left|y_1-y_2\right|\le {\delta }_2\mathrm{,}{\Vert {\tilde{z}}_1\ominus {\tilde{z}}_2\Vert }_{{ℝ}_F}\le {\delta }_3\}\mathrm{.}\end{array}$

where

$\begin{array}{c}{\Vert {\tilde{z}}_1\ominus {\tilde{z}}_2\Vert }_{{ℝ}_F}=\Vert \left[{\tilde{z}}_{1_{\alpha }}^L\mathrm{,}{\tilde{z}}_{1_{\alpha }}^U\right]{\ominus }_{int}\left[{\tilde{z}}_{2_{\alpha }}^L\mathrm{,}{\tilde{z}}_{2_{\alpha }}^U\right]\Vert \\ =\inf \left\{\left|t-u\right|\mathrm{,}t\in \left[{\tilde{z}}_{1_{\alpha }}^L\mathrm{,}{\tilde{z}}_{1_{\alpha }}^U\right]\mathrm{,}u\in \left[{\tilde{z}}_{2_{\alpha }}^L\mathrm{,}{\tilde{z}}_{2_{\alpha }}^U\right]\mathrm{<mo>\; ;\; </mo>}\forall \alpha \in \left[\mathrm{0,1}\right]\right\}\mathrm{.}\end{array}$

For this reason,

$\begin{array}{l}{\omega }_{\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]}\left(\tilde{f}\mathrm{,}{\delta }_1\mathrm{,}{\delta }_2\mathrm{,}{\delta }_3\right)\\ =\sup \{D\left(\tilde{f}\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_{1_{\alpha }}^L\right)\mathrm{,}\tilde{f}\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_{1_{\alpha }}^U\right)\mathrm{,}\tilde{f}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_{2_{\alpha }}^L\right)\mathrm{,}\tilde{f}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_{2_{\alpha }}^U\right)\right)\mathrm{<mo>\; ;\; </mo>}\\ \left|x_1-x_2\right|\le {\delta }_1\mathrm{,}\left|y_1-y_2\right|\le {\delta }_2\mathrm{,}\inf \mathrm{<mo>\; \{\; </mo>}\left|t-u\right|\mathrm{,}t\in \left[{\tilde{z}}_{1_{\alpha }}^L\mathrm{,}{\tilde{z}}_{1_{\alpha }}^U\right]\mathrm{,}u\in \left[{\tilde{z}}_{2_{\alpha }}^L\mathrm{,}{\tilde{z}}_{2_{\alpha }}^U\right]\mathrm{<mo>\; ;\; </mo>}\forall \alpha \in \left[\mathrm{0,1}\right]\}\mathrm{.}\end{array}$

$\forall \left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)\mathrm{,}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\in \left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right].$

$\omega$ is the modulus of oscillation of the function $\tilde{f}$ on $\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]$ . If furthermore, $\tilde{f}$ is continuous on $\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]$ , then $\omega$ is called the uniform continuity modulus of $\tilde{f}$ .

We can prove the following theorem from the definition 3.1.

Theorem 3.2. The following statements, concerning the modulus of oscillation are true.

1)

$D\left(\tilde{f}\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)\mathrm{,}\tilde{f}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\right)\le {\omega }_{\left(\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\right)}\left(\tilde{f}\mathrm{,}\left|x_1-x_2\right|\mathrm{,}\left|y_1-y_2\right|\mathrm{,}{\Vert {\tilde{z}}_1\ominus {\tilde{z}}_2\Vert }_{{ℝ}_F}\right)$

$\forall \left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)\mathrm{,}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\in \left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]$ ,

2) ${\omega }_{\left(\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\right)}\left(\tilde{f}\mathrm{,}{\delta }_1\mathrm{,}{\delta }_2\mathrm{,}{\delta }_3\right)$ is a non-decreasing mapping in ${\delta }_1\mathrm{,}{\delta }_2$ and ${\delta }_3$ ,

3) ${\omega }_{\left(\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\right)}\left(\tilde{f}\mathrm{,0,0,0}\right)=0$ ,

4) ${\omega }_{\left(\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\right)}\left(\tilde{f}\mathrm{,}m{\delta }_1\mathrm{,}n{\delta }_2\mathrm{,}s{\delta }_3\right)\le mns{\omega }_{\left(\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\right)}\left(\tilde{f}\mathrm{,}{\delta }_1\mathrm{,}{\delta }_2\mathrm{,}{\delta }_3\right)$ $\forall {\delta }_1\mathrm{,}{\delta }_2\mathrm{,}{\delta }_3\ge 0$ and $m\mathrm{,}n\mathrm{,}s\in \mathbb{N}$ ,

5) $\begin{array}{l}{\omega }_{\left(\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\right)}\left(\tilde{f}\mathrm{,}{\lambda }_1{\delta }_1\mathrm{,}{\lambda }_2{\delta }_2\mathrm{,}{\lambda }_3{\delta }_3\right)\\ \le \left({\lambda }_1+1\right)\left({\lambda }_2+1\right)\left({\lambda }_3+1\right){\omega }_{\left(\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\right)}\left(\tilde{f}\mathrm{,}{\delta }_1\mathrm{,}{\delta }_2\mathrm{,}{\delta }_3\right)\end{array}$ for any ${\delta }_1\mathrm{,}{\delta }_2\mathrm{,}{\delta }_3\mathrm{,}{\lambda }_1\mathrm{,}{\lambda }_2\mathrm{,}{\lambda }_3\ge 0$ .

6) If $\left[e\mathrm{,}f\right]\subseteq \left[a\mathrm{,}b\right]\mathrm{<mo>\; ;\; </mo>}\left[g\mathrm{,}h\right]\subseteq \left[c\mathrm{,}d\right]$ and $\left[\tilde{i}\mathrm{,}\tilde{j}\right]\preccurlyeq \left[\tilde{p}\mathrm{,}\tilde{q}\right]$ , then

${\omega }_{\left(\left[e\mathrm{,}f\right]\times \left[g\mathrm{,}h\right]\otimes \left[\tilde{i}\mathrm{,}\tilde{j}\right]\right)}\left(\tilde{f}\mathrm{,}{\delta }_1\mathrm{,}{\delta }_2\mathrm{,}{\delta }_3\right)\le {\omega }_{\left(\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\right)}\left(\tilde{f}\mathrm{,}{\delta }_1\mathrm{,}{\delta }_2\mathrm{,}{\delta }_3\right)$ .

Proof (6) According to the hypothesis,

$\begin{array}{l}\sup \{D\left(\tilde{f}\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)\mathrm{,}\tilde{f}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\right)\mathrm{<mo>\; ;\; </mo>}\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)\mathrm{,}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\in \left[e\mathrm{,}f\right]\\ \times \left[g\mathrm{,}h\right]\otimes \left[\tilde{i}\mathrm{,}\tilde{j}\right]\mathrm{,}\left|x_1-x_2\right|\le {\delta }_1\mathrm{,}\left|y_1-y_2\right|\le {\delta }_2\mathrm{,}{\Vert {\tilde{z}}_1\ominus {\tilde{z}}_2\Vert }_{{ℝ}_F}\le {\delta }_3\}\\ \le \sup \{D\left(\tilde{f}\left(x_1,y_1,{\tilde{z}}_1\right),\tilde{f}\left(x_2,y_2,{\tilde{z}}_2\right)\right);\left(x_1,y_1,{\tilde{z}}_1\right),\left(x_2,y_2,{\tilde{z}}_2\right)\in \left[a,b\right]\\ \times \left[c,d\right]\otimes \left[\tilde{p},\tilde{q}\right],\left|x_1-x_2\right|\le {\delta }_1,\left|y_1-y_2\right|\le {\delta }_2,{\Vert {\tilde{z}}_1\ominus {\tilde{z}}_2\Vert }_{{ℝ}_F}\le {\delta }_3\}\end{array}$

which is prove the relation.

We can prove similarly the other statements.

Lemma 3.3.

1) If $\tilde{f}$ and $\tilde{h}$ are Henstock triple integrable mappings and if $D\left(\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\mathrm{,}\tilde{h}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\right)$ is Lebesgue integrable, then

$\begin{array}{l}D\left(\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{h}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\right)\\ \le \left(L\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}D\left(\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\mathrm{,}\tilde{h}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{.}\end{array}$

2) Let $\tilde{f}\mathrm{<mo>\; :\; </mo>}\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\to {ℝ}_F$ be a Henstock triple integrabe bounded mapping.

Then, $\forall \left(u\mathrm{,}v\mathrm{,}\tilde{w}\right)\in \left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]$ , the function ${\phi }_{\left(u\mathrm{,}v\mathrm{,}\tilde{w}\right)}\mathrm{<mo>\; :\; </mo>}\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\to {ℝ}^{+}$ defined by ${\phi }_{\left(u\mathrm{,}v\mathrm{,}\tilde{w}\right)}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)=D\left(\tilde{f}\left(u\mathrm{,}v\mathrm{,}\tilde{w}\right)\mathrm{,}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\right)$ is Lebesgue integrable on $\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]$ .

Proof (2) If $\tilde{f}$ is Henstock integrable and bounded on $\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]$ , then it follows that ${\tilde{f}}_{-}^{\alpha }\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)$ and ${\tilde{f}}_{+}^{\alpha }\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)$ are Henstock triple integrable with $\alpha \in \left[\mathrm{0,1}\right]$ . Therefore, ${\tilde{f}}_{-}^{\alpha }\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)$ and ${\tilde{f}}_{+}^{\alpha }\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)$ are Lebesgue measurable and uniformly bounded $\forall \alpha \in \left[\mathrm{0,1}\right]$ , [7] . Moreover,

$\begin{array}{l}\phi \left(x,y,\tilde{z}\right)=D\left(\tilde{f}\left(x_1,y_1,{\tilde{z}}_1\right),\tilde{f}\left(x_2,y_2,{\tilde{z}}_2\right)\right)\\ =\underset{\alpha \in \left[0,1\right]}{\sup }\max \left\{\left|{\tilde{f}}_{-}^{\alpha }\left(x_1,y_1,{\tilde{z}}_1\right)-{\tilde{f}}_{-}^{\alpha }\left(x_2,y_2,{\tilde{z}}_2\right)\right|,\left|{\tilde{f}}_{+}^{\alpha }\left(x_1,y_1,{\tilde{z}}_1\right)-{\tilde{f}}_{+}^{\alpha }\left(x_2,y_2,{\tilde{z}}_2\right)\right|\right\}\\ =\underset{{\alpha }_n\in \left[0,1\right]}{\sup }\max \left\{\left|{\tilde{f}}_{-}^{{\alpha }_n}\left(x_1,y_1,{\tilde{z}}_1\right)-{\tilde{f}}_{-}^{{\alpha }_n}\left(x_2,y_2,{\tilde{z}}_2\right)\right|,\left|{\tilde{f}}_{+}^{{\alpha }_n}\left(x_1,y_1,{\tilde{z}}_1\right)-{\tilde{f}}_{+}^{r_n}\left(x_2,y_2,{\tilde{z}}_2\right)\right|\right\},\end{array}$

where the ${\alpha }_n\left(n\in \mathbb{N}\right)$ are the rational numbers in $\left[\mathrm{0,1}\right]$ . According to Lebesgue’s dominated convergence theorem, it follows that $\phi \left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)$ is Lebesgue integrable over $\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]$ and what completes the proof.

Keeping now three integrals we reach the following definitions.

Definition 3.4. A function $\tilde{f}\mathrm{<mo>\; :\; </mo>}\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\to {ℝ}_F$ is said to be $\left(L_1\mathrm{,}{L}_2\mathrm{,}{L}_3\right)$ Lipschitz if for any $\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)\mathrm{,}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\in \left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]$ ,

$D\left(\tilde{f}\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)\mathrm{,}\tilde{f}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\right)\le L_1\left|x_1-x_2\right|+L_2\left|y_1-y_2\right|+L_3{\Vert {\tilde{z}}_1\ominus {\tilde{z}}_2\Vert }_{{ℝ}_F}\mathrm{.}$

In order to introduce triple Simpson’s rule for evaluating FHTI, firstly we prove the following theorem.

Theorem 3.5. Let $\tilde{f}\mathrm{<mo>\; :\; </mo>}\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\to {ℝ}_F$ be a fuzzy Henstock integrable, bounded mapping. Then, for any subdivision $a=x_0<x_1<x_2<\cdots <x_m=b$ , $c=y_0<y_1<y_2<\cdots <y_n=d$ , $\tilde{p}={\tilde{z}}_0\prec {\tilde{z}}_1\prec {\tilde{z}}_2\prec \cdots \prec {\tilde{z}}_s=\tilde{q}$ and any points ${\xi }_i\in \left[x_{i-1}\mathrm{,}{x}_i\right]$ , ${\eta }_j\in \left[y_{j-1}\mathrm{,}{y}_j\right]$ , ${\tilde{\zeta }}_k\in \left[{\tilde{z}}_{k-1}\mathrm{,}{\tilde{z}}_k\right]$ we have

$\begin{array}{l}D(\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\stackrel{}{}\\ \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\odot \tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right))\\ \le {\displaystyle \underset{i=1}{\overset{m}{\sum }}}{\displaystyle \underset{j=1}{\overset{n}{\sum }}}{\displaystyle \underset{k=1}{\overset{s}{\sum }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\\ {\omega }_{\left(\left[x_{i-1}\mathrm{,}{x}_i\right]\times \left[y_{j-1}\mathrm{,}{y}_j\right]\otimes \left[{\tilde{z}}_{k-1}\mathrm{,}{\tilde{z}}_k\right]\right)}\left(\tilde{f}\mathrm{,}\left(x_i-x_{i-1}\right)\mathrm{,}\left(y_j-y_{j-1}\right)\mathrm{,}{\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\right)\mathrm{.}\end{array}$

Proof Since that the Henstock integral is additive related to interval [13] , hence,

$\begin{array}{l}D(\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\stackrel{}{}\\ \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\odot \tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right))\\ =D(\underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(FHTI\right){\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\\ \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\odot \tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right))\mathrm{.}\end{array}$

Since it’s clear that $\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{k}\text{d}\tilde{z}\text{d}y\text{d}x=\left(b-a\right)\left(d-c\right)\odot \left(\tilde{q}\ominus \tilde{p}\right)\otimes \tilde{k}$ for any fuzzy constant $\tilde{k}\in {ℝ}_F$ , we obtain

$\begin{array}{l}D(\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\stackrel{}{}\\ \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\odot \tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right))\\ =D(\underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(FHTI\right){\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\\ \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(FHTI\right){\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}\tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\text{d}\tilde{z}\text{d}y\text{d}x)\mathrm{.}\end{array}$

By the fourth property of the theorem 3.2, we have

$\begin{array}{l}D(\underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(FHTI\right){\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\\ \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(FHTI\right){\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}\tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\text{d}\tilde{z}\text{d}y\text{d}x)\\ \le \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}D(\left(FHTI\right){\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}\tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\\ \left(FHTI\right){\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}\tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\text{d}\tilde{z}\text{d}y\text{d}x)\end{array}$

Since the functions $D\left(\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\mathrm{,}\tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\right)$ are Lebesgue integrable for $i=1,\cdots ,m;j=1,\cdots ,n$ and $k=\mathrm{1,}\cdots \mathrm{,}s$ from lemma 3.3 we have

$\begin{array}{l}D(\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\stackrel{}{}\\ \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\odot \tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right))\\ \le \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(L\right){\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}D\left(\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\mathrm{,}\tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{.}\end{array}$

From the first property of the theorem 3.2 applied to each of the above integrals we have

$\begin{array}{l}D(\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\stackrel{}{}\\ \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\odot \tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right))\\ \le \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(L\right){\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}{\omega }_{\left(\left[x_{i-1}\mathrm{,}{x}_i\right]\times \left[y_{j-1}\mathrm{,}{y}_j\right]\otimes \left[{\tilde{z}}_{k-1}\mathrm{,}{\tilde{z}}_k\right]\right)}\\ \left(\tilde{f}\mathrm{,}\left(x_i-x_{i-1}\right)\mathrm{,}\left(y_j-y_{j-1}\right)\mathrm{,}{\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\right)\text{d}\tilde{z}\text{d}y\text{d}x\\ ={\displaystyle \underset{i=1}{\overset{m}{\sum }}}{\displaystyle \underset{j=1}{\overset{n}{\sum }}}{\displaystyle \underset{k=1}{\overset{s}{\sum }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\\ {\omega }_{\left(\left[x_{i-1}\mathrm{,}{x}_i\right]\times \left[y_{j-1}\mathrm{,}{y}_j\right]\otimes \left[{\tilde{z}}_{k-1}\mathrm{,}{\tilde{z}}_k\right]\right)}\left(\tilde{f}\mathrm{,}\left(x_i-x_{i-1}\right)\mathrm{,}\left(y_j-y_{j-1}\right)\mathrm{,}{\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\right)\mathrm{,}\end{array}$

which completes the proof.

Corollary 3.6. Let $\tilde{f}\mathrm{<mo>\; :\; </mo>}\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\to {ℝ}_F$ be a Henstock triple integrable, bounded mapping. Then, for $m=n=s=2$ we have

$\begin{array}{l}D(\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\stackrel{}{}\\ \underset{i=1}{\overset{2}{{\displaystyle \oplus }}}\underset{j=1}{\overset{2}{{\displaystyle \oplus }}}\underset{k=1}{\overset{2}{{\displaystyle \oplus }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\odot \tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right))\\ \le \left(\alpha -a\right)\left(\beta -c\right){\Vert \tilde{\gamma }\ominus \tilde{p}\Vert }_{{ℝ}_F}{\omega }_{\left[a\mathrm{,}\alpha \right]\times \left[c\mathrm{,}\beta \right]\otimes \left[\tilde{p}\mathrm{,}\tilde{\gamma }\right]}\left(\tilde{f}\mathrm{,}\left(\alpha -a\right)\mathrm{,}\left(\beta -c\right)\mathrm{,}{\Vert \tilde{\gamma }\ominus \tilde{p}\Vert }_{{ℝ}_F}\right)\\ +\left(\alpha -a\right)\left(\beta -c\right){\Vert \tilde{q}\ominus \tilde{\gamma }\Vert }_{{ℝ}_F}{\omega }_{\left[a\mathrm{,}\alpha \right]\times \left[c\mathrm{,}\beta \right]\otimes \left[\tilde{\gamma }\mathrm{,}\tilde{q}\right]}\left(\tilde{f}\mathrm{,}\left(\alpha -a\right)\mathrm{,}\left(\beta -c\right)\mathrm{,}{\Vert \tilde{q}\ominus \tilde{\gamma }\Vert }_{{ℝ}_F}\right)\\ +\left(\alpha -a\right)\left(d-\beta \right){\Vert \tilde{\gamma }\ominus \tilde{p}\Vert }_{{ℝ}_F}{\omega }_{\left[a\mathrm{,}\alpha \right]\times \left[\beta \mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{\gamma }\right]}\left(\tilde{f}\mathrm{,}\left(\alpha -a\right)\mathrm{,}\left(d-\beta \right)\mathrm{,}{\Vert \tilde{\gamma }\ominus \tilde{p}\Vert }_{{ℝ}_F}\right)\end{array}$

$\begin{array}{l}+\left(\alpha -a\right)\left(d-\beta \right){\Vert \tilde{q}\ominus \tilde{\gamma }\Vert }_{{ℝ}_F}{\omega }_{\left[a\mathrm{,}\alpha \right]\times \left[\beta \mathrm{,}d\right]\otimes \left[\tilde{\gamma }\mathrm{,}\tilde{q}\right]}\left(\tilde{f}\mathrm{,}\left(\alpha -a\right)\mathrm{,}\left(d-\beta \right)\mathrm{,}{\Vert \tilde{q}\ominus \tilde{\gamma }\Vert }_{{ℝ}_F}\right)\\ +\left(b-\alpha \right)\left(\beta -c\right){\Vert \tilde{\gamma }\ominus \tilde{p}\Vert }_{{ℝ}_F}{\omega }_{\left[\alpha \mathrm{,}b\right]\times \left[c\mathrm{,}\beta \right]\otimes \left[\tilde{p}\mathrm{,}\tilde{\gamma }\right]}\left(\tilde{f}\mathrm{,}\left(b-\alpha \right)\mathrm{,}\left(\beta -c\right){\Vert \tilde{\gamma }\ominus \tilde{p}\Vert }_{{ℝ}_F}\right)\\ +\left(b-\alpha \right)\left(\beta -c\right){\Vert \tilde{q}\ominus \tilde{\gamma }\Vert }_{{ℝ}_F}{\omega }_{\left[\alpha \mathrm{,}b\right]\times \left[c\mathrm{,}\beta \right]\otimes \left[\tilde{\gamma }\mathrm{,}\tilde{q}\right]}\left(\tilde{f}\mathrm{,}\left(b-\alpha \right)\mathrm{,}\left(\beta -c\right){\Vert \tilde{q}\ominus \tilde{\gamma }\Vert }_{{ℝ}_F}\right)\\ +\left(b-\alpha \right)\left(d-\beta \right){\Vert \tilde{\gamma }\ominus \tilde{p}\Vert }_{{ℝ}_F}{\omega }_{\left[\alpha \mathrm{,}b\right]\times \left[\beta \mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{\gamma }\right]}\left(\tilde{f}\mathrm{,}\left(b-\alpha \right)\mathrm{,}\left(d-\beta \right)\mathrm{,}{\Vert \tilde{\gamma }\ominus \tilde{p}\Vert }_{{ℝ}_F}\right)\\ +\left(b-\alpha \right)\left(d-\beta \right){\Vert \tilde{q}\ominus \tilde{\gamma }\Vert }_{{ℝ}_F}{\omega }_{\left[\alpha \mathrm{,}b\right]\times \left[\beta \mathrm{,}d\right]\otimes \left[\tilde{\gamma }\mathrm{,}\tilde{q}\right]}\left(\tilde{f}\mathrm{,}\left(b-\alpha \right)\mathrm{,}\left(d-\beta \right)\mathrm{,}{\Vert \tilde{q}\ominus \tilde{\gamma }\Vert }_{{ℝ}_F}\right)\mathrm{,}\end{array}$

for any $\alpha \in \left[a\mathrm{,}b\right]\mathrm{,}\beta \in \left[c\mathrm{,}d\right]$ and $\tilde{\gamma }\in \left[\tilde{p}\mathrm{,}\tilde{q}\right]$ , $\left(u\mathrm{,}v\mathrm{,}\tilde{w}\right)\in \left[a\mathrm{,}\alpha \right]\times \left[c\mathrm{,}\beta \right]\otimes \left[\tilde{p}\mathrm{,}\tilde{\gamma }\right]$ and $\left(u^{\prime }\mathrm{,}{v}^{\prime }\mathrm{,}{\tilde{w}}^{\prime }\right)\in \left[\alpha \mathrm{,}b\right]\times \left[\beta \mathrm{,}d\right]\otimes \left[\tilde{\gamma }\mathrm{,}\tilde{q}\right]$ where ${\xi }_1=u$ , ${\xi }_2=u^{\prime }$ ; ${\eta }_1=v$ , ${\eta }_2=v^{\prime }$ ; ${\tilde{\zeta }}_1=\tilde{w}$ , ${\tilde{\zeta }}_2={\tilde{w}}^{\prime }$ .

Proof It’s clear that for $m=2,n=2$ and $s=2$ in the theorem 3.5 the inequality stated above is obtained.

The corollary below gives the fuzzy variant of Simpson’s threefold rule with a new error bound.

Corollary 3.7. Let $\tilde{f}\mathrm{<mo>\; :\; </mo>}\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\mapsto {ℝ}_F$ be a Henstock triple integrable, bounded mapping. Then,

$\begin{array}{l}D(\left(ITFHK\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\frac{\left(b-a\right)\left(d-c\right){\Vert \tilde{q}\ominus \tilde{p}\Vert }_{{ℝ}_F}}{216}\\ \odot \{\left[\left(\tilde{f}\left(a\mathrm{,}c\mathrm{,}\tilde{p}\right)\oplus \tilde{f}\left(a\mathrm{,}c\mathrm{,}\tilde{q}\right)\right)\oplus 4\odot \left(\tilde{f}\left(\frac{a+b}{2}\mathrm{,}c\mathrm{,}\tilde{p}\right)\oplus \tilde{f}\left(\frac{a+b}{2}\mathrm{,}c\mathrm{,}\tilde{q}\right)\right)\oplus \left(\tilde{f}\left(b\mathrm{,}c\mathrm{,}\tilde{p}\right)\oplus \tilde{f}\left(b\mathrm{,}c\mathrm{,}\tilde{q}\right)\right)\right]\\ \oplus 4\odot \{\left[\tilde{f}\left(a\mathrm{,}c\mathrm{,}\frac{\tilde{p}\oplus \tilde{q}}{2}\right)\oplus 4\odot \tilde{f}\left(\frac{a+b}{2}\mathrm{,}c\mathrm{,}\frac{\tilde{p}\oplus \tilde{q}}{2}\right)\oplus \tilde{f}\left(b\mathrm{,}c\mathrm{,}\frac{\tilde{p}\oplus \tilde{q}}{2}\right)\right]\\ \oplus \left[\tilde{f}\left(a\mathrm{,}\frac{c+d}{2}\mathrm{,}\tilde{p}\right)\oplus 4\odot \tilde{f}\left(\frac{a+b}{2}\mathrm{,}\frac{c+d}{2}\mathrm{,}\tilde{p}\right)\oplus \tilde{f}\left(b\mathrm{,}\frac{c+d}{2}\mathrm{,}\tilde{p}\right)\right]\\ \oplus \left[\tilde{f}\left(a\mathrm{,}\frac{c+d}{2}\mathrm{,}\frac{\tilde{p}\oplus \tilde{q}}{2}\right)\oplus 4\odot \tilde{f}\left(\frac{a+b}{2}\mathrm{,}\frac{c+d}{2}\mathrm{,}\frac{\tilde{p}\oplus \tilde{q}}{2}\right)\oplus \tilde{f}\left(b\mathrm{,}\frac{c+d}{2}\mathrm{,}\frac{\tilde{p}\oplus \tilde{q}}{2}\right)\right]\end{array}$

$\begin{array}{l}\oplus \left[\tilde{f}\left(a\mathrm{,}\frac{c+d}{2}\mathrm{,}\tilde{q}\right)\oplus 4\odot \tilde{f}\left(\frac{a+b}{2}\mathrm{,}\frac{c+d}{2}\mathrm{,}\tilde{q}\right)\oplus \tilde{f}\left(b\mathrm{,}\frac{c+d}{2}\mathrm{,}\tilde{q}\right)\right]\\ \oplus \left[\tilde{f}\left(a\mathrm{,}d\mathrm{,}\frac{\tilde{p}\oplus \tilde{q}}{2}\right)\oplus 4\odot \tilde{f}\left(\frac{a+b}{2}\mathrm{,}d\mathrm{,}\frac{\tilde{p}\oplus \tilde{q}}{2}\right)\oplus \tilde{f}\left(b\mathrm{,}d\mathrm{,}\frac{\tilde{p}\oplus \tilde{q}}{2}\right)\right]\}\\ \oplus \left[\left(\tilde{f}\left(a\mathrm{,}d\mathrm{,}\tilde{p}\right)\oplus \tilde{f}\left(a\mathrm{,}d\mathrm{,}\tilde{q}\right)\right)\oplus 4\odot \left(\tilde{f}\left(\frac{a+b}{2}\mathrm{,}d\mathrm{,}\tilde{p}\right)\oplus \tilde{f}\left(\frac{a+b}{2}\mathrm{,}d\mathrm{,}\tilde{q}\right)\right)\oplus \left(\tilde{f}\left(b\mathrm{,}d\mathrm{,}\tilde{p}\right)\oplus \tilde{f}\left(b\mathrm{,}d\mathrm{,}\tilde{q}\right)\right)\right]\})\\ \le 27\left(b-a\right)\left(d-c\right){\Vert \tilde{q}\ominus \tilde{p}\Vert }_{{ℝ}_F}{\omega }_{\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]}\left(\tilde{f}\mathrm{,}\frac{b-a}{6}\mathrm{,}\frac{d-c}{6}\mathrm{,}\frac{{\Vert \tilde{q}\ominus \tilde{p}\Vert }_{{ℝ}_F}}{6}\right)\end{array}$

Proof

This inequality follows from the previous corollary by setting $\alpha =\frac{5a+b}{6}$ , $\beta =\frac{5c+d}{6}$ , $\tilde{\gamma }=\frac{5\odot \tilde{p}\oplus \tilde{q}}{6}$ , $u=a$ , $v=\frac{a+b}{2}$ , $w=b$ , $u^{\prime }=c$ , $v^{\prime }=\frac{c+d}{2}$ , $w^{\prime }=d$ , $\tilde{u}=\tilde{p}$ , $\tilde{v}=\frac{\tilde{p}\oplus \tilde{q}}{2}$ and $\tilde{w}=\tilde{q}$ .

Theorem 3.8. Let $\tilde{f}\mathrm{<mo>\; :\; </mo>}\left[a\mathrm{,}b\right]\times \left[c\mathrm{,}d\right]\otimes \left[\tilde{p}\mathrm{,}\tilde{q}\right]\to {ℝ}_F$ be a Lipschitz mapping with the constants $L_1\mathrm{,}{L}_2$ and $L_3$ . Then, for any subdivision

${\Delta }_m:a=x_0<x_1<x_2<\cdots <x_m=b$ , ${\Delta }_n:c=y_0<y_1<y_2<\cdots <y_n=d$

and

${\Delta }_s:\tilde{p}={\tilde{z}}_0\prec {\tilde{z}}_1\prec {\tilde{z}}_2\prec \cdots \prec {\tilde{z}}_s=\tilde{q}$ . $\forall {\xi }_i\in \left[x_{i-1}\mathrm{,}{x}_i\right]$ , $i=\mathrm{1,2,}\cdots \mathrm{,}m$ ;

${\eta }_j\in \left[y_{j-1}\mathrm{,}{y}_j\right]$ , $j=\mathrm{1,2,}\cdots \mathrm{,}n$ and ${\tilde{\zeta }}_k\in \left[{\tilde{z}}_{k-1}\mathrm{,}{\tilde{z}}_k\right]$ , $k=\mathrm{1,2,}\cdots \mathrm{,}s$ ;

we have

$\begin{array}{l}D(\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\stackrel{}{}\\ \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\odot \tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right))\\ \le {\displaystyle \underset{i=1}{\overset{m}{\sum }}}{\displaystyle \underset{j=1}{\overset{n}{\sum }}}{\displaystyle \underset{k=1}{\overset{s}{\sum }}}(L_1\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}{\left(x_i-x_{i-1}\right)}^2\\ +L_2\left(x_i-x_{i-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}{\left(y_j-y_{j-1}\right)}^2\\ +L_3\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}^2)\mathrm{.}\end{array}$

Proof Similar to the proof of theorem 3.5 we have

$\begin{array}{l}D(\left(FHTI\right){\displaystyle {\int }_a^b}{\displaystyle {\int }_c^d}{\displaystyle {\int }_{\tilde{p}}^{\tilde{q}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{,}\stackrel{}{}\\ \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\odot \tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right))\\ \le \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(L\right){\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}D\left(\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\mathrm{,}\tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{.}\end{array}$

We obtain by the definition of a Lipschitz mapping

$\begin{array}{l}\underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}\left(L\right){\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}D\left(\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\mathrm{,}\tilde{f}\left({\xi }_i\mathrm{,}{\eta }_j\mathrm{,}{\tilde{\zeta }}_k\right)\right)\text{d}\tilde{z}\text{d}y\text{d}x\mathrm{.}\\ \le \underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}(L_1{\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}\left|x-{\xi }_i\right|\text{d}\tilde{z}\text{d}y\text{d}x\\ \oplus L_2{\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}\left|y-{\eta }_j\right|\text{d}\tilde{z}\text{d}y\text{d}x\oplus L_3{\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}{\Vert \tilde{z}\ominus {\tilde{\zeta }}_k\Vert }_{{ℝ}_F}\text{d}\tilde{z}\text{d}y\text{d}x)\end{array}$

It follows by direct computation that

$\begin{array}{l}\underset{i=1}{\overset{m}{{\displaystyle \oplus }}}\underset{j=1}{\overset{n}{{\displaystyle \oplus }}}\underset{k=1}{\overset{s}{{\displaystyle \oplus }}}(L_1{\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}\left|x-{\xi }_i\right|\text{d}\tilde{z}\text{d}y\text{d}x\\ \oplus L_2{\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}\left|y-{\eta }_j\right|\text{d}\tilde{z}\text{d}y\text{d}x\oplus L_3{\displaystyle {\int }_{x_{i-1}}^{x_i}}{\displaystyle {\int }_{y_{j-1}}^{y_j}}{\displaystyle {\int }_{{\tilde{z}}_{k-1}}^{{\tilde{z}}_k}}{\Vert \tilde{z}\ominus {\tilde{\zeta }}_k\Vert }_{{ℝ}_F}\text{d}\tilde{z}\text{d}y\text{d}x)\\ =\frac{1}{2}{\displaystyle \underset{i=1}{\overset{m}{\sum }}}{\displaystyle \underset{j=1}{\overset{n}{\sum }}}{\displaystyle \underset{k=1}{\overset{s}{\sum }}}\mathrm{<mo\; stretchy="false">\; (\; </mo>}{L}_1\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\left[{\left(x_i-{\xi }_i\right)}^2-{\left(x_{i-1}-{\xi }_i\right)}^2\right]\\ +L_2\left(x_i-x_{i-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}\left[{\left(y_j-{\eta }_j\right)}^2-{\left(y_{j-1}-{\eta }_j\right)}^2\right]\\ +L_3\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right)\left[{\Vert {\tilde{z}}_k\ominus {\tilde{\zeta }}_k\Vert }_{{ℝ}_F}^2-{\Vert {\tilde{z}}_{k-1}\ominus {\tilde{\zeta }}_k\Vert }_{{ℝ}_F}^2\right])\end{array}$

$\begin{array}{l}\le \frac{1}{2}{\displaystyle \underset{i=1}{\overset{m}{\sum }}}{\displaystyle \underset{j=1}{\overset{n}{\sum }}}{\displaystyle \underset{k=1}{\overset{s}{\sum }}}\mathrm{<mo\; stretchy="false">\; (\; </mo>}{L}_1\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}{\left(x_i-x_{i-1}\right)}^2\\ +L_2\left(x_i-x_{i-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}{\left(y_j-y_{j-1}\right)}^2\\ +L_3\left(x_i-x_{i-1}\right)\left(y_j-y_{j-1}\right){\Vert {\tilde{z}}_k\ominus {\tilde{z}}_{k-1}\Vert }_{{ℝ}_F}^2)\mathrm{.}\end{array}$

4. Numerical Example

Let $\tilde{f}\mathrm{<mo>\; :\; </mo>}\left[\mathrm{0,1}\right]\times \left[\mathrm{1,2}\right]\otimes \left[\tilde{1}\mathrm{,}\tilde{2}\right]\to {ℝ}_F$ , $\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)=x^2+3y\oplus \tilde{z}$ where $\tilde{1}=\left(0,1,2\right)$ ; $\tilde{2}=\left(1,2,3\right)$ ; $\tilde{z}=\left(z-1,z,z+1\right)$ , and where $\left(a_1\mathrm{,}{a}_2\mathrm{,}{a}_3\right)$ is a triangular fuzzy number such that

$\mu \left(x\right)=\{\begin{array}{ll}\frac{x-a_1}{a_2-a_1}\hfill & a_1\le x\le a_2\hfill \\ \frac{a_3-x}{a_3-a_2}\hfill & a_2\le x\le a_3\hfill \\ 0\hfill & \text{otherwise}\hfill \end{array}$

We must compute the integral

$\left(FHTI\right){\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\tilde{1}}^{\tilde{2}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x$

numerically.

Firstly let’s define the different α-cuts below:

${\left(\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\right)}^{\alpha }=\left[x^2+3y+z-1+\alpha \mathrm{,}{x}^2+3y+z+1-\alpha \right]\mathrm{,}$

${\tilde{z}}^{\alpha }=\left[z-1+\alpha \mathrm{,}z+1-\alpha \right]\mathrm{,}$

${\tilde{1}}^{\alpha }=\left[\alpha \mathrm{,2}-\alpha \right]\mathrm{,}$

${\tilde{2}}^{\alpha }=\left[1+\alpha \mathrm{,3}-\alpha \right]\mathrm{,}\text{and}$

${\tilde{I}}^{\alpha }={\left(\left(FHTI\right){\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\tilde{1}}^{\tilde{2}}}\tilde{f}\left(x\mathrm{,}y\mathrm{,}\tilde{z}\right)\text{d}\tilde{z}\text{d}y\text{d}x\right)}^{\alpha }=\left[\min S\mathrm{,}\max S\right]\mathrm{.}$

where

$\begin{array}{l}S=\{{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{1+\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x;{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{1+\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x;\\ {\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{3-\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x;{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{3-\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x;\\ {\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{1+\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x;{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{1+\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x;\\ {\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{3-\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x;{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{3-\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x;\end{array}$

$\begin{array}{l}{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{1+\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x;{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{1+\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x;\\ {\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{3-\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x;{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{3-\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x;\\ {\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{1+\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x;{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{1+\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x;\\ {\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{3-\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x;{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{3-\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x\}.\end{array}$

We must calculate these sixteen integrals while noting that they are two by two equal.

We remark that

$\begin{array}{l}D\left(\tilde{f}\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)\mathrm{,}\tilde{f}\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\right)\\ =\underset{\alpha \in \left[\mathrm{0,1}\right]}{\sup }\max \{\left|{\tilde{f}}_{-}^{\alpha }\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)-{\tilde{f}}_{-}^{\alpha }\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\right|\mathrm{,}\left|{\tilde{f}}_{+}^{\alpha }\left(x_1\mathrm{,}{y}_1\mathrm{,}{\tilde{z}}_1\right)-{\tilde{f}}_{+}^{\alpha }\left(x_2\mathrm{,}{y}_2\mathrm{,}{\tilde{z}}_2\right)\right|\\ =\underset{\alpha \in \left[\mathrm{0,1}\right]}{\sup }\max \{\left|\left(x_1^2-x_2^2\right)+3\left(y_1-y_2\right)+\left(z_1-z_2\right)+\left(\alpha -1\right)\right|\mathrm{,}\\ \left|\left(x_1^2-x_2^2\right)+3\left(y_1-y_2\right)+\left(z_1-z_2\right)+\left(1-\alpha \right)\right|\}\end{array}$

$\begin{array}{l}\le \underset{\alpha \in \left[\mathrm{0,1}\right]}{\sup }\max \left\{\left|x_1-x_2\right|\left(\left|x_1+x_2\right|\mathrm{,}\left|x_1+x_2\right|\right)\right\}+\underset{\alpha \in \left[\mathrm{0,1}\right]}{\sup }\max \left\{\left|y_1-y_2\right|\left(\mathrm{3,3}\right)\right\}\\ +\underset{\alpha \in \left[\mathrm{0,1}\right]}{\sup }\max \left\{\left|z_1-z_2\right|\left(\mathrm{1,1}\right)\right\}\\ \le 2\left|x_1-x_2\right|+3\left|y_1-y_2\right|+\left|z_1-z_2\right|\mathrm{.}\end{array}$

i.e. $\tilde{f}$ is a Lipschitz mapping with $L_1=2,L_2=3$ and $L_3=1$ .

Indeed,

${\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{1+\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{1+\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x,$ (4.1)

because $\text{d}{\tilde{z}}_{-}^{\alpha }=\text{d}{\tilde{z}}_{+}^{\alpha }=\text{d}z$

From where, 4.1 becomes

$\begin{array}{c}{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{1+\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}z\text{d}y\text{d}x={\displaystyle {\int }_{\alpha }^{1+\alpha }}{\displaystyle {\int }_1^2}{\displaystyle {\int }_0^1}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}x\text{d}y\text{d}z\\ ={\displaystyle {\int }_{\alpha }^{1+\alpha }}{\displaystyle {\int }_1^2}{\displaystyle {\int }_0^1}\left(x^2+3y+z-1+\alpha \right)\text{d}x\text{d}y\text{d}z\\ ={\displaystyle {\int }_{\alpha }^{1+\alpha }}{\displaystyle {\int }_1^2}{\left[\frac{x^3}{3}+\left(3y+z-1+\alpha \right)\right]}_{x=0}^1\text{d}y\text{d}z\\ ={\displaystyle {\int }_{\alpha }^{1+\alpha }}{\displaystyle {\int }_1^2}\left(3y-\frac{2}{3}+\alpha +z\right)\text{d}y\text{d}z\end{array}$

$\begin{array}{l}={\displaystyle {\int }_{\alpha }^{1+\alpha }}{\left[\frac{3}{2}{y}^2+\left(-\frac{2}{3}+\alpha +z\right)y\right]}_{y=1}^2\text{d}z\\ ={\displaystyle {\int }_{\alpha }^{1+\alpha }}\left[\left(\frac{23}{6}+\alpha \right)+z\right]\text{d}z\\ ={\left[\left(\frac{23}{6}+\alpha \right)z+\frac{1}{2}{z}^2\right]}_{z=\alpha }^{1+\alpha }\\ =\frac{13}{3}+2\alpha .\end{array}$

Proceeding in the same way for the other integrals we obtain

$\begin{array}{c}{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{3-\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{3-\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x\\ ={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{3-\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}z\text{d}y\text{d}x\\ =-2{\alpha }^2-\frac{23}{3}\alpha +16\end{array}$

$\begin{array}{c}{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{1+\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{1+\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x\\ ={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{1+\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}z\text{d}y\text{d}x\\ =\frac{6{\alpha }^2+29\alpha -16}{3}\end{array}$

$\begin{array}{c}{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{3-\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{3-\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x\\ ={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{3-\alpha }}{\tilde{f}}_{-}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}z\text{d}y\text{d}x\\ =\frac{19}{3}\end{array}$

$\begin{array}{c}{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{1+\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{1+\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x\\ ={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{1+\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}z\text{d}y\text{d}x\\ =\frac{19}{3}\end{array}$

$\begin{array}{c}{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{3-\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{3-\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x\\ ={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{\alpha }^{3-\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}z\text{d}y\text{d}x\\ =2{\alpha }^2-\frac{53}{3}\alpha +22\end{array}$

$\begin{array}{c}{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{1+\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{1+\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x\\ ={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{1+\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}z\text{d}y\text{d}x\\ =\frac{-6{\alpha }^2+47\alpha -22}{3}\end{array}$

$\begin{array}{c}{\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{3-\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{-}^{\alpha }\text{d}y\text{d}x={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{3-\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}{\tilde{z}}_{+}^{\alpha }\text{d}y\text{d}x\\ ={\displaystyle {\int }_0^1}{\displaystyle {\int }_1^2}{\displaystyle {\int }_{2-\alpha }^{3-\alpha }}{\tilde{f}}_{+}^{\alpha }\left(x,y,\tilde{z}\right)\text{d}z\text{d}y\text{d}x\\ =\frac{25}{3}-2\alpha .\end{array}$

Finally S becomes

$\begin{array}{l}S=\{\frac{13}{3}+2\alpha \mathrm{,}-2{\alpha }^2-\frac{23}{3}\alpha +\mathrm{16,}\frac{6{\alpha }^2+29\alpha -16}{3}\mathrm{,}\frac{19}{3}\mathrm{,2}{\alpha }^2-\frac{53}{3}\alpha +\mathrm{22,}\\ \frac{-6{\alpha }^2+47\alpha -22}{3}\mathrm{,}\frac{25}{3}-2\alpha \}\mathrm{.}\end{array}$

We have for $\alpha =1$ that ${\tilde{I}}^1=\left[\mathrm{6,3333\; <mo>\; ;\; </mo>6,3333}\right]$ , for $\alpha =0,9$

$S=\left\{\mathrm{6,1333\; <mo>\; ;\; </mo>7,4800\; <mo>\; ;\; </mo>4,9867\; <mo>\; ;\; </mo>6,3333\; <mo>\; ;\; </mo>7,7200\; <mo>\; ;\; </mo>5,1467\; <mo>\; ;\; </mo>6,5333}\right\}$

this is shows that ${\tilde{I}}^{\mathrm{0,9}}=\left[\mathrm{4,9867\; <mo>\; ;\; </mo>7,7200}\right]$ .

In Table 1 below we summarize the results for the different values of $\alpha$ .

5. Conclusion

In this article, it was a question of introducing and evaluating the fuzzy triple integral of Henstock-Kurzweil whose one of the variables is fuzzy, using Simpson’s rule and Hausdorff distance. Therefore, we compute the integral on a quasi-fuzzy parallelepiped domain. In this direction, we have established and demonstrated a theorem which shows the upper limit of the distance between the exact values and approximate. For the rest, it would be possible to do the same analyzes for the integral of a fuzzy function with three variables of which two (or three) are fuzzy on a fuzzy parallelepipedal domain.