• <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.

    精品亚洲成国产av| 亚洲欧洲精品一区二区精品久久久| 国产高清国产精品国产三级| 欧美国产精品va在线观看不卡| 亚洲国产欧美日韩在线播放| 亚洲视频免费观看视频| 久久精品国产亚洲av高清一级| 一区在线观看完整版| 国产精品欧美亚洲77777| 色播在线永久视频| 老司机靠b影院| 久久久精品免费免费高清| 叶爱在线成人免费视频播放| 一区二区av电影网| 91精品伊人久久大香线蕉| 一级黄色大片毛片| 午夜福利,免费看| 狠狠婷婷综合久久久久久88av| 日本色播在线视频| 9色porny在线观看| 久久久久国产精品人妻一区二区| 国产真人三级小视频在线观看| 岛国毛片在线播放| 国产极品粉嫩免费观看在线| 在线观看免费高清a一片| 啦啦啦啦在线视频资源| 午夜免费观看性视频| 久久人妻福利社区极品人妻图片 | 日日夜夜操网爽| 国产精品国产三级国产专区5o| 日本欧美视频一区| 亚洲七黄色美女视频| 日本vs欧美在线观看视频| 日日摸夜夜添夜夜爱| 亚洲精品第二区| 欧美97在线视频| 欧美日韩综合久久久久久| 麻豆国产av国片精品| 亚洲欧洲精品一区二区精品久久久| 热99国产精品久久久久久7| 亚洲国产精品国产精品| 久久中文字幕一级| e午夜精品久久久久久久| 久久午夜综合久久蜜桃| 久久久精品区二区三区| 亚洲男人天堂网一区| 中国美女看黄片| 熟女少妇亚洲综合色aaa.| av不卡在线播放| 丁香六月天网| 久久国产精品大桥未久av| 国产成人一区二区三区免费视频网站 | 一区二区三区四区激情视频| 老司机靠b影院| 亚洲av美国av| 亚洲熟女精品中文字幕| 免费一级毛片在线播放高清视频 | 欧美精品啪啪一区二区三区 | 宅男免费午夜| 亚洲av欧美aⅴ国产| 欧美黄色淫秽网站| 色婷婷久久久亚洲欧美| 亚洲人成网站在线观看播放| 亚洲av欧美aⅴ国产| 又粗又硬又长又爽又黄的视频| 高清欧美精品videossex| 嫩草影视91久久| 性高湖久久久久久久久免费观看| 亚洲欧洲国产日韩| 高清不卡的av网站| 韩国精品一区二区三区| 欧美日韩亚洲综合一区二区三区_| 精品熟女少妇八av免费久了| 国产成人a∨麻豆精品| 飞空精品影院首页| 欧美黄色片欧美黄色片| bbb黄色大片| 午夜福利视频精品| 精品人妻一区二区三区麻豆| 亚洲国产精品一区二区三区在线| 一级片'在线观看视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产毛片av蜜桃av| 亚洲国产欧美一区二区综合| 中国国产av一级| 一级黄片播放器| 久久精品aⅴ一区二区三区四区| 在线看a的网站| 狂野欧美激情性xxxx| 99九九在线精品视频| 精品人妻一区二区三区麻豆| 亚洲欧美一区二区三区黑人| 高清av免费在线| 老司机深夜福利视频在线观看 | 精品国产一区二区三区久久久樱花| 久久久精品国产亚洲av高清涩受| 亚洲中文日韩欧美视频| 久久国产精品男人的天堂亚洲| 欧美激情高清一区二区三区| 可以免费在线观看a视频的电影网站| 国产成人精品久久二区二区免费| 女人高潮潮喷娇喘18禁视频| 国产免费福利视频在线观看| 亚洲精品国产色婷婷电影| 国产精品一区二区精品视频观看| 两个人免费观看高清视频| 欧美精品人与动牲交sv欧美| 国产激情久久老熟女| 丝袜脚勾引网站| 久久久欧美国产精品| 久久天躁狠狠躁夜夜2o2o | 亚洲五月色婷婷综合| 国产又爽黄色视频| 性色av乱码一区二区三区2| 国产一区亚洲一区在线观看| 九色亚洲精品在线播放| 国产日韩欧美在线精品| 一级毛片黄色毛片免费观看视频| 国产麻豆69| 王馨瑶露胸无遮挡在线观看| 国产一区有黄有色的免费视频| 免费女性裸体啪啪无遮挡网站| 日韩大片免费观看网站| 黄色视频在线播放观看不卡| 日韩欧美一区视频在线观看| 久久精品成人免费网站| 欧美久久黑人一区二区| 亚洲专区中文字幕在线| 日本a在线网址| 国产亚洲欧美精品永久| 亚洲欧美中文字幕日韩二区| 伊人久久大香线蕉亚洲五| 亚洲av电影在线进入| 少妇人妻 视频| 精品久久久精品久久久| 欧美+亚洲+日韩+国产| 国产伦理片在线播放av一区| 亚洲欧美成人综合另类久久久| 在线观看免费高清a一片| 97在线人人人人妻| 亚洲黑人精品在线| 黄色a级毛片大全视频| 丝袜在线中文字幕| 国产有黄有色有爽视频| av一本久久久久| 亚洲精品在线美女| 大陆偷拍与自拍| 久久久久视频综合| 亚洲精品日本国产第一区| 美女午夜性视频免费| 国产精品 国内视频| 丰满迷人的少妇在线观看| 亚洲精品国产av蜜桃| 大陆偷拍与自拍| 大码成人一级视频| 97精品久久久久久久久久精品| 亚洲第一青青草原| 亚洲欧洲国产日韩| a 毛片基地| 欧美老熟妇乱子伦牲交| 午夜av观看不卡| 少妇猛男粗大的猛烈进出视频| 日韩制服丝袜自拍偷拍| 亚洲国产精品一区二区三区在线| 国产成人av激情在线播放| 欧美日韩亚洲高清精品| 国产伦人伦偷精品视频| 精品第一国产精品| 久久人人97超碰香蕉20202| 欧美激情高清一区二区三区| 欧美日韩视频高清一区二区三区二| 91字幕亚洲| 国产av精品麻豆| 一级片'在线观看视频| 一级黄片播放器| 亚洲国产日韩一区二区| 久久久久视频综合| 女人被躁到高潮嗷嗷叫费观| 国产成人精品久久久久久| 欧美国产精品一级二级三级| 999精品在线视频| 国产成人精品在线电影| 国产片内射在线| 999精品在线视频| 香蕉丝袜av| 夫妻午夜视频| 亚洲国产欧美日韩在线播放| 人成视频在线观看免费观看| 亚洲,一卡二卡三卡| 丝袜喷水一区| 久久精品亚洲av国产电影网| 好男人视频免费观看在线| 国产又爽黄色视频| 69精品国产乱码久久久| www日本在线高清视频| 久久毛片免费看一区二区三区| 欧美亚洲 丝袜 人妻 在线| 免费在线观看日本一区| 国产精品麻豆人妻色哟哟久久| 久久久亚洲精品成人影院| 亚洲国产精品一区三区| tube8黄色片| 色精品久久人妻99蜜桃| 丝袜美足系列| 1024香蕉在线观看| 一本—道久久a久久精品蜜桃钙片| 日韩 亚洲 欧美在线| 国产一区二区激情短视频 | 一区二区三区乱码不卡18| 男女无遮挡免费网站观看| 91国产中文字幕| 王馨瑶露胸无遮挡在线观看| 亚洲成人国产一区在线观看 | 欧美av亚洲av综合av国产av| 日本色播在线视频| 亚洲欧美精品综合一区二区三区| 十八禁网站网址无遮挡| 精品久久蜜臀av无| 久久人人97超碰香蕉20202| 大型av网站在线播放| 国产精品一区二区免费欧美 | 丝袜美腿诱惑在线| 国产日韩一区二区三区精品不卡| 日本五十路高清| 一级,二级,三级黄色视频| 精品卡一卡二卡四卡免费| 80岁老熟妇乱子伦牲交| 一级毛片 在线播放| 国产精品一区二区精品视频观看| av有码第一页| 亚洲情色 制服丝袜| 国产爽快片一区二区三区| 性色av乱码一区二区三区2| 国产日韩一区二区三区精品不卡| 国产精品久久久av美女十八| 一区二区三区精品91| 亚洲国产最新在线播放| 捣出白浆h1v1| 黑人欧美特级aaaaaa片| 青春草视频在线免费观看| 欧美日韩一级在线毛片| 亚洲熟女毛片儿| 人人妻人人澡人人爽人人夜夜| 国产人伦9x9x在线观看| 成年动漫av网址| bbb黄色大片| 欧美成人午夜精品| 一级片免费观看大全| 精品少妇内射三级| 成年人免费黄色播放视频| 美女扒开内裤让男人捅视频| 91成人精品电影| 两性夫妻黄色片| 国产精品久久久久久精品电影小说| 国产欧美亚洲国产| 日本vs欧美在线观看视频| 脱女人内裤的视频| 亚洲一区二区三区欧美精品| 欧美在线一区亚洲| 脱女人内裤的视频| 免费久久久久久久精品成人欧美视频| 一区在线观看完整版| 久久国产亚洲av麻豆专区| 少妇粗大呻吟视频| 成在线人永久免费视频| 国产精品免费大片| 亚洲欧美精品综合一区二区三区| 久久毛片免费看一区二区三区| 在线观看免费午夜福利视频| 啦啦啦在线免费观看视频4| 午夜福利影视在线免费观看| 丝瓜视频免费看黄片| 波野结衣二区三区在线| 各种免费的搞黄视频| 日本欧美国产在线视频| 美女中出高潮动态图| 91国产中文字幕| 亚洲av片天天在线观看| 丁香六月天网| 久久天躁狠狠躁夜夜2o2o | 欧美国产精品va在线观看不卡| 欧美日韩黄片免| 99国产精品免费福利视频| a级毛片在线看网站| 国产极品粉嫩免费观看在线| 中文字幕最新亚洲高清| 多毛熟女@视频| 美女扒开内裤让男人捅视频| 国产精品一区二区免费欧美 | 极品人妻少妇av视频| 天堂8中文在线网| 亚洲欧美激情在线| avwww免费| 亚洲激情五月婷婷啪啪| 日韩免费高清中文字幕av| 久久天堂一区二区三区四区| 亚洲综合色网址| 免费在线观看影片大全网站 | 亚洲精品国产av成人精品| 日韩,欧美,国产一区二区三区| 男女国产视频网站| 欧美在线黄色| 一本—道久久a久久精品蜜桃钙片| 天天躁夜夜躁狠狠躁躁| 国产欧美日韩综合在线一区二区| 免费看av在线观看网站| 欧美日韩成人在线一区二区| 一级黄片播放器| 亚洲av日韩在线播放| 亚洲成人手机| 国产成人精品在线电影| 精品卡一卡二卡四卡免费| 在线天堂中文资源库| av天堂久久9| 国产日韩欧美亚洲二区| 国产精品久久久久久人妻精品电影 | 亚洲一码二码三码区别大吗| 亚洲精品日本国产第一区| 久久人妻福利社区极品人妻图片 | 男的添女的下面高潮视频| netflix在线观看网站| www.999成人在线观看| 多毛熟女@视频| 精品一区二区三卡| 国产成人免费观看mmmm| 午夜免费成人在线视频| 免费在线观看黄色视频的| 久久国产精品影院| 母亲3免费完整高清在线观看| 国产精品久久久久成人av| 亚洲av电影在线进入| 国精品久久久久久国模美| 亚洲国产精品999| 中文精品一卡2卡3卡4更新| 一个人免费看片子| xxx大片免费视频| 新久久久久国产一级毛片| 欧美日韩一级在线毛片| 91字幕亚洲| 亚洲天堂av无毛| 日韩一区二区三区影片| 免费看十八禁软件| a 毛片基地| 男女下面插进去视频免费观看| 亚洲自偷自拍图片 自拍| 成在线人永久免费视频| 日韩免费高清中文字幕av| 免费高清在线观看日韩| 午夜福利视频精品| 夫妻性生交免费视频一级片| 好男人视频免费观看在线| 看十八女毛片水多多多| 亚洲国产欧美网| 女性生殖器流出的白浆| 亚洲中文日韩欧美视频| 最近中文字幕2019免费版| netflix在线观看网站| 久久女婷五月综合色啪小说| 王馨瑶露胸无遮挡在线观看| 90打野战视频偷拍视频| 人妻人人澡人人爽人人| 亚洲七黄色美女视频| 日本欧美国产在线视频| 国产精品久久久久久人妻精品电影 | 日韩大片免费观看网站| 日本欧美视频一区| 亚洲精品国产一区二区精华液| 亚洲国产看品久久| www.精华液| 精品一品国产午夜福利视频| 丰满人妻熟妇乱又伦精品不卡| 一区福利在线观看| 中文字幕最新亚洲高清| 91精品三级在线观看| 校园人妻丝袜中文字幕| 女性生殖器流出的白浆| 一区二区三区乱码不卡18| 香蕉丝袜av| 高清欧美精品videossex| 丝袜美腿诱惑在线| 国产老妇伦熟女老妇高清| 午夜福利免费观看在线| 男女高潮啪啪啪动态图| av线在线观看网站| 大型av网站在线播放| 亚洲欧美日韩高清在线视频 | 色婷婷av一区二区三区视频| 中文字幕制服av| 如日韩欧美国产精品一区二区三区| 一二三四在线观看免费中文在| 满18在线观看网站| 免费日韩欧美在线观看| 高清黄色对白视频在线免费看| 亚洲一区中文字幕在线| 亚洲激情五月婷婷啪啪| 热re99久久精品国产66热6| 国产精品一国产av| 久久久久久久国产电影| 亚洲伊人久久精品综合| 国产又爽黄色视频| 丰满人妻熟妇乱又伦精品不卡| 久久天堂一区二区三区四区| 国产成人av激情在线播放| 中文字幕人妻丝袜一区二区| xxxhd国产人妻xxx| 成人亚洲精品一区在线观看| 又黄又粗又硬又大视频| av网站免费在线观看视频| 一级黄色大片毛片| 久久天躁狠狠躁夜夜2o2o | 欧美成狂野欧美在线观看| 悠悠久久av| 日本av免费视频播放| 久久久久久人人人人人| 99久久精品国产亚洲精品| 精品亚洲成a人片在线观看| 成人国语在线视频| 在线观看免费视频网站a站| 亚洲熟女毛片儿| 久久久亚洲精品成人影院| 日韩一区二区三区影片| 亚洲熟女精品中文字幕| 亚洲精品一区蜜桃| 午夜激情久久久久久久| 另类亚洲欧美激情| 一级毛片电影观看| 黄片小视频在线播放| 午夜免费鲁丝| 黑丝袜美女国产一区| 亚洲国产最新在线播放| av网站免费在线观看视频| 99热全是精品| 麻豆国产av国片精品| 亚洲精品国产av蜜桃| 免费不卡黄色视频| 日韩视频在线欧美| 久久久精品94久久精品| 久久亚洲国产成人精品v| av在线播放精品| 久久久久精品人妻al黑| 久久久久久久精品精品| 精品国产一区二区三区久久久樱花| 国产xxxxx性猛交| 免费日韩欧美在线观看| 国产精品av久久久久免费| 精品高清国产在线一区| 免费女性裸体啪啪无遮挡网站| 欧美日韩精品网址| 嫁个100分男人电影在线观看 | 性少妇av在线| 美女午夜性视频免费| 国产高清国产精品国产三级| av有码第一页| 免费日韩欧美在线观看| 又大又爽又粗| 999久久久国产精品视频| 最新在线观看一区二区三区 | 18禁黄网站禁片午夜丰满| 日韩制服丝袜自拍偷拍| av网站在线播放免费| 18禁国产床啪视频网站| 免费一级毛片在线播放高清视频 | 蜜桃在线观看..| av网站在线播放免费| 熟女av电影| 久久 成人 亚洲| 国产亚洲精品久久久久5区| 黄色视频在线播放观看不卡| 日本色播在线视频| 老熟女久久久| 免费久久久久久久精品成人欧美视频| 欧美激情 高清一区二区三区| 啦啦啦中文免费视频观看日本| 老司机影院毛片| 9热在线视频观看99| 成人手机av| 国产在线一区二区三区精| 操美女的视频在线观看| 精品亚洲成a人片在线观看| 欧美中文综合在线视频| www.999成人在线观看| 精品久久久精品久久久| 久久精品成人免费网站| 日韩制服丝袜自拍偷拍| 国产成人啪精品午夜网站| 一本一本久久a久久精品综合妖精| 日日摸夜夜添夜夜爱| 天堂中文最新版在线下载| 麻豆乱淫一区二区| 人人妻人人澡人人看| 国产亚洲av片在线观看秒播厂| 欧美另类一区| 亚洲中文av在线| 色婷婷久久久亚洲欧美| 亚洲熟女精品中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 精品第一国产精品| 丝袜美腿诱惑在线| 亚洲欧美激情在线| av片东京热男人的天堂| 国产精品久久久人人做人人爽| 亚洲 欧美一区二区三区| 国产精品av久久久久免费| 久久久久国产精品人妻一区二区| bbb黄色大片| 国产精品香港三级国产av潘金莲 | 成年人免费黄色播放视频| 高清欧美精品videossex| 日韩大码丰满熟妇| 丝瓜视频免费看黄片| 欧美中文综合在线视频| 女人久久www免费人成看片| 中文乱码字字幕精品一区二区三区| 久久99精品国语久久久| 国产野战对白在线观看| 中国国产av一级| 成在线人永久免费视频| 麻豆乱淫一区二区| 国产一区二区在线观看av| avwww免费| 一区二区三区四区激情视频| cao死你这个sao货| 午夜精品国产一区二区电影| 久久久国产一区二区| 亚洲成av片中文字幕在线观看| 国产在视频线精品| 亚洲精品久久久久久婷婷小说| 国产在线一区二区三区精| 男女国产视频网站| 成在线人永久免费视频| 建设人人有责人人尽责人人享有的| 亚洲国产av新网站| 91麻豆精品激情在线观看国产 | 亚洲午夜精品一区,二区,三区| 日本av免费视频播放| 国产精品麻豆人妻色哟哟久久| 青青草视频在线视频观看| 99精国产麻豆久久婷婷| 天堂俺去俺来也www色官网| 极品人妻少妇av视频| 国产不卡av网站在线观看| 精品卡一卡二卡四卡免费| 久久精品成人免费网站| 国产日韩欧美亚洲二区| 国产精品一区二区在线不卡| 国产精品.久久久| 精品福利观看| 国产精品秋霞免费鲁丝片| 热99久久久久精品小说推荐| 美女福利国产在线| 亚洲精品日韩在线中文字幕| av天堂在线播放| 成人18禁高潮啪啪吃奶动态图| 精品国产一区二区三区久久久樱花| 麻豆乱淫一区二区| 各种免费的搞黄视频| www.999成人在线观看| 中文乱码字字幕精品一区二区三区| 99精品久久久久人妻精品| h视频一区二区三区| 久久 成人 亚洲| 欧美精品人与动牲交sv欧美| 99久久人妻综合| 国产视频首页在线观看| 视频区欧美日本亚洲| 欧美精品亚洲一区二区| 午夜老司机福利片| 免费看不卡的av| 午夜91福利影院| 中文欧美无线码| 天天躁狠狠躁夜夜躁狠狠躁| 久久青草综合色| 国产精品久久久久久人妻精品电影 | 国产高清videossex| 国产精品一区二区在线观看99| 国产日韩欧美视频二区| 精品久久蜜臀av无| 亚洲中文av在线| 国产亚洲av高清不卡| 两人在一起打扑克的视频| 黄片小视频在线播放| 欧美性长视频在线观看| 久9热在线精品视频| 精品人妻在线不人妻| 午夜日韩欧美国产| 18在线观看网站| 精品亚洲成国产av| 国产精品一区二区免费欧美 | 国产亚洲午夜精品一区二区久久| 女性生殖器流出的白浆| 亚洲成人免费电影在线观看 | 中文字幕人妻丝袜一区二区| 曰老女人黄片| 色婷婷av一区二区三区视频| 男女床上黄色一级片免费看| 热99国产精品久久久久久7| 又大又黄又爽视频免费| 另类亚洲欧美激情| 99国产精品一区二区三区| 久久精品久久久久久噜噜老黄| 1024香蕉在线观看| 免费在线观看视频国产中文字幕亚洲 | 操美女的视频在线观看| 99国产精品免费福利视频| 国产精品一区二区精品视频观看| 精品少妇久久久久久888优播| 久久久久久久久久久久大奶| 精品人妻1区二区| 黄色怎么调成土黄色| 国产伦人伦偷精品视频| 亚洲 国产 在线| 欧美精品av麻豆av| 国产有黄有色有爽视频| 中文字幕制服av| 亚洲av电影在线观看一区二区三区|