• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A Nonstandard Computational Investigation of SEIR Model with Fuzzy Transmission,Recovery and Death Rates

    2023-12-15 03:58:58AhmedMsmaliFazalDayanMuhammadRafiqNaumanAhmedAbdullahAliAhmadiniandHassanHamali
    Computers Materials&Continua 2023年11期

    Ahmed H.Msmali,Fazal Dayan,Muhammad Rafiq,Nauman Ahmed,Abdullah Ali H.Ahmadini and Hassan A.Hamali

    1Department of Mathematics,College of Science,Jazan University,Jazan,45142,Saudi Arabia

    2Department of Mathematics,School of Science,University of Management and Technology,Lahore,Pakistan

    3Department of Mathematics,Faculty of Science&Technology,University of Central Punjab,Lahore,Pakistan

    4Department of Mathematics and Statistics,University of Lahore,Lahore,Pakistan

    5Department of Medical Laboratory Technology,Faculty of Applied Medical Sciences,Jazan University,P.O.Box 1906,Jazan,45142,Saudi Arabia

    ABSTRACT In this article,a Susceptible-Exposed-Infectious-Recovered(SEIR)epidemic model is considered.The equilibrium analysis and reproduction number are studied.The conventional models have made assumptions of homogeneity in disease transmission that contradict the actual reality.However,it is crucial to consider the heterogeneity of the transmission rate when modeling disease dynamics.Describing the heterogeneity of disease transmission mathematically can be achieved by incorporating fuzzy theory.A numerical scheme nonstandard,finite difference(NSFD) approach is developed for the studied model and the results of numerical simulations are presented.Simulations of the constructed scheme are presented.The positivity,convergence and consistency of the developed technique are investigated using mathematical induction,Jacobean matrix and Taylor series expansions respectively.The suggested scheme preserves all these essential characteristics of the disease dynamical models.The numerical and simulation results reveal that the proposed NSFD method provides an adequate representation of the dynamics of the disease.Moreover,the obtained method generates plausible predictions that can be used by regulators to support the decision-making process to design and develop control strategies.Effects of the natural immunity on the infected class are studied which reveals that an increase in natural immunity can decrease the infection and vice versa.

    KEYWORDS Epidemic model;fuzzy parameters;NSFD scheme;convergence;positivity;consistency

    1 Introduction

    Mathematical modeling has played a key role in better understanding the behavior of infectious diseases.Infectious disease modeling is a powerful technique that has been used to research the dynamic transmission of infections,forecast the trajectory of an outbreak,and assess the most effective epidemic disease control measures.The transmission of infectious diseases and their prevention are frequently better understood through the use of mathematical models.To find solutions to manage diseases in the population,numerous studies have been carried out using the mathematical model method.To learn more about the dynamics of infectious disease transmission and control,numerous models on the disease’s dynamics have been created and studied.Over the past few years,there has been a collective worldwide initiative to expedite advancements in creating a global surveillance network aimed at addressing the challenges posed by emerging and re-emerging infectious diseases during pandemics.To quickly assess potentially urgent situations,scientists from medicine and molecular biology to computer science and applied mathematics have collaborated.Consequently,the utilization of mathematical modeling plays a vital role in predicting,evaluating,and managing potential outbreaks.In order to gain a deeper comprehension of contagious dynamics,it is necessary to extensively investigate the influence of various variables.These variables encompass aspects ranging from the microscopic level of host-pathogen interactions to the dynamics of transmission between hosts,alongside the broader ecological,social,economic,and demographic factors prevalent worldwide.Thorough studies focusing on these factors are essential[1,2].

    Mathematical models have been critical in understanding the spread of directly transmissible infectious diseases.Since the mid-twentieth century,mathematical epidemiology appears to have grown at an exponential rate.Many different models for infectious diseases have been developed,mathematically analyzed,and applied.The concept underlying compartmental models involve classifying the entire population according to their epidemiological condition.The type of epidemic model is determined by the number of compartments considered.A mathematical model is a condensed representation of existing clinical and biological knowledge that incorporates elements that seem pertinent to the question being researched regarding disease dynamics and progression.As a result,a model is a simplified representation of a complex reality,and its structure is determined by the processes under consideration and extrapolated[3].

    Mathematical and statistical models,as well as computational engineering,are helping to shed light on the problem and aid in decision making.In 1766,the first publication on mathematical modeling of epidemics was published.To study smallpox mortality in England,which was one in 14 of total mortality at the time,a mathematical model developed by Daniel Bernoulli.Bernoulli’s model demonstrated that immunization anti-virus treatment would add about three years to life expectancy at birth.This work’s English translation can be found in[4].The deterministic compartmental epidemic modeling was developed by Kermack et al.[5-7].Enko created a probabilistic model to describe the measles epidemic in discrete time and calculated the number of contacts between infected and susceptible people in the population[8].After that,a huge number of models were created,examined,and used both qualitatively and quantitatively to study a variety of infectious diseases.In addition to being used to compare,plan,implement,evaluate,and optimize various detection,preventive,therapeutic,and control programs,mathematical models have become useful tools for studying the transmission and management of infectious diseases.

    The susceptible,exposed,infectious,and recovered (SEIR) model is used to imitate the spread of an infectious illness within a population.It is frequently used to investigate the dynamics of infectious diseases and to assess the possible impact of various interventions such as vaccination campaigns,social distancing measures,and quarantine regimes.The compartment structure may differ depending on the nature of the disease [9].Ngwa et al.[10] developed a malaria model in an endemic scenario with a changeable human and mosquito population.Sheikh investigated an SEIR model that took into account limited treatment resources[11].Zhang et al.investigated an SEIR model that included continuous inflows of new susceptible,exposed,infective,and recovered.This model includes a disease-related death rate[12].Zhou et al.conducted a study on an SEIR epidemic model that incorporated a saturated recovery process[13].Yi et al.conducted studied the dynamic behaviors of an SEIR model that incorporated seasonal transmission rate forcing.In order to assess the dynamics of the system,the study considered scenarios with one,two,and three variable parameters[14].Sun et al.examined an SEIR model with varying population sizes and vaccination strategies.They identified three threshold parameters crucial for disease eradication,including the total number of infectives and their proportion in the population[15].Kolawole et al.investigated an SEIR model to determine the significance of saturation.The Runge-Kutta(RK)method was used to solve the model numerically.The susceptible and infected compartments of the model are significantly impacted by the saturation term,according to the results of the numerical simulation.The saturation phrase,on the other hand,affects susceptible individuals more so than infected individuals[16].Hukmah et al.developed a COVID-19 disease model by modifying the SEIR model to include quarantine and vaccination compartments [17].Chen et al.created an SEIR model for information propagation in complex networks with the effect of a hot search.Mathematical and numerical simulations were also carried out[18].Yang et al.[19]studied an SEIR epidemic model that included both sexual and non-sexual transmission routes.Mohajan investigated a COVID-19 SEIR model[20].

    Many infectious diseases can be stopped from spreading by susceptible people’s immunity,and certain infections provide recovered patients temporary or permanent protection from re-infection.The dynamic spread of the susceptible,infectious and recovered(SIR),SEIR,susceptible,vaccinated,exposed,infectious and recovered(SVEIR),and many other epidemic models had been investigated thoroughly by a number of researchers.For a better comprehension of the disease’s transmission,Olopade et al.investigated an SEIR mathematical model including natural immunity and treatment rates[21].Immunity can be obtained through targeted immunization;it can also be obtained naturally after a person has fully recovered from infection;and,in some cases,the mother’s antibodies may be transferred to the newborn to confer a certain level of immunity.Immunity duration varies depending on the disease,with some diseases providing almost lifelong immunity and others only providing a brief period of insensitivity.Vaccine-induced immunity frequently requires a boost after a period of time,such as when the vaccine’s effectiveness decreases due to a lack of disease exposure.In the case of measles,for example,vaccinated individuals have less immunity than those who have natural immunity[22,23].Some mass vaccination models were presented by Scherer and McLean,and they can be used to provide crucial insights into the anticipated effects of immunization campaigns.They demonstrate how to determine the vaccination coverage rate necessary to eradicate an infection,examine the effects of vaccine-induced immunity that deteriorates with time,and examine the interactions between strains of infectious agents that are both vaccine vulnerable and vaccine resistant[24].A mathematical model for cholera epidemics was given by Sanches et al.[25]and takes seasonality,host immunity decline,and cholera control measures into account.Innate immunity was the subject of research by Jain et al.in disease modeling[26].

    Zadeh first proposed the fuzzy theory in 1965 [27].Fuzzy theory plays an important role in mathematical modeling by providing a way to handle uncertain or ambiguous information in a mathematical framework.Traditional mathematical models assume that all variables are precise and can be accurately measured or calculated.In many real-world circumstances,however,variables may be inaccurate or impossible to assess,leading to model uncertainty.Fuzzy theory provides a mathematical framework for dealing with uncertainty by allowing variables to take on values that are not explicitly stated,but rather are characterized by degrees of membership in a set.This allows for more flexible and realistic variable representation,allowing mathematical models to reflect the complex and uncertain nature of real-world systems more accurately.In mathematical modeling,fuzzy parameters are used to indicate the degree of uncertainty or imprecision in parameter values.When the interactions between variables are complex and exact data is lacking,fuzzy parameters are extremely useful.They allow for the quantitative representation of uncertainty and imprecision,as well as the incorporation of subjective knowledge and expert opinion into mathematical models.Fuzzy theory has been applied in a wide range of mathematical modeling applications.Researchers can construct more accurate and robust models that can handle uncertainty and imprecision by combining fuzzy theory into mathematical models,resulting in better predictions and judgments.The fuzzy theory has been applied in various ways to develop and enhance epidemic models.Barros et al.[28] and Mondal et al.[29] investigated epidemic models with fuzzy transmission coefficients.Mishra et al.suggested a Susceptible Infectious-Recovered-Susceptible(SIRS)model for the fuzzy transmission of worms in a computer network[30].

    Many traditional schemes,such as forward Euler,Runge-Kutta,and others,occasionally fail by generating oscillations,chaos,and false steady states.One option for avoiding all of these numerical instabilities is to build schemes using the nonstandard finite-difference method.This technique,pioneered by Mickens[31],has resulted in the development of new numerical schemes that preserve physical qualities such as stability,positivity,and boundedness,etc.Using the NSFD approach,Verma et al.investigated the generalized Burgers-Huxley (GBH) equation [32].To assess influenza disease,Khalsaraei et al.presented an NSFD technique.Some intriguing aspects of the developed scheme are investigated,and numerical results are presented[33].Razza et al.investigated a stochastic model of HIV/AIDS based on NSFD theory in a bisexual population,taking counseling and antiretroviral therapy into account,and concluded that the stochastic model of the HIV/AIDS epidemic is more pragmatic than the deterministic model.All of the important properties of the disease dynamical model were retained by the NSFD stochastic scheme[34].Jawaz et al.investigated a delayed reactiondiffusion epidemic model and developed an NSFD-based scheme to study HIV/AIDS [35].The NSFD theory proposed by Mickens is extensively used in disease mathematical and numerical modeling[36,37],just to mention a few.The susceptible-infected-treatment-recovered(SITR)model was developed by Aslefallah et al.to explore the dynamics of the new coronavirus (COVID-19).Error analysis is performed to assess the model’s accuracy[38].Singular boundary approach is used to solve a modified anomalous diffusion process in two dimensions with a nonlinear source term,initial and Dirichlet-type boundary conditions.In terms of Riemann-Liouville fractional derivatives,they modeled the process as a two-dimensional nonlinear time-fractional sub-diffusion equation[39].The singular boundary method (SBM) is used to solve the telegraph equation in two dimensions,regardless of domain shape[40].A study focuses on the spectral meshless radial point interpolation of two-dimensional pseudo parabolic equations.To approximate time derivatives,the finite difference time-stepping approach is used.This meshless approach’s stability and convergence are thoroughly discussed and theoretically verified[41].Allehiany et al.investigated a Covid-19 model with fuzziness and numerically solved it using the NSFD scheme[42].Alhebshi et al.looked into a computer virus model that used fuzzy criteria[43].

    The novelty of the current work is the development,accomplishment,and numerical analysis of the first order explicit numerical technique with NSFD in a fuzzy environment,particularly with fuzzy parameters.The use of fuzzy theory in mathematical modeling has the potential to improve our understanding of disease transmission and to inform the development of more effective control strategies.To our knowledge,the model under study has never been studied in the literature in the NSFD and fuzzy senses,and this is the first study of this model in this sense.The primary goal of this study is to investigate an effective implicit numerical integration method for the solution of the SEIR model natural immunity and treatment rates.The convergence,positivity and consistency of the developed scheme are investigated in this work.We claim that the tool for the suggested numerical integration scheme is trustworthy,effective,and adopts all the dynamical characteristics for long-term disease behavior.

    2 SEIR Mathematical Model with Fuzzy Parameters

    In this research,the model proposed by Olopade et al.[21]is took into account.

    The fuzzy model that corresponds to the model mentioned above can be expressed as

    The disease transmission,recovery from the disease,and death due to disease rates of the infected individuals are considered fuzzy numbers due to their uncertain natures.These parameters are denoted byβ(?),θ(?)andδ(?),respectively,and are defined below:

    The death rateδ(?)will be higher at the higher virus level,i.e.,?0<?and 1-ξ,(ξ≥0)is the maximum death whileθ0>0 is the minimum recovery rate.

    2.1 Equilibrium Analysis

    Case 1:If If? <?min,we obtain

    C0(S0,E0,I0,R0)=

    Case 2: If?min <?≤?M,then we haveβ(?)=andC*(S*,E*,I*,R*) is obtained,where

    2.2 The Fuzzy Basic Reproductive Number(BRN)

    The BRNR0is given by

    SinceR0is a direct function of the virus?,it can be analyzed as follows:

    Case 1:If? <?min,then we haveβ(?)=0 and we obtainR0(?)=0.

    Case 2:If?min <?≤?M,then we haveβ(?)and we obtain

    R0(?)=

    2.3 Stability Analysis

    At disease free equilibrium(DFE),the Jacobean matrix of the system(2)is given by

    whereA1=μ+σ+k,A2=μ+δ+ωandA3=μ+θ.

    Eigenvalues of the above matrix areλ1=-μ,λ2=-A1,λ3=-A2andλ4=-A3.All eigenvalues of the Jacobean of the system (2) at DFE are less than one which shows that the system is locally asymptotically stable atC0.

    3 Numerical Modeling

    3.1 NSFD Scheme

    The NSFD scheme consists of a general set of numerical analysis methods that provide numerical solutions for differential equations by building a discrete model.R.E.Mickens was the one who first proposed it as a method to increase the accuracy and stability of numerical solutions.He used nonstandard analysis concepts,a type of mathematical logic,to create difference equations that more precisely mimic the continuous differential equations.The creation of the NSFD schemes for ordinary differential equations (ODEs) was one of Mickens’early contributions.He demonstrated that by employing nonstandard difference formulas that integrate higher-order derivatives,it is possible to improve accuracy and properly capture the oscillatory nature of solutions.Mickens extended his work to NSFD methods for partial differential equations (PDEs) in the late 1990s and early 2000s.The advantages of NSFD approaches over typical finite difference systems in terms of accuracy,stability,and efficiency are established.NSFD scheme for the system(2)is

    3.2 Positivity of the NSFD Scheme

    Given that all the state variables within the model represent segments of a population,it is necessary for at least one of them to be positive,while the other subpopulations must remain nonnegative at every instance in time.Ensuring the positivity of state variables is a crucial aspect in compartmental epidemic models.This fundamental property is thoroughly examined and upheld through the application of the principle of mathematical induction within an implicit numerical integration scheme.The following theorem is particularly advantageous in effectively addressing this concern.

    Theorem:LetS,E,IandRinvolved in the NSFD scheme are positive att=0;moreover,ifθ,π,μ,?,β,σ,k,ωandδare also positive,thenSn+1≥0,En+1≥0,In+1≥0 andRn+1≥0.

    Proof:Consider the system(7),we have

    Puttingn=1,in the above scheme,we have

    Next,assume that the above system of equations ensures that the value of has the attribute of positivity forn=2,3,4,...,n-1,i.e.,Sn≥0;En≥0;In≥0;andRn≥0;forn=2,3,4,...,n-1.

    The positivity will now be examined for a random positivity integern∈Z,we observe that

    Therefore,for any positive integer values ofn,the proposed NSFD scheme guarantees positivity for the state variables,S,E,IandR.

    3.3 Convergence Analysis

    If the Eigenvalues of the system (7) lie within the unit circle the suggested implicit numerical integration scheme demonstrates convergence.To prove this,assume that

    The Jacobean matrix corresponding to the system(7)is

    3.4 Consistency Analysis

    The Taylor’s series is used to check the consistency of the proposed scheme.From first equation of the system(7),

    For the consistency,the Taylor’s series expansion is considered as follows:

    Substituting the value ofSn+1in the above equation and after some simplifications,we get

    by applying Taylor’s series on the last equation of system (7).It is therefore concluded that our proposed scheme is consistent of order 1.

    4 Numerical Simulations

    Table 1 displays the values of parameters used for numerical simulations.

    Table 1: Values of parameters

    The graphical behavior of exposed populations at two different values of the step sizes are shown in Fig.1 at DFE and endemic equilibrium(EE)points for all three cases.The graphs show smooth,stable and converging behavior at both the step sizes and the increase in the values of the step size does not affect it.Convergence and positivity are the two major characteristics of these types of models as these models consist of the population which cannot be negative.The developed scheme preserve these feature which can be seen in these graphs.

    Fig.2 depicts the graphical behavior of infected populations at two different values of the step sizes at DFE and EE points for all three cases.The graphs exhibit convergent and stable behavior for both step sizes,and an increase in step size values has no impact on them.Since the behavior of the method is not affected with an increase in the values of the step size.It can be concluded that the NSFD theory can be used to study the long-term behavior of the model.This is an interesting feature of the developed method which many other classical methods such as Euler and many other do not keep at increasing step sizes.This behavior shows the superiority of the NSFD scheme over Euler’s method in fuzzy conditions.

    The graphs of the recovered population are shown in Fig.3 at different values of the step size for case 1,case 2 and case 3.The graphs behave similarly at different values of the step size and converge smoothly.We can conclude from this behavior that the NSFD method is capable of reflecting the dynamics of the studied model in fuzzy conditions.The typical standard schemes that exist in the literature can cause chaos and misleading variations for some passions of the discretization constraints[44].

    The effects of the natural immunity are presented in Fig.4.Low immunity causes greater infection and the infection decreases with an increase in the natural immunity.Immunity boosting is crucial in epidemiology,particularly in the context of infectious diseases.When a pathogen enters the body,the immune system mobilizes to defend against the invader.This response can include a variety of mechanisms such as antibody production,immune cell activation,and inflammation.The severity and duration of the illness can be determined by the strength of the immune response.Immune enhancement can help prevent or lessen the severity of infections.Vaccines,for example,stimulate the immune system,causing it to produce a protective response against specific pathogens.This can help prevent infection or lessen the severity of an infection if it occurs.Aside from vaccines,other ways to boost immunity include eating a healthy diet,exercising regularly,getting enough sleep,and reducing stress.These measures can help improve the immune system’s overall function,making it more capable of fighting infections.

    Figure 1:Exposed populations at different step sizes

    Figure 2:Infected populations at different step sizes

    Figure 3:Infected populations at different step sizes

    Figure 4:Effects of natural immunity on infected population

    5 Conclusion

    Epidemic modeling is a useful tool for understanding and controlling the spread of infectious diseases.In real life,however,the parameters used to model the disease may not be precisely known.Fuzzy theory can be used in such cases to incorporate uncertainty into models.In this article,an SEIR epidemic model is studied.The positivity and boundedness of the model are investigated.Equilibrium points and basic reproduction number are studied in fuzzy senses.For the numerical solution of the model,an NSFD method is developed.Positivity,convergence and consistency of the proposed method are investigated.Positivity is one of the major characteristic of these epidemic models since non positive values are meaningless due to the fact that these models deals with population which cannot be negative.Many classical standard finite difference techniques do not hold this property for some values of the discretization parameter.The proposed technique preserves positivity for all values of the discretization parameter.Convergence is the other important property in these type of models and the suggested method preserves this.

    Natural immunity is an important aspect of disease dynamics,and failing to account for it can result in a variety of modeling disasters.We have shown mathematically and numerically that increasing natural immunity can lead to a decrease in infection.Immunity boosting plays a critical role in epidemiology by preventing and controlling infectious diseases,protecting vulnerable populations,achieving herd immunity,preventing future pandemics,improving the treatment of infections,and promoting individual responsibility for health.This study could help in the modeling of other similar diseases.Natural immunity can be applied to other compartmental models.Future work can consider saturated incidence,treatment,delayed,fractional,and stochastic directions.

    Acknowledgement:The authors extend their appreciation to the Deputyship for Research and Innovation,Ministry of Education in Saudi Arabia for funding this research work through Project Number ISP22-6.

    Funding Statement:This research was funded by the Ministry of Education in Saudi Arabia of funder Grant Number ISP22-6 and the APC was funded by the Ministry of Education in Saudi Arabia.

    Author Contributions:Study conception and design: Fazal Dayan and Muhammad Rafiq;data collection:Fazal Dayan and Nauman Ahmed;analysis and interpretation of results:Abdullah Ali H.Ahmadini and Hassan A.Hamali;draft manuscript preparation: Fazal Dayan,Ahmed H.Msmali.All authors have reviewed the results and approved the final version of the manuscript.

    Availability of Data and Materials:Data will be made available on request.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    一夜夜www| 国产三级在线视频| 亚洲欧美日韩高清在线视频| 久久伊人香网站| 99精品在免费线老司机午夜| 50天的宝宝边吃奶边哭怎么回事| 精品欧美一区二区三区在线| 女人被狂操c到高潮| 真人一进一出gif抽搐免费| 91精品三级在线观看| 亚洲国产精品sss在线观看 | av国产精品久久久久影院| 九色亚洲精品在线播放| 午夜福利影视在线免费观看| 久久人人爽av亚洲精品天堂| 国产精品美女特级片免费视频播放器 | 国产一区在线观看成人免费| 午夜精品国产一区二区电影| 少妇 在线观看| 亚洲国产精品合色在线| 欧美 亚洲 国产 日韩一| 亚洲 欧美一区二区三区| 久久人人爽av亚洲精品天堂| 校园春色视频在线观看| 久久香蕉激情| 日本免费a在线| 成人免费观看视频高清| 免费高清视频大片| 亚洲国产精品999在线| 国产视频一区二区在线看| 91在线观看av| 女人爽到高潮嗷嗷叫在线视频| 一区二区三区国产精品乱码| av天堂久久9| 亚洲成a人片在线一区二区| 大型av网站在线播放| 久久精品人人爽人人爽视色| 欧美激情极品国产一区二区三区| 88av欧美| 男女床上黄色一级片免费看| 99精品欧美一区二区三区四区| 搡老岳熟女国产| 人人妻,人人澡人人爽秒播| 国产深夜福利视频在线观看| 婷婷精品国产亚洲av在线| 亚洲精品国产区一区二| 欧美日韩国产mv在线观看视频| 精品福利观看| 欧美激情高清一区二区三区| 男人舔女人下体高潮全视频| 最新美女视频免费是黄的| 久9热在线精品视频| av网站免费在线观看视频| 啦啦啦 在线观看视频| 99国产综合亚洲精品| 女生性感内裤真人,穿戴方法视频| 天天影视国产精品| 亚洲人成网站在线播放欧美日韩| 国产精品美女特级片免费视频播放器 | 高清黄色对白视频在线免费看| 女性生殖器流出的白浆| 亚洲成人国产一区在线观看| 黑人巨大精品欧美一区二区蜜桃| 美女国产高潮福利片在线看| 精品一区二区三区av网在线观看| bbb黄色大片| 亚洲av第一区精品v没综合| 亚洲成人精品中文字幕电影 | 国产成人系列免费观看| 麻豆久久精品国产亚洲av | 色婷婷av一区二区三区视频| 大型av网站在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 女人被狂操c到高潮| 熟女少妇亚洲综合色aaa.| 欧美黄色淫秽网站| 欧美中文日本在线观看视频| 婷婷丁香在线五月| 法律面前人人平等表现在哪些方面| 国产av一区在线观看免费| 国产精品一区二区在线不卡| 丁香六月欧美| 精品高清国产在线一区| 一区二区三区激情视频| av有码第一页| 久久草成人影院| 超碰成人久久| 国产免费av片在线观看野外av| 久久天躁狠狠躁夜夜2o2o| 多毛熟女@视频| 亚洲欧美精品综合一区二区三区| 国产成年人精品一区二区 | 亚洲成人免费av在线播放| 香蕉丝袜av| 久久久国产成人精品二区 | 18禁裸乳无遮挡免费网站照片 | 久久九九热精品免费| 丰满人妻熟妇乱又伦精品不卡| 亚洲色图av天堂| 久久精品国产亚洲av香蕉五月| av网站免费在线观看视频| 国内久久婷婷六月综合欲色啪| 国产精品一区二区精品视频观看| 男女床上黄色一级片免费看| 日韩一卡2卡3卡4卡2021年| 亚洲av成人一区二区三| 欧美日韩黄片免| 一进一出抽搐gif免费好疼 | 亚洲人成电影免费在线| 一级作爱视频免费观看| www.www免费av| 大香蕉久久成人网| 精品久久久久久电影网| 婷婷丁香在线五月| 中文字幕最新亚洲高清| 啦啦啦在线免费观看视频4| 欧美乱色亚洲激情| 一二三四在线观看免费中文在| 免费在线观看影片大全网站| 欧美日本亚洲视频在线播放| 男女下面进入的视频免费午夜 | 国产激情久久老熟女| 日日摸夜夜添夜夜添小说| 欧美人与性动交α欧美软件| 亚洲自拍偷在线| 亚洲国产精品一区二区三区在线| 黄片播放在线免费| 99国产极品粉嫩在线观看| 欧美日韩乱码在线| a在线观看视频网站| 国产黄色免费在线视频| 精品乱码久久久久久99久播| 狠狠狠狠99中文字幕| 一进一出抽搐动态| 精品国产一区二区久久| 日韩大码丰满熟妇| 欧美日韩一级在线毛片| 久久人妻av系列| 国产精品一区二区精品视频观看| 久久国产精品男人的天堂亚洲| 国产亚洲欧美在线一区二区| 91国产中文字幕| √禁漫天堂资源中文www| 香蕉久久夜色| 99久久人妻综合| 91精品三级在线观看| 男人的好看免费观看在线视频 | 国产高清videossex| 热99国产精品久久久久久7| 99久久人妻综合| 黄色视频不卡| 久久性视频一级片| 国产高清videossex| 激情在线观看视频在线高清| 国产野战对白在线观看| 国产精品日韩av在线免费观看 | e午夜精品久久久久久久| 久久国产亚洲av麻豆专区| 精品国产亚洲在线| 国产精品国产av在线观看| 亚洲黑人精品在线| 色综合欧美亚洲国产小说| 亚洲专区字幕在线| 亚洲国产毛片av蜜桃av| 天天添夜夜摸| 欧美日韩福利视频一区二区| 亚洲欧美日韩高清在线视频| 一边摸一边做爽爽视频免费| 在线观看免费视频网站a站| 久久人妻熟女aⅴ| 操美女的视频在线观看| 级片在线观看| 韩国av一区二区三区四区| 午夜两性在线视频| 日韩精品免费视频一区二区三区| 99精品欧美一区二区三区四区| 亚洲黑人精品在线| 夜夜夜夜夜久久久久| 99精品在免费线老司机午夜| 夜夜看夜夜爽夜夜摸 | 18美女黄网站色大片免费观看| 人妻久久中文字幕网| 国产熟女xx| 午夜精品国产一区二区电影| 一进一出好大好爽视频| 亚洲黑人精品在线| 法律面前人人平等表现在哪些方面| 欧美激情久久久久久爽电影 | 久久中文看片网| 日本wwww免费看| 亚洲午夜精品一区,二区,三区| 亚洲av成人av| 国产片内射在线| 好男人电影高清在线观看| 交换朋友夫妻互换小说| 好看av亚洲va欧美ⅴa在| 80岁老熟妇乱子伦牲交| 亚洲一区高清亚洲精品| 午夜精品国产一区二区电影| 黑丝袜美女国产一区| 日韩一卡2卡3卡4卡2021年| 高清黄色对白视频在线免费看| 五月开心婷婷网| bbb黄色大片| 亚洲色图综合在线观看| 精品少妇一区二区三区视频日本电影| 精品一区二区三卡| 亚洲一码二码三码区别大吗| 国产激情久久老熟女| 色尼玛亚洲综合影院| 国产国语露脸激情在线看| 国产精品一区二区免费欧美| 亚洲精品成人av观看孕妇| 18禁国产床啪视频网站| 亚洲五月婷婷丁香| 黄频高清免费视频| 欧美人与性动交α欧美精品济南到| 丰满人妻熟妇乱又伦精品不卡| 无限看片的www在线观看| 亚洲精品中文字幕在线视频| 一a级毛片在线观看| 嫩草影院精品99| 18禁黄网站禁片午夜丰满| 国产aⅴ精品一区二区三区波| 黄色成人免费大全| 久久亚洲精品不卡| 老鸭窝网址在线观看| 高清在线国产一区| 少妇 在线观看| 黄片小视频在线播放| 嫁个100分男人电影在线观看| 高清在线国产一区| 国产亚洲av高清不卡| 88av欧美| 一级作爱视频免费观看| 免费一级毛片在线播放高清视频 | 变态另类成人亚洲欧美熟女 | 亚洲精品国产精品久久久不卡| 中亚洲国语对白在线视频| 美女午夜性视频免费| 亚洲av电影在线进入| 很黄的视频免费| 亚洲欧美一区二区三区久久| 妹子高潮喷水视频| 日韩三级视频一区二区三区| 日韩一卡2卡3卡4卡2021年| 一区二区三区国产精品乱码| 亚洲欧美日韩无卡精品| 久久人妻熟女aⅴ| 美女 人体艺术 gogo| 精品久久久久久成人av| 美女国产高潮福利片在线看| 一区福利在线观看| 国产不卡一卡二| 亚洲九九香蕉| 女人被狂操c到高潮| 色老头精品视频在线观看| 久久精品aⅴ一区二区三区四区| 夜夜爽天天搞| 国产高清国产精品国产三级| 日本黄色视频三级网站网址| 老汉色av国产亚洲站长工具| 老司机午夜十八禁免费视频| 又黄又粗又硬又大视频| 国产不卡一卡二| www.www免费av| 伦理电影免费视频| 亚洲 欧美 日韩 在线 免费| 宅男免费午夜| 亚洲精华国产精华精| 正在播放国产对白刺激| 欧美日韩黄片免| 90打野战视频偷拍视频| 很黄的视频免费| 久久久久久久久中文| 午夜影院日韩av| 久久久久久久精品吃奶| 成人18禁在线播放| 亚洲成人免费电影在线观看| 88av欧美| bbb黄色大片| 超碰成人久久| 一a级毛片在线观看| 最近最新中文字幕大全免费视频| 国产亚洲欧美精品永久| 侵犯人妻中文字幕一二三四区| 免费高清视频大片| 两人在一起打扑克的视频| 黑人操中国人逼视频| 久久狼人影院| 国产成人欧美在线观看| 久久人人精品亚洲av| 99在线视频只有这里精品首页| 国产精品免费视频内射| 99国产精品一区二区蜜桃av| 超碰成人久久| 一夜夜www| 精品卡一卡二卡四卡免费| www.999成人在线观看| 日本三级黄在线观看| 99国产极品粉嫩在线观看| 人成视频在线观看免费观看| 日本撒尿小便嘘嘘汇集6| 日韩大码丰满熟妇| 男女高潮啪啪啪动态图| 一本综合久久免费| 久久亚洲精品不卡| 国产国语露脸激情在线看| 欧美 亚洲 国产 日韩一| 男女做爰动态图高潮gif福利片 | 女人被躁到高潮嗷嗷叫费观| 国产成人精品在线电影| 久久久久久久久免费视频了| 人人妻人人爽人人添夜夜欢视频| 久久影院123| 国产一区在线观看成人免费| 成人三级做爰电影| 久久人人精品亚洲av| 黄色片一级片一级黄色片| 色哟哟哟哟哟哟| 在线观看午夜福利视频| 精品少妇一区二区三区视频日本电影| 视频在线观看一区二区三区| 久热爱精品视频在线9| 国产成人精品久久二区二区91| 精品久久久精品久久久| 中文字幕另类日韩欧美亚洲嫩草| 又黄又粗又硬又大视频| 亚洲欧美一区二区三区黑人| 欧美性长视频在线观看| av在线天堂中文字幕 | 亚洲男人天堂网一区| svipshipincom国产片| 日本欧美视频一区| 欧美另类亚洲清纯唯美| 国产又色又爽无遮挡免费看| 国产亚洲欧美精品永久| 一区二区日韩欧美中文字幕| 丰满人妻熟妇乱又伦精品不卡| 国产精品香港三级国产av潘金莲| 欧洲精品卡2卡3卡4卡5卡区| 天天影视国产精品| 国产熟女午夜一区二区三区| 我的亚洲天堂| 老司机福利观看| 老熟妇乱子伦视频在线观看| 亚洲成av片中文字幕在线观看| 一进一出抽搐动态| 久久久久国产精品人妻aⅴ院| 午夜亚洲福利在线播放| 午夜福利影视在线免费观看| 日本a在线网址| 国产成+人综合+亚洲专区| 国产精品自产拍在线观看55亚洲| 女人爽到高潮嗷嗷叫在线视频| 嫩草影院精品99| 精品卡一卡二卡四卡免费| 久久天堂一区二区三区四区| 18禁裸乳无遮挡免费网站照片 | 国产99久久九九免费精品| 亚洲一区二区三区色噜噜 | 久99久视频精品免费| 成人永久免费在线观看视频| 精品乱码久久久久久99久播| 国产区一区二久久| 亚洲精品成人av观看孕妇| 淫秽高清视频在线观看| 美女扒开内裤让男人捅视频| 一区二区三区激情视频| 精品乱码久久久久久99久播| 国产野战对白在线观看| 国产亚洲精品久久久久5区| 国产av一区在线观看免费| 最好的美女福利视频网| 母亲3免费完整高清在线观看| 日韩大尺度精品在线看网址 | av免费在线观看网站| 在线天堂中文资源库| 久久久久国产一级毛片高清牌| 超碰97精品在线观看| а√天堂www在线а√下载| 国产激情欧美一区二区| 亚洲av日韩精品久久久久久密| 日日摸夜夜添夜夜添小说| 国产亚洲欧美在线一区二区| 欧美乱码精品一区二区三区| 亚洲av第一区精品v没综合| 欧美日韩瑟瑟在线播放| 18美女黄网站色大片免费观看| 精品少妇一区二区三区视频日本电影| 69av精品久久久久久| 成人精品一区二区免费| 欧洲精品卡2卡3卡4卡5卡区| 一二三四在线观看免费中文在| 成人免费观看视频高清| 最新在线观看一区二区三区| 夜夜看夜夜爽夜夜摸 | 久久精品亚洲av国产电影网| 757午夜福利合集在线观看| 淫秽高清视频在线观看| 大码成人一级视频| 女人被狂操c到高潮| 亚洲欧美激情在线| 国产精品久久视频播放| 在线观看66精品国产| 黄片大片在线免费观看| 一进一出抽搐gif免费好疼 | 国产又色又爽无遮挡免费看| 色尼玛亚洲综合影院| 亚洲成人国产一区在线观看| 精品国内亚洲2022精品成人| 国产精品自产拍在线观看55亚洲| 色婷婷久久久亚洲欧美| 日韩高清综合在线| 丰满饥渴人妻一区二区三| 啪啪无遮挡十八禁网站| 日韩免费高清中文字幕av| 午夜91福利影院| 精品一品国产午夜福利视频| 免费在线观看完整版高清| 黑人猛操日本美女一级片| 波多野结衣av一区二区av| 97超级碰碰碰精品色视频在线观看| 999久久久精品免费观看国产| 一a级毛片在线观看| x7x7x7水蜜桃| 热re99久久精品国产66热6| 久久天躁狠狠躁夜夜2o2o| 欧美日韩精品网址| 国产免费现黄频在线看| 法律面前人人平等表现在哪些方面| 18禁观看日本| 欧美日韩福利视频一区二区| 激情视频va一区二区三区| 国产亚洲欧美在线一区二区| 90打野战视频偷拍视频| 性色av乱码一区二区三区2| 欧美日韩瑟瑟在线播放| 久久久久九九精品影院| 午夜精品国产一区二区电影| 性欧美人与动物交配| 国产又色又爽无遮挡免费看| 两个人看的免费小视频| 久久香蕉国产精品| 国产精品久久久久成人av| 大型av网站在线播放| 欧美成人午夜精品| 无人区码免费观看不卡| 精品国产一区二区三区四区第35| 成人18禁高潮啪啪吃奶动态图| 欧美日本中文国产一区发布| 日韩免费高清中文字幕av| 一个人观看的视频www高清免费观看 | 丝袜在线中文字幕| 18禁黄网站禁片午夜丰满| 亚洲成人免费av在线播放| 精品久久久久久成人av| 国产精品一区二区三区四区久久 | 国产成人一区二区三区免费视频网站| 免费少妇av软件| 亚洲精品国产色婷婷电影| 亚洲av成人av| 成人永久免费在线观看视频| 亚洲一区二区三区欧美精品| 一进一出抽搐gif免费好疼 | 国产成年人精品一区二区 | 久久亚洲真实| 亚洲精品中文字幕在线视频| 欧美丝袜亚洲另类 | 亚洲七黄色美女视频| 黑人巨大精品欧美一区二区mp4| 男人操女人黄网站| 亚洲自拍偷在线| 夫妻午夜视频| 国产午夜精品久久久久久| 热re99久久国产66热| 性色av乱码一区二区三区2| 久热爱精品视频在线9| 久9热在线精品视频| 好男人电影高清在线观看| 免费在线观看黄色视频的| 久久婷婷成人综合色麻豆| 久久久久久人人人人人| 电影成人av| 超碰成人久久| 日韩视频一区二区在线观看| 久久久国产成人精品二区 | 人成视频在线观看免费观看| 成人永久免费在线观看视频| 亚洲av电影在线进入| 亚洲 欧美一区二区三区| 亚洲久久久国产精品| 精品国产国语对白av| 免费搜索国产男女视频| 每晚都被弄得嗷嗷叫到高潮| 成年女人毛片免费观看观看9| 久久天堂一区二区三区四区| 丝袜美足系列| 亚洲五月天丁香| 一进一出抽搐gif免费好疼 | www.999成人在线观看| 在线观看舔阴道视频| 国产精品 国内视频| 国产精品综合久久久久久久免费 | 国产精品一区二区免费欧美| 亚洲九九香蕉| 丰满饥渴人妻一区二区三| 性少妇av在线| 国产一卡二卡三卡精品| 亚洲欧美一区二区三区久久| 日韩高清综合在线| 亚洲精品一区av在线观看| 久久久久精品国产欧美久久久| 久久香蕉国产精品| 欧美日韩av久久| 99精品久久久久人妻精品| 国产一区在线观看成人免费| 精品国内亚洲2022精品成人| 国产成人精品无人区| 一区二区三区激情视频| 欧美黄色淫秽网站| 人人妻人人爽人人添夜夜欢视频| 一本大道久久a久久精品| 国产精华一区二区三区| 亚洲成人免费av在线播放| 久久国产精品人妻蜜桃| 国产精品av久久久久免费| 亚洲欧美激情在线| 最新在线观看一区二区三区| 亚洲男人天堂网一区| 高清欧美精品videossex| 桃红色精品国产亚洲av| 午夜影院日韩av| 精品第一国产精品| 黄频高清免费视频| 老司机靠b影院| 欧美一级毛片孕妇| 国产一区二区三区综合在线观看| 欧美激情久久久久久爽电影 | 看黄色毛片网站| 一个人观看的视频www高清免费观看 | 亚洲精品久久午夜乱码| 叶爱在线成人免费视频播放| 精品国产超薄肉色丝袜足j| 十分钟在线观看高清视频www| 亚洲精品一二三| 国产亚洲精品一区二区www| 高清毛片免费观看视频网站 | 女警被强在线播放| 婷婷六月久久综合丁香| 一区二区三区国产精品乱码| 国产黄色免费在线视频| 91成人精品电影| 久久人妻av系列| 国产伦人伦偷精品视频| 12—13女人毛片做爰片一| 免费久久久久久久精品成人欧美视频| 亚洲精品粉嫩美女一区| 精品福利观看| 日韩一卡2卡3卡4卡2021年| 日韩精品免费视频一区二区三区| √禁漫天堂资源中文www| 国产一卡二卡三卡精品| 日本一区二区免费在线视频| 精品久久久精品久久久| 久久人妻福利社区极品人妻图片| 成人永久免费在线观看视频| 在线国产一区二区在线| 国产一区二区三区综合在线观看| 国产精品野战在线观看 | 极品人妻少妇av视频| 老司机靠b影院| 视频区图区小说| 欧美在线黄色| 国产高清videossex| 他把我摸到了高潮在线观看| 国产激情欧美一区二区| 亚洲av第一区精品v没综合| 成人特级黄色片久久久久久久| 制服诱惑二区| av中文乱码字幕在线| 黑人猛操日本美女一级片| 欧美黑人欧美精品刺激| 国产精品电影一区二区三区| 国产一区二区在线av高清观看| 国产熟女xx| 亚洲aⅴ乱码一区二区在线播放 | 午夜福利欧美成人| 制服人妻中文乱码| 国产三级黄色录像| 一本大道久久a久久精品| 日韩人妻精品一区2区三区| 精品一区二区三区视频在线观看免费 | 久久午夜综合久久蜜桃| 欧美日韩亚洲综合一区二区三区_| 欧美日韩黄片免| www.精华液| 丝袜美腿诱惑在线| 欧美亚洲日本最大视频资源| 亚洲男人的天堂狠狠| 亚洲人成77777在线视频| 欧美日韩黄片免| 夫妻午夜视频| 免费观看人在逋| 亚洲专区中文字幕在线| 中文字幕最新亚洲高清| 亚洲国产精品一区二区三区在线| 97碰自拍视频| 精品一区二区三区av网在线观看| 夜夜爽天天搞| 无遮挡黄片免费观看| 国产成人免费无遮挡视频| a级毛片在线看网站| 亚洲成人精品中文字幕电影 | 日韩欧美三级三区|