Loading [MathJax]/jax/output/SVG/jax.js
Research article

A higher-order numerical scheme for system of two-dimensional nonlinear fractional Volterra integral equations with uniform accuracy

  • Received: 31 January 2023 Revised: 10 March 2023 Accepted: 15 March 2023 Published: 03 April 2023
  • MSC : 65R20, 65D30, 65L12

  • We give a modified block-by-block method for the nonlinear fractional order Volterra integral equation system by using quadratic Lagrangian interpolation based on the classical block-by-block method. The core of the method is that we divide its domain into a series of subdomains, that is, block it, and use piecewise quadratic Lagrangian interpolation on each subdomain to approximate κ(x,y,s,r,u(s,r)). Our proposed method has uniform accuracy and its convergence order is O(h4αx+h4βy). We give a strict proof for the error analysis of the method, and give several numerical examples to verify the correctness of the theoretical analysis.

    Citation: Ziqiang Wang, Qin Liu, Junying Cao. A higher-order numerical scheme for system of two-dimensional nonlinear fractional Volterra integral equations with uniform accuracy[J]. AIMS Mathematics, 2023, 8(6): 13096-13122. doi: 10.3934/math.2023661

    Related Papers:

    [1] Boina Anil Kumar, S. K. Paikray, Hemen Dutta . Cost optimization model for items having fuzzy demand and deterioration with two-warehouse facility under the trade credit financing. AIMS Mathematics, 2020, 5(2): 1603-1620. doi: 10.3934/math.2020109
    [2] Messaoud Bounkhel, Lotfi Tadj . Model predictive control based integration of pricing and production planning. AIMS Mathematics, 2024, 9(1): 2282-2307. doi: 10.3934/math.2024113
    [3] Yuanlin Ding . Existence and stability analysis of solutions for periodic conformable differential systems with non-instantaneous impulses. AIMS Mathematics, 2025, 10(2): 4040-4066. doi: 10.3934/math.2025188
    [4] Gunduz Caginalp . Fat tails arise endogenously from supply/demand, with or without jump processes. AIMS Mathematics, 2021, 6(5): 4811-4846. doi: 10.3934/math.2021283
    [5] Velusamy Kavitha, Mani Mallika Arjunan, Dumitru Baleanu . Non-instantaneous impulsive fractional-order delay differential systems with Mittag-Leffler kernel. AIMS Mathematics, 2022, 7(5): 9353-9372. doi: 10.3934/math.2022519
    [6] Kui Liu . Stability analysis for $ (\omega, c) $-periodic non-instantaneous impulsive differential equations. AIMS Mathematics, 2022, 7(2): 1758-1774. doi: 10.3934/math.2022101
    [7] Mohamed Adel, M. Elsaid Ramadan, Hijaz Ahmad, Thongchai Botmart . Sobolev-type nonlinear Hilfer fractional stochastic differential equations with noninstantaneous impulsive. AIMS Mathematics, 2022, 7(11): 20105-20125. doi: 10.3934/math.20221100
    [8] Huanhuan Zhang, Jia Mu . Periodic problem for non-instantaneous impulsive partial differential equations. AIMS Mathematics, 2022, 7(3): 3345-3359. doi: 10.3934/math.2022186
    [9] Thabet Abdeljawad, Pshtiwan Othman Mohammed, Hari Mohan Srivastava, Eman Al-Sarairah, Artion Kashuri, Kamsing Nonlaopon . Some novel existence and uniqueness results for the Hilfer fractional integro-differential equations with non-instantaneous impulsive multi-point boundary conditions and their application. AIMS Mathematics, 2023, 8(2): 3469-3483. doi: 10.3934/math.2023177
    [10] Ahmed Abouelregal, Meshari Alesemi, Husam Alfadil . Thermoelastic reactions in a long and thin flexible viscoelastic cylinder due to non-uniform heat flow under the non-Fourier model with fractional derivative of two different orders. AIMS Mathematics, 2022, 7(5): 8510-8533. doi: 10.3934/math.2022474
  • We give a modified block-by-block method for the nonlinear fractional order Volterra integral equation system by using quadratic Lagrangian interpolation based on the classical block-by-block method. The core of the method is that we divide its domain into a series of subdomains, that is, block it, and use piecewise quadratic Lagrangian interpolation on each subdomain to approximate κ(x,y,s,r,u(s,r)). Our proposed method has uniform accuracy and its convergence order is O(h4αx+h4βy). We give a strict proof for the error analysis of the method, and give several numerical examples to verify the correctness of the theoretical analysis.



    Inventory control is a critical activity for a company organization in the domain of inventory management. Various realistic assumptions about inventory control characteristics like client demand, degradation rate, preservation, backordering (partial or complete) and so on are essential for successful inventory management. In addition, numerous business schemes (for example, advance payment, trade credit and price discount) are necessary to construct an accurate and realistic inventory model.

    Proper control of inventory is an essential task for a business organization. For proper inventory management, different realistic assumptions regarding customer demand, deterioration rate, preservation, stock-out scenario, etc., need to be properly adopted during inventory analysis. Furthermore, to ensure the applicability of inventory modeling in the practical competitive businesses, many practical business strategies (viz., advance payment, permissible delay in payment and price discount) are needed to be taken into account when developing the inventory model. In inventory modeling, client demand is quite crucial, as it controls the main characteristics for any type of product management scheme. Customer demand can be influenced by various factors, such as time, stock, product quality and price. Product price has a considerable effect on demand because the demand of a product always changes as the selling price of the product increases or decreases. Consequently, to introduce a product to the market, a producer must be aware of the client's needs, as well as the product's price.

    The success or failure of any practitioner in business largely depends on how they manage the stock of items. This activity becomes indispensable, especially when the items are perishable, because deteriorated items do not bring any revenue, as they increase the deterioration cost/opportunity cost for the practitioner. Moreover, not all items (for example, vegetables, fish and fruits) perish from the moment of storage, but after some time from the moment of storage, so they are known as non-immediate/non-instantaneous perishable/decaying items in stock analysis. As a result, the ordering decision for perishable products ultimately measures the success or failure of any practitioner to a great extent. On the other hand, suppliers/manufacturers of perishable items are very aware of the orders of their customers (retailers), as any order cancellation can cause a dramatic loss to their business. Apparently, they impose a prepayment strategy for their customers (retailers) to control the order cancellations. Under this prepayment strategy, the practitioner pays the supplier/manufacturer a percentage of the total acquisition price in advance, and the rest of the acquisition price is paid at the time of receiving the lot. However, the practitioner incurs an opportunity cost against the prepayment amount because they cannot earn any revenue before receiving the goods. Consequently, the following research question is raised: Is there any effect of prepayment strategy on the practitioner's order policy for non-instantaneous decay items? Under this prepayment strategy, the inventory decisions for the perishable items become much more complicated when the client demand varies with respect to a variable selling price. To the best of the authors' knowledge, there is no research work on inventory control to help practitioners by providing the optimal order policy for non-instantaneous decay items under the conditions of prepayment strategies and non-linear price-sensitive client demand. With the ambition of proper stock management, this study, for the first time, investigates the consequences of a prepayment scheme on the practitioner's ordering decision for non-instantaneous decay items while the customer demand changes non-linearly with respect to selling price.

    Realizing the consequences of selling price, several researchers in the stock management area have developed inventory models that take into account price-dependent demand rates. In general, consumer demand shows an upward trend when the unit selling price falls. Nagaraju et al. [20] developed a two-tier inventory system for price-dependent demand that uses both centralized and decentralized tactics to optimize lot size and inventory selections. Ghoreshi et al. [9] explored the impact of client demand for linear pricing on a practitioner's profit within a trade credit scheme for deteriorated items. Banerjee and Agrawal [2] investigated price-dependent demand and proposed joint discounting and ordering decisions for degraded products with a short shelf life. For deteriorating items, Chen et al. [3] developed an optimal pricing and replenishment method for when client demand varies due to the resultant impacts of stock amount, price and product storage time. Adopting both perfect and imperfect products from a production process, Ruidas et al. [28] formulated a production model focusing on price- and stock-related customer demand. Considering the impact of carbon emissions on the environment, Ruidas et al. [29] developed another production model for client demand that is related to linear pricing uses the selling price as a decision variable. Recently, Ruidas et al. [30,31] studied two other production models, adopting a price-sensitive customer demand. The non-linear price-related demand structures reflect the behavior of the clients more appropriately than the linear form. Pando et al. [21] considered the consequences of price in an exponential structure on client demand; they found the best pricing strategy for a retailer to maximize returns on stock management. By considering a non-linear price-sensitive demand structure instead of the linear form, San-José et al. [33] described the optimal pricing policy for a practitioner to minimize profit per unit time. To provide optimal pricing and inventory decisions as related to perishable items for practitioners, this study also investigated client demand structure as related to non-linear pricing.

    During stock-in, every perishable product/commodity deteriorates or loses usefulness due to many factors, such as temperature, dryness, perishability, spoilage, vaporization and time. As a result, the degradation effects play a significant role in inventory control system analysis. Ghare and Schreder [8], to our knowledge, were the first to use this approach to design an inventory modeling problem. Covert and Philip [4] improved this method by including two-parameter Weibull distributed degradation. Subsequently, several scholars conducted several studies based on static or variable degradation. Khan et al. [15] obtained closed-form optimal inventory policies for decay items, with the deterioration rate adopted as constant. Replacing the constant degradation rate by a time-varying degradation rate, Khan et al. [13] explored the best pricing and inventory decisions for a retailer, assuming expiration dates for perishable items. However, not all items (for example, vegetables, fish and fruits) perish from the moment of storage, but after some time from the moment of storage; they are are known as non-immediate/non-instantaneous perishable/decaying items in stock analysis. Ghoreshi et al. [9] explored the impact of client demand for a non-instantaneous decay rate on ordering decisions within a trade credit scheme. Yadav et al. [38] developed an inventory model for deteriorating electronic components within warehousing by using a genetic algorithm. Khan et al. [14] further investigated the effects of delayed decay initiation in a two-warehouse storage environment based on different decay initiation moments in different warehouses. Recently, Khan et al. [16] described a stock management scheme for non-instantaneous decay items under the conditions of an advance payment business scheme. However, these aforementioned authors assumed a constant or linear price-sensitive client demand structure while formulating the model. This study explores the non-instantaneous decay feature for perishable items under the conditions of non-linear price-sensitive demand for the first time.

    Due to volatility in the business markets, sometimes practitioners cannot fulfill the clients' desires on time. As a result, practitioners encounter a stock-out situation. During this stock-out period, a portion of customers either waits for the arrival of new products or moves elsewhere to fulfill their desires. When all clients wait for the new products to arrive, the situation is called "full backordering", which is not compatible for many perishable items. In addition, when only a portion of the customers wait to have their desires fulfilled, the situation is called "partial backordering", which is more compatible with many practical business situations. For degrading products with time-varying demand and holding costs, Dutta and Kumar [6] described a partial backordering environment in an inventory procedure. Jaggi et al. [10] investigated the impact of inflation and the time value of money on an inventory management system by considering perishable goods and partially backlog shortages. Under the conditions of price-sensitive demand and inventory age-associated holding costs, Rastogi and Singh [25] established an inventory model for fluctuating degrading and partial backlog pharmaceutical products. Rahman et al. [23] established the concept of a preservation facility and created an inventory model to account for price-sensitive demand and partial backordering. Udayakumar et al. [36] investigated a partial backlog shortage scenario in a practitioner's stock management scheme, focusing on deteriorating products under constant demand.

    For security of payment, advance payment is an essential method in inventory management. It is a business plan wherein a merchant orders a product by paying a part of the entire purchase price in equal installments before receiving the product, and the remaining amount is paid after receiving the product. An advance payment strategy ensures payment, as well as timely delivery of products. However, the practitioner incurs an opportunity cost against the prepayment amount because they cannot earn any revenue before receiving the goods. Zhang [40] presented this notion in his model with a set per-payment cost. By combining advance payment with delayed payment, Zia and Taleizadeh [42] created an inventory model to achieve the best ordering strategy. Zhang et al. [41] investigated a two-stage supply chain system with an advance payment method. Li et al. [17] proposed a hybrid advance, cash and delay payment strategy, looking at cash flow analysis, for a degrading inventory system. Taleizadeh [35] further investigated the impact of advance payment conditions on the inventory systems. In addition, Khan et al. [11] proposed a two-warehouse inventory model under the conditions of an advance payment scheme, whereas Shaikh et al. [34] studied an advance payment inventory model with interval-valued inventory costs. Furthermore, Manna et al. [18] explored the prepayment method in an interval environment to handle uncertainty. Khan et al. [15] proposed a flexible prepayment scheme in accordance with the purchased quantity. Alshanbari et al. [1] investigated an advance payment system for degrading goods under the conditions of a capital cost reduction facility utilizing multiple installments. Duary et al. [5] proposed an inventory model with advance payment for a discount facility with a deteriorating item. To improve the financing efficiency through the supply chain, Pfohl and Gomm [22], Gelsomino et al. [7] and Xu et al. [37] conducted separate systematic literature review studies on supply chain financing. Recently, Marchi et al. [19] studied reverse factoring as a technique to improve cash flow in a supply chain. In order to clarify the research gaps of a prepayment mechanism under stock management, Table 1 is presented.

    Table 1.  Comparative assessment of prepayment regulation under non-linear price dependent demand.
    Authors Demand Payment method for acquisition cost Deterioration Shortage Backordering Objective function
    Constant/ others Price-dependent in advance cash Instantaneous Non-instantaneous No Fully Partially with a Cost Minimization Profit maximization
    Ghoreshi et al. [20] Linear form Variable rate
    Zia and Taleizadeh [42] Constant Fixed rate
    Banerjee and Aggarwal [2] Non-linear form
    Li et al. [17] Exponential form
    Taleizadeh [35] Constant Fixed rate
    San-José et al. [32] Logit form
    Khan et al. [11] Constant Fixed rate
    Shaikh et al. [34] Constant Fixed rate
    Khan et al. [14] Constant Fixed rate
    Manna et al. [18] Linear form Variable rate
    Khan et al. [15] Constant Fixed rate
    Alshanbari et al. [1] Linear form Variable rate
    San-Jose et al. [32] Non-linear form
    Rahman et al. [24] Stock sensitive Variable rate
    Khan et al. [16] Time sensetive Variable rate
    Duary et al. [5] Time sensetive Variable rate
    Rezagholifam
    et al. [26]
    Inventory level sensitive Linear form
    This study Non-linear form Variable rate

     | Show Table
    DownLoad: CSV

    Table 1 reveals that only two studies [14,16] investigated the effects of a prepayment business strategy on ordering decisions for non-instantaneous decay items given a constant or time-sensitive client demand. However, they did not take into account the impacts of selling price on the clients' behavior. To the authors' best knowledge, there is no research work on inventory control to help practitioners by providing optimal ordering policies for non-instantaneous decay items under the conditions of prepayment strategies and price-sensitive client demand. The fundamental contribution of this study is to fill this research gap by designing an inventory model for a non-instantaneous decay item with non-linear price-dependent demand given an advance payment business policy. For the first time, this work investigates the consequences of a prepayment scheme on the practitioner's ordering decision for non-instantaneous decay items given that the customer demand changes non-linearly with respect to selling price. Moreover, the backlog rate varies based on how long consumers have been waiting, whereas a multiple installment facility is adopted to reduce the capital cost in the prepayment amount. To maximize the average profit, the existence of the optimal replenishing and pricing policies for the practitioner are explored theoretically. In addition, MATLAB software has been used to provide the graphic illustration of the average profit's concavity.

    Finally, the results of varying various parameters in the proposed model are shown, leading to a useful conclusion and guidelines for the practitioner.

    The remainder of the study is delineated as follows. The problem definition, along with all of the necessary hypotheses and notations, is delivered in Section 3. Section 4 exposes the problem formulation of the described inventory procedure. All analytical outcomes are derived for the optimal polices in Section 5, while Section 6 delivers numerical illustrations to verify the described inventory procedure. Following the sensitivity analysis, several managerial insights are derived in Section 7. Finally, Section 8 points out the conclusions of the study, along with several future research lines.

    Consider a business situation in which a practitioner who stocks a non-instantaneous decay item from a supplier to satisfy a portion of the future client desires given that the clients' behavior varies non-linearly with respect to selling price. Due to the perishable nature of the item, the supplier does not want to face cancellation or postponement of any order from the practitioner, as any order cancellation or postponement would cause a dramatic loss to their business. To prevent order cancellation or postponement, the supplier imposes a business scheme known as a prepayment mechanism. With this arrangement, the practitioner has to confirm the order by paying a percentage of the acquisition cost in advance, and the remaining percentage is to be paid at the time of receipt of the order. In addition, the practitioner incurs an opportunity cost against the prepayment amount because they cannot earn any revenue before receiving the goods. The main purpose of this study is to provide a complete decision framework for the practitioner by establishing optimal pricing and ordering decisions for a non-instantaneous decay item under a prepayment business contract to maximize the average profit.

    The following suppositions and notations are used to build the inventory model.

    (ⅰ) Nonlinear price-dependent demand is taken into account while developing the model. More precisely, the mathematical form of the demand structure is taken as D(p)=apb, where a>0 and b>1 [36]. The parameter a reveals the highest possible customer demand without any effect of price, and the parameter b measures the reduced customer demand after taking into account the consequences of price.

    (ⅱ) The deterioration starts after the td unit of time from the moment of storage, while the rate of deterioration θ(0<θ<<1) is constant during the deterioration occurring period. The deterioration starting moment, td, can be assessed from past experiences and by using the maximum likelihood technique. In this study, for simplicity, the value of td is assumed as a known constant such that 0tdt1, where t1 denotes the time point when stock becomes zero level [36].

    (ⅲ) Any deteriorated product is neither replaceable nor reparable [20].

    (ⅳ) The horizon of inventory planning is unlimited and the lead time is M unit of time, where the value of M is known and offered by the supplier [16].

    (ⅴ) The practitioner has to confirm the order by paying a percentage K of the acquisition cost in advance with the help of n equal installments, and the remaining percentage is to be paid at the time of receipt of the order [16].

    (ⅵ) A stock-out situation is encountered by the practitioner and a variable backlogging rate (in exponential form) is taken into consideration based on the customer waiting times. Mathematically, this rate is represented as eδx, where x denotes the waiting time of the consumers and δ>0 is the backlogging parameter [39].

    Notation Unit Description
    A $/order Replenishment cost
    a Constant Demand parameter (a>0)
    b Constant Demand parameter (b>0)
    Cp $/unit Purchase cost
    h $/unit/unit of time Holding cost
    Cs $/unit/unit of time Shortage cost
    Cd $/unit/unit of time Deterioration cost
    C1 $/unit/unit of time Lost sale cost
    I(t) Units Level of inventory any time t
    θ Constant Deterioration rate (0<θ<<1)
    δ Constant Backlogging parameter (δ>0)
    M Unit of time Lead time
    n Constant Installment number to prepay
    K Constant Fraction of the acquisition price to prepay (K[0,1])
    td Unit of time Time at which deterioration starts
    R Units Backlogged units
    S Units Maximum inventory level
    ˆX $/cycle Total cyclic cost
    X $/cycle Total cyclic profit
    TP(p,t1,t2) $/unit of time Total profit per unit of time
    Decision variables
    p $/unit Pre-unit selling price
    t1 Unit of time Time point when stock becomes zero level
    t2 Unit of time Shortage period

    Initially, an order is placed by the retailer with (S+R) units of a sole item by reimbursing a certain part K of the purchase cost with n equal several part payment facility within the receiving time M, and the products are received after the payment of all of the remaining amount of the purchase price at time t=0. Then, R units are utilized instantly to satisfy the backlogged demand and the level of inventory becomes S (see Figure 1).

    Figure 1.  Graphic representation of inventory and paying mechanism.

    The governing differential equations are as follows:

    dI(t)dt=D(p),0<ttd, (1)
    dI(t)dt+αI(t)=D(p),td<tt1, (2)
    dI(t)dt=D(p)eδ(t1+t2t),t1<tt1+t2, (3)

    with the initial conditions I(t)=S at t=0, I(t)=0 at t=t1 and I(t)=Ratt=t1+t2.

    Using the initial condition of inventory level at t=0, from Eq (1), one has

    I(t)=D(p)t+S. (4)

    The inventory level for t(td,t1], given Eq (2) with I(t)=0 at t=t1, is

    I(t)=D(p)α{eα(t1t)1}. (5)

    The shortage level for t(t1,t1+t2], given Eq (3) with I(t)=R at t=t1+t2, is

    I(t)=D(p)δ{1eδ(t1+t2t)}R. (6)

    The continuity property of the inventory level at t=td, as according to Eqs (4) and (5), provides

    S=D(p)td+D(p)α{eα(t1td)1}. (7)

    Again, the continuity property of the inventory level at t=t1, given Eqs (5) and (6), delivers

    R=D(p)δ(1eδt2). (8)

    For a single business period, the practitioner bears the following costs:

    Ordering cost =A:

    Carrying cost of inventory

    =htd0I(t)dt+ht1tdI(t)dt=h{D(p)t2d2+tdD(p)α{eα(t1td)1}}+hD(p)α2{eα(t1td)α(t1td)1}.

    Shortage cost

    =Cbt1+t2t1I(t)dt=CbD(p)δ2(1eδt2δt2eδt2).

    Purchase cost

    =CP(S+R)=CPD(p)[td+1α{eα(t1td)1}+1δ(1eδt2)].

    Capital cost

    =IC[(1+2+...+n)MnKCP(S+R)n]
    =ICMKCPD(p)(n+1)2n[td+1α{eα(t1td)1}+1δ(1eδt2)].

    Lost sale cost

    =C1t1+t2t1(1eδ(t1+t2t))Ddt=C1D(p){t21δ(1eδt2)}.

    As a result, the total cyclic cost for the practitioner is

    ˆX=[ordering cost+carrying cost+shortage cost+purchasecost+capital cost+lostsalecost]
    =A+h{D(p)t2d2+tdD(p)α{eα(t1td)1}}+hD(p)α2{eα(t1td)α(t1td)1}+CbD(p)δ2(1eδt2δt2eδt2)+{1+(n+1)2nICMK}CPD(p)[td+1α{eα(t1td)1}+1δ(1eδt2)]+C1D(p){t21δ(1eδt2)}.

    The sales revenue for the practitioner within a single business period is

    p[t10D(p)dt+t1+t2t1D(p)eδ(t1+t2t)dt]=pD(p){t1+1δ(1eδt2)}.

    Consequently, the profit for the practitioner within a single business period is

    X=pD(p){t1+1δ(1eδt2)}Ah{D(p)t2d2+tdD(p)α{eα(t1td)1}}
    hD(p)α2{eα(t1td)α(t1td)1}CbD(p)δ2(1eδt2δt2eδt2)
    {1+(n+1)2nICMK}CPD(p)[td+1α{eα(t1td)1}+1δ(1eδt2)]C1D(p){t21δ(1eδt2)}.

    Finally, the optimization problem for the entire system takes the form

    MaximizeTP(p,t1,t2)=Xt1+t2, (9)

    which is subject to 0tdt1,t1>0 and t2>0.

    To compute the optimum solution of Problem (9), the first-order partial derivatives of TP(p,t1,t2) with respect to p,t1 and t2 should be calculated as

    TP(p,t1,t2)t1=X(t1+t2)2+D(p)(t1+t2)[phtdeα(t1td)hα{eα(t1td)1}{1+(n+1)2nICMK}CPeα(t1td)], (10)
    TP(p,t1,t2)t2=X(t1+t2)2+D(p)(t1+t2)[peδt2Cbt2eδt2C1(1eδt2){1+(n+1)2nICMK}CPeδt2] (11)

    and

    TP(p,t1,t2)p=abp(b+1)(t1+t2)[p(1b1){t1+1δ(1eδt2)}+h{t2d2+tdα{eα(t1td)1}}+hα2{eα(t1td)α(t1td)1}+Cbδ2(1eδt2δt2eδt2)+{1+(n+1)2nICMK}CP[td+1α{eα(t1td)1}+1δ(1eδt2)]+C1{t21δ(1eδt2)}]. (12)

    The optimality is explored in two ways. First, the optimum replenishing policy is explored for a given p>0. For this purpose, the necessary conditions from Eqs (10) and (11) are

    D(p)(t1+t2)[phtdeα(t1td)hα{eα(t1td)1}{1+(n+1)2nICMK}CPeα(t1td)]=X (13)

    and

    D(p)(t1+t2)[peδt2Cbt2eδt2C1(1eδt2){1+(n+1)2nICMK}CPeδt2]=X. (14)

    Since the right-hand sides of both Eqs (13) and (14) are the same, one can find

    Cbt2eδt2+(p+C1)(1eδt2)+{1+(n+1)2nICMK}CPeδt2=htdeα(t1td)+hα{eα(t1td)1}+{1+(n+1)2nICMK}CPeα(t1td). (15)

    From Eq (14), one can write

    D(p)(t1+t2)[peδt2Cbt2eδt2C1(1eδt2){1+(n+1)2nICMK}CPeδt2]pD(p){t1+1δ(1eδt2)}Ah{D(p)t2d2+tdD(p)α{eα(t1td)1}}hD(p)α2{eα(t1td)α(t1td)1}CbD(p)δ2(1eδt2δt2eδt2){1+(n+1)2nICMK}CPD(p)[td+1α{eα(t1td)1}+1δ(1eδt2)]C1D(p){t21δ(1eδt2)}. (16)

    The first-order derivative of the left-hand side of Eq (15) with respect to t2 is

    [Cb(1δt2)+(p+C1)δ{1+(n+1)2nICMK}CPδ]eδt2.

    Consequently, the left-hand side of Eq (15) is continuous; it increases for all t2[0,˜t2] and decreases for all t2, where

    ˜t2=Cb+(p+C1)δ{1+(n+1)2nICMK}CPδCbδ.

    Consequently, the left-hand side of Eq (15) has a maximum at ˜t2, and the maximum is

    (p+C1)δ{CP+(n+1)2nICMKCPpC1}+(Cbδ+p+C1)eCb+(p+C1)δ{1+(n+1)2nICMK}CPδCb.

    Regarding the right-hand side, the right-hand side of Eq (15) is increasing for all t1[td,) and tends to infinity when t1. As a result, there exists a unique ˜t1 such that

    htdeα(˜t1td)+hα{eα(˜t1td)1}+{1+(n+1)2nICMK}CPeα(˜t1td)(p+C1)δ{CP+(n+1)2nICMKCPpC1}
    +(Cbδ+p+C1)eCb+(p+C1)δ{1+(n+1)2nICMK}CPδCb.

    Moreover, there exists a unique t1[td,) for any t'2(0,˜t2) such that Eq (15) is maintained. Similarly, for any given t'2(˜t2,), there exists a unique t1[td,) such that Eq (15) is maintained. As a result, t1 can be computed in terms of t2 uniquely, and the expression of t1 in terms of t2, given Eq (15), is

    t1=td+1αln[Cbt2eδt2+(p+C1)(1eδt2)+{1+(n+1)2nICMK}CPeδt2+hαhtd+hα+{1+(n+1)2nICMK}CP]. (17)

    Let us define the auxiliary function from Eq (16) as

    Ψ(t2)=D(p)(t1+t2)[peδt2Cbt2eδt2C1(1eδt2){1+(n+1)2nICMK}CPeδt2]pD(p){t1+1δ(1eδt2)}+A+h{D(p)t2d2+tdD(p)α{eα(t1td)1}}+hD(p)α2{eα(t1td)α(t1td)1}+C1D(p){t21δ(1eδt2)}+CbD(p)δ2(1eδt2δt2eδt2)+{1+(n+1)2nICMK}CPD(p)[td+1α{eα(t1td)1}+1δ(1eδt2)]. (18)

    The first-order derivative of Ψ(t2) with respect to t2 is

    dΨ(t2)dt2=D(p)(dt1dt2+1)[peδt2Cbt2eδt2C1(1eδt2){1+(n+1)2nICMK}CPeδt2]+D(p)(t1+t2)[δpeδt2Cb(eδt2δt2eδt2)C1δeδt2+δ{1+(n+1)2nICMK}CPeδt2]pD(p){dt1dt2+eδt2}+htdD(p)eα(t1td)dt1dt2+hD(p)α{eα(t1td)1}dt1dt2+CbD(p)t2eδt2+{1+(n+1)2nICMK}CPD(p)[eα(t1td)dt1dt2+eδt2]+C1D(p){1eδt2},

    i.e.,

    dΨ(t2)dt2=D(p)(t1+t2)[Cb(1δt2)+δ(p+C1)δ{1+(n+1)2nICMK}CP]eδt2. (19)

    Now, the existence of the optimal replenishing policy can be observed from the following theorem:

    Theorem. For any p>0,

    (a) if Ψ(˜t2)<0, then the optimal point (t1,t2) is unique and maximizes TP(t1,t2|p) given t2(0,˜t2).

    (b) if Ψ(˜t2)0, the optimal value of t2 is t2.

    Proof. See Appendix.

    For any p>0, Theorem (a) shows that, if Ψ(˜t2)<0, then TP(t1,t2|p) holds its global maximum value at the single point (t1,t2), where t2(0,˜t2). Now, at the point (t1,t2), from Eq (14), one can find

    TP(t1,t2|p)=D(p)[peδt2Cbt2eδt2C1(1eδt2){1+(n+1)2nICMK}CPeδt2].

    Again, Theorem (b) shows that, if Ψ(˜t2)0, then t2. Then, from Eq (17), one can compute

    t1=td+1αln[p+C1+hαhtd+hα+{1+(n+1)2nICMK}CP]

    and limt2TP(t1,t2|p)=C1D(p). This yields that the given p>0 provides a negative profit; hence, this p>0 is unstable and it should be improved.

    Now, the unique optimal selling price per unit is computed after analyzing its existence. For any t1>td and t2>0, taking the first-order derivative with respect to p and performing some simplifications, the necessary condition to compute p is

    dTP(p|t1,t2)dp=[D(p)+D(p){pCPCP(n+1)2nICMK}]{1+1δt2eδt2δ(t1+t2)}D(p)(t1+t2)[Cbδ2(1eδt2δt2eδt2)C1δ(1δt2eδt2)h{t2d2+tdα{eα(t1td)1}}+hα2{eα(t1td)α(t1td)1}+{1+(n+1)2nICMK}CP[td+1α{eα(t1td)1}]{1+(n+1)2nICMK}CPt1]=0. (20)

    Since

    {1+(n+1)2nICMK}CP[td+1α{eα(t1td)1}]{1+(n+1)2nICMK}CPt1
    ={1+(n+1)2nICMK}CP[td+1α{1+α(t1td)+12α2(t1td)2+...1}t1]>0,
    1eδt2δt2eδt2>0,1δt2eδt2<0,eα(t1td)1

    and

    eα(t1td)α(t1td)1>0,

    the coefficient of D'(p)(t1+t2) is always positive. In addition, D'(p)=abp(b+1) is negative always. Consequently, the necessary condition (i.e., Eq (20)) is solvable for p if

    D(p)+D(p){pCPCP(n+1)2nICMK}<0,

    that is,

    p(1b)+bCP{1+(n+1)2nICMK}<0

    as

    1+1δt2eδt2δ(t1+t2)=t2(t1+t2)+1eδt2δ(t1+t2)>0.

    Furthermore, if the per unit gross profit including the capital cost D(p){pCPCP(n+1)2nICMK} is strictly concave in p, then its derivative D(p)+D(p){pCPCP(n+1)2nICMK} is strictly decreasing; hence,

    2D(p)+D(p){pCPCP(n+1)2nICMK}<0.

    As a result, the second-order derivative with respect to p is

    d2TP(p|t1,t2)dp2={2D(p)+D(p)[pCP{1+(n+1)2nICMK}]}{1+1δt2eδt2δ(t1+t2)}D(p)(t1+t2)[Cbδ2(1eδt2δt2eδt2)C1δ(1δt2eδt2)h{t2d2+tdα{eα(t1td)1}}+hα2{eα(t1td)α(t1td)1}+{1+(n+1)2nICMK}CP[td+1α{eα(t1td)1}]{1+(n+1)2nICMK}CPt1]<0. (21)

    Therefore, a single selling price p can be found which maximizes the total profit of the practitioner. In addition, the selling price p which satisfies the equation

    p(1b)+bCP{1+(n+1)2nICMK}=0

    is the lower bound of the optimal selling price (say, pl). As a result,

    pl=bCP{1+(n+1)2nICMK}b1.

    Summarizing the theoretical outcomes, the following algorithm is established to achieve the optimal solution for the retailer.

    Algorithm

    Step 1: Input the values of known parameters.
    Step 2: Start with setting j=1 and p=p(j)=bCP{1+(n+1)2nICMK}b1, which is the solution of p(1b)+bCP{1+(n+1)2nICMK}=0.
    Step 3: Put p(j) and ˜t2=Cb+(p(j)+C1)δ{1+(n+1)2nICMK}CPδCbδ into Eq (15) to compute the corresponding value of t1, that is, ˜t1. Utilizing these values, given Eq (18), find the value of Ψ(˜t2).
    Step 4: If Ψ(˜t2)<0, then go to Step 5; otherwise, set p(j+1)=p(j)+κ, where κ is a sufficiently small positive number, and go to Step 3 after updating j by j+1.
    Step 5: From Eq (18), find the value of t2 by solving Ψ(t2)=0. Determine the corresponding t1 from Eq (15) and then go to Step 5.
    Step 6: Compute the value of p from Eq (20). Set p(j+1)=p.
    Step 7: If |p(j+1)p(j+1)| is sufficiently small, go to Step 8; otherwise, set j=j+1 and go to Step 3.
    Step 8: Report the optimal solution: p=p(j+1), t1 and t2 with the maximum profit TP(p,t1,t2).
    Step 9: End

    In the next section, we illustrate the described model by solving three numerical examples.

    Utilizing the solution algorithm in the previous section, three numerical examples were solved while LINGO 18 software was used to solve all the necessary equations during the implementation of the algorithm. In addition, MATLAB software was used to illustrate the concavity of the objective function by using the same numerical case to evaluate the applicability of the suggested model.

    Example 1. Let A=$200 /order, a=3500, b=1.5, Cp=$30 /unit, C1=$10 /unit, Ch=$1 /unit, α=0.05, δ=0.4, M=0.25 months, n=20, Ic=0.01/$/month, K=0.4, td=0.2 month and Cs=$15 /unit. Therefore, adopting the algorithm, the optimal solutions for this example are p=115.8991, t1=6.5999 months, t2=0.3964 months, S=21.7184 units, R=1.0284 units, Q=22.7468 units and TP(p,t1,t2)=$187.2284.

    The concavity of the objective function for Example 1 is shown pairwise with respect to the independent decision variables (p, t1 and t2) in Figures 24.

    Figure 2.  Concavity of the average profit function with respect to p and t1.
    Figure 3.  Concavity of the average profit function with respect to p and t2.
    Figure 4.  Concavity of the average profit function with respect to t1 and t2.

    Example 2. Let A=$250/order, a=2500, b=1.2, Cp=$35/unit, C1=$10/unit, Ch=$1/unit, α=0.05, δ=0.4, M=0.5 months, n=20, Ic=0.01/$/month, K=0.5, td=0.5 month and Cs=$15/unit. Hence, adopting the algorithm, the optimal solutions for this example are p=266.3658, t1=6.8282 months, t2=0.1985 months, S=24.4005 units, R=0.5862 units, Q=24.9867 units and TP(p,t1,t2)=$645.4862.

    Example 3. Let A=$250/order, a=2200, b=1.4, Cp=$55/unit, C1=$15/unit, Ch=$1.5/unit, α=0.07, δ=0.5, M=0.8 months, n=25, Ic=0.05/$/month, K=0.45, td=0.4 months and Cs=$17/unit. Then, adopting the algorithm, the optimal solutions for this example are p=283.5804, t1=8.5979 months, t2=0.5030 months, S=9.2970 units, R=0.3604 units, Q=9.6574 units and TP(p,t1,t2)=$135.6230.

    This section deals with the variation patterns of optimal policies p,t1,t2,S,R, along with the profit per unit of time due to the changes of the original value of the key parameters in Example 1. This process is conducted by altering only a single parameter at a time from -20% to +20% and keeping the values of the rest of the parameters as the same. The corresponding computational results are illustrated graphically in Figures 58.

    Figure 5.  Percent change with respect to a.
    Figure 6.  Percent change with respect to b.
    Figure 7.  Percent change with respect to Cp.
    Figure 8.  Percent change with respect to A.

    From the sensitivity analysis in Figure 5, it is concluded that the profit per unit time (TP) is equally sensitive, whereas the positive stock period (t1) and stock-out period (t2) are relatively responsive and the initial stock (S) and highest shortage (R) are less responsive for changing the value of a. The optimal selling price for the retailer falls gradually when the demand parameter a increases. Though the optimal selling price decreases, the maximum profit per unit of time increases. This observation suggests that the decision manager should implement some effective marketing strategies (for example, attractive advertising campaigns) to increase potential market demand and then reduce the selling price slightly.

    According to Figure 6, the profit per unit time (TP) and stock-out period t2 are highly responsive, whereas the positive stock period (t1), initial stock (S) and highest shortage level (R) are equally responsive in terms of altering the value of the price-sensitive parameter b. The variation of the price-sensitive parameter (b) reveals the dramatic consequences on the stock-out period. When the parameter b increases, the client demand decreases; hence, the manager stocks a relatively small amount. As a result, the stock-out period increases significantly. On the other hand, the profit of the retailer falls as client demand falls considerably when the parameter b increases.

    Figure 7 exposes that profit per unit of time (TP) is equally sensitive, while both the initial stock (S) and maximum shortage level (R) are highly sensitive with respect to the variation of the per unit acquisition price (Cp). On the another hand, the positive stock period (t1) and stock-out period (t2) are relatively sensitive when Cp increases or decreases. The optimal selling price for the retailer increases when the per unit acquisition price increases. However, the profit per unit time decreases. In this case, increasing the per unit selling price does not ensure an increase of the profit. This outcome suggests that the decision manager should negotiate the per unit acquisition price with the supplier. Moreover, when the per unit acquisition price (Cp) rises, the selling price (p) also increases; hence, the client demand D(p) falls considerably. As a result, both the initial stock (S) and maximum shortage level (R) falls when the per unit acquisition price increases.

    Figure 8 indicates that the profit per unit of time (TP) has only negative consequences due to the variation of the ordering cost per order (A), while all of the rest (independent and dependent) of the decision variables have positive impacts with respect to the changes of A. When the ordering cost increases, the decision manager purchases more quantities so that the ordering cost against each unit falls and, therefore, the stock-in period (t1) and entire cycle length increase. To face a higher ordering cost per order, the decision manager sets up a higher selling price for each unit. In this case, though the unit selling price increases, the profit per unit time decreases.

    Taking into account the behavior with respect to the variations of the inventory parameters, the following management insights have been drawn for the inventory manager.

    (a) As both the demand parameters (a and b) have the greatest consequences for the manager's profit per unit of time, the decision manager should implement some effective marketing strategies (for example, attractive advertising campaigns) to increase potential market demand and then reduce the selling price slightly.

    (b) Though the optimal selling price for the retailer increases, the profit per unit time decreases when the per unit acquisition price increases. In this case, increasing the per unit selling price does not ensure an increase of the profit. This outcome suggests that the decision manager should negotiate the per unit acquisition price with the supplier.

    (c) When the ordering cost for placing every order increases, the decision manager should purchase more quantities so that the ordering cost against each unit decreases.

    (d) As

    dTPdM={(n+1)2nICK}CPD(p)[td+1α{eα(t1td)1}+1δ(1eδt2)]<0,

    the inventory manager's profit is strictly decreasing with respect to the allowed lead time from the supplier. When the lead time (M) increases, the inventory manger's capital cost or opportunity cost from the prepayment amount increases; hence, the profit per unit of time falls. This observation suggests that the decision manager should negotiate the lead time with the supplier.

    (e) Since

    dTPdK={(n+1)2nICM}CPD(p)[td+1α{eα(t1td)1}+1δ(1eδt2)]<0,

    the profit per unit of time is strictly decreasing when the fraction of the acquisition price to prepay increases. When the fraction of the acquisition price to prepay (K) increases, the inventory manger's capital cost or opportunity cost from the prepayment amount increases; hence, the profit per unit of time decreases. Therefore, the decision manager should negotiate with the supplier to reduce the fraction of the acquisition price for prepayment.

    (f) As

    dTPdn=ICKCPD(p)2n2[td+1α{eα(t1td)1}+1δ(1eδt2)]>0,

    the profit per unit of time shows an upward trend when the installment frequency increases. A higher installment number for prepayment helps the manager to reduce the capital cost or opportunity cost against the prepayment amount. Consequently, the inventory manger should motivate the supplier to allow a higher installment frequency to implement the prepayment regulation.

    In this study, a stock control system was investigated for degrading products with non-linear price-dependent client demand and advanced payment regulations with a multiple installment plan given exponential partial backordering. To achieve the optimal joint pricing and replenishment schedule for the retailer, some important characteristics were explored. The nature of the decision variables in the formulated inventory problem is complex; hence, the optimal values of these variables cannot be obtained via a simple algebraic approach. As a result, the optimal pricing and replenishment scheduling are achieved by employing a numerical technique in an algorithmic manner. Three numerical examples were used to ensure the efficacy of the algorithm; the concavity of the average profit function was then visually shown pairwise with respect to the decision variables by using MATLAB software. To obtain management insights into optimal policies, the impacts of key parameters were explored by altering the values of these parameters in an inventory problem. Precisely, the decision manager should implement some effective marketing strategies (for example, attractive advertising campaigns) to increase potential market demand and then reduce the selling price slightly in order to increase the profit per unit of time. When the ordering cost for placing every order increases, the decision manager should purchase more quantities so that the ordering cost against each unit decreases. The inventory manager should motivate the supplier to allow a higher installment frequency to implement the prepayment regulation so that the capital cost or opportunity cost against the prepayment amount is reduced and, hence, the profit per unit of time is increased. Moreover, to increase the profit, the decision manager should negotiate with the supplier to reduce both the lead time and fraction of the acquisition price for accomplishing the prepayment.

    This work considers the installment frequency as an input parameter to complete the prepayment without incurring any cost against each installment, which is a limitation of this study. Therefore, an interesting extension of this work would be to investigate the optimal installment frequency for the retailer.

    No effect of product storage time on client demand was considered while developing the model. However, the behavior of clients depends not only on the unit selling price, but also on the storage period of the deteriorating products. Hence, another immediate future line of this study is to investigate the combined effects of price and time on the demand structure, either in additive or multiplicative form. It might also be an interesting expansion of a model without shortage to loosen the zero-ending situation by using a non-ending inventory model.

    We would like to thank King Saud University, Riyadh, Saudi Arabia for the Research Supporting Project grant (RSP-2021/323).

    We have no conflicts with competing interests to disclose.

    Proof of the theorem.

    (a) Since Cb(1δt2)+δ(p+C1)δ{1+(n+1)2nICMK}CP is positive for all t2 and negative for all t2(˜t2,), Eq (19) reveals that Ψ(t2) is strictly decreasing for all t2 and increasing for all t2(˜t2,). As a result, Ψ(t2) has a minimum at t2=˜t2, and the minimum is Ψ(˜t2). First, the interval [0,˜t2] is considered. When t2=0, Eq (17) shows that t1<td, which contradicts the restriction t1td. As a result, when t2=0, t1=td. Thus, one has

    Ψ(0)=A+hD(p)t2d2+{1+(n+1)2nICMK}CPD(p)td>0.

    Because Ψ(t2) is strictly decreasing for all t2, under the condition of Ψ(˜t2)<0, the intermediate value theorem ensures a unique solution t2(0,˜t2) such that Ψ(t2)=0. Now, exploiting

    Cb(1δt2)+(p+C1)δ{1+(n+1)2nICMK}CPδ>0

    for all t2<˜t2, one has

    2TP(t1,t2|p)t21|(t1,t2)=(t1,t2)=D(p)(t1+t2)[αhtdeα(t1td)+heα(t1td)α{1+(n+1)2nICMK}CPeα(t1td)]<0,
    2TP(t1,t2|p)t22|(t1,t2)=(t1,t2)=D(p)(t1+t2)[Cb(1δt2)+(p+C1)δ{1+(n+1)2nICMK}CPδ]eδt2<0.

    And,

    2TP(t1,t2|p)t1t2|(t1,t2)=(t1,t2)=2TP(t1,t2|p)t2t1|(t1,t2)=(t1,t2)=0.

    Therefore, the determinant of the Hessian matrix at (t1,t2) is

    {D(p)}2eδt2(t1+t2)2[αhtdeα(t1td)+heα(t1td)α{1+(n+1)2nICMK}CPeα(t1td)][Cb(1δt2)+(p+C1)δ{1+(n+1)2nICMK}CPδ]>0.

    Since the first and second principal minors of the Hessian matrix are negative and positive, respectively, (t1,t2) is the optimal solution.

    Again, the interval is considered now. Because Ψ(t2) is increasing for all t2(˜t2,), there are two possible cases corresponding to the condition Ψ(˜t2)<0, namely, (i) there is no t2(˜t2,) such that Ψ(t2)=0, and (ii) there is a single point ˆt2 in the interval (˜t2,) such that Ψ(ˆt2)=0. For the first possibility, the optimal solution for t2 does not exist. On the another hand, for the second possibility, the determinant of the Hessian matrix at (t1,t2) becomes negative as

    Cb(1δt2)+(p+C1)δ{1+(n+1)2nICMK}CPδ<0

    for all t2>˜t2. As a result, the point (t1,t2) is not the optimal solution.

    (b) Because Ψ(t2) has a global minimum at t2=˜t2, Ψ(t2) is always positive for all t2 given the condition Ψ(˜t2)>0. Now, exploiting Eqs (11) and (18), one has TP(t1,t2|p)t2=Ψ(t2)(t1+t2)2>0. This highlights that TP(t1,t2|p) has a higher value for higher values of t2. Therefore, TP(t1,t2|p) achieves the maximum value at the point t2. Again, when Ψ(˜t2)=0, TP(t1,t2|p)t2|t2=˜t2=0. Since TP(t1,t2|p) is increasing in both the intervals (0,˜t2) and (˜t2,) in this case, the point t2=˜t2 is an inflection point; hence, TP(t1,t2|p) achieves the maximum value at the point t2.

    This completes the proof.



    [1] K. A. Ahmad, R. Ezzati, K. M. Afshar, Solving systems of fractional two-dimensional nonlinear partial Volterra integral equations by using Haar wavelets, J. Appl. Anal., 27 (2021), 239–257. https://doi.org/10.1515/JAA-2021-2050 doi: 10.1515/JAA-2021-2050
    [2] A. Karimi, K. Maleknejad, R. Ezzati, Numerical solutions of system of two-dimensional Volterra integral equations via Legendre wavelets and convergence, Appl. Numer. Math., 156 (2020), 228–241. https://doi.org/10.1016/j.apnum.2020.05.003 doi: 10.1016/j.apnum.2020.05.003
    [3] H. Liu, J. Huang, W. Zhang, Y. Ma, Meshfree approach for solving multi-dimensional systems of Fredholm integral equations via barycentric Lagrange interpolation, Appl. Math. Comput., 346 (2019), 295–304. https://doi.org/10.1016/j.amc.2018.10.024 doi: 10.1016/j.amc.2018.10.024
    [4] J. Xie, M. Yi, Numerical research of nonlinear system of fractional Volterra-Fredholm integral-differential equations via Block-Pulse functions and error analysis, J. Comput. Appl. Math., 345 (2019), 159–167. https://doi.org/10.1016/j.cam.2018.06.008 doi: 10.1016/j.cam.2018.06.008
    [5] P. Gonzˊalez-Rodelas, M. Pasadas, A. Kouibia, B. Mustafa, Numerical solution of linear Volterra integral equation systems of second kind by radial basis functions, Mathematics, 10 (2022), 223. https://doi.org/10.3390/MATH10020223 doi: 10.3390/MATH10020223
    [6] A. R. Yaghoobnia, R. Ezzati, Using Bernstein multi-scaling polynomials to obtain numerical solution of Volterra integral equations system, Comput. Appl. Math., 39 (2020), 608–616. https://doi.org/10.1007/s40314-020-01198-4 doi: 10.1007/s40314-020-01198-4
    [7] A. Jafarian, S. Measoomy, S. Abbasbandy, Artificial neural networks based modeling for solving Volterra integral equations system, Appl. Soft Comput., 27 (2015), 391–398. https://doi.org/ 10.1016/j.asoc.2014.10.036 doi: 10.1016/j.asoc.2014.10.036
    [8] J. Cao, C. Xu, A high order schema for the numercial solution of the fractional ordinary differential equations, J. Comput. Phys., 238 (2013), 154–168. https://doi.org/10.1016/j.jcp.2012.12.013 doi: 10.1016/j.jcp.2012.12.013
    [9] R. Katani, S. Shahmorad, Block by block method for the systems of nonlinear Volterra integral equtions, Appl. Math. Model., 34 (2010), 400–406. https://doi.org/10.1016/j.apm.2009.04.013 doi: 10.1016/j.apm.2009.04.013
    [10] H. H. Sorkun, S. Yalçinbaş, Approximate solutions of linear Volterra integral equation systems with variable coefficients, Appl. Math. Model., 34 (2010), 3451–3464. https://doi.org/10.1016/j.apm.2010.02.034 doi: 10.1016/j.apm.2010.02.034
    [11] M. I. Berenguer, D. Gámez, A. I. G. Guillem, M. R. Galán, M. C. S. Pérez, Biorthogonal systems for solving Volterra integral equation systems of the second kind, J. Comput. Appl. Math., 235 (2010), 1875–1883. https://doi.org/10.1016/j.cam.2010.07.011 doi: 10.1016/j.cam.2010.07.011
    [12] K. Maleknejad, A. S. Shamloo, Numerical solution of singular Volterra integral equations system of convolution type by using operational matrices, Appl. Math. Comput., 195 (2007), 500–505. https://doi.org/10.1016/j.amc.2007.05.001 doi: 10.1016/j.amc.2007.05.001
    [13] A. Tahmasbi, O. S. Fard, Numerical solution of linear Volterra integral equations system of the second kind, Appl. Math. Comput., 201 (2008), 547–552. https://doi.org/10.1016/j.amc.2007.12.041 doi: 10.1016/j.amc.2007.12.041
    [14] M. Rabbani, K. Maleknejad, N. Aghazadeh, Numerical computational solution of the Volterra integral equations system of the second kind by using an expansion method, Appl. Math. Comput., 187 (2006), 1143–1146. https://doi.org/10.1016/j.amc.2006.09.012 doi: 10.1016/j.amc.2006.09.012
    [15] S. Yalçinbaş, K. Erdem, Approximate solutions of nonlinear Volterra equation systems, Internat. J. Modern Phys. B., 24 (2010), 6235–6258. https://doi.org/10.1142/S0217979210055524 doi: 10.1142/S0217979210055524
    [16] M. A. Zaky, I. G. Ameen, N. A. Elkot, E. H. Doha, A unified spectral collocation method for nonlinear systems of multi-dimensional integral equations with convergence analysis, Appl. Numer. Math., 161 (2021), 27–45. https://doi.org/10.1016/j.apnum.2020.10.028 doi: 10.1016/j.apnum.2020.10.028
    [17] E. Babolian, M. Mordad, A numerical method for solving systems of linear and nonlinear integral equations of the second kind by hat basis functions, Comput. Math. Appl., 62 (2011), 187–198. https://doi.org/10.1016/j.camwa.2011.04.066 doi: 10.1016/j.camwa.2011.04.066
    [18] F. Mirzaee, E. Hadadiyan, Solving system of linear Stratonovich Volterra integral equations via modification of hat functions, Appl. Math. Comput., 293 (2017), 254–264. https://doi.org/10.1016/j.amc.2016.08.016 doi: 10.1016/j.amc.2016.08.016
    [19] E. Babolian, J. Biazar, A. R. Vahidi, On the decomposition method for system of linear equations and system of linear Volterra integral equations, Appl. Math. Comput., 147 (2004), 19–27. https://doi.org/10.1016/S0096-3003(02)00644-6 doi: 10.1016/S0096-3003(02)00644-6
    [20] K. Maleknejad, M. Shahrezaee, Using Runge-Kutta method for numerical solution of the system of Volterra integral equation, Appl. Math. Comput., 149 (2004), 399–410. https://doi.org/10.1016/s0096-3003(03)00148-6 doi: 10.1016/s0096-3003(03)00148-6
    [21] W. Jiang, Z. Chen, Solving a system of linear Volterra integral equations using the new reproducing kernel method, Appl. Math. Comput., 219 (2013), 10225–10230. https://doi.org/10.1016/j.amc.2013.03.123 doi: 10.1016/j.amc.2013.03.123
    [22] J. Biazar, H. Ebrahimi, Chebyshev wavelets approach for nonlinear systems of Volterra integral equations, Comput. Math. Appl., 63 (2012), 608–616. https://doi.org/10.1016/j.camwa.2011.09.059 doi: 10.1016/j.camwa.2011.09.059
    [23] F. Mirzaee, S. Hoseini, A new collocation approach for solving systems of high-order linear Volterra integro-differential equations with variable coefficients, Appl. Math. Comput., 311 (2017), 272–282. https://doi.org/10.1016/j.amc.2017.05.031 doi: 10.1016/j.amc.2017.05.031
    [24] D. Conte, S. Shahmorad, Y. Talaei, New fractional Lanczos vector polynomials and their application to system of Abel-Volterra integral equations and fractional differential equations, J. Comput. Appl. Math., 366 (2020), 112409. https://doi.org/10.1016/j.cam.2019.112409 doi: 10.1016/j.cam.2019.112409
    [25] K. Sadri, K. Hosseini, D. Baleanu, S. Salahshour, A high-accuracy Vieta-Fibonacci collocation scheme to solve linear time-fractional telegraph equations, Wave. Rand. Complex Media, 2022, 2135789. https://doi.org/10.1080/17455030.2022.2135789
    [26] S. Khadijeh, H. Kamyar, B. Dumitru, S. Soheil, P. Choonkil, Designing a matrix collocation method for fractional delay intgro-differential equations with weakly singular kernels based on Vieta-Fibonacci polynomials, Fractal Fract., 6 (2022). https://doi.org/10.3390/fractalfract6010002
    [27] J. Dixon, S. Mxkee, Weakly singular discrete Gronwall inequalities, Z. Angew. Math. Mech., 66 (1978), 535–544. https://doi.org/10.1002/zamm.19860661107 doi: 10.1002/zamm.19860661107
  • This article has been cited by:

    1. Md. Al-Amin Khan, Leopoldo Eduardo Cárdenas-Barrón, Neale R. Smith, Rafael Ernesto Bourguet-Díaz, Imelda de Jesús Loera-Hernández, Rodrigo Ernesto Peimbert-García, Gerardo Treviño-Garza, Armando Céspedes-Mota, Effects of all-units discount on pricing and replenishment policies of an inventory model under power demand pattern, 2023, 10, 2330-2674, 10.1080/23302674.2022.2161712
    2. Md. Al-Amin Khan, Leopoldo Eduardo Cárdenas-Barrón, Gerardo Treviño-Garza, Armando Céspedes-Mota, A prepayment installment decision support framework in an inventory system with all-units discount against link-to-order prepayment under power demand pattern, 2023, 213, 09574174, 119247, 10.1016/j.eswa.2022.119247
    3. Jui-Jung Liao, H. M. Srivastava, Shy-Der Lin, Inventory models for non-instantaneous deteriorating items with expiration dates under the joined effect of preservation technology and linearly time-dependent holding cost when order-size linked to advance payment, 2024, 337, 0254-5330, 197, 10.1007/s10479-024-05909-6
    4. Valentín Pando, Luis A. San-José, Joaquín Sicilia, David Alcaide-López-de-Pablo, Pricing decision in a newsvendor model with partial backorders under normal probability distribution for the demand, 2024, 132, 0307904X, 57, 10.1016/j.apm.2024.04.037
    5. Valentín Pando, Luis A. San-José, Joaquín Sicilia, David Alcaide-López-de-Pablo, An inventory model with price- and stock-dependent demand and time- and stock quantity-dependent holding cost under profitability maximization, 2024, 164, 03050548, 106520, 10.1016/j.cor.2023.106520
    6. Md. Al-Amin Khan, Ali Akbar Shaikh, Aminur Rahman Khan, Adel Fahad Alrasheedi, Advertising and pricing strategies of an inventory model with product freshness-related demand and expiration date-related deterioration, 2023, 73, 11100168, 353, 10.1016/j.aej.2023.04.059
    7. Md. Rukonuzzaman, Md. Al-Amin Khan, Aminur Rahman Khan, Ali AlArjani, Md. Sharif Uddin, El-Awady Attia, Effects of a quantity-based discount frame in inventory planning under time-dependent demand: A case study of mango businesses in Bangladesh, 2023, 35, 10183647, 102840, 10.1016/j.jksus.2023.102840
    8. Md. Al-Amin Khan, Leopoldo Eduardo Cárdenas-Barrón, Gerardo Treviño-Garza, Armando Céspedes-Mota, Installment for advance payment and pricing decisions for an inventory system under power demand pattern and all-units discount, 2023, 265, 09255273, 108951, 10.1016/j.ijpe.2023.108951
    9. Fleming Akhtar, Md. Al-Amin Khan, Ali Akbar Shaikh, Adel Fahad Alrasheedi, Interval valued inventory model for deterioration, carbon emissions and selling price dependent demand considering buy now and pay later facility, 2024, 15, 20904479, 102563, 10.1016/j.asej.2023.102563
    10. P. Saranyaa, C. Sugapriya, D. Nagarajan, Hexagonal fuzzy programming techniques for sustainable inventory management and precredit payment in refractory systems, 2025, 2330-7706, 1, 10.1080/23307706.2024.2443510
    11. Anuz Kumar Chakrabarty, Md. Al-Amin Khan, Adel Fahad Alrasheedi, Aminur Rahman Khan, Modeling for sustainable greenhouse flower-plant production with trapezoidal type demand and partial backordering under environmental regulations, 2025, 11, 2363-6203, 10.1007/s40808-024-02194-x
    12. Fleming Akhtar, Hachen Ali, Md. Al-Amin Khan, Ali Akbar Shaikh, Optimal control approach for solving a supply chain problem under variable demand and emissions tax regulation with an unknown production rate, 2025, 2095-7513, 10.1007/s42524-025-4110-6
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1358) PDF downloads(58) Cited by(0)

Figures and Tables

Tables(11)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog