DDNS Discrete Dynamics in Nature and Society 1607-887X 1026-0226 Hindawi Publishing Corporation 10.1155/2015/842792 842792 Research Article Mathematical Modelling, Simulation, and Optimal Control of the 2014 Ebola Outbreak in West Africa http://orcid.org/0000-0002-9643-5054 Rachah Amira 1 http://orcid.org/0000-0001-8641-2505 Torres Delfim F. M. 2 Yuan Sanling 1 Mathématiques pour l’Industrie et la Physique Institut de Mathématiques de Toulouse Université Paul Sabatier 31062 Toulouse Cedex 9 France ups-tlse.fr 2 Center for Research and Development in Mathematics and Applications (CIDMA) Department of Mathematics University of Aveiro 3810-193 Aveiro Portugal ua.pt 2015 1452015 2015 26 12 2014 28 02 2015 1452015 2015 Copyright © 2015 Amira Rachah and Delfim F. M. Torres. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The Ebola virus is currently one of the most virulent pathogens for humans. The latest major outbreak occurred in Guinea, Sierra Leone, and Liberia in 2014. With the aim of understanding the spread of infection in the affected countries, it is crucial to modelize the virus and simulate it. In this paper, we begin by studying a simple mathematical model that describes the 2014 Ebola outbreak in Liberia. Then, we use numerical simulations and available data provided by the World Health Organization to validate the obtained mathematical model. Moreover, we develop a new mathematical model including vaccination of individuals. We discuss different cases of vaccination in order to predict the effect of vaccination on the infected individuals over time. Finally, we apply optimal control to study the impact of vaccination on the spread of the Ebola virus. The optimal control problem is solved numerically by using a direct multiple shooting method.

1. Introduction

Ebola is a lethal virus for humans. It was previously confined to Central Africa but recently was also identified in West Africa . As of October 8, 2014, the World Health Organization (WHO) reported 4656 cases of Ebola virus deaths, with most cases occurring in Liberia . The extremely rapid increase of the disease and the high mortality rate make this virus a major problem for public health. Typically, patients present fever, malaise, abdominal, headache, and asthenia. One week after the onset of symptoms, a rash often appears followed by haemorrhagic complications, leading to death after an average of 10 days in 50%–90% of infections . Ebola is transmitted through direct contact with blood, bodily secretions and tissues of infected ill or dead humans and nonhuman primates .

The main aim of this work is to understand the dynamics of the Liberian population infected by Ebola virus in 2014, by using an appropriate mathematical model. More precisely, we begin by considering a system of ordinary differential equations to describe the 2014 Ebola outbreak, which is nothing else than an epidemic SIR model, that is, a model based on the division of the population into three groups: the susceptible, the infected, and the recovered . We simulate the model obtained by parameters estimated on November 4, 2014, by Kaurov . Our objectives are to better understand the outbreak and to predict the effect of vaccination on the infected individuals over time.

In order to deal with the epidemic of Ebola, governments have decided to implement tough measures. For example, some have applied quarantine procedures while others have opted for mass vaccination plans as a precaution to the epidemic . In this context, we consider an optimal control problem to study the effect of vaccination on the spread of virus.

The text is organized as follows. In Section 2 we introduce a basic mathematical model to describe the dynamics of the Ebola virus during the 2014 outbreak in Liberia. After the mathematical modelling, we use the obtained model to simulate it in Section 3 with the parameters estimated from recent statistical data based on the WHO report of the 2014 Ebola outbreak . Following , in Section 4 we study different cases of vaccination of individuals by adding a vaccination term to the model. Section 5 presents an optimal control problem, which we use to study the impact of a vaccination campaign on the spread of the virus. We end with Section 6 of conclusions and future work, where the results are summarized and some research perspectives presented.

2. The Basic Mathematical Model

The objective of this section is to describe mathematically the dynamics of the population infected by the Ebola virus. The dynamics is described by a system of differential equations. This system is based on the common SIR (susceptible-infectious-recovery) epidemic model, where the population is divided into three groups: the susceptible group, denoted by S(t), the infected group, denoted by I(t), and the recovered group, denoted by R(t). The total population, assumed constant during the short period of time under study, is given by N=S(t)+I(t)+R(t).

