﻿<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD with MathML3 v1.2 20190208//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd">
<article
    xmlns:mml="http://www.w3.org/1998/Math/MathML"
    xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="article">
  <front>
    <journal-meta>
      <journal-id journal-id-type="publisher-id">JML</journal-id>
      <journal-title-group>
        <journal-title>Journal of Mathematics Letters</journal-title>
      </journal-title-group>
      <issn pub-type="epub"></issn>
      <issn pub-type="ppub"></issn>
      <publisher>
        <publisher-name>Science Publications</publisher-name>
      </publisher>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.31586/jml.2024.846</article-id>
      <article-id pub-id-type="publisher-id">JML-846</article-id>
      <article-categories>
        <subj-group subj-group-type="heading">
          <subject>Article</subject>
        </subj-group>
      </article-categories>
      <title-group>
        <article-title>
          
        </article-title>
      </title-group>
      <contrib-group>
<contrib contrib-type="author">
<name>
<surname>Prognosis</surname>
<given-names>Modelling Population Growth</given-names>
</name>
<xref rid="af1" ref-type="aff">1</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Akaligwo</surname>
<given-names>Emmanuel</given-names>
</name>
<xref rid="af1" ref-type="aff">1</xref>
<xref rid="cr1" ref-type="corresp">*</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Aharanwa</surname>
<given-names>Boniface</given-names>
</name>
<xref rid="af2" ref-type="aff">2</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Aderotimi</surname>
<given-names>Joshua</given-names>
</name>
<xref rid="af2" ref-type="aff">2</xref>
</contrib>
      </contrib-group>
<aff id="af1"><label>1</label> Department of Mathematics, Federal University Lokoja, Kogi State, Nigeria</aff>
<aff id="af2"><label>2</label> Department of Mathematics/Statistics, Imo State Polytechnic Omuma, Imo State, Nigeria</aff>
<author-notes>
<corresp id="c1">
<label>*</label>Corresponding author at: Department of Mathematics, Federal University Lokoja, Kogi State, Nigeria
</corresp>
</author-notes>
      <pub-date pub-type="epub">
        <day>19</day>
        <month>01</month>
        <year>2024</year>
      </pub-date>
      <volume>2</volume>
      <issue>1</issue>
      <history>
        <date date-type="received">
          <day>09</day>
          <month>12</month>
          <year>2023</year>
        </date>
        <date date-type="rev-recd">
          <day>10</day>
          <month>01</month>
          <year>2024</year>
        </date>
        <date date-type="accepted">
          <day>18</day>
          <month>01</month>
          <year>2024</year>
        </date>
        <date date-type="pub">
          <day>19</day>
          <month>01</month>
          <year>2024</year>
        </date>
      </history>
      <permissions>
        <copyright-statement>&#xa9; Copyright 2024 by authors and Trend Research Publishing Inc. </copyright-statement>
        <copyright-year>2024</copyright-year>
        <license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/4.0/">
          <license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p>
        </license>
      </permissions>
      <abstract>
        Logistic growth model and its variants have been adjudged to be the most appropriate model for forecasting human population. However, in this article, we estimated the carrying capacity of Abuja using the logistic model. Then, we presented the parameters used to ascertain that the logistic model has the best fit in modelling population growth of Abuja over time. Meanwhile, a population growth sensitivity analysis is presented for the year 1962 to 2200.The result shows that by the year 2050, Abuja population growth rate will be out of control, if nothing substantial is implemented. Similarly, from the year 2150, the results show that stability will return again. Furthermore, the result of the error analysis conducted on the logistic model shows that Abuja has a growing population and that logistic growth model with MAPE and RMSE values of 0.98% and 7,817.07 respectively is the most accurate. The study concludes that logistic growth model with R&#x02212;squared value of 0.776 has the best fit for population growth projection of Abuja. With approximate growth rate at 9.3% per annum, the projected population of Abuja will hit 30,220,701 million by the year 2039 all things being equal. Therefore, we recommend that the government should invest in massive agricultural reforms to accommodate the growing population, expand Abuja by developing its suburbs, and engage in massive reorientation of the populace on the dangers of uncontrolled births and the education of the girl child.
      </abstract>
      <kwd-group>
        <kwd-group><kwd>Census; Population Growth; Carrying Capacity; Population Size</kwd>
</kwd-group>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec1">
<title>Introduction</title><p>Nigeria is one of the fastest growing countries in the world, with an estimated population of over 180 million and an annual population growth rate of 2.9%. Nigeria is the most populous black nation in Sub-Sahara Africa and the tenth most populous in the world according to [
<xref ref-type="bibr" rid="R1">1</xref>,<xref ref-type="bibr" rid="R2">2</xref>,<xref ref-type="bibr" rid="R3">3</xref>]. This sounds interesting, but the question is, do the government have appropriate data and policy direction for effective monitoring and control of the population growth rate in such a way that the growing population maintains stability? No doubt it is the duty of government to monitor and manage population and demography related issues. In 1988, the armed forces ruling council headed by General Ibrahim Babangida, approved a policy document titled &#x26;#x0201c;Nigeria Policy on Population for Unity, Progress and Self-reliance&#x26;#x0201d; [
<xref ref-type="bibr" rid="R4">4</xref>]. This policy was a proof of government seriousness and concern for family planning as part of the overall socio-economic development of the country. In 2004, another policy of government was introduced by President Olusegun Obasanjo Administration. This document was called&#x26;#x0201d; Nigeria Policy on Population for Sustainable Development&#x26;#x0201d; [
<xref ref-type="bibr" rid="R5">5</xref>]. The policy recognizes that population factors, social and economic development, and environmental issues are irrevocably inter-related and are therefore critical to the achievement of sustainable development in Nigeria. Since then, several administrations have had one birth control measure or the other, but all of these policies have failed to categorically state the result of the population sensitivity analysis which informed such hastily decision seven though it has good intentions. Presently, there is no population sensitivity analysis data across the nation and in particular Abuja. In an attempt to fill the gap, [
<xref ref-type="bibr" rid="R6">6</xref>] used a mathematical approach and concluded that by 2040, the population growth rate of Nigeria will be out of control stating that everything should be done by the government both in policy and enforcement to ensure that the population growth is been checked early enough and controlled. Conversely, in this article, a logistic population growth model that incorporates a theoretical carrying capacity is investigated. The stability of the growing population is analyzed using the population growth sensitivity analysis. The expected time of instability is presented and recommendations given. </p>
<p></p>
<p><bold>Study Area:</bold> Abuja also known as center of unity is the capital of Nigeria is one of the fastest growing cities in the world with a projected population of over 3 million in 2023 with an annual growth rate of 9.3% [
<xref ref-type="bibr" rid="R7">7</xref>]. Geographically, Abuja is located on coordinates 9<sup>0</sup> 4&#x26;#x02032;N and 7<sup>0</sup> 29&#x26;#x02032;E covering an area of 570 sq miles and a population density of Abuja is located north of the confluence of River Niger and River Benue. It is bordered by the states of Niger to the West and North, Kaduna to the northeast, Nasarawa to the east and south, and Kogi to the southwest. The indigenous tribe in Abuja includes; Koro, Gbari, Gade, Nupe, Gwandara, Dibo, Bassa, Ganagana, and Ebira, and other settler groups like the Hausa, Igbo and Yoruba [
<xref ref-type="bibr" rid="R8">8</xref>]. Abuja is currently made up of six local government areas, namely; Abaji, Abuja Municipal, Bwari, Gwagwalada, Kuje and Kwali. Abuja was the centre of extensive trading in the days preceding British colonization. A wide range of products were involved in the intense trading activity, attesting to the degree of diversification in the economy which, though largely agrarian, had a vibrant manufacturing sector whereby sheanut oil, honey, benniseed, locust bean cakes, ginger, peppers, kola nuts and rice were complemented by iron products, textile, mats leather items and pottery. Abuja was carved out from Niger, Plateau and Kwara states in 1976.</p>
</sec><sec id="sec2">
<title>Materials and Methods</title><p><italic>Exponential Model</italic>: The exponential model is represented by the differential equation</p>

<inline-formula><math><semantics><mrow><mfrac><mrow><mi>d</mi><mi>N</mi></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mi>r</mi><mi>N</mi></mrow></semantics></math></inline-formula><p>with initial condition <math><semantics><mrow><mi>N</mi><mfenced separators="|"><mrow><mn>0</mn></mrow></mfenced><mo>=</mo><msub><mrow><mi>N</mi></mrow><mrow><mn>0</mn></mrow></msub></mrow></semantics></math>. Integrating and multiplying both sides by <math><semantics><mrow><mfrac><mrow><mi>d</mi><mi>t</mi></mrow><mrow><mi>N</mi></mrow></mfrac></mrow></semantics></math> gives</p>

<inline-formula><math><semantics><mrow><mrow><munderover><mo stretchy="false">∫</mo><mrow><mn>0</mn></mrow><mrow><mi>t</mi></mrow></munderover><mrow><mfrac><mrow><mn>1</mn></mrow><mrow><mi>N</mi></mrow></mfrac></mrow></mrow><mo>=</mo><mrow><munderover><mo stretchy="false">∫</mo><mrow><mn>0</mn></mrow><mrow><mi>t</mi></mrow></munderover><mrow><mi>r</mi><mi>s</mi><mi>d</mi><mi>s</mi></mrow></mrow></mrow></semantics></math></inline-formula><p>Therefore,</p>

