Symbiosis Evolution Model and Behavior of Multiple Resource Agents in the Smart Elderly Care Service Ecosystem

: Population aging has become an important factor restricting China (cid:48) s social and economic development. The smart health and elderly care industry has developed rapidly in the past ﬁve years. However, the service resources among various elderly service providers are relatively isolated and scattered. In other words, the core management problem in the components of the smart elderly care service ecosystem is how to deal with the relationships of interest among multiple resource agents. Thus, the main contribution of this study is to employ symbiosis theory and the logistic growth model to construct a model of the evolution of the symbiosis of multiple resource agents in the smart elderly care service ecosystem. Then, we carry out a stability analysis, and analyze the evolutionary model of two resource agents (cid:48) symbiosis under different values of interdependence coefﬁcients. Finally, we use computer simulations to dynamically simulate the model and comparatively analyze the population density of the hospital–nursing home symbiotic relationship using real cases in China. According to the study, we ﬁnd that the enterprise goal in the smart elderly care service ecosystem should be to maximize the overall value of the multiple resource agents, and the result of the symbiotic evolution between different resource agents depends on the symbiotic interdependence coefﬁcient, while the resource agent uses different strategies under different symbiosis models. Therefore, regulation is needed to ensure the relative fairness of the distribution of value co-creation in the smart elderly care service ecosystem when the resource agent takes actions that beneﬁt itself. Of course, when the ecosystem is in a reciprocal symbiosis model, each resource agent beneﬁts from the activities of the other resource agents, which is ideal in reality; in other words, the best symbiosis model between the two resource agents should be the similar reciprocal symbiosis model.