To create the equations that describe the population of each group along time, let us start with the fact that the population of the susceptible group will be reduced as the infected come into contact with them with a rate of infection β. This means that the change in the population of susceptible is equal to the negative product of β with S(t) and I(t): (1)dStdt=-βStItN.Now, let us create the equation that describes the infected group over time, knowing that the population of this group changes in two ways:

people leave the susceptible group and join the infected group, thus adding to the total population of infected a term βS(t)I(t);

people leave the infected group and join the recovered group, reducing the infected population by -μI(t).

This is written as (2)dItdt=βStItN-μIt.Finally, the equation that describes the recovery population is based on the individuals recovered from the virus at rate μ. This means that the recovery group is increased by μ multiplied by I(t):(3)dRtdt=μIt.Figure 1 shows the relationship between the three variables of our SIR model. As shown in the next section, this simple model describes well the 2014 Ebola outbreak in Liberia for suitable chosen parameters β and μ and appropriate initial conditions S(0), I(0), and R(0).

The scheme of the susceptible-infected-recovered model (1)–(3).

3. Numerical Simulations of the Basic Model and Discussion

With the aim of approximating the 2014 Ebola outbreak in Liberia, we now simulate the complete set of (1)–(3) that describe the SIR model: (4)dStdt=-βStItN,dItdt=βStItN-μIt,dRtdt=μIt,where β is the rate of infection and μ the rate of recovery. The data for the outbreak in Liberia is based on the cumulative numbers of total reported cases (confirmed and probable) taken from the official data reported by the Ministry of Health of Liberia up to October 26, 2014, as provided by the World Health Organisation (WHO)  (in particular, the Ebola virus disease cases reported each week from Liberia are given in Figure 2 of ). The model fits well the reported data of confirmed (means infected) cases in Liberia. Note that the exact population size N does not need to be known to estimate the model parameters as long as the number of cases is small compared to the total population size . To estimate the parameters of the model, we adapted the initialisation of S and I with the reported data of  by fitting the actual data of confirmed cases in Liberia. The data that we are studying is the number of individuals with confirmed infection between July 14 and October 12, 2014. The result of fitting is shown in Figure 2. The comparison between the curve of infected obtained by our simulation and the reported data of confirmed cases by WHO shows that the mathematical model (4) fits well the real data by using β=0.000318 as the rate of infection and μ=0.0175 as the rate of recovery. The initial numbers of susceptible, infected, and recovered groups, are given by S(0)=460, I(0)=12, and R(0)=0, respectively. The choice of S(0) is in agreement with the data shown in Figure 2 of , as the biggest value in the histogram is 460. To interpret the simulation results, we describe the evolution of each group of individuals over time and the link between them. Figure 3 shows that the susceptible group begins to plummet and, simultaneously, the number of infected begins to rise. This is due to how infectious the Ebola virus is.

Graph of infected I(t) obtained from (4) with β=0.000318, μ=0.0175, S(0)=460, I(0)=12, and R(0)=0 versus the real data of confirmed cases for the 2014 Ebola outbreak occurred in Liberia from July 14 (day t=0) and October 12, 2014 (day t=90).

The evolution of susceptible (S(t)), infected (I(t)), and recovered (R(t)) for the period July 14 (t=0) and October 12, 2014 (t=90), given by the Ebola Liberia model (4) with β=0.000318, μ=0.0175, S(0)=460, I(0)=12, and R(0)=0. The real data of confirmed cases for the 2014 Ebola outbreak occurred in Liberia, taken from , is shown with circles.

When we compare our results with those of Astacio et al. , when they modeled the 1976 outbreak of Ebola in Yambuku, a small village in Mongala District in northern Democratic Republic of Congo (previously Zaire), and the 1995 outbreak of Ebola in Kikwit, in the southwestern part of the Democratic Republic of Congo, we note that the evolution of each group obtained in our simulations is reasonable. The main difference between the 2014 outbreak of Ebola in Liberia and the previous outbreaks of Ebola in Yambuku and Kikwit studied in  is at the end of the evolution of infected: while in the previous 1976 and 1995 outbreaks of Ebola in Yambuku and Kikwit the infected group goes to zero, now in the 2014 outbreak of Ebola in Liberia the infected group remains important with a high number of infected individuals by the virus, which describes the reality of what is currently happening in Liberia. One can conclude that the 2014 Ebola outbreak of Liberia is much more dangerous than the previous Ebola outbreaks.