<inline-formula><math><semantics><mrow><mi>N</mi><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi>N</mi></mrow><mrow><mn>0</mn></mrow></msub><msup><mrow><mi>e</mi></mrow><mrow><mi>r</mi><mi>t</mi></mrow></msup></mrow></semantics></math></inline-formula><p>where the number of population at time <math><semantics><mrow><mi>t</mi></mrow></semantics></math> is <math><semantics><mrow><mi>N</mi><mo>(</mo><mi>t</mi><mo>)</mo></mrow></semantics></math>, <math><semantics><mrow><mi>N</mi><mfenced separators="|"><mrow><mn>0</mn></mrow></mfenced><mo>=</mo><msub><mrow><mi>N</mi></mrow><mrow><mn>0</mn></mrow></msub></mrow></semantics></math> is the initial population size, <math><semantics><mrow><mi>r</mi></mrow></semantics></math> is the population growth rate and <math><semantics><mrow><mi>t</mi><mo>=</mo><mi>t</mi><mi>i</mi><mi>m</mi><mi>e</mi></mrow></semantics></math>. Malthus proposed the exponential model. It is assumed that there are unlimited resources to sustain growth and there is no competition within the population.</p>
<p><italic>Logistic Growth Model</italic>: The logistic growth model is represented by the differential equation</p>

<inline-formula><math><semantics><mrow><mfrac><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">N</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">t</mi></mrow></mfrac><mo>=</mo><mi mathvariant="normal">r</mi><mi mathvariant="normal">N</mi><mfenced separators="|"><mrow><mn>1</mn><mo>-</mo><mfrac><mrow><mi mathvariant="normal">N</mi></mrow><mrow><mi mathvariant="normal">K</mi></mrow></mfrac></mrow></mfenced></mrow></semantics></math></inline-formula><p>where r is the population growth rate and K is the carrying capacity for <math><semantics><mrow><mi>N</mi><mo>(</mo><mi>t</mi><mo>)</mo></mrow></semantics></math> as <math><semantics><mrow><mi>t</mi><mo>→</mo><mi>∞</mi></mrow></semantics></math>.</p>
<p>In the logistic model, a population grows until it attains a maximum capacity [
<xref ref-type="bibr" rid="R6">6</xref>]. Verhulst proposed the logistic model and it does not assume unlimited resources. Instead, it assumes there is a carrying capacity <math><semantics><mrow><mi>K</mi></mrow></semantics></math> for the population. If the population is above K, then the population will decrease, but if below, then it will increase. The solution to the logistic equation can be obtained by applying the separation of variables method. From the equation above, we have:</p>

<inline-formula><math><semantics><mrow><mi mathvariant="normal">N</mi><mfenced separators="|"><mrow><mi mathvariant="normal">t</mi></mrow></mfenced><mo>=</mo><mfrac><mrow><mi mathvariant="normal">K</mi></mrow><mrow><mn>1</mn><mo>+</mo><mi mathvariant="normal">A</mi><msup><mrow><mi mathvariant="normal">e</mi></mrow><mrow><mo>-</mo><mi mathvariant="normal">r</mi><mi mathvariant="normal">t</mi></mrow></msup></mrow></mfrac><mi mathvariant="normal"> </mi></mrow></semantics></math></inline-formula><p>Where</p>

<inline-formula><math><semantics><mrow><mi mathvariant="normal">A</mi><mo>=</mo><mfrac><mrow><mi mathvariant="normal">K</mi><mo>-</mo><msub><mrow><mi mathvariant="normal">N</mi></mrow><mrow><mn>0</mn></mrow></msub></mrow><mrow><msub><mrow><mi mathvariant="normal">N</mi></mrow><mrow><mn>0</mn></mrow></msub></mrow></mfrac></mrow></semantics></math></inline-formula><p>Then, the equation above is the classical logistic growth model. Furthermore, the logistic model can be modified and written as </p>

<inline-formula><math><semantics><mrow><mi mathvariant="normal">N</mi><mfenced separators="|"><mrow><mi mathvariant="normal">t</mi></mrow></mfenced><mo>=</mo><mfrac><mrow><mi mathvariant="normal">K</mi></mrow><mrow><mn>1</mn><mo>+</mo><mi mathvariant="normal">A</mi><msup><mrow><mi mathvariant="normal">k</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msup></mrow></mfrac></mrow></semantics></math></inline-formula><p>Where <math><semantics><mrow><mi>k</mi><mo>=</mo><msup><mrow><mi>e</mi></mrow><mrow><mo>-</mo><mi>r</mi></mrow></msup><mo>.</mo></mrow></semantics></math></p>
<p><italic>Theorem 2.1</italic>. The least upper bound of the logistic model</p>

<inline-formula><math><semantics><mrow><mi>N</mi><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mo>=</mo><mfrac><mrow><mi>K</mi></mrow><mrow><mn>1</mn><mo>+</mo><msup><mrow><mi>A</mi><mi>k</mi></mrow><mrow><mi>t</mi></mrow></msup></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>K</mi></mrow><mrow><mn>2</mn></mrow></mfrac><mo>.</mo><mi mathvariant="normal"> </mi></mrow></semantics></math></inline-formula><p><italic>Proof</italic><italic>.</italic><italic> </italic>From the definition of the least upper bound, we have</p>

<inline-formula><math><semantics><mrow><mn>2</mn><mi mathvariant="normal">K</mi><mo>&gt;</mo><mo>(</mo><mi mathvariant="normal">K</mi><mo>-</mo><mn>2</mn><mi mathvariant="normal">ϵ</mi><mo>)</mo><mo>(</mo><mn>1</mn><mo>+</mo><mi mathvariant="normal">A</mi><mi mathvariant="normal">k</mi><mi mathvariant="normal">t</mi><mo>)</mo></mrow></semantics></math></inline-formula><p>Hence, </p>

<inline-formula><math><semantics><mrow><msup><mrow><mi>k</mi></mrow><mrow><mi>t</mi></mrow></msup><mo>&lt;</mo><mfrac><mrow><mi>K</mi><mo>+</mo><mn>2</mn><mi>ϵ</mi></mrow><mrow><mi>A</mi><mo>(</mo><mi>K</mi><mo>-</mo><mn>2</mn><mi>ϵ</mi><mo>)</mo></mrow></mfrac><mi mathvariant="normal"> </mi></mrow></semantics></math></inline-formula><p>The application of a little algebra gives us</p>

<inline-formula><math><semantics><mrow><mi>t</mi><mo>&lt;</mo><mrow><mrow><mi mathvariant="italic">ln</mi></mrow><mo>⁡</mo><mrow><mfenced open="|" close="|" separators="|"><mrow><mfrac><mrow><mi>K</mi><mo>+</mo><mn>2</mn><mi>ϵ</mi></mrow><mrow><mi>A</mi><mo>(</mo><mi>K</mi><mo>-</mo><mn>2</mn><mi>ϵ</mi><mo>)</mo></mrow></mfrac></mrow></mfenced><mo>-</mo><mrow><mrow><mi mathvariant="italic">ln</mi></mrow><mo>⁡</mo><mrow><mi>k</mi><mo>.</mo></mrow></mrow></mrow></mrow><mi mathvariant="normal"> </mi></mrow></semantics></math></inline-formula><p><italic>Population Growth Sensitivity</italic>: To calculate the growth sensitivity of the predicted population, we differentiate the sensitivity model using the quotient rule. Hence</p>

<inline-formula><math><semantics><mrow><mfrac><mrow><mi>d</mi><mi>N</mi><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mo>-</mo><mi>K</mi><mi>A</mi><msup><mrow><mi>k</mi></mrow><mrow><mi>t</mi></mrow></msup><mi>r</mi></mrow><mrow><msup><mrow><mo>[</mo><mn>1</mn><mo>+</mo><mi>A</mi><msup><mrow><mi>k</mi></mrow><mrow><mi>t</mi></mrow></msup><mo>]</mo></mrow><mrow><mn>2</mn></mrow></msup></mrow></mfrac><mi mathvariant="normal"> </mi><mo>.</mo></mrow></semantics></math></inline-formula><p>Then the equation above is known as population growth sensitivity model. </p>
</sec><sec id="sec3">
<title>Results</title><p><italic>Population Data Visualization</italic>: The population data visualization for Abuja from 1962&#x26;#x02013;2019 is plotted as demonstrated in Fig. 2 using Matlab R2023a;</p>
<fig id="fig1">
<label>Figure 1</label>
<caption>
<p>Population raw data simulation for Abuja</p>
</caption>
<graphic xlink:href="846.fig.001" />
</fig><fig id="fig2">
<label>Figure 2</label>
<caption>
<p>Growth Sensitivity simulation for Abuja</p>
</caption>
<graphic xlink:href="846.fig.002" />
</fig><fig id="fig3">
<label>Figure 3</label>
<caption>
<p><b> </b>Logistic Model simulation for Abuja for 1962 to 2200</p>
</caption>
<graphic xlink:href="846.fig.003" />
</fig><fig id="fig4">
<label>Figure 4</label>
<caption>
<p>Exponential Model simulation for 1962 to 2020</p>
</caption>
<graphic xlink:href="846.fig.004" />
</fig><fig id="fig5">
<label>Figure 5</label>
<caption>
<p>Logistic Model simulation for 1962 to 2020</p>
</caption>
<graphic xlink:href="846.fig.005" />
</fig><title>3.1. Estimation of Parameters</title><p>We estimate the parameters involved for the logistic model and obtained the following values below. All codes were written in Matlab R2012b, and installed on a personal computer with Intel &#x26;#x000ae; Core (TM) i5-2600 CPU@2.30 GHz and 8.00GB (7.78GB usable) RAM running on Windows 8.1 operating system.</p>
<p>Example 3.1: From [
<xref ref-type="bibr" rid="R6">6</xref>], the population data published and gazette for Abuja is <math><semantics><mrow><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mn>0</mn></mrow></msub><mo>=</mo><mn>1,543,293</mn></mrow></semantics></math>, <math><semantics><mrow><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>=</mo><mn>1,693,706</mn></mrow></semantics></math> and <math><semantics><mrow><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mn>2</mn></mrow></msub><mo>=</mo><mn>1,858,777</mn><mi mathvariant="bold-italic"> </mi></mrow></semantics></math>for 2007, 2008, and 2009 respectively. Hence, by applying [
<xref ref-type="bibr" rid="R6">6</xref>], we have</p>

