the economic impact of dementia is already enormous , with total costs greater than  177 billion in europe alone in 2008.1 in 2010 , the total estimated worldwide costs of dementia were 604 billion . about 70% of the costs occurred in western europe and north america.2 the annual cost per person with dementia is estimated at almost  20,000 . the economic impacts of ad are shown in several studies.311 mounting an effective response to the ad epidemic will require international research collaboration and innovative approaches to sharing , integration , and analysis of diverse types of data from multiple sources.12 among the subtypes of dementia , ad is generally considered the most common , accounting for between 60% and 80% of all dementia cases.13 there are various studies predicting the number of patients with ad or dementia in different regions . some of them are based on expert opinions , and some of them are based on mathematical models . a detailed literature review on the global occurrence of dementia is shown in prince et al.14 the literature containing information about the rates of dementia from 1980 to 2009 is studied here as well . the prediction of the number of patients with ad based on mathematical models for the years 2030 and 2050 can also be found in this review . this historical review as well as the prediction contains information about regions all around the world for males and females and for different age groups . another interesting review and prediction for 2030 and 2050 is shown in sosa - ortiz et al.15 the numbers of patients with ad are from all of the world but the patients are not divided by sex . a multistate probabilistic model for all continents has been introduced in brookmeyer et al,16 predicting the number of patients with ad until the year 2050 . the patients were divided into two groups based on the progression of ad . the number of patients with ad has been predicted for each continent separately . in this study , one of the other predictions is based on the method of delphi consensus.17 this method has been used to predict the number of patients with dementia until 2040 based on the opinions of experts without using simulation tools . a comparison of the results of these predictions can be seen in table 1 . our research is focused on predicting the spread of ad in the european union ( eu ) ; thus , we use data from related studies about this region . there are also studies about other regions , such as the us,13,18 people s republic of china,19 or predictions of the spread of ad in developing countries.20 ad is a slowly progressing neurodegenerative brain disease with irreversible brain effects;21 it is the most common cause of dementia22,23 and is responsible for an estimated 60% to 80% of the cases.13 dementia is characterized by a decline in memory , language , problem - solving capability , and other cognitive skills that affect a person s ability to perform everyday activities . ad was first identified more than 100 years ago , but 70 years passed before it was recognized as the most common cause of dementia and a major killer.24 criteria for all - cause dementia and ad dementia are presented in mckhann et al.25 another study26 provides a systematic review of the global literature on the prevalence of dementia ( 19802009 ) and a meta - analysis estimates the prevalence and numbers of those who are affected.27 ad is a neurodegenerative disease characterized by a progressive cognitive decline , leading to the need for complete care within several years after clinical diagnosis . people younger than 65 can also develop the disease , although this is much more rare . the global prevalence of dementia in the 65 + age group was estimated to be 24 million in 2001 , and is predicted to double every 20 years , affecting more than 80 million people worldwide by 2040.28,29 predictions of the population with ad and dementia in europe from different sources are in table 1 . unfortunately , the sources chose different years for the data collection and predictions ; thus , for some years there are no data . models help us to understand dynamics of real - world processes by mimicking actual but simplified forces that are assumed to result in the system s behavior . modeling dynamic systems is essential to our understanding of real - world phenomena.30 as pegels noted,31 the computer modeling process is to the mind what the telescope and the microscope are to the eye . dynamic models try to reflect changes in real or simulated time and take into account that the model components are constantly changing as a result of previous conditions and current influences . system dynamics was introduced by forrester32 as a method for modeling and analyzing the behavior of complex social systems , particularly in an industrial context . it has been used to examine various social , economical , and environmental systems33 in which a holistic view is important and feedback loops are critical for understanding the interrelationships . the approach has attracted attention in recent years since computer software has become more available to help managers with the modeling of systems and decision making . system dynamics represents a set of techniques to describe , observe , and examine the behavior of dynamic systems . the process of model creation and simulation is computer aided , providing a graphical representation of a model and enabling the inspection of the system state at any time . the computer assisted modeling of a system s dynamics provides the ability to run a simulation several times with different values of the variables , allowing further investigation of the influence of the different factors on the system . sterman34 writes :  system dynamics is a perspective and set of conceptual tools that enable us to understand the structure and dynamics of complex systems . system dynamics is also a rigorous modeling method that enables us to build formal computer simulations of complex systems and use them to design more effective policies and organizations . models can show how basic demographic progresses shape the size and composition of a population , revealing the underlying dynamics of growth and attrition , concentration and dispersal , and aging and renewal.35 population dynamics models serve a wide variety of uses , including the explanation of previous population fluctuations as well as projecting future population progression under hypothetical scenarios.36 the following formula has been used for the population modeling . no detailed model of migration was implemented due to the lack of reliable data regarding migration , although both emigration and immigration are considered in the population size . n(t+1)=b(t)n(t)q(t)n(t)=r(t)n(t),where b(t ) and q(t ) are positive functions . in the following , we will denote 
 nx(t )  size of a cohort of the general population at age xmx(t )  size of a cohort of the population with ad at the age xtqx(t )  specific death rate at age x for the total populationqx(t )  specific death rate at age x for the population with adtpx(t )  specific survival rate at age x for the general populationpx(t )  specific survival rate at age x for the population with adrx  prevalence rate at age x. nx(t )  size of a cohort of the general population at age x mx(t )  size of a cohort of the population with ad at the age x tqx(t )  specific death rate at age x for the total population qx(t )  specific death rate at age x for the population with ad tpx(t )  specific survival rate at age x for the general population px(t )  specific survival rate at age x for the population with ad rx  prevalence rate at age x. our theoretical model consists of two systems of recurrence relations . ,n100(t));n1(t+1)=tp0n0(t);nx+1(t+1)=tpxnx(t);n100(t+1)=tp99n99(t).(1 ) the function f computes the births from the total population ; each function nx(t + 1 ) computes the number of people in the cohort x that  survived  from the previous year when 1 year younger . system 2 describes the population with ad : 
 m65(t+1)=r65n65(t);m66(t+1)=p65m65(t)+k66n66(t);mx+1(t+1)=pxmx(t)+kx+1nx+1(t);m100(t+1)=p99m99(t)+k100n100(t).(2 ) the number mx(t ) is the number of people x - years old with ad at time t. the coefficient r65 denotes the prevalence rate of 65-year - olds and is used for finding how many patients with ad are in the population . we have to combine the proportion of surviving patients from the previous year with those newly - developed in the recent year . our estimations are based on data from eurostat ( eu28 , the year 2013).37 but this only includes some of the data required for our model ; therefore , we have to estimate or calculate the missing data . in ad , the prevalence rate is estimated to be from around 1.5%2.5% in the 6569-year bracket , reaching almost 40% in the 9094-year age group , and it increases exponentially.38 the age - standardized prevalence for those who are older than 59 years varied in a narrow band , 5%7% in the majority of regions of the world , with a slightly higher prevalence in latin america ( 8.5% ) , and a significantly lower prevalence in the four sub - saharan african regions ( 2%4% ) . the estimation was that 35.6 million people lived with dementia worldwide in 2010 , but according to an estimation , these numbers are predicted to almost double every 20 years in the future : to 65.7 million in 2030 and 115.4 million in 2050 . in 2010 , 58% of all people suffering from dementia lived in countries with low or middle incomes , with this proportion anticipated to rise to 63% in 2030 and 71% in 2050.14 an exponential function for the prevalence rate depending on age has been created based on the data from:38 
 rx=0.0142e0.1161x . we also need to find an expression for the proportion kx(t ) of new patients with ad in each year - cohort . for the prevalence rate , one has 
 rx = mx(t)nx(t),(3)and for the cohort x of patients with ad one has 
 mx(t)=px1(t1)mx1(t1)+kx(t)nx(t).(4 ) the whole cohort is composed from the patients who survived from the previous year ( and who were 1 year younger ) ; we denote this part by 
 mx(t ) and from new patients who are diagnosed this year . from equation 3 , we get 
 mx(t)=rxnx(t ) . if we put equations 4 and 5 together , we can write 
 rxnx(t)=mx(t)+kx(t)nx(t)(5)from which we get 
 kx(t)=rxmx(t)nx(t).(6 ) firstly , we estimate the evolution of the probability of the death of members of the individual age cohorts ( x = 0 ,  , n ) based on the available data from the eurostat database from 2003 to 2012 . we used the linear regression model with unknown parameters ax and bx , 
 tqx(t)=axt+bx . a simple linear function would not be suitable because after a couple of years the death rate would become negative . we used a linear rational function for modeling the deceleration of the death rate with a zero limit with increasing time . the parameters ax and bx are calculated coefficients and t is the corresponding time of simulation . another challenge is to find the specific death rate ( qx ) for patients with ad . we know the corresponding life expectancy ( the function ex ) of patients with ad . but we need the values of qx , the death rate according to the age x , for our model . let us have a look at some basic formulas from demography : 
 qx=1px , lx+1=lxpx , dx = lxlx+1,lx = lx+1+axdx , tx = tx+1+lx , ex = txlx , where lx stands for the number of people who survive to age x , dx stands for the number of deaths during the year , lx gives the number of person - years lived between ages x , and x+1 , ax is the coefficient of the distribution of deaths during the year , usually for a newborn cohort a0 = 0.08 is used , and for the rest ax = 0.5 , and tx gives the total number of person - years lived beyond age x , by all members of the cohort . we write n for the last year ( the year that none survive ) . then 
 tn = ln , ln = andn , dn = ln , from which we get 
 en = tnln = lnln = anlnlnan . for the previous age cohort , we have 
 en1=tn1ln1=ln+ln1ln1=andn+ln+an1dn1ln1=anln+ln+an1ln1an1lnln1=ln(an+1an1)an1ln1ln1 if we write 
 bn1=an+1an1,we get 
 en1=ln1pn1bn1ln1+an1=pn1bn1+an1and finally we can evaluate pn1 : 
 pn1=en1an1bn1.(7 ) for now on , we will write 
 bn1=ani+1 + 1ani . for the next step , we have 
 en2=tn2ln2=ln+ln1+ln2ln2=andn+ln+an1dn1+ln1+an2dn2ln2=anln+ln+an1(ln1ln)+ln1+an2(ln2ln1)ln2=ln(an+1an1)+ln1(an1 + 1an2)+an2ln2ln2=lnbn1+ln1bn2+an2ln2ln2=ln1pn1bn1+ln2pn2bn2ln2+an2=pn1pn2ln2bn1ln2+pn2bn2+an2=pn1pn2bn1+pn2bn2+an2=pn2(pn1bn1+bn2)+an2 . then , 
 pn2=en2an2bn1pn1+bn2.(8 ) similarly , 
 pn3=en3an3bn1pn1pn2+bn2pn2+bn3.(9 ) now , one can see the pattern from formulas 79 and we can formalize the general formula for the birth rate at the age nx : 
 pnx = enxanxi=1xbnij=1x1pnj(10)for x = 1 ,  , n. and from this survival rate , we can easily compute the required death rate to be qx = 1 px . a population model based on the formulas from  estimation of parameters  section the input to the model is the current size of the eu population based on the eurostat data . we used the simulation software stella to organize the formulas and implement them into formal models that can be run on a computer . stella was chosen as the software for this research because it is a powerful tool for modeling complex systems and their development in time . because of the complexity of dynamic systems , the use of formal models and numbers is essential : they help us dispel the complexity of many real - world processes and force us to be specific . the results of the model are not affected by randomness . the population model and the model for the population with ad consist of basic elements of the stock and flow diagram , such as stocks , flows , converters , and action connectors.33 as we can see in figures 1 and 2 , we use all the elements . the elements named  age x  are stocks , which means the place where the people in a given age are waiting to move to the next state ( age + 1 ) . the only way to enter this stock is through the inflow named  populationbirth .  this inflow calculates the number of incoming ( birth ) people from the sum of the current total population and birth rate . the outflows from the states and from the model are through the outflows named  death x ,  which is calculated from the  agespecificdeathrate . the structure of the model of the population with ad is shown in figure 2 . the main variables in this model are the numbers of patients with ad at each age . the inflow to the first stock of patients at age 65 is given by the relationship between  prevalencerate  and the size of cohort of people of age 65 . other inflows are given by the relationships of the coefficient  k  calculated from equation 6 and the size of the corresponding cohort of people at the given age . the other outflow is the aging of patients to the next age - cohort s stock . the last stock of patients at age 100 has only one outflow ,  death ,  because it is the last stock in this model . in ad , the prevalence rate is estimated to be from around 1.5%2.5% in the 6569-year bracket , reaching almost 40% in the 9094-year age group , and it increases exponentially.38 the age - standardized prevalence for those who are older than 59 years varied in a narrow band , 5%7% in the majority of regions of the world , with a slightly higher prevalence in latin america ( 8.5% ) , and a significantly lower prevalence in the four sub - saharan african regions ( 2%4% ) . the estimation was that 35.6 million people lived with dementia worldwide in 2010 , but according to an estimation , these numbers are predicted to almost double every 20 years in the future : to 65.7 million in 2030 and 115.4 million in 2050 . in 2010 , 58% of all people suffering from dementia lived in countries with low or middle incomes , with this proportion anticipated to rise to 63% in 2030 and 71% in 2050.14 an exponential function for the prevalence rate depending on age has been created based on the data from:38 
 rx=0.0142e0.1161x . we also need to find an expression for the proportion kx(t ) of new patients with ad in each year - cohort . for the prevalence rate , one has 
 rx = mx(t)nx(t),(3)and for the cohort x of patients with ad one has 
 mx(t)=px1(t1)mx1(t1)+kx(t)nx(t).(4 ) the whole cohort is composed from the patients who survived from the previous year ( and who were 1 year younger ) ; we denote this part by 
 mx(t ) and from new patients who are diagnosed this year . from equation 3 , we get 
 mx(t)=rxnx(t ) . if we put equations 4 and 5 together , we can write 
 rxnx(t)=mx(t)+kx(t)nx(t)(5)from which we get 
 kx(t)=rxmx(t)nx(t).(6 ) firstly , we estimate the evolution of the probability of the death of members of the individual age cohorts ( x = 0 ,  , n ) based on the available data from the eurostat database from 2003 to 2012 . we used the linear regression model with unknown parameters ax and bx , 
 tqx(t)=axt+bx . a simple linear function would not be suitable because after a couple of years the death rate would become negative . we used a linear rational function for modeling the deceleration of the death rate with a zero limit with increasing time . the parameters ax and bx another challenge is to find the specific death rate ( qx ) for patients with ad . we know the corresponding life expectancy ( the function ex ) of patients with ad . but we need the values of qx , the death rate according to the age x , for our model . let us have a look at some basic formulas from demography : 
 qx=1px , lx+1=lxpx , dx = lxlx+1,lx = lx+1+axdx , tx = tx+1+lx , ex = txlx , where lx stands for the number of people who survive to age x , dx stands for the number of deaths during the year , lx gives the number of person - years lived between ages x , and x+1 , ax is the coefficient of the distribution of deaths during the year , usually for a newborn cohort a0 = 0.08 is used , and for the rest ax = 0.5 , and tx gives the total number of person - years lived beyond age x , by all members of the cohort . . then 
 tn = ln , ln = andn , dn = ln , from which we get 
 en = tnln = lnln = anlnlnan . for the previous age cohort , we have 
 en1=tn1ln1=ln+ln1ln1=andn+ln+an1dn1ln1=anln+ln+an1ln1an1lnln1=ln(an+1an1)an1ln1ln1 if we write 
 bn1=an+1an1,we get 
 en1=ln1pn1bn1ln1+an1=pn1bn1+an1and finally we can evaluate pn1 : 
 pn1=en1an1bn1.(7 ) for now on , we will write 
 bn1=ani+1 + 1ani . for the next step , we have 
 en2=tn2ln2=ln+ln1+ln2ln2=andn+ln+an1dn1+ln1+an2dn2ln2=anln+ln+an1(ln1ln)+ln1+an2(ln2ln1)ln2=ln(an+1an1)+ln1(an1 + 1an2)+an2ln2ln2=lnbn1+ln1bn2+an2ln2ln2=ln1pn1bn1+ln2pn2bn2ln2+an2=pn1pn2ln2bn1ln2+pn2bn2+an2=pn1pn2bn1+pn2bn2+an2=pn2(pn1bn1+bn2)+an2 . then , 
 pn2=en2an2bn1pn1+bn2.(8 ) similarly , 
 pn3=en3an3bn1pn1pn2+bn2pn2+bn3.(9 ) now , one can see the pattern from formulas 79 and we can formalize the general formula for the birth rate at the age nx : 
 pnx = enxanxi=1xbnij=1x1pnj(10)for x = 1 ,  , n. and from this survival rate , we can easily compute the required death rate to be qx = 1 px . in ad , the prevalence rate is estimated to be from around 1.5%2.5% in the 6569-year bracket , reaching almost 40% in the 9094-year age group , and it increases exponentially.38 the age - standardized prevalence for those who are older than 59 years varied in a narrow band , 5%7% in the majority of regions of the world , with a slightly higher prevalence in latin america ( 8.5% ) , and a significantly lower prevalence in the four sub - saharan african regions ( 2%4% ) . the estimation was that 35.6 million people lived with dementia worldwide in 2010 , but according to an estimation , these numbers are predicted to almost double every 20 years in the future : to 65.7 million in 2030 and 115.4 million in 2050 . in 2010 , 58% of all people suffering from dementia lived in countries with low or middle incomes , with this proportion anticipated to rise to 63% in 2030 and 71% in 2050.14 an exponential function for the prevalence rate depending on age has been created based on the data from:38 
 rx=0.0142e0.1161x . we also need to find an expression for the proportion kx(t ) of new patients with ad in each year - cohort . for the prevalence rate , one has 
 rx = mx(t)nx(t),(3)and for the cohort x of patients with ad one has 
 mx(t)=px1(t1)mx1(t1)+kx(t)nx(t).(4 ) the whole cohort is composed from the patients who survived from the previous year ( and who were 1 year younger ) ; we denote this part by 
 mx(t ) and from new patients who are diagnosed this year . from equation 3 , we get 
 mx(t)=rxnx(t ) . if we put equations 4 and 5 together , we can write 
 rxnx(t)=mx(t)+kx(t)nx(t)(5)from which we get 
 kx(t)=rxmx(t)nx(t).(6 ) firstly , we estimate the evolution of the probability of the death of members of the individual age cohorts ( x = 0 ,  , n ) based on the available data from the eurostat database from 2003 to 2012 . we used the linear regression model with unknown parameters ax and bx , 
 tqx(t)=axt+bx . a simple linear function would not be suitable because after a couple of years the death rate would become negative . we used a linear rational function for modeling the deceleration of the death rate with a zero limit with increasing time . the parameters ax and bx are calculated coefficients and t is the corresponding time of simulation . another challenge is to find the specific death rate ( qx ) for patients with ad . we know the corresponding life expectancy ( the function ex ) of patients with ad . but we need the values of qx , the death rate according to the age x , for our model . let us have a look at some basic formulas from demography : 
 qx=1px , lx+1=lxpx , dx = lxlx+1,lx = lx+1+axdx , tx = tx+1+lx , ex = txlx , where lx stands for the number of people who survive to age x , dx stands for the number of deaths during the year , lx gives the number of person - years lived between ages x , and x+1 , ax is the coefficient of the distribution of deaths during the year , usually for a newborn cohort a0 = 0.08 is used , and for the rest ax = 0.5 , and tx gives the total number of person - years lived beyond age x , by all members of the cohort . we write n for the last year ( the year that none survive ) . then 
 tn = ln , ln = andn , dn = ln , from which we get 
 en = tnln = lnln = anlnlnan . for the previous age cohort , we have 
 en1=tn1ln1=ln+ln1ln1=andn+ln+an1dn1ln1=anln+ln+an1ln1an1lnln1=ln(an+1an1)an1ln1ln1 if we write 
 bn1=an+1an1,we get 
 en1=ln1pn1bn1ln1+an1=pn1bn1+an1and finally we can evaluate pn1 : 
 pn1=en1an1bn1.(7 ) for now on , we will write 
 bn1=ani+1 + 1ani . for the next step , we have 
 en2=tn2ln2=ln+ln1+ln2ln2=andn+ln+an1dn1+ln1+an2dn2ln2=anln+ln+an1(ln1ln)+ln1+an2(ln2ln1)ln2=ln(an+1an1)+ln1(an1 + 1an2)+an2ln2ln2=lnbn1+ln1bn2+an2ln2ln2=ln1pn1bn1+ln2pn2bn2ln2+an2=pn1pn2ln2bn1ln2+pn2bn2+an2=pn1pn2bn1+pn2bn2+an2=pn2(pn1bn1+bn2)+an2 . then , 
 pn2=en2an2bn1pn1+bn2.(8 ) similarly , 
 pn3=en3an3bn1pn1pn2+bn2pn2+bn3.(9 ) now , one can see the pattern from formulas 79 and we can formalize the general formula for the birth rate at the age nx : 
 pnx = enxanxi=1xbnij=1x1pnj(10)for x = 1 ,  , n. and from this survival rate , we can easily compute the required death rate to be qx = 1 px . a population model based on the formulas from  estimation of parameters  section was created using stock and flow diagrams ( figures 1 and 2 ) . the input to the model is the current size of the eu population based on the eurostat data . we used the simulation software stella to organize the formulas and implement them into formal models that can be run on a computer . stella was chosen as the software for this research because it is a powerful tool for modeling complex systems and their development in time . because of the complexity of dynamic systems , the use of formal models and numbers is essential : they help us dispel the complexity of many real - world processes and force us to be specific . the results of the model are not affected by randomness . the population model and the model for the population with ad consist of basic elements of the stock and flow diagram , such as stocks , flows , converters , and action connectors.33 as we can see in figures 1 and 2 , we use all the elements . the elements named  age x  are stocks , which means the place where the people in a given age are waiting to move to the next state ( age + 1 ) . the only way to enter this stock is through the inflow named  populationbirth . this inflow calculates the number of incoming ( birth ) people from the sum of the current total population and birth rate . the outflows from the states and from the model are through the outflows named  death x ,  which is calculated from the  agespecificdeathrate . the structure of the model of the population with ad is shown in figure 2 . the main variables in this model are the numbers of patients with ad at each age . the inflow to the first stock of patients at age 65 is given by the relationship between  prevalencerate  and the size of cohort of people of age 65 . other inflows are given by the relationships of the coefficient  k  calculated from equation 6 and the size of the corresponding cohort of people at the given age . the other outflow is the aging of patients to the next age - cohort s stock . the last stock of patients at age 100 has only one outflow ,  death ,  because it is the last stock in this model . there was only one run of the model because there is no randomness in the system dynamics model of the population growth ; therefore , the results are always the same . in the result of the simulation , the prediction of the total population is similar to the prediction from eurostat ( table 2 and figure 3 ) . moreover , our model provides information about the number of people with ad and also the characteristics of the population with ad . most of other predictions are until the year 2050 ( table 1 ) . the number of patients with ad in those predictions rises . we have simulated 70 years to the future and it can be seen that number of patients with ad starts to slowly decrease from the year 2060 because of changes in eu population size and its structure . our population structure has been verified by data from eurostat , which provides population prediction until the year 2080 . simulating more than 70 years would make our model less accurate because there could be major changes in population structure . the weaknesses of the model is that population structure can change rapidly and may not correspond with eurostat prediction . our model is easily alterable and new data about population or disease characteristic can be simply implemented . the data from eurostat are indicated by the dotted line for better comparison with our results in figure 3 . rapid growth can be seen at the beginning of the simulation ; it slows down after ~30 years , followed by a slow decrease . if we look only at the elderly population ( people over 80 years old , figure 4 ) , we can see that the model also fits the prediction from the organisation for economic co - operation and development ( oecd),39 especially at the beginning of the simulation . our prediction is quite optimistic compared to the oecd in the following years : the proportion of the elderly population will decrease , and therefore the population will get younger . moreover , we obtained a prediction for the population of patients with ad ( figure 5 ) . it can be seen that the population of patients with ad will increase at first , but at later times will decrease . it looks like the population growth of patients with ad follows the growth of the general population with a delay . the right y - axis shows the numbers of the ad population , and the left y - axis shows the numbers of the general population . the maximal value of the population of patients with ad is ~15 million around the year 2060 . this proportion will increase from 1.5% to almost 3% in 50 years and then it will decrease to 2.5% after another 20 years . the behavior of this proportion is similar to the behavior of the population of patients with ad . after a closer look at the chosen characteristics , such as median or average , a population that has been in the process of aging starts to become younger after a couple of decades ( figure 7 ) . this will not apply to patients with ad . a system dynamics model simulating the population growth of the eu and the number of patients with ad has been introduced . based on the characteristics of ad , it was possible to simulate the spread of ad till the year 2080 and explore different population characteristics . the modeling results of the population growth equal the predictions of eurostat.37 we were able to create a model of population growth similar to the eurostat prediction where it is , moreover , possible to observe the number of people with ad and different population characteristics of those patients . the data from eurostat do not provide information about the number of people with ad at all . also , compared to other models and predictions , the population of the eu was 508 million and the number of patients with ad was 7.5 million in the year 2013 . based on eurostat and our predictions , in the year 2040 the population of the eu will be 524 million . based on our predictions , the number of patients with ad in 2040 will be 13.1 million ; the eu population in 2080 will be 520 million with 13.7 million patients with ad . a system dynamics modeling approach has been used for the prediction of the number of patients with ad in the eu population till the year 2080 . these results can be used to determine the economic burden of the treatment of these patients . with different input data , such a simulation can also be used for different regions as well as for different noncontagious disease predictions .