Introduction
With the continuous improvement of social development and living standards, human life expectancy has significantly improved. Low mortality and low birth rates have become universal phenomena around the world [1]. The proportion of the elderly population has gradually increased, leading to related developments in developed and developing countries around the world [2]. China is one of the countries with a relatively high degree of population aging and faces severe challenges. According to the latest data from the National Bureau of Statistics of China, as of the end of 2019, the number of people aged 60 and over in the country was 253.88 million, accounting for 18.1% of the population. Among them, the population aged 65 and over was 176.03 million people, or 12.6% (http://www.stats.gov.cn/tjsj/zxfb/202002/t20200228_1728913.html, accessed on 28 February 2020), and this number is expected to reach nearly 500 million in China by 2050.
improve in regards to the selection of resource agent objective function and the demand of global resource subject allocation [45,46]. In the field of elderly care services, the city of Oxford in the United Kingdom integrates services for the elderly in management and governance via needs and results goals, service user paths, infrastructure construction, human resources, and workplace management (Oxfordshire County Council, 21 December 2015. Older People's Joint Commissioning Strategy 2013-2016. http://www.oxfordshire.gov. uk/cms/news/2013/jan/older-peoples-joint-commissioning-strategy-2013-2016 accessed on accessed on 28 February 2020).
To make up for the deficiency, many scholars have introduced the concept of "symbiosis" from biology into research on resource agents and have constructed the common governance model of resource agents from the perspective of symbiosis [47]. The key to symbiotic management of resource agents is to clarify the mutual relationship between the symbiotic units (each resource agent) and reveal its essential relation [48], that is, to study the symbiotic behavior model of resource agents. This model should be dynamic or symbiotic evolution [49]. Within the existing literature, the theoretical basis of the symbiotic evolution model is still derived from biology, that is, the logistic growth model [50] used to study the growth, development, or reproduction of plants and animals [51,52]. It is widely used in socio-economic research and is a population dynamic evolution model based on private resource consumption. For example, scholars use the logistic growth model to describe the optimal management of agricultural commodities and fishery production resources [53]; to determine the dynamical mechanism of urbanization [54]; to describe the market penetration of many new products and technologies [55]; to give a summary of world energy usage and source substitution [56]; to describe the industrial revolution as an evolutionary process [57]; and so on. Furthermore, the stable and sustainable development of the service ecosystem requires a set of reasonable income distribution mechanisms to regulate the benefits relationship, enhance the cooperative ability among enterprises, and enhance the efficiency of value creation of the whole service ecosystem.
Under the smart elderly care model, reasonable and effective technical means can improve the efficiency of support for the elderly [58,59]. The ultimate goal of smart elderly care is to combine emerging information technology to establish a comprehensive elderly care service system platform [60,61] at the system level to meet the needs of the elderly and improve the quality of elderly care [62][63][64]. The construction of a smart elderly care service ecosystem can not only realize information sharing between multiple resource agents, but also greatly enhance the coordination capabilities between multiple resource agents, and provide the possibility to establish a mutually beneficial and symbiotic elderly care service ecosystem. The management science problem of how to use symbiosis theory to study the symbiosis mode between two resource agents and even multiple resource agents, and realize the symbiotic evolution of interests among multiple resource agents, should be solved at present.
A review of the existing literature shows that the construction of a multiple resource symbiosis governance model in the smart elderly care service ecosystem has important theoretical research significance and practical application value. The starting point of multiple resource symbiosis governance should not only realize the maximization of profits of a single agent, but also achieve the maximization of the value of multiple resource agents in the service ecosystem, and then a harmonious symbiosis of multiple resource subjects should be built. The key lies in recognizing the symbiotic evolution of the multiple resource agent to promote the sustainable development of enterprises.
Therefore, the main contribution of this study is to employ symbiosis theory and the logistic growth model to construct a model of the evolution of the symbiosis of multiple resource agents in the smart elderly care service ecosystem. This study then verifies the rationality of the constructed theoretical model from two aspects: model simulation and real cases. The rest of this paper is organized as follows: In Section 2, we construct the model of the evolution of the symbiosis of multiple resource agents based on the logistic growth model. Then, in Section 3 we carry out a stability analysis, and then analyze the evolutionary model of two resource agents symbiosis under different values of interdependence coefficients. Next, Section 4 provides the model simulation and a case study, respectively. Finally, Section 5 concludes the paper.

Logistic Growth Model
A multiple resource agent in the smart elderly care service ecosystem, denoted by i, can be viewed as a different population, the density and number of which are in dynamic change. This paper attempts to construct a model of the evolutionary dynamic of multiple resource agents' symbiosis using the logistic growth model. The logistic growth model is a model of population evolution based on exclusive resource consumption. It is used to describe the evolutionary trajectory of species' populations over time in the biological world. When only one biological population exists in a natural environment, the model can be used to describe its number. The differential equation of the logistic growth model can be written as follows [50]: where x(t) is the population density or the number of people at a time t; r is the inherent growth rate; N is the maximum number of people allowed by the environmental resources; and 1 − x N reflects the retardation effect of the population on its own growth caused by the consumption of limited resources [65].
When the logistic growth model is introduced into the field of economics, according to the relevant theories and mathematical models of modern evolutionary economics, the multiple resource agent symbiosis evolutionary dynamic model construction process can be shown as follows.

Model Construction
In this section we describe the model construction. First, the following assumptions are made: Hypothesis 1 (H1). This paper selects revenue v i as the measure of the population density of the resource agent i; the expansion of population size includes the increase of both population number and income, which can be denoted as a function of time, v i (t). It includes not only the meaning in the daily sense, but also the impact of various exogenous influencing factors such as technological innovation, information, specialization and collaboration, institutional arrangements, and market environment [66]. Due to social and market constraints, under certain environmental constraints, the limited nature of resources makes the revenue created by the industry limited. In this case, the resource agent's benefits obtained from the value creation process are limited.

Hypothesis 2 (H2).
The type space and strategy space of the multiple resource agent are assumed to be exogenously given and unchanging; that is, no new species or strategy is involved in this model. Second, the model is constructed as follows: where the revenue of the enterprise v i is the population density of the resource agent; u denotes the strategy of the population, which, in economics, can represent a certain genetic trait of the participant i (e.g., investment routines of shareholders); and f i (u, x) denotes the fitness of the species i, representing an average variation of the density of the species i from one period to another. When only one stakeholder i exists alone in the symbiotic circle, the fitness function of this population can be denoted as follows: where r i is the endogenous growth rate of resource agent i; k i is the resource carrying capacity of the environment where species i resides, that is, the maximum revenue in the symbiotic environment where the stakeholder i resides; v i /k i is the ratio of the size level of stakeholder i against the maximum size, that is, the natural saturation; and 1 − v i /k i denotes the drag effect caused by stakeholder i's consumption of limited resources. We assume that the number of multiple resource agents in the symbiotic circle is s, and that a positive interdependence relation (j) dominates the symbiotic relations among multiple resource agents; that is, j exerts facilitating effects on the growth of revenue v i of i; hence, the fitness function of each agent can be denoted as follows: . . . . . .
Therefore, the model of the evolutionary dynamics of multiple resource agent symbiosis can be denoted as follows: where α ij is the interdependence coefficient, which indicates a facilitating effect of j on i if α ij > 0; no influence if α ij = 0; and a hindering effect when α ij < 0, where α ij is called the competition coefficient.

Model Construction with Two Resource Agents
Assume that in the circle of two symbiotic resource agents, v i (t) indicates revenue generated by i and v j (t) indicates revenue by j. When i or j exist alone, the model of evolutionary symbiosis can be separately denoted as follows: Because it has an interdependence relationship between the two resource agents, the model of evolutionary symbiosis can be separately denoted as follows: where α i , the interdependence coefficient of i, indicates a facilitating effect of j on i; and α j , the interdependence coefficient of j, indicates a facilitating effect of i on j. The combination of Equations (7) and (8) can be considered the model of evolutionary symbiosis:

Evolutionary Model of Multiple Resource Agents' Strategy Learning Rules
When resource agent i changes its strategy u i , the original dynamic evolution model will undergo drastic changes, which can be analyzed with the replicator dynamic model from the discipline of ecology. Equation (11) is a generalized fitness function, from which a generalized multiple resource agent evolution model is obtained. Through appropriate transformation, the above model can be transformed into a dynamic model of strategy frequency (or strategy learning).
where K = n s ∑ i=1 v i is the total number of the population; p i = v i K is the number of species using strategy; u i is a proportion of the population; In the model, the type space and strategy space of species are both exogenous given invariants; the model does not involve the generation of new species or strategies. Therefore, such an evolution model is a local evolution analysis. The above evolution model can describe the selection process of strategies during participants interaction; participants in such models are less conscious and usually adopt unconscious learning rules, and changes in strategy are mainly influenced by natural selection (market competition).

Analysis of the Two Resource Agents' Symbiotic Behavior
According to the results of building the model, we can analyze the two resource agents' symbiotic behaviors. When two resource agents have reached a stable symbiotic state, it can be denoted by the following differential equations: When α i = 0 and α j = 0, then, We have the equilibrium point (k i ,k j ); that is, the incomes for parties of i and j are k i and k j , respectively. According to the service ecosystem perspective, this is an independent symbiosis mode between i and j; that is, no symbiotic relation exists and no party has made contributions to the counterparty; in disregard of other stakeholders, the corporation has not made any new value-added, which is obviously contradictory to the corporation s vision. Therefore, the independent symbiotic mode is inappropriate in contemporary corporations.
When α i = 0 and α j = 0, or α i = 0 and α j = 0, then, We have the equilibrium points (k i , k j (1 + α j )) or (k i (1 + α i ), k j ). The former scenario denotes that i, whose income of equilibrium is k i , cannot obtain benefits from j, but j, whose income of equilibrium is k j (1 + α j ), can gain benefits from i. The latter scenario denotes that i, whose income of equilibrium is k j (1 + α j ), can retain benefits from j, but j, whose income of equilibrium is k i , cannot gain benefits from i. This mode is called partial benefit symbiosis; that is, the resource agent who has not received an effective facilitating effect from its counterparty can reach the original equilibrium value of k i or k j under the existing environmental carrying capacities, while the other resource agent, who has received effective facilitating effects from its counterparty, can break through the existing environmental carrying capacities to form new equilibrium value.
When α i = 0 and α j = 0, then, Equation (16) gives the equilibrium points ( , where the following conditions should be fulfilled depending on realistic meanings: On this occasion, α i α j < 1. From the perspective of a stability analysis on the differential equation, the Taylor expansion of differential equations should be carried out at the equilibrium points: Substituting ( ) into the expansion equation gives the coefficient matrix A: According to the stability theory of differential equations, the condition for a stable node is α i α j < 1; the fulfillment of this condition also falls into the following scenarios: Scenario 1: when α i < 1, α j > 1, and α i α j < 1; or when α i > 1, α j < 1, and α i α j < 1. When α i < 1, α j > 1, and α i α j < 1, the effect of i on j is significantly larger, and vice versa; as a result, the evolutionary process of dynamic growth for i and j will gradually converge toward an equilibrium point. This dynamic evolutionary trend is dissymmetrical, so we call it an asymmetric reciprocal symbiosis model. From the perspective of economics and management, both parties can obtain benefits from the counterparty, that is, share the created corporate value, but such a sharing mechanism is dissymmetrical. In this model, both symbiotic parties can produce new value, which is distributed as per the dissymmetrical mechanism, giving rise to the asynchrony of both parties evolutions. Scenario 2: when α i < 1, α j < 1, and α i α j < 1. When α i < 1, α j < 1, and α i α j < 1, the effect of i on j is not much different from the effect of j on i, with nearly the same or equal effect. The dynamic evolutionary trend is close or similar to the state of symmetric evolution. Therefore, when α i = α j , we call it a symmetric reciprocal symbiosis model, while when α i = α j , we call it a similar reciprocal symbiosis model. From the perspective of economics and management, under the symmetric reciprocal symbiosis model, each party can benefit from the behaviors of the counterparty, and both symbiotic parties can produce new value, which is distributed as per a symmetric mechanism, and the evolutions of both parties are synchronized. However, the symmetric reciprocal symbiosis model is ideal in reality, and the best symbiosis model between the two resource agents should be the similar reciprocal symbiosis model.
When the Interdependence Coefficient is Negative When the dependency coefficient is less than zero, it can be specifically analyzed as follows: When α i = 0 and α i < 0, or α i < 0 and α j = 0, then i and j are in the symbiotic mode, one side is damaged, and the other party is neither profitable nor harmful. When α i > 0 and α j < 0, or α i < 0 and α j > 0, then i and j are in the parasitic mode, although it is not necessarily harmful to the host; but there is a transfer of the host's interest to the parasite, and no new value is generated. When α i < 0 and α j < 0, i and j are in a competitive symbiotic relationship. According to the degree of influence of the competing parties on each other, i and j are in a symmetric competitive symbiotic relationship or an asymmetric competitive symbiotic relationship.
To summarize, two resource agents symbiotic relations and equilibrium points can be determined according to difference values of the interdependence coefficient. The above scenarios are summarized in Table 1. Table 1. Two resource agents symbiosis behavior modes under different values of interdependence coefficients.

Interdependence Coefficient Equilibrium Points Symbiosis Behavior Mode Remarks
Cannot exist in the long term Asymmetric reciprocal symbiosis Asynchronous evolutions; instability Similar reciprocal symbiosis Similar to equilibrium Non-cooperative win-win model Competitive symbiosis

Model Simulation
The numerical simulation of the symbiosis model assigns the parameters of the model based on the experience of large enterprises in China [66]. Based on previous research [67], this paper assumes that the natural growth rates of resource agents i and j, respectively, are 0.03 and 0.01. The variable K represents the natural load and can represent the total value in an economic sense, assuming that the value load of i and j each is 100 units, and the initial population density of resource subjects i and j is 1 unit.
This part gives the simulation results of 8 kinds of symbiosis behavior mode in Table 1. The left plot shows the population density of the resource agent, V, over several thousand iterations, while the right plot shows the strategy frequency of each resource agent, P, over several thousand iterations. At last, we take three resource agents as the case of the simulation of the evolution curve of multiple resource symbiotic dynamics.

Simulation of the Independent Symbiosis Model
The independent symbiosis (α i = 0, α j = 0) represents two resource agents that have no relationship with each other. In Figure 1a, after nearly 1000 iterations, resource agent i grows faster than resource agent j due to its larger internal growth rate, and it reaches the environmental load first. In Figure 1b, the two resource agents reach the frequency of about 0.5 from opposite directions after about 1000 iterations. The difference is that i first rises and then falls, whereas j first falls and then rises. Figure 1 shows that the two resource agents are using the same strategy; that is, the strategy of development independent of each other is equal to the development opportunities of the two resource subjects.

Simulation of the Parasitic Symbiosis Model
One of the parasitic symbiosis models is 0 i α > and 0 j α < ; that is, j has a positive effect on i while i has a negative effect on j. In other words, i is a parasite, which is conducive to the development of i but not to the development of host j. Figure 2 presents in the dynamic evolution of population density, i quickly increases to the environmental carrying capacity, whereas j first increases and then decreases and basically does not share the rights and interests. In the choice of strategy, the frequency of the parasitic i continues to rise to 1, indicating that i prefers the parasitic strategy, while j tends to give up or go to a p value of 0 under the parasitic strategy. One of the parasitic symbiosis models is 0 i α < and 0 j α > ; that is, i has a positive effect on j, while j has a negative effect on i. In other words, j is a parasite, which is conducive to the development of j but not to the development of host i. Figure 3 presents in the dynamic evolution of population density, due to its own internal growth, i first increases and then declines to basically share no rights and interests, whereas j quickly grows nearly to the environmental load. In the selection of dynamic strategy, the frequency of the parasitic j tends to be 1 after decreasing first and then increasing, indicating that j prefers the parasitic strategy, while i tends to give up or go to a P value of 0 under the parasitic strategy.

Simulation of the Parasitic Symbiosis Model
One of the parasitic symbiosis models is α j < 0 and α i > 0; that is, j has a positive effect on i while i has a negative effect on j. In other words, i is a parasite, which is conducive to the development of i but not to the development of host j. Figure 2 presents in the dynamic evolution of population density, i quickly increases to the environmental carrying capacity, whereas j first increases and then decreases and basically does not share the rights and interests. In the choice of strategy, the frequency of the parasitic i continues to rise to 1, indicating that i prefers the parasitic strategy, while j tends to give up or go to a p value of 0 under the parasitic strategy.

Simulation of the Parasitic Symbiosis Model
One of the parasitic symbiosis models is 0 i   and 0 j   ; that is, j has a positive effect on i while i has a negative effect on j. In other words, i is a parasite, which is conducive to the development of i but not to the development of host j. Figure 2 presents in the dynamic evolution of population density, i quickly increases to the environmental carrying capacity, whereas j first increases and then decreases and basically does not share the rights and interests. In the choice of strategy, the frequency of the parasitic i continues to rise to 1, indicating that i prefers the parasitic strategy, while j tends to give up or go to a p value of 0 under the parasitic strategy. One of the parasitic symbiosis models is 0 i   and 0 j   ; that is, i has a positive effect on j, while j has a negative effect on i. In other words, j is a parasite, which is conducive to the development of j but not to the development of host i. Figure 3 presents in the dynamic evolution of population density, due to its own internal growth, i first increases and then declines to basically share no rights and interests, whereas j quickly grows nearly to the environmental load. In the selection of dynamic strategy, the frequency of the parasitic j tends to be 1 after decreasing first and then increasing, indicating that j prefers the parasitic strategy, while i tends to give up or go to a P value of 0 under the parasitic strategy.  One of the parasitic symbiosis models is α i < 0 and α j > 0; that is, i has a positive effect on j, while j has a negative effect on i. In other words, j is a parasite, which is conducive to the development of j but not to the development of host i. Figure 3 presents in the dynamic evolution of population density, due to its own internal growth, i first increases and then declines to basically share no rights and interests, whereas j quickly grows nearly to the environmental load. In the selection of dynamic strategy, the frequency of the parasitic j tends to be 1 after decreasing first and then increasing, indicating that j prefers the parasitic strategy, while i tends to give up or go to a P value of 0 under the parasitic strategy.

Simulation of the Parasitic Symbiosis Model
One of the parasitic symbiosis models is 0 i   and 0 j   ; that is, j has a positive effect on i while i has a negative effect on j. In other words, i is a parasite, which is conducive to the development of i but not to the development of host j. Figure 2 presents in the dynamic evolution of population density, i quickly increases to the environmental carrying capacity, whereas j first increases and then decreases and basically does not share the rights and interests. In the choice of strategy, the frequency of the parasitic i continues to rise to 1, indicating that i prefers the parasitic strategy, while j tends to give up or go to a p value of 0 under the parasitic strategy. One of the parasitic symbiosis models is 0 i   and 0 j   ; that is, i has a positive effect on j, while j has a negative effect on i. In other words, j is a parasite, which is conducive to the development of j but not to the development of host i. Figure 3 presents in the dynamic evolution of population density, due to its own internal growth, i first increases and then declines to basically share no rights and interests, whereas j quickly grows nearly to the environmental load. In the selection of dynamic strategy, the frequency of the parasitic j tends to be 1 after decreasing first and then increasing, indicating that j prefers the parasitic strategy, while i tends to give up or go to a P value of 0 under the parasitic strategy.

Simulation of the Competitive Symbiosis Model
The competitive symbiosis mode is α i < 0 and α j < 0. The dynamic simulation curve of the competitive symbiosis model shows that although i and j are competitive, the absolute value of the dependence coefficient is less than 1, and both the population densities increase and then stabilize, which is reflected in the change of strategy of the dynamic selection curve. Figure 4 presents the state of this trade-off, forming a symmetrical symbiotic state of competition.

Simulation of the Competitive Symbiosis Model
The competitive symbiosis mode is 0 i   and < 0 . The dynamic simulation curve of the competitive symbiosis model shows that although i and j are competitive, the absolute value of the dependence coefficient is less than 1, and both the population densities increase and then stabilize, which is reflected in the change of strategy of the dynamic selection curve. Figure 4 presents the state of this trade-off, forming a symmetrical symbiotic state of competition.

Simulation of Partial Harm Symbiosis Model
One of the partial harm symbiosis models is that the resource agent i is harmful to j, losing the benefit to j, which makes its population density increase slowly. The dynamic selection curve of the strategy changes in the opposite direction. As shown in Figure 5, the frequency of resource agent i tends toward 1, and it is willing to continue to execute the strategy, while j tends to give up the strategy. One of the partial harm symbiosis models is that the resource agent j is harmful to i, losing the benefit of i, which makes its population density increase change from quickly to slowly, while the population density of j tends to increase quickly. The dynamic selection curve of the strategy changes in the opposite direction. As shown in Figure 6, the resource agent i first rises and then falls to 0; thus, i chooses to give up the strategy. Meanwhile, resource agent j first falls and then rises to 1, showing that it is willing to continue to execute the strategy.

Simulation of Partial Harm Symbiosis Model
One of the partial harm symbiosis models is that the resource agent i is harmful to j, losing the benefit to j, which makes its population density increase slowly. The dynamic selection curve of the strategy changes in the opposite direction. As shown in Figure 5, the frequency of resource agent i tends toward 1, and it is willing to continue to execute the strategy, while j tends to give up the strategy.
R PEER REVIEW 11 of 22

Simulation of the Competitive Symbiosis Model
The competitive symbiosis mode is 0 i   and < 0 . The dynamic simulation curve of the competitive symbiosis model shows that although i and j are competitive, the absolute value of the dependence coefficient is less than 1, and both the population densities increase and then stabilize, which is reflected in the change of strategy of the dynamic selection curve. Figure 4 presents the state of this trade-off, forming a symmetrical symbiotic state of competition.

Simulation of Partial Harm Symbiosis Model
One of the partial harm symbiosis models is that the resource agent i is harmful to j, losing the benefit to j, which makes its population density increase slowly. The dynamic selection curve of the strategy changes in the opposite direction. As shown in Figure 5, the frequency of resource agent i tends toward 1, and it is willing to continue to execute the strategy, while j tends to give up the strategy. One of the partial harm symbiosis models is that the resource agent j is harmful to i, losing the benefit of i, which makes its population density increase change from quickly to slowly, while the population density of j tends to increase quickly. The dynamic selection curve of the strategy changes in the opposite direction. As shown in Figure 6, the resource agent i first rises and then falls to 0; thus, i chooses to give up the strategy. Meanwhile, resource agent j first falls and then rises to 1, showing that it is willing to continue to execute the strategy.  One of the partial harm symbiosis models is that the resource agent j is harmful to i, losing the benefit of i, which makes its population density increase change from quickly to slowly, while the population density of j tends to increase quickly. The dynamic selection curve of the strategy changes in the opposite direction. As shown in Figure 6, the resource agent i first rises and then falls to 0; thus, i chooses to give up the strategy. Meanwhile, resource agent j first falls and then rises to 1, showing that it is willing to continue to execute the strategy.
losing the benefit of i, which makes its population density increase change from quickly to slowly, while the population density of j tends to increase quickly. The dynamic selection curve of the strategy changes in the opposite direction. As shown in Figure 6, the resource agent i first rises and then falls to 0; thus, i chooses to give up the strategy. Meanwhile, resource agent j first falls and then rises to 1, showing that it is willing to continue to execute the strategy.

Simulation of Partial Benefit Symbiosis Model
One of the partial benefit symbiosis models is that resource agent j benefits from the behavior of i, which makes the population density of j rise. While i is not affected by j, its population growth rate first rises and then declines, depending on its internal growth rate, indicating that its dynamic benefits are still affected under the interaction with j. The effect is reflected in the change of the strategy's dynamic selection curve due to its own internal growth; the curve of i appears to rise first in this strategy selection and then tends to give up the strategy. Meanwhile, the curve of j appears to fall first and then rise. It can be seen in Figure 7.
R PEER REVIEW 12 of 22

Simulation of Partial Benefit Symbiosis Model
One of the partial benefit symbiosis models is that resource agent j benefits from the behavior of i, which makes the population density of j rise. While i is not affected by j, its population growth rate first rises and then declines, depending on its internal growth rate, indicating that its dynamic benefits are still affected under the interaction with j. The effect is reflected in the change of the strategy's dynamic selection curve due to its own internal growth; the curve of i appears to rise first in this strategy selection and then tends to give up the strategy. Meanwhile, the curve of j appears to fall first and then rise. It can be seen in Figure 7. The second partial benefit symbiosis model is that resource agent i benefits from the behavior of j. Because i's own internal growth rate is higher than that of j, its population density rises faster and then tends to be stable. However, j is not affected by i and depends on its internal growth rate, making its population density rise slowly and tend to be stable. As reflected in the change of the strategy's dynamic selection curve, i still tends to rise and then decline in the strategy selection, but the downward trend is significantly less than the degree not affected by j, tending to be stable at 0.6; that is, in the long run, the frequency of the selection of this strategy by i is 0.6. Meanwhile, j can benefit from its own growth, and the strategy selection curve shows a trend of falling first and then rising and tends to be stable at 0.4; that is, in the long run, the frequency of the selection of this strategy by j is 0.4. It can be seen in Figure 8.  The second partial benefit symbiosis model is that resource agent i benefits from the behavior of j. Because i's own internal growth rate is higher than that of j, its population density rises faster and then tends to be stable. However, j is not affected by i and depends on its internal growth rate, making its population density rise slowly and tend to be stable. As reflected in the change of the strategy's dynamic selection curve, i still tends to rise and then decline in the strategy selection, but the downward trend is significantly less than the degree not affected by j, tending to be stable at 0.6; that is, in the long run, the frequency of the selection of this strategy by i is 0.6. Meanwhile, j can benefit from its own growth, and the strategy selection curve shows a trend of falling first and then rising and tends to be stable at 0.4; that is, in the long run, the frequency of the selection of this strategy by j is 0.4. It can be seen in Figure 8. then decline in the strategy selection, but the downward trend is significantly less than the degree not affected by j, tending to be stable at 0.6; that is, in the long run, the frequency of the selection of this strategy by i is 0.6. Meanwhile, j can benefit from its own growth, and the strategy selection curve shows a trend of falling first and then rising and tends to be stable at 0.4; that is, in the long run, the frequency of the selection of this strategy by j is 0.4. It can be seen in Figure 8.

Simulation of Asymmetric Reciprocal Symbiosis Model
Asymmetrical reciprocal symbiosis is actually an asynchronous evolution, which brings great instability. One of the patterns is that the influence on resource agent i by resource agent j is higher than that of resource agent i on resource agent j, because the internal growth rate of i is greater than that of j, so the population density of i grows faster, and the population density of j first rises and then falls. These patterns are reflected in the changes of the strategy's dynamic selection curve: in Figure 9, for i, the curve of the strategy selection continues to rise and tends to 1, so i sticks to the strategy, while j shows a continuous decline in the strategy selection, tending to 0, indicating that j gives up the strategy.

Simulation of Asymmetric Reciprocal Symbiosis Model
Asymmetrical reciprocal symbiosis is actually an asynchronous evolution, which brings great instability. One of the patterns is that the influence on resource agent i by resource agent j is higher than that of resource agent i on resource agent j, because the internal growth rate of i is greater than that of j, so the population density of i grows faster, and the population density of j first rises and then falls. These patterns are reflected in the changes of the strategy's dynamic selection curve: in Figure 9, for i, the curve of the strategy selection continues to rise and tends to 1, so i sticks to the strategy, while j shows a continuous decline in the strategy selection, tending to 0, indicating that j gives up the strategy. The second of the patterns is that the influence of resource agent i on resource agent j is higher than that of resource agent j on resource agent i, because the internal growth rate of i is greater than that of j, so the population density of i first rises and then falls, and the population density of j first rises slowly and then grows more quickly. These patterns are reflected in the changes of the strategy's dynamic selection curve: as shown in Figure  10, for i, the curve of the strategy selection tends to first rise and then fall, so i tends to give up the strategy, while for j, the curve of the strategy selection tends to first fall and then rise, indicating that j tends to choose the strategy.  The second of the patterns is that the influence of resource agent i on resource agent j is higher than that of resource agent j on resource agent i, because the internal growth rate of i is greater than that of j, so the population density of i first rises and then falls, and the population density of j first rises slowly and then grows more quickly. These patterns are reflected in the changes of the strategy's dynamic selection curve: as shown in Figure 10, for i, the curve of the strategy selection tends to first rise and then fall, so i tends to give up the strategy, while for j, the curve of the strategy selection tends to first fall and then rise, indicating that j tends to choose the strategy. The second of the patterns is that the influence of resource agent i on resource agent j is higher than that of resource agent j on resource agent i, because the internal growth rate of i is greater than that of j, so the population density of i first rises and then falls, and the population density of j first rises slowly and then grows more quickly. These patterns are reflected in the changes of the strategy's dynamic selection curve: as shown in Figure  10, for i, the curve of the strategy selection tends to first rise and then fall, so i tends to give up the strategy, while for j, the curve of the strategy selection tends to first fall and then rise, indicating that j tends to choose the strategy.

Simulation of a Similar Reciprocal Symbiosis Model
Similar to the reciprocal symbiosis model, both i and j will influence each other, and the degree of influence is weak; that is, the dependency coefficient is less than 1. In this model, the population density curve and the dynamic strategy selection curve of the two

Simulation of a Similar Reciprocal Symbiosis Model
Similar to the reciprocal symbiosis model, both i and j will influence each other, and the degree of influence is weak; that is, the dependency coefficient is less than 1. In this model, the population density curve and the dynamic strategy selection curve of the two resource subjects are basically synchronized, but they are slightly different because of the different degree of influence between i and j. The i population density curve shows the trend of rapid growth and then decline, while the j population density curve shows a steady growth trend; but the trends of both i and j are synchronous, and the j population density curve is below that of i. These patterns are reflected in the changes of the strategy's dynamic selection curve; as shown in Figure 11, there are fluctuations in i and early in j that show that the frequency of i increases rapidly, whereas the frequency of j decreases rapidly; then i decreases steadily to 0.6, while j rises to 0.4. Overall, the strategy is slightly beneficial to i. Figure 10. One of the asymmetric reciprocal symbiosis models; the value is (0.6, 1.2).

Simulation of a Similar Reciprocal Symbiosis Model
Similar to the reciprocal symbiosis model, both i and j will influence each other, and the degree of influence is weak; that is, the dependency coefficient is less than 1. In this model, the population density curve and the dynamic strategy selection curve of the two resource subjects are basically synchronized, but they are slightly different because of the different degree of influence between i and j. The i population density curve shows the trend of rapid growth and then decline, while the j population density curve shows a steady growth trend; but the trends of both i and j are synchronous, and the j population density curve is below that of i. These patterns are reflected in the changes of the strategy's dynamic selection curve; as shown in Figure 11, there are fluctuations in i and early in j that show that the frequency of i increases rapidly, whereas the frequency of j decreases rapidly; then i decreases steadily to 0.6, while j rises to 0.4. Overall, the strategy is slightly beneficial to i. When resource agent j is slightly affected by i, the i population density curve shows a trend of rapid growth first and then tends to decrease, while the j population density When resource agent j is slightly affected by i, the i population density curve shows a trend of rapid growth first and then tends to decrease, while the j population density curve shows a rapid growth trend and exceeds the i curve at about 1000 iterations, but the trends of both i and j are synchronous. These patterns are reflected in the changes of the strategy's dynamic selection curves; as shown in Figure 12, the frequency of i rises first and then decreases, while the frequency of j decreases first and then rises and exceeds the i curve at about 1000 iterations. Overall, the strategy is slightly beneficial to j.

R PEER REVIEW 14 of 22
curve shows a rapid growth trend and exceeds the i curve at about 1000 iterations, but the trends of both i and j are synchronous. These patterns are reflected in the changes of the strategy's dynamic selection curves; as shown in Figure 12, the frequency of i rises first and then decreases, while the frequency of j decreases first and then rises and exceeds the i curve at about 1000 iterations. Overall, the strategy is slightly beneficial to j.

Simulation of Symmetric Reciprocal Symbiosis Model
The symmetric reciprocal symbiosis model describes two resource agents with the same influence on each other and the interdependence coefficient of less than 1, which is an ideal state. In this state, the population density of resource agent i increases rapidly and then decreases, and the population density curve of j increases first slowly and then rapidly, and then the population density curves of the two overlap for thousands of iterations. These patterns are reflected in the changes of the strategy's dynamic selection curves; as shown in Figure 13, the frequency of i first rises and then falls, while the frequency of j first falls and then rises, and it overlaps with the i curve at a frequency of 0.5 from about 1000 iterations onward, indicating that the selection probabilities of the two resource agents in the strategy are the same. That is, the development opportunities of the

Simulation of Symmetric Reciprocal Symbiosis Model
The symmetric reciprocal symbiosis model describes two resource agents with the same influence on each other and the interdependence coefficient of less than 1, which is an ideal state. In this state, the population density of resource agent i increases rapidly and then decreases, and the population density curve of j increases first slowly and then rapidly, and then the population density curves of the two overlap for thousands of iterations. These patterns are reflected in the changes of the strategy's dynamic selection curves; as shown in Figure 13, the frequency of i first rises and then falls, while the frequency of j first falls and then rises, and it overlaps with the i curve at a frequency of 0.5 from about 1000 iterations onward, indicating that the selection probabilities of the two resource agents in the strategy are the same. That is, the development opportunities of the two resource entities are equal. an ideal state. In this state, the population density of resource agent i increases rapidly and then decreases, and the population density curve of j increases first slowly and then rapidly, and then the population density curves of the two overlap for thousands of iterations. These patterns are reflected in the changes of the strategy's dynamic selection curves; as shown in Figure 13, the frequency of i first rises and then falls, while the frequency of j first falls and then rises, and it overlaps with the i curve at a frequency of 0.5 from about 1000 iterations onward, indicating that the selection probabilities of the two resource agents in the strategy are the same. That is, the development opportunities of the two resource entities are equal.

Simulation of Multiple Resource Agents Model
Simulation of the evolution curve of multiple resource agents' symbiotic dynamics (in this case, three resource agents) can be established on the basis of that of two resource agents, assuming that the three resource agents have an inherent growth rate of 0.2, the environmental carrying rate is 100, the interdependence coefficient of j on i is 0.1, the interdependence coefficient of m on i is 0.1, the interdependence coefficient of i on j is 0.5, the interdependence coefficient of m on j is 0.5, the interdependence coefficient of i on m is 0.9, and the interdependence coefficient of j on m is 0.9. It can be seen in the Figure 14 that the interdependence coefficient of i is 0.2, the interdependence coefficient of j is 1, and the interdependence coefficient of m is 1.8. That is, the three resource agents are in a symbiotic state; their population density evolution curves and strategy selection curves are shown in Figure 14. This model is obviously beneficial to the resource agent m.

Simulation of Multiple Resource Agents Model
Simulation of the evolution curve of multiple resource agents' symbiotic dynamics (in this case, three resource agents) can be established on the basis of that of two resource agents, assuming that the three resource agents have an inherent growth rate of 0.2, the environmental carrying rate is 100, the interdependence coefficient of j on i is 0.1, the interdependence coefficient of m on i is 0.1, the interdependence coefficient of i on j is 0.5, the interdependence coefficient of m on j is 0.5, the interdependence coefficient of i on m is 0.9, and the interdependence coefficient of j on m is 0.9. It can be seen in the Figure 14 that the interdependence coefficient of i is 0.2, the interdependence coefficient of j is 1, and the interdependence coefficient of m is 1.8. That is, the three resource agents are in a symbiotic state; their population density evolution curves and strategy selection curves are shown in Figure 14. This model is obviously beneficial to the resource agent m.

Case Study
American nursing homes can provide medical treatment, long-term care, and endof-life care, which can be seen as a service complex of medical and elderly care [68][69][70]. Moreover, medical care and elderly care in China are separate [71,72]; that is, the division of labor between hospitals and nursing homes is clearly different from that in other countries. Many scholars have borrowed from the American nursing home model and proposed an old-age care model with Chinese characteristics. Scholars have suggested that medical resources and elderly care resources should be organically integrated to provide professional and continuous health care services for the elderly to form a service model that integrates medical care and elderly care [2,73]. The government has also issued a series of policy documents to promote medical care integration (http://www.gov.cn/xinwen/2019-10/26/content_5445271.htm accessed on 25 September 2020). The integration of medical care and elderly care is of great significance to cope with the challenges of aging. How to handle the relationship between hospitals and nursing homes remains to be solved in China. It is also the core duality in the construction of a smart elderly care service ecosystem. Chinese companies represented by JAG(a diversified investment enterprise

Case Study
American nursing homes can provide medical treatment, long-term care, and endof-life care, which can be seen as a service complex of medical and elderly care [68][69][70]. Moreover, medical care and elderly care in China are separate [71,72]; that is, the division of labor between hospitals and nursing homes is clearly different from that in other countries. Many scholars have borrowed from the American nursing home model and proposed an old-age care model with Chinese characteristics. Scholars have suggested that medical resources and elderly care resources should be organically integrated to provide professional and continuous health care services for the elderly to form a service model that integrates medical care and elderly care [2,73]. The government has also issued a series of policy documents to promote medical care integration (http://www.gov.cn/xinwen/2019-10/26/content_5445271.htm accessed on 25 Septem-ber 2020). The integration of medical care and elderly care is of great significance to cope with the challenges of aging. How to handle the relationship between hospitals and nursing homes remains to be solved in China. It is also the core duality in the construction of a smart elderly care service ecosystem. Chinese companies represented by JAG (a diversified investment enterprise group, Anhui, Jingan have explored the integration of medical and elderly care and the construction of a smart elderly care service ecosystem.
JAG was established in 2001. It has integrated health care, elderly care services, international trade, financial services, urban public services, real estate development, property management, construction and building materials, and so on. The Health Industry Management Group was one of three major industry sectors in JAG, which was established in 2010. It started with the establishment of JAG Nursing Home and JAG Hospital. By providing more comprehensive and professional elderly care services, it builds the core layer of a smart elderly care service ecosystem based on the combination of medical and elderly care. On this basis, JAG will combine elderly care with the real estate industry and implement the strategy of building a health care complex for the elderly. The health care complex will include the use of the Internet, big data, and other information technology, on site, to implement the integration of institutions, communities, and homes; in terms of services, the integration of medical care, nursing care, and housekeeping will be implemented. In order to achieve this strategy, JAG has successively invested in institutions such as vocational training schools, nursing homes and hospitals in many cities, and other institutions. This investment will help JAG cover the whole field, all the elements, and all the functions of elderly care, and implement the strategies of basic informatization, management informatization, and service informatization to establish the internal ecosystem network of the enterprise. Finally, by establishing a development model of production, education, and research, the social network of the enterprise will be extended.
In the smart elderly care service ecosystem, the relationship between the nursing home and hospital is the most basic dual relationship or symbiosis. While providing medical services to society, JAG Hospital also provides medical security for the elderly care services of JAG Nursing Home. JAG Nursing Home is the geriatrics department of JAG Hospital. It has a geriatrics department stationed on every floor, realizing the integrated medical care and elderly care service model, and it has achieved good performance.
This article intends to introduce the hospital-nursing home symbiosis into the evolution model constructed in the previous sections, using a case analysis of JAG, which fits the model to a certain extent and meets the two basic assumptions. First, hospitals and nursing homes are the two main resource entities in the smart elderly care service ecosystem, and the acquisition and changes of their operating income are related to factors such as resource input, the symbiotic relationship of nursing homes or hospitals, information technology innovations, institutional arrangements, market competition, and so on. Both hospital and nursing home incomes are functions of time t and have an upper limit of allowable resources under certain conditions, which meets the first assumption. Second, the JAG medical and elderly care integration model is actually the core link in the construction of the smart elderly care service ecosystem. That is to say, since 2011, the type space and strategic space in the JAG medical and elderly care integration service ecosystem have not fundamentally changed, which meets the second assumption. Third, similar to the original logistic growth model, the population density curve of the evolution model of multiple resource symbiosis shows an S-shaped growth over time; that is, it generally shows an upward trend that accelerates first and then decelerates within a certain period of time. However, after a certain point in time, it will decelerate and then accelerate. In this paper, the operating incomes of JAG Nursing Home and JAG Hospital are used as the population density parameters to study the trend and direction of the change curve for the strategy of combining medical care and elderly care, with Figure 15 showing the details. The figure shows that the population density curve changes of nursing homes and hospitals are in line with the original model's trends. of time. However, after a certain point in time, it will decelerate and then accelerate. this paper, the operating incomes of JAG Nursing Home and JAG Hospital are used the population density parameters to study the trend and direction of the change cur for the strategy of combining medical care and elderly care, with Figure 15 showing t details. The figure shows that the population density curve changes of nursing homes a hospitals are in line with the original model's trends. Next, this paper carries out the model calculation and analysis of the hospital-nu ing home symbiosis model.

Estimation of Hospital-Nursing Home Dependence Coefficient
In the model of the evolution of symbiosis, the determination of the dependence efficient is the key to the model calculation, but there is no unified measurement standa in the existing literature. The dependence coefficient reflects the degree of influence o species on the population density of another species. The degree of symbiosis in biologi symbiosis theory is similar, so this paper uses the degree of symbiosis to estimate the d pendence coefficient.
The degree of symbiosis refers to the degree of correlation between the changes the quality parameters of the two symbiosis units, reflecting the degree of mutual inf ence of the quality parameters of the two symbiosis units; the quality parameters refer the parameters that reflect the intrinsic properties of the symbiosis unit. This article sele operating income as the main quality parameter of hospitals and nursing homes. The sy biosis degree ij  reflects the influence of the change rate of the main quality paramet Next, this paper carries out the model calculation and analysis of the hospital-nursing home symbiosis model.

Estimation of Hospital-Nursing Home Dependence Coefficient
In the model of the evolution of symbiosis, the determination of the dependence coefficient is the key to the model calculation, but there is no unified measurement standard in the existing literature. The dependence coefficient reflects the degree of influence of a species on the population density of another species. The degree of symbiosis in biological symbiosis theory is similar, so this paper uses the degree of symbiosis to estimate the dependence coefficient.
The degree of symbiosis refers to the degree of correlation between the changes of the quality parameters of the two symbiosis units, reflecting the degree of mutual influence of the quality parameters of the two symbiosis units; the quality parameters refer to the parameters that reflect the intrinsic properties of the symbiosis unit. This article selects operating income as the main quality parameter of hospitals and nursing homes. The symbiosis degree δ ij reflects the influence of the change rate of the main quality parameters of nursing homes v i on the change rate of hospitals' main quality parameters v j . Therefore, the symbiosis degree δ ij is used to estimate the dependence coefficient of hospitals, that is, α j = δ ij (see Equation (20)). Similarly, the symbiosis degree δ ji reflects the influence of the change rate of the main quality parameters of the hospital v j on the change rate of nursing homes' main quality parameters v i . Therefore, the symbiosis degree δ ji is used to estimate the dependency coefficient of nursing homes, that is, α i = δ ji (see Equation (21)).
where ∆ is the difference, which is equal to the difference between the amount of this year and the previous year, and d is the differential, which represents the derivative relationship. The difference is converted into a derivative relationship in the formula, which is convenient for solving by the linear regression method, below.

Hospital-Nursing Home Population Density Function Relationship and Calculation of Dependence Coefficient
Equations (20) and (21) show that the solution of the dependence coefficient requires the derivative or functional relationship between the operating income of the nursing home v i and the operating income of the hospital v j to establish a linear relationship between the two variables; we then use IBM SPSS STATISTICS 22.0 to perform the linear regression analysis. The pair of functional relations (22) and (23) are obtained, respectively. The coefficient of determination of linear fit or linear regression R 2 is 0.948, which is close to 1, indicating that the two have a good linear fit; the linear regression coefficient is significant at the confidence level of 0.05, which also indicates that the linear fit is reasonable. The standardized residual ZRE on the vertical axis and the operating income of the nursing home or the hospital on the horizontal axis are all around e = 0 random fluctuations, indicating that the fitting effect is good and the linear regression equation is satisfactory.
In this case study, α i = δ ji = 3.268v i /v j , α j = δ ij = 0.275v j /v i , and α i α j = 0.899 < 1. This paper calculated and analyzed the JAG hospital-nursing home symbiosis behavior pattern, as shown in Table 2.  Table 2 shows that the hospital-nursing home symbiosis behavior pattern from 2012 to 2019 basically belongs to the asymmetric reciprocal symbiosis pattern because of α i α j < 1. In 2012-2016 and 2019, nursing homes had a greater degree of beneficial influence on hospitals than vice versa and were more beneficial to the hospitals. In 2018, the hospitals had a greater degree of beneficial influence on nursing homes than vice versa and were more beneficial to the nursing homes, and in 2017, the degree of mutual influence between the two parties was basically the same, making the symbiotic relationship between the two similar to the symmetrical and mutually beneficial symbiosis. Expressing the dependence coefficient curves of the two, as shown in Figure 16, shows that the impact on nursing homes under the combined medical care model tends to be positive and shows an upward trend, and the dependence coefficient curves of the two tend to be close; that is, the two population density (income) evolutions tend to be synchronized. In addition, Figure 16 shows that under the asymmetric reciprocal symbiosis model from 2012 to 2019, the incomes of nursing homes and hospitals, and the total income of the medical-care integration, all showed an increasing trend, indicating that a good medical-care symbiosis relationship is beneficial to income growth, and that nursing home and hospitals tend to evolve simultaneously. an upward trend, and the dependence coefficient curves of the two tend to be close; that is, the two population density (income) evolutions tend to be synchronized. In addition, Figure 16 shows that under the asymmetric reciprocal symbiosis model from 2012 to 2019, the incomes of nursing homes and hospitals, and the total income of the medical-care integration, all showed an increasing trend, indicating that a good medical-care symbiosis relationship is beneficial to income growth, and that nursing home and hospitals tend to evolve simultaneously. Figure 16. Change of incomes under the dynamic evolution of the hospital-nursing home symbiotic relationship. Figure 16. Change of incomes under the dynamic evolution of the hospital-nursing home symbiotic relationship.

Conclusions
All participants or resource agents in the service ecosystem have their own independent market goal, which makes it difficult to integrate the goals in the ecosystem. Individuals driven by profit maximization excessively pursue short-term profit maximization, affecting the entire service ecosystem's realization of the goal of value co-creation. From the perspective of an evolutionary game, this paper uses the logistic growth model to build a model of the symbiotic evolution of multiple resources in the smart elderly care service ecosystem, uses computer simulations to dynamically simulate the model, and comparatively analyzes the population density of the hospital-nursing home symbiotic relationship using real cases in China.
The conclusions and suggestions are as follows: (1) The result of the symbiotic evolution between different resource agents in the smart elderly care service ecosystem depends on the symbiotic interdependence coefficient. (2) The enterprise goal in this ecosystem should be to maximize the overall value of the multiple resource agents. (3) When this ecosystem is viewed using a reciprocal symbiosis model, each resource agent benefits from the activities of the other resource agents, resulting in an increase in environmental capacity. (4) When the symbiotic model of the ecosystem is a competitive model, different resource agents are damaged by each other's activities, and the environmental capacity is less than the environmental capacity of independent development. (5) When the symbiotic model of the smart elderly care service ecosystem is a parasitic model, the population size of the parasitic resource agent is bigger than that of an independent development, and the size of the other resource agent is smaller than that of an independent development. (6) The strategic choice of resource agent has an important impact on the evolution of the value creation of the smart elderly care service ecosystem. When there is no revenue regulator in the system that can adjust the corresponding revenue distribution rules, the resource agent (in the interest of marketization) will take actions that benefit itself. Therefore, regulation is needed to ensure the relative fairness of the distribution of value co-creation in the smart elderly care service ecosystem, especially for resource agents that have short-term returns. (7) The case of JAG shows that it is reasonable to use the symbiotic evolution model constructed in this article to analyze the symbiotic relationship between the key resource agents of the smart elderly care service ecosystem.
The symbiosis of multiple resources is an inevitable product of the development of the smart elderly care service ecosystem, up to a certain stage. This ecosystem is a complex system of economic factors and social factors interacting together. Under the action of many factors, the evolution of this ecosystem appears as a complex behavior and trajectory. Therefore, in the next step of this research, we plan to study the following aspects: (1) We will use game theory and system dynamics theory to study the benefit allocation criteria of multiple resource agents in the smart elderly care service ecosystem, and verify the effectiveness of the model of dynamic symbiotic evolution constructed in this paper. (2) In the future, we will study the impact of strategy changes by introducing a mutation mechanism and then obtain a continuous replicator-dynamic mutant model. (3) The smart elderly care service ecosystem needs certain institutional support to maximize its system benefits; that is, the institutions and institutional arrangements in the smart elderly care service ecosystem are of great significance for studying the value co-creation and symbiosis of multiple resource agents. Social network analysis can be used in the future to study how the relationship between different social situations and reward and punishment mechanisms affect the behavior of multiple resource agents, and how the behavior of resource agents further affects the results of value co-creation.