<inline-formula><math><semantics><mrow><mi>K</mi><mo>=</mo><mfrac><mrow><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>P</mi></mrow><mrow><mn>0</mn></mrow></msub><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>-</mo><mn>2</mn><msub><mrow><mi>P</mi></mrow><mrow><mn>0</mn></mrow></msub><msub><mrow><mi>P</mi></mrow><mrow><mn>2</mn></mrow></msub><mo>+</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub><msub><mrow><mi>P</mi></mrow><mrow><mn>2</mn></mrow></msub></mrow></mfenced></mrow><mrow><msubsup><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow><mrow><mn>2</mn></mrow></msubsup><mo>-</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>0</mn></mrow></msub><msub><mrow><mi>P</mi></mrow><mrow><mn>2</mn></mrow></msub></mrow></mfrac><mo>,</mo><mi mathvariant="normal"> </mi><mi mathvariant="normal"> </mi><mi mathvariant="normal"> </mi><mi mathvariant="normal"> </mi><mi>A</mi><mo>=</mo><mfrac><mrow><mi>K</mi><mo>-</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>0</mn></mrow></msub></mrow><mrow><msub><mrow><mi>P</mi></mrow><mrow><mn>0</mn></mrow></msub></mrow></mfrac><mo>,</mo><mi mathvariant="normal"> </mi><mi mathvariant="normal"> </mi><mi mathvariant="normal"> </mi><mi>k</mi><mo>=</mo><mfrac><mrow><msub><mrow><mi>P</mi></mrow><mrow><mn>0</mn></mrow></msub><mo>(</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>2</mn></mrow></msub><mo>-</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>)</mo></mrow><mrow><msub><mrow><mi>P</mi></mrow><mrow><mn>2</mn></mrow></msub><mo>(</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>-</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>0</mn></mrow></msub><mo>)</mo></mrow></mfrac></mrow></semantics></math></inline-formula><p>we can obtain the carrying capacity <math><semantics><mrow><mi mathvariant="bold-italic">K</mi></mrow></semantics></math> and population growth rate <math><semantics><mrow><mi mathvariant="bold-italic">r</mi></mrow></semantics></math> from solving the above system of equations as</p>

<inline-formula><math><semantics><mrow><mi mathvariant="bold-italic">K</mi><mo>=</mo><mn>16</mn><mo>,</mo><mn>946</mn><mo>,</mo><mn>311</mn><mo>,</mo><mn>461</mn><mo>,</mo><mi mathvariant="normal"> </mi><mi mathvariant="normal"> </mi><mi mathvariant="normal"> </mi><mi mathvariant="normal"> </mi><mi mathvariant="normal"> </mi><mi mathvariant="bold-italic">r</mi><mo>=</mo><mn>0</mn><mo>.</mo><mn>093009</mn><mi mathvariant="normal"> </mi><mi mathvariant="normal"> </mi><mi mathvariant="bold-italic">A</mi><mo>=</mo><mn>10979</mn><mo>.</mo><mn>62</mn><mo>,</mo><mi mathvariant="normal"> </mi><mi mathvariant="bold-italic">k</mi><mo>=</mo><mn>0</mn><mo>.</mo><mn>911184881</mn><mi mathvariant="normal"> </mi></mrow></semantics></math></inline-formula><title>3.2. Error Analysis</title><p><italic>Mean Absolute Percentage Error (MAPE):</italic> The Mean Absolute Percentage Error calculation is given by the formulae</p>

<inline-formula><math><semantics><mrow><mi mathvariant="normal">M</mi><mi mathvariant="normal">A</mi><mi mathvariant="normal">P</mi><mi mathvariant="normal">E</mi><mo>=</mo><mfrac><mrow><mn>1</mn></mrow><mrow><mi mathvariant="normal">N</mi></mrow></mfrac><mrow><munderover><mo stretchy="false">∑</mo><mrow><mi mathvariant="normal">t</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi mathvariant="normal">N</mi></mrow></munderover><mrow><mfenced open="|" close="|" separators="|"><mrow><mfrac><mrow><msub><mrow><mi mathvariant="normal">A</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub><mo>-</mo><msub><mrow><mi mathvariant="normal">F</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub></mrow><mrow><msub><mrow><mi mathvariant="normal">A</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub></mrow></mfrac></mrow></mfenced><mo>×</mo><mn>100</mn><mi mathvariant="normal">%</mi></mrow></mrow><mi> </mi></mrow></semantics></math></inline-formula><p>where <math><semantics><mrow><msub><mrow><mi>A</mi></mrow><mrow><mi>t</mi></mrow></msub></mrow></semantics></math> denotes actual value, <math><semantics><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi>t</mi></mrow></msub></mrow></semantics></math>  denotes forecast value and <math><semantics><mrow><mi>N</mi></mrow></semantics></math> is the number of observations.</p>
<p><italic>Root Mean Square Error (RMSE)</italic><italic>: </italic>The Root Mean Square Error calculation is given by the formulae</p>