introductionalzheimer s disease ( ad ) is a slowly progressing neurodegenerative brain disease with irreversible brain effects ; it is the most common cause of dementia . with increasing age , the probability of suffering from ad increases . 
 in this research , population growth of the european union ( eu ) until the year 2080 and the number of patients with ad are modeled.aimthe aim of this research is to predict the spread of ad in the eu population until year 2080 using a computer simulation.methodsfor the simulation of the eu population and the occurrence of ad in this population , a system dynamics modeling approach has been used . 
 system dynamics is a useful and effective method for the investigation of complex social systems . over the past decades 
 , its applicability has been demonstrated in a wide variety of applications . in this research , this method has been used to investigate the growth of the eu population and predict the number of patients with ad . 
 the model has been calibrated on the population prediction data created by eurostat.resultsbased on data from eurostat , the eu population until year 2080 has been modeled . in 2013 , the population of the eu was 508 million and the number of patients with ad was 7.5 million . based on the prediction , in 2040 , the population of the eu will be 524 million and the number of patients with ad will be 13.1 million . by the year 2080 , the eu population will be 520 million and the number of patients with ad will be 13.7 million.conclusionsystem dynamics modeling approach has been used for the prediction of the number of patients with ad in the eu population till the year 2080 . 
 these results can be used to determine the economic burden of the treatment of these patients . with different input data , 
 the simulation can be used also for the different regions as well as for different noncontagious disease predictions .