4. Modelling and Simulation of Ebola Outbreak with Vaccination

In this section, we simulate the model by using parameters estimated on November 4, 2014, by Kaurov , who studied statistically the recent data of the World Health Organisation (WHO) . In his statistical study, Kaurov studied the outbreak by modeling it with Wolfram’s Mathematica language. Here we use his SIR model for modelling and time-discrete equations for resolution of the set of equations through Matlab.

4.1. The Basic Model with Kaurov’s Estimation of Parameters

Let us recall the complete set of equations considered by  (5)dStdt=-βStIt,dItdt=βStIt-μIt,dRtdt=μIt,where, as before, β is the rate of infection and μ is the rate of recovery. The parameters obtained by Kaurov, knowing that 95% of population is susceptible and 5% of population is infected, are β=0.2 and μ=0.1 . In agreement, the initial susceptible, infected, and recovered populations are given, respectively, by S(0)=0.95, I(0)=0.05, and R(0)=0 (note that in contrast with Sections 2 and 3, the variables are now dimensionless).

Our numerical simulations were done with the ODE solver of Matlab. The results are shown in Figure 4. To interpret the results, we describe the evolution of each group over time and the links between them. Figure 4 shows that the susceptible group begins to plummet due to how infectious the virus is and, simultaneously, the infected group’s number begins to rise.

Solution to the 2014 outbreak Ebola model (5) with β=0.2, μ=0.1, S(0)=0.95, I(0)=0.05, and R(0)=0.

An important dimensionless quantity that plays a key role in any SIR model is the well-known basic reproduction number R0 . This number tells us how fast the disease will spread at the epidemic. For model (5), the basic reproduction number is simply given by β/μ. Let us now discuss the comparison between the different values of R0 obtained recently. On September 2, 2014, Althaus studied the case of the 2014 Ebola outbreak occurred in Guinea, in Sierra Leone, and in Liberia . In his study, Althaus used the WHO data occurred in the three countries. He found different values for the basic reproduction number, depending on the country: in Guinea R0=1.51, in Sierra Leone the value of R0 is given by 2.53, while in Liberia R0 is equal to 1.59. Let us now compare between these values and the value for the basic reproduction number studied by Kaurov, knowing that his study, which was done on November 4, 2014, is more recent than the one of Althaus, done on September 2, 2014. As the parameters estimated statistically by Kaurov are for the three countries, our comparison will be between R0,A, the mean of the three basic reproduction numbers found by Althaus, and R0,K, the number found more recently by Kaurov. The value of R0,A is given by 1.876, where the more recent R0,K is given by 2. We claim that the basic reproduction number increased between September and November 2014.

4.2. Dynamics in Case of Vaccination

Infectious diseases have tremendous influence on human life. In last decades, controlling infectious diseases has been an increasingly complex issue. A strategy to control infectious diseases is through vaccination. Now our idea is to study the effect of vaccination in practical Ebola situations. The case of the French nurse cured of Ebola with the help of an experimental vacine is a proof of the possibility of treatment by vaccinating infected individuals .

In this section, we improve the SIR model described in Section 4.1. The system of equations that describe the susceptible-infectious-recovery (SIR) model with vaccination is (6)dStdt=-βStIt-υSt,dItdt=βStIt-μIt,dRtdt=μIt+υSt,where β is the rate of infection, μ is the rate of recovery, and υ is the percentage of individuals vaccinated every day . Figure 5 shows the relationship between the variables of the SIR model with vaccination (6).

The scheme of the susceptible-infected-recovered model with vaccination (6).

4.3. Numerical Simulation in Case of Vaccination

We simulate the SIR model with vaccination (6) in order to predict the evolution of every group of individuals in case of vaccination. Our study is based in the test of different rates of vaccinations and their effect on the curve of every group. The results are shown in Figures 6 and 7. These figures present different rates of vaccination and their effect on the peak of the curve of infected individuals along time. The time-dependent curve of infected individuals shows that when the percentage of vaccinated individuals increases, the peak of the curve of infected individuals is less important and the period of infection (the corresponding number of days) is shorter while the number of recovered gets more important. In fact, in case of infection without vaccination, the curve of the infected group goes to zero in about 90 days (see Figure 6), where in the presence of vaccination the curve of the infected group goes to zero in about 45 days (see Figure 7). This shows the efficiency of using vaccination in controlling Ebola. Table 1 presents the maximum percentage of infected individuals, corresponding to different rates of vaccinations, and the respective number of days for the peak of the infected.