<inline-formula><math><semantics><mrow><mi mathvariant="normal">R</mi><mi mathvariant="normal">M</mi><mi mathvariant="normal">S</mi><mi mathvariant="normal">E</mi><mo>=</mo><msqrt><mrow><munderover><mo stretchy="false">∑</mo><mrow><mi mathvariant="normal">t</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi mathvariant="normal">N</mi></mrow></munderover><mrow><mfrac><mrow><msup><mrow><mo>(</mo><msub><mrow><mi mathvariant="normal">A</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub><mo>-</mo><msub><mrow><mi mathvariant="normal">F</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub><mo>)</mo></mrow><mrow><mn>2</mn></mrow></msup></mrow><mrow><mi mathvariant="normal">N</mi></mrow></mfrac></mrow></mrow></msqrt><mi> </mi></mrow></semantics></math></inline-formula><p>Where <math><semantics><mrow><msub><mrow><mi>A</mi></mrow><mrow><mi>t</mi></mrow></msub></mrow></semantics></math> denotes actual values, <math><semantics><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi>t</mi></mrow></msub></mrow></semantics></math> denotes forecast values and <math><semantics><mrow><mi>N</mi></mrow></semantics></math> is the number of observations. According to [
<xref ref-type="bibr" rid="R9">9</xref>] MAPE&lt;10% is highly accurate, 10%&#x26;#x02212;20% is good and 21%&#x26;#x02212;50% is reasonable forecasting, while>50% is in accurate forecasting. However, the lower MAPE values are better because they signify smaller percentage errors. </p>
<p></p>
<p><italic>R-Squared</italic>: The R-squared calculation is given by the formulae</p>

<inline-formula><math><semantics><mrow><mi mathvariant="normal">R</mi><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">q</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">d</mi><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mrow><mrow><mo stretchy="false">∑</mo><mrow><msup><mrow><mo>(</mo><msub><mrow><mi mathvariant="normal">y</mi></mrow><mrow><mi mathvariant="normal">i</mi></mrow></msub><mo>-</mo><mover accent="true"><mrow><mi mathvariant="normal">y</mi></mrow><mo>^</mo></mover><mo>)</mo></mrow><mrow><mn>2</mn></mrow></msup></mrow></mrow></mrow><mrow><mrow><mo stretchy="false">∑</mo><mrow><msup><mrow><mo>(</mo><msub><mrow><mi mathvariant="normal">y</mi></mrow><mrow><mi mathvariant="normal">i</mi></mrow></msub><mo>-</mo><mover accent="true"><mrow><mi mathvariant="normal">y</mi></mrow><mo>-</mo></mover><mo>)</mo></mrow><mrow><mn>2</mn></mrow></msup></mrow></mrow></mrow></mfrac><mi> </mi></mrow></semantics></math></inline-formula><p>where <math><semantics><mrow><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></semantics></math> is the actual value and <math><semantics><mrow><mover accent="true"><mrow><mi>y</mi></mrow><mo>^</mo></mover></mrow></semantics></math> is the predicted value. While <math><semantics><mrow><mover accent="true"><mrow><mi>y</mi></mrow><mo>-</mo></mover></mrow></semantics></math> is the mean value. The closer the <math><semantics><mrow><msup><mrow><mi>R</mi></mrow><mrow><mn>2</mn></mrow></msup></mrow></semantics></math> value is to 1, the better the correlation. TheTable <xref ref-type="table" rid="tabtable 1below"> table 1below</xref> shows the results of the error analysis conducted on exponential and logistic population growth models using the simulated data of Abuja from1962 to 2016 shows the following results.</p>
<table-wrap id="tab1">
<label>Table 1</label>
<caption>
<p><b>Table 1</b><b>.</b> Error Analysis of different Growth Models</p>
</caption>

<table>
<thead>
<tr>
<th align="center"><bold>Error Analysis</bold></th>
<th align="center"><bold>Exponential Model</bold></th>
<th align="center"><bold>Logistic Model</bold></th>
<th align="center"></th>
</tr>
</thead>
<tbody>
<tr>
<td align="center">%MAPE</td>
<td align="center">0.98</td>
<td align="center">0.98</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">RMSE</td>
<td align="center">7,811.07</td>
<td align="center">7,817.07</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">R&#x00026;#8722;Squared</td>
<td align="center">0.776</td>
<td align="center">0.776</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">Linear regression</td>
<td align="center">y=60131x&#x00026;#8722;1E+08</td>
<td align="center">y=60129x&#x00026;#8722;1E+08</td>
<td align="center"></td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>

</fn>
</table-wrap-foot>
</table-wrap><p></p>
<title>3.3. Simulations</title><p>We present the population size simulation for Abuja using exponential model and logistic model. A sensitivity analysis model from 2007&#x26;#x02013;2199 is also presented.</p>
<table-wrap id="tab2">
<label>Table 2</label>
<caption>
<p><b></b> Population simulation showing different models.</p>
</caption>

