روش برون‌یابی برای حل عددی یک مدل بیماری‌های عفونی بومی

در این مقاله یک مدل بیماری عفونی را که به شکل دستگاه معادلات انتگرال ولترای نوع دوم غیرخطی است در نظر می‌گیریم. این مدل از نوع مدل sir است. روش برون‌یابی به‌سوی حد ریچاردسون را برای حل این سیستم طوری طراحی می‌کنیم که با یک فرایند تکرار، سیستم غیرخطی با درجۀ دقت خوب قابل حل باشد. الگوریتم حل چنین سیستم‌هایی به‌طور کامل تشریح می‌شود. این الگوریتم دارای نوعی ساختار تو در تو است که باعث می‌شود از اطلاعات پیشین در زمان‌های بعدی بتوان استفاده کرد و همین امر برنامه‌نویسی این الگوریتم را جالب کرده است. این الگوریتم با هر زبان سطح بالا قابل برنامه‌نویسی است، که ما این فرآیند را با زبان برنامه‌نویسی mathematica انجام داده‌ایم. تحلیل خطا هم از دیدگاه تئوری آنالیز عددی و هم با استفاده از چند مثال نمونه به‌طور شفاف نشان داده شده است، برای این منظور طیفی از نمونه مسائلی را با استفاده از تبدیلات لاپلاس طراحی کردیم که جواب تحلیلی داشته باشند.

extrapolation method for numerical solution of a model for endemic infectious diseases

introductionmany infectious diseases are endemic in a population. in other words they present for several years. suppose that the population size is constant and the population is uniform. in the sir model the population is divided into three disjoint classes which change with time t and let ,  and be the fractions of the population that susceptible, infectious and removed, respectively. this model formulated as the following system of nonlinear volterra integral equation.where ,  and are unknown functions and other constants and functions are known.  the susceptibles are transferred at a rate equal to  times the number of infectives, where  is a constant.  is the nonincreasing probablity function of remaining infectious  units after becoming infectious, with  and  and is dominated by a decaying exponential, such as gamma distributed. since the population size is constant, the birth rate must be equal to the death rate . the death rate is the same for susceptibles, infectives and removed individuals. the fraction  of newborns are immunized so that the flow rate of immunized newborns into the removed class is . the initial susceptible and removed fractions be  and and  be the fraction of  the population that was initially infectious and is still alive and infectious at time .material and methodswe apply the richardson extrapolation method for numerical solution of this model, so that the nonlinear system is solvable by an iterative process with a good accuracy. the algorithm of such systems completely described. this algorithm has a kind of nested structure, which cause we use the lag data in the future times, and it is the interesting section of programing of the algorithm. this algorithm is ready for programing with every program language, which we do this process by mathematica programing software. convergence and accuracy of the method is illustrated by either theoretical and numerical analysis, and some benchmark sample problems. for this aim by using laplase transform, we sketch a spectrum of sample problems. these problems have analytical solution and appropriate for comparison with numerical solutions.results and discussionwe solve some test examples by using present technique to demonstrate the efficiency, high accuracy and the simplicity of the present method. the main advantage of the method is the applicability of  method for a large interval of time, as the algorithm shows. numerical results shows the accuracy of the method for a long time interval.conclusionthe following conclusions were drawn from this research.the proposed algorithm is very suitable for mathematical programing.many cancer problems have such structure and the method is applicable for them.this method has two characteristics, solve a nonlinear problem and use of previous solution in new interval. so the method is applicable for various kind of problems with little additional works../files/site1/files/51/%d8%a8%d8%a7%d8%a8%d8%a7%db%8c%d8%a7%d8%b1.pdf