Maximum percentage of infected individuals I corresponding to different rates υ of vaccination and respective number of days for the peak of I.

Rate of vaccination Maximum percentage of infected Days to the peak of I(t)
0 17% 26.52
0.005 15% 23.47
0.01 13% 22.07
0.02 10% 17.09
0.03 9% 13.68
0.06 7% 9.28

A solution to the SIR model (6) with β=0.2, μ=0.1, S(0)=0.95, I(0)=0.05, R(0)=0, and different rates of vaccination: υ=0 (a),   υ=0.005 (b), and υ=0.01 (c).

A solution to the SIR model (6) with β=0.2, μ=0.1, S(0)=0.95, I(0)=0.05, R(0)=0, and different rates of vaccination: υ=0.02 (a), υ=0.03 (b), and υ=0.06 (c).

5. The Optimal Control Problem

Optimal control has been recently used with success in a series of epidemiological problems: see [17, 2527] and references therein. In this section, we present an optimal control problem by introducing into model (5) a control u(t) representing the vaccination rate at time t. The control u(t) is the fraction of susceptible individuals being vaccinated per unit of time. We assume that all susceptible individuals that are vaccinated are transferred directly to the removed class. Then, the mathematical model with control is given by the following system of nonlinear differential equations: (7)dS(t)dt=-βS(t)I(t)-u(t)S(t),dI(t)dt=βS(t)I(t)-μI(t),dR(t)dt=μI(t)+u(t)S(t),where S(0)0, I(0)0, and R(0)0 are given. Our goal is to reduce the infected individuals and the cost of vaccination. Precisely, our optimal control problem consists of minimizing the objective functional (8)Ju=0tfIt+A2u2tdt,subject to 0u(t)0.9 for t[0,tf], where the parameters A0 and tf denote, respectively, the weight on cost and the duration of the vaccination program. We present the numerical solution to this optimal control problem, comparing the results with the numerical simulation results obtained in Section 4 via Matlab. As in Section 4, the rate of infection is given by β=0.2, the recovered rate is given by μ=0.1, and we use S(0)=0.95, I(0)=0.05, and R(0)=0 for the initial number of susceptible, infected, and recovered populations, respectively (i.e., at the beginning 95% of population is susceptible and 5% is already infected with Ebola). Note that when u(t)0 model (7) reduces to (5); and when u(t)υ model (7) reduces to (6).

The numerical simulation of the system without control (u(t)0) was done with Matlab in Section 4.1. For our numerical solution of the optimal control problem, we have used the ACADO solver , which is based on a multiple shooting method, including automatic differentiation and based ultimately on the semidirect multiple shooting algorithm of Bock and Pitt . ACADO is a self-contained public domain software environment written in C++ for automatic control and dynamic optimization. Figures 8 and 9 show the significant difference in the number of infected individuals with and without control. In fact, in case of optimal control, the function I decreases greatly: the maximum value on the infected curve I under optimal control is 5.6% against 17.9% without any control (see Figure 9).

The function I(t) of infected individuals with and without control.

Comparison between the peaks of infected individuals I(t) in case of optimal control versus without control.

In Figure 10, we see that the number of susceptible S decreases more rapidly during the vaccination campaign. It reaches 11% at the end of the campaign against 19% in the absence of optimal control.

Susceptible individuals S(t) in case of optimal control versus without control.

Figure 11 shows that the number of recovered individuals increases rapidly. The number R(tf) at the end of the optimal control vaccination period is 88.7% instead of 80.5% without control. Figure 12 gives a representation of the optimal control u(t).

Recovered individuals R(t) in case of optimal control versus without control.

The optimal control u for problem (7)-(8).

6. Conclusion