<table>
<thead>
<tr>
<th align="center"><bold>Year</bold></th>
<th align="center"><bold>Time(t)</bold></th>
<th align="center"><bold>NPC Data</bold></th>
<th align="center"><bold>Exponential Model</bold></th>
<th align="center"><bold>Logistic Model</bold></th>
<th align="center"><bold>Sensitivity Model</bold></th>
<th align="center"></th>
</tr>
</thead>
<tbody>
<tr>
<td align="center">2007</td>
<td align="center">0</td>
<td align="center">1,543,293</td>
<td align="center">1,543,293</td>
<td align="center">1,543,293</td>
<td align="center">143527.0452</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2008</td>
<td align="center">1</td>
<td align="center">1,693,706</td>
<td align="center">1,693,707</td>
<td align="center">1,693,706</td>
<td align="center">157514.1336</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2009</td>
<td align="center">2</td>
<td align="center">1,858,777</td>
<td align="center">1,858,780</td>
<td align="center">1,858,777</td>
<td align="center">172864.0015</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2010</td>
<td align="center">3</td>
<td align="center">2,039,937</td>
<td align="center">2,039,942</td>
<td align="center">2,039,934</td>
<td align="center">189709.3651</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2011</td>
<td align="center">4</td>
<td align="center">2,238,752</td>
<td align="center">2,238,761</td>
<td align="center">2,238,744</td>
<td align="center">208195.8528</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2012</td>
<td align="center">5</td>
<td align="center">2,456,945</td>
<td align="center">2,456,957</td>
<td align="center">2,456,928</td>
<td align="center">228483.2589</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2013</td>
<td align="center">6</td>
<td align="center">2,696,403</td>
<td align="center">2,696,419</td>
<td align="center">2,696,372</td>
<td align="center">250746.9188</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2014</td>
<td align="center">7</td>
<td align="center">2,959,199</td>
<td align="center">2,959,219</td>
<td align="center">2,959,147</td>
<td align="center">275179.2161</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2015</td>
<td align="center">8</td>
<td align="center">3,247,608</td>
<td align="center">3,247,633</td>
<td align="center">3,247,526</td>
<td align="center">301991.2366</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2016</td>
<td align="center">9</td>
<td align="center">3,564,126</td>
<td align="center">3,564,157</td>
<td align="center">3,564,002</td>
<td align="center">331414.5804</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2039</td>
<td align="center">32</td>
<td align="center">&#x00026;nbsp;</td>
<td align="center">30,263,745</td>
<td align="center">30,220,701</td>
<td align="center">2805784.642</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2060</td>
<td align="center">53</td>
<td align="center">&#x00026;nbsp;</td>
<td align="center">213,357,987</td>
<td align="center">210,817,405</td>
<td align="center">19363987.38</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2090</td>
<td align="center">83</td>
<td align="center">&#x00026;nbsp;</td>
<td align="center">3,473,789,823</td>
<td align="center">2,884,739,754</td>
<td align="center">222633388.4</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2107</td>
<td align="center">100</td>
<td align="center">&#x00026;nbsp;</td>
<td align="center">16,882,256,619</td>
<td align="center">8,461,080,785</td>
<td align="center">394039070.4</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2110</td>
<td align="center">103</td>
<td align="center">&#x00026;nbsp;</td>
<td align="center">22,315,157,725</td>
<td align="center">9,635,825,443</td>
<td align="center">386620577.9</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2150</td>
<td align="center">143</td>
<td align="center">&#x00026;nbsp;</td>
<td align="center">920,858,296,579</td>
<td align="center">16,640,479,158</td>
<td align="center">27931803.21</td>
<td align="center"></td>
</tr>
<tr>
<td align="center">2199</td>
<td align="center">192</td>
<td align="center">&#x00026;nbsp;</td>
<td align="center">87,759,479,694,267</td>
<td align="center">16,943,045,374</td>
<td align="center">303716.9377</td>
<td align="center"></td>
</tr>
</tbody>
</table>
</table-wrap><p></p>
</sec><sec id="sec4">
<title>Discussion</title><p>It is clear that within the first few days of population increase, the exponential and logistic models predict identical numbers. Gradually, however, the numbers predicted by the logistic model fall below those from the exponential model. By 2107, when the population is about 50% of carrying capacity, the value predicted for the exponential model has doubled that of the logistic model.</p>
<p>Table 2 shows that the Nigerian population from 1962 to 2016 increases from 24,314 people to 3,564,126 million. InFigure <xref ref-type="fig" rid="figfigure 1"> figure 1</xref>, it is observed that the populations over six decades from 1962 till 2016 increases. InFigure <xref ref-type="fig" rid="figfigure 5"> figure 5</xref> the growth increases exponentially as the year increases in the first century from 2017 to 2039 with population increase from 3,911,313 million to 30,220,701 million people. It is observed that the population for the next century from 2107 to 2200 increases slightly and then tends to stability. Hence the population growth becomes slightly stable and then remains stable after 2200. InFigure <xref ref-type="fig" rid="figfigure 3"> figure 3</xref> the growth sensitivity increases for the next three and half decades from 2015 to 2107.This means that the growth rate is still under some form of control or bound. After 2107 the growth sensitivity result shows that the population growth will become stable. From 2050 to 2150 the growth rate is out of control showing an increase in growth sensitivity which implies that everything should be done by the government both in policy and enforcement to ensure that the population growth is been checked early enough and controlled. This will bring the required stability of the population growth to an early time of the first century and not a later time of the second century. Finally, zero growth sensitivity is stable implying that in the second century the population growth of Nigeria beyond 2200 will attain stability.</p>
</sec><sec id="sec5">
<title>Conclusions</title><p>In general, it can be seen that the total population of Abuja increases every year and the estimation result using logistic model has a relatively small error so that the estimation result of the population approaches the actual census result. Besides, the logistic population growth model was used to simulate the future population of Abuja up to 2199. As the population growth moves closer and closer to the carrying capacity, and there are no developmental plans and policies to engage and accommodate the population growth. In that case, there will be a crowding effect on the population density, which will bring about competition, resulting in hardship and hunger. Eventually, this will give birth to crises, insecurity, kidnapping, and other societal crimes. Regarding these factors mentioned, it is expedient that the government invest in massive agricultural reforms to accommodate the growing population, expand Abuja by developing its suburbs, and engage in massive reorientation of the populace on the dangers of uncontrolled births and the education of the girl child. A growing population when strategically harnessed can build a strong economy and the manpower that can build a strong army to defend its territory.</p>
<p></p>
<p><bold>Acknowledgments:</bold><bold> </bold>We want to acknowledge the support received from Simons Foundation for sub-sahara Africa nationals based at the Department of Mathematics and Statistical Sciences, Botswana International University of Science and Technology, Botswana.</p>
<p><bold>Conflicts of Interest</bold>: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.</p>
</sec>
  </body>
  <back>
    <ref-list>
      <title>References</title>
      
<ref id="R1">
<label>[1]</label>
<mixed-citation publication-type="other">National Population Commission (NPC) Nigeria. (2009). Final results of 2006 Census Official Gazette of 2nd February, 2009. Abuja, Nigeria: National Population Commission.
</mixed-citation>
</ref>
<ref id="R2">
<label>[2]</label>
<mixed-citation publication-type="other">Nigeria Demographic and Health Survey (2008). Report of the National Population Commission Federal Republic of Nigeria, Abuja, Nigeria.
</mixed-citation>
</ref>
<ref id="R3">
<label>[3]</label>
<mixed-citation publication-type="other">National Bureau of Statistics (NBS) Nigeria (2018). Demographic Statistics Bulletin. Official Gazette of May, 2018. Abuja, Nigeria: National Bureau of Statistics.
</mixed-citation>
</ref>
<ref id="R4">
<label>[4]</label>
<mixed-citation publication-type="other">Abu,A.S., Akpan,G.U., Okon,U.A., &#x00026; Marvelous,A.I. (2018). Impact of Population Growth Rate on Unemployment in Nigeria. BINGHAM JOURNAL OF ECONOMICS AND ALLIED STUDIES (BJEAS) VOL. 1 NO. 2.
</mixed-citation>
</ref>
<ref id="R5">
<label>[5]</label>
<mixed-citation publication-type="other">Ingiabuna,T.E.,&#x00026;Uzobo,E.(2016). Population and development in Nigeria: An assessment of the National Policy on Population and Sustainable Development. International Journal of Development and Management Review, 11(1), 80-102.
</mixed-citation>
</ref>
<ref id="R6">
<label>[6]</label>
<mixed-citation publication-type="other">Ekakitie,O.,&#x00026; Ekereke,L.(2019).Growth sensitivity of the Nigerian Population and a prediction for the future. International Journal of Applied Science and Mathematical Theory, 5(3).
</mixed-citation>
</ref>
<ref id="R7">
<label>[7]</label>
<mixed-citation publication-type="other">Manu, S. L., &#x00026; Shikaa, S. (2023). Mathematical modeling of Taraba State population growth using exponential and logistic models. Results in Control and Optimization, 100265.
</mixed-citation>
</ref>
<ref id="R8">
<label>[8]</label>
<mixed-citation publication-type="other">Thomas-Emeagwali, G. (1989). Notes on the history of Abuja, central Nigeria. African Study Monographs, 9(4), 191-196.
</mixed-citation>
</ref>
<ref id="R9">
<label>[9]</label>
<mixed-citation publication-type="other">NASIR,M.O.,&#x00026;IBINAYIN,S.(2020) MODELING FOR PROGNOSIS OF NASARAWA STATE POPULATION GROWTH. International Journal of Pure and Applied Science, Vol.19, No.9
</mixed-citation>
</ref>
<ref id="R0">
<label>[0]</label>
<mixed-citation publication-type="other">
</mixed-citation>
</ref>
<ref id="R0">
<label>[0]</label>
<mixed-citation publication-type="other">
</mixed-citation>
</ref>
    </ref-list>
  </back>
</article>