We investigated a mathematical model that provides a good description of the 2014 Ebola outbreak in Liberia: the model fits well the data of confirmed cases in Liberia provided by the World Health Organization. The evolution of infected individuals shows that the infected population remains important, describing the actual evolution of the virus in Liberia. This evolution is different than Ebola’s outbreak in the last decades. Our study of the SIR model with vaccination shows that vaccination is a very efficient factor in reducing the number of infected individuals in a short period of time and increasing the number of recovered individuals. We also discussed an optimal control problem related to the impact of a vaccination strategy on the spread of the 2014 outbreak of Ebola in West Africa. The numerical simulation of both systems, with and without control, shows that an optimal control strategy greatly helps to reduce the number of infected and susceptible individuals and increases the number of recovered individuals.

As future work, we plan to include in our study other factors. In particular, we intend to include in the mathematical model the treatment of infected individuals with a quarantine procedure. Another interesting line of research is to investigate the effect of impulsive vaccination.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research was initiated while Rachah was visiting the Department of Mathematics of University of Aveiro, Portugal. The hospitality of the host institution and the financial support provided by the Institut de Mathématiques de Toulouse, France, are here gratefully acknowledged. Torres was supported by funds through the Portuguese Foundation for Science and Technology (FCT), within CIDMA project UID/MAT/04106/2013 and OCHERA project PTDC/EEI-AUT/1450/2012, cofinanced by FEDER under POFC-QREN with COMPETE reference FCOMP-01-0124-FEDER-028894.

Barry M. Traoré F. A. Sako F. B. Kpamy D. O. Bah E. I. Poncin M. Keita S. Cisse M. Touré A. Ebola outbreak in Conakry, Guinea: epidemiological, clinical, and outcome features Médecine et Maladies Infectieuses 2014 44 11-12 491 494 10.1016/j.medmal.2014.09.009 Lewnard J. A. Ndeffo Mbah M. L. Alfaro-Murillo J. A. Altice F. L. Bawo L. Nyenswah T. G. Galvani A. P. Dynamics and control of Ebola virus transmission in Montserrado, Liberia: a mathematical modelling analysis The Lancet Infectious Diseases 2014 14 12 1189 1195 10.1016/s1473-3099(14)70995-8 WHO Report of an International Study Team. Ebola haemorrhagic fever in Sudan 1976 Bulletin of the World Health Organization 1978 56 2 247 270 Report of an International Commission Ebola haemorrhagic fever in Zaire, 1976 Bulletin of the World Health Organization 1978 56 2 271 293 Okware S. I. Omaswa F. G. Zaramba S. Opio A. Lutwama J. J. Kamugisha J. Rwaguma E. B. Kagwa P. Lamunu M. An outbreak of Ebola in Uganda Tropical Medicine and International Health 2002 7 12 1068 1075 10.1046/j.1365-3156.2002.00944.x 2-s2.0-0036885737 Legrand J. Grais R. F. Boelle P. Y. Valleron A. J. Flahault A. Understanding the dynamics of Ebola epidemics Epidemiology and Infection 2007 135 4 610 621 2-s2.0-34247475255 10.1017/S0950268806007217 Dowell S. F. Mukunu R. Ksiazek T. G. Khan A. S. Rollin P. E. Peters C. J. Transmission of Ebola hemorrhagic fever: a study of risk factors in family members, Kikwit, Democratic Republic of the Congo, 1995. Commission de Lutte contre les Epidémies a Kikwit Journal of Infectious Diseases 1999 179 supplement 1 S87 S91 10.1086/514284 2-s2.0-0033061837 Borio L. Inglesby T. Peters C. J. Schmaljohn A. L. Hughes J. M. Jahrling P. B. Ksiazek T. Johnson K. M. Meyerhoff A. O'Toole T. Ascher M. S. Bartlett J. Breman J. G. Eitzen E. M. Jr. Hamburg M. Hauer J. Henderson D. A. Johnson R. T. Kwik G. Layton M. Lillibridge S. Nabel G. J. Osterholm M. T. Perl T. M. Russell P. Tonat K. Hemorrhagic fever viruses as biological weapons: Medical and public health management Journal of the American Medical Association 2002 287 18 2391 2405 10.1001/jama.287.18.2391 2-s2.0-0037042244 Peters C. J. LeDuc J. W. An introduction to Ebola: the virus and the disease Journal of Infectious Diseases 1999 179 supplement 1 2-s2.0-0033054445 Chowell G. Hayman J. M. Bettencourt L. M. A. Castillo-Chavez C. Mathematical and Statistical Estimation Approaches in Epidemiology 2009 Dordrecht, The Netherlands Springer Zeng D. Chen H. Castillo-Chavez C. Lober W. B. Thurmond M. Infectious Disease Informatics and Biosurveillance 2011 27 New York, NY, USA Springer Integrated Series in Information Systems Chowell G. Hengartner N. W. Castillo-Chavez C. Fenimore P. W. Hyman J. M. The basic reproductive number of Ebola and the effects of public health measures: the cases of Congo and Uganda Journal of Theoretical Biology 2004 229 1 119 126 10.1016/j.jtbi.2004.03.006 MR2082775 2-s2.0-2642573296 Kaurov V. Modeling a pandemic like ebola with the Wolfram language Technical Communication & Strategy 2014 http://blog.wolfram.com/2014/11/04/modeling-a-pandemic-like-ebola-with-the-wolfram-language/ Farrington C. P. Kanaan M. N. Gay N. J. Branching process models for surveillance of infectious diseases controlled by mass vaccination Biostatistics 2003 4 2 279 295 10.1093/biostatistics/4.2.279 2-s2.0-0142089895 Forum on Medical and Public Health Preparedness for Catastrophic Events The 2009 H1N1 Influenza Vaccination Campaign: Summary of a Workshop Series 2010 Washington, DC, USA The National Academies Press Robert W. Tindle Vaccines for Human Papillomavirus Infection and Disease Medical Intelligence Unit 1999 Austin, Tex, USA Medical Intelligence Unit, R. G. Landes Company Rodrigues H. S. Monteiro M. T. T. Torres D. F. M. Vaccination models and optimal control strategies to dengue Mathematical Biosciences 2014 247 1 12 10.1016/j.mbs.2013.10.006 MR3148910 2-s2.0-84892390417 WHO Ebola Response Roadmap—Situation Report Update 2014 World Health Organization http://www.who.int/csr/disease/ebola/situation-reports/en World Health Organisation (WHO) Ebola Response Roadmap Situation Report 2014 Geneva, Switzerland World Health Organisation http://apps.who.int/iris/bitstream/10665/137376/1/roadmapsitrep_29Oct2014_eng.pdf?ua=1 Althaus C. L. Estimating the reproduction number of Ebola virus (EBOV) during the 2014 outbreak in West Africa PLoS Currents Outbreaks 2014 10.1371/currents.outbreaks.91afb5e0f279e7f29e7056095255b288 Astacio J. Briere D. Guilléon M. Martinez J. Rodriguez F. Valenzuela-Campos N. Mathematical models to study the outbreaks of Ebola Report 1996 BU-1365-M Lewis M. A. Dietel M. Scriba P. Raff W. K. Biology and Epidemiology of Hormone Replacement Therapy 2006 Berlin, Germany Springer Valleron A. J. Schwartz D. Goldberg M. Salamon R. Collectif Lépidémiologie Humaine, Conditions de Son Développement en France, et Rôle des Mathématiques 2006 462 Institut de France Académie des Sciences Yon M. J. Modélisation de la propagation d'un virus 2010 INSA Rouen Rodrigues H. S. Monteiro M. T. T. Torres D. F. M. Dynamics of dengue epidemics when using optimal control Mathematical and Computer Modelling 2010 52 9-10 1667 1673 2-s2.0-77956012968 MR2719552 10.1016/j.mcm.2010.06.034 Rodrigues P. Silva C. J. Torres D. F. M. Cost-effectiveness analysis of optimal control measures for tuberculosis Bulletin of Mathematical Biology 2014 76 10 2627 2645 10.1007/s11538-014-0028-6 MR3266821 Silva C. J. Torres D. F. M. Optimal control for a tuberculosis model with reinfection and post-exposure interventions Mathematical Biosciences 2013 244 2 154 164 10.1016/j.mbs.2013.05.005 MR3101449 2-s2.0-84881227998 Ariens D. Houska B. Ferreau H. J. ACADO Toolkit User's Manual, Toolkit for Automatic Control and Dynamic Optimization 2010 http://www.acadotoolkit.org Bock H. G. Pitt K. J. A multiple shooting algorithm for direct solution of optimal control problems Proceedings of the 9th IFAC World Congress 1984 Budapest, Hungary Pergamon Press 243 247