Applications of a constrained mechanics methodology in economics
Abstract
The paper presents instructive interdisciplinary applications of constrained mechanics calculus in economics on a level appropriate for the undergraduate physics education. The aim of the paper is:

to meet the demand for illustrative examples suitable for presenting the background of the highly expanding research field of econophysics even on the undergraduate level and

to enable the students to understand deeper the principles and methods routinely used in mechanics by looking at the well known methodology from the different perspective of economics.
Two constrained dynamic economic problems are presented using the economic terminology in an intuitive way. First, the Phillips model of business cycle is presented as a system of forced oscillations and the general problem of two interacting economies is solved by the nonholonomic dynamics approach. Second, the CassKoopmansRamsey model of economical growth is solved as a variational problem with a velocity dependent constraint using the vakonomic approach. The specifics of the solution interpretation in economics compared to mechanics is discussed in detail, a discussion of the nonholonomic and vakonomic approaches to constrained problems in mechanics and economics is provided and an economic interpretation of the Lagrange multipliers (possibly surprising for the students of physics) is carefully explained.
The paper can be used by the undergraduate students of physics interested in interdisciplinary physics applications to get in touch with current scientific approach to economics based on a physical background or by university teachers as an attractive supplement to the classical mechanics lessons.
1 Introduction
The economic science has been influenced by physical concepts from its very beginning.^{1}^{1}1 Let us mention Daniel Bernoulli, as an example, who was the originator of utilitybased preferences, or one of the founders of neoclassical economic theory, Irving Fisher, who was originally trained in physics under J. W. Gibbs. . Although it is for contributions to physics and mathematics that Newton is celebrated, the Newtonian principles, formulated at the end of the seventeenth century, have powerfully influenced most branches of science  economics undoubtedly. Newton’s mechanics brought the doctrine of scientific determinism, the principle that all events are the inescapable results of preceding causes, for which (until the work of Planck and Einstein in twentieth century) the scientists tend to think of nature as a mechanical device whose behavior could be revealed by observation, experimentation, measurement and calculation. The idea of nature governed by natural laws dominated the new world order and many scholars have presumed human behavior and economics to be governed by such laws as well (see [1]).
As the economic thinking were developing through the centuries, the possibility that economic science can be inspired by physics was continuously debated. Nevertheless, the newly built physical theories and methodologies were repeatedly applied to economic problems by scholars trying to capture observed economic behavior. As a well developed branch of physics perfectly equipped with mathematical apparatus, mechanics has permanently served an inspiration for theoretical constructions in economics (for detailed discussion see [2]).
Looking for recently arisen intersections of physics and economics one arrives to econophysics which describes the phenomena of development and dynamics of economic systems by using a strictly physically motivated methodology. The official birth of the term ”econophysics” dates back to a paper by H. Stanley in 1996 [3]. Currently, mainly the applications of statistical physics and nonlinear dynamics are considered a core of econophysics (see [4]), but in broader context, econophysics can be considered an interdisciplinary research field applying theories and methods originally developed by physicists in order to solve problems in economics (for approach employing solely the classical mechanics methods see [5]).
Econophysics is a modern, quickly expanding interdisciplinary branch of science, which has already been transferred from the area of purely scholar interest into the real world including the establishment of graduate and postgraduate university study programs. Although econophysics has become a part of university physics education, the teaching aids are still being developed. Since the core of econophysics lessons requires passing advanced physics courses, very little is available for undergraduate students to satisfy their curiosity about this ”econophysics fashion”. In this paper we present instructive examples of how the methodology of constrained dynamic systems, commonly used in classical mechanics, can be used for solving economic problems and, in this way, we supply instruments to present and demonstrate this interdisciplinary topic on the undergraduate level of physics education. For recent educational contributions dealing with econophysics (but lacking concrete examples or inappropriate for undergraduates) see [6], [7], [8], [9]. Being intended for the students of physics without preceding knowledge of economic theories, the examples given in the paper are based on intuitive economic terminology and supplemented by easytoread explanations and interpretations of the economical background.
Apart from the illustrative potential of the examples, the look at standard methodology procedures used in mechanics (such as formulation of the equations of motion or handling the constrained systems) from a completely different point of view enables the reader to reach a better understanding of the physical background of these principles  a point which is often replaced by the calculation routine otherwise.
The paper can be used by undergraduate students of physics interested in solving interdisciplinary applications to get an idea about employing physical apparatus in economic problems. Moreover, the examples provided in this paper can be utilized by university teachers as an attractive supplement of traditional mechanics courses on the undergraduate level or as a motivation inviting the students to enroll for advanced interdisciplinary courses.
In Sec. 2 the classical Phillips model of business cycle is solved via the physical model of forced oscillations under friction: in Sec. 2.1 the economic model is briefly developed in an intuitive way, in Sec 2.2 the problem of two interacting economies is stated and an appropriate, physically motivated model of a nonholonomic system is used for its solution and in Sec. 2.3 we demonstrate that the results obtained are in good qualitative correspondence with the observed behavior of two economies. In Sec. 3, the CassKoopmansRamsey model of economic growth is investigated by the calculus of variations. It is shown, that the Lagrangian has the meaning of the overall (current value) utility in economy and the existing connection among the economic quantities is modeled by a velocitydependent constraint. The vakonomic approach is used for solving the system and the examples of a typical economic solution procedure and the interpretation of the solution are presented and discussed with respect to the standard procedures used in mechanics. A discussion of the difference between the nonholonomic and vakonomic approach together with the possibilities of their application in economics and mechanics is provided in Appendix 1. A special paragraph is devoted to an interesting economic interpretation of the Lagrange multipliers and an additional example and comments on this topic are given in Appendix 2.
2 Physically motivated business cycle description
Business cycle (or economic cycle) refers to economywide fluctuations in economic activity over several months or years. These fluctuations occur around a longterm growth trend, and typically involve periods of rapid economic growth (an expansion), and periods of relative stagnation or decline (a contraction or recession, see e.g. [10]). Business cycles exist in economy of any country or region (e.g. we can detect business cycles in Spanish economy or business cycles of European Union economy). Understanding how this cycles come into existence and what are their determinants, is crucial for making good economical policies in a particular country/region since the level of economic activity is linked with the standard of living in a country. Therefore mathematical models of business cycles form an important and intensively studied topic in the science of economics^{2}^{2}2Note that currently there is a number of different mathematical models of business cycle and still new are being developed. Generally in economic science, for each economic phenomenon (such as business cycle), there exist several competing approaches that use different assumptions and mathematical apparatus for building the model. Whatever the details of these approaches, all of them aim to follow the current and forecast the future behavior of the economic quantities as precisely as possible. Interesting from the physical point of view is that there is no unique ”correct” model in almost any economic discipline. There are many models for every phenomenon, each having its advantages and disadvantages and each corresponding with reality to an own extent..
The business cycle is described via the time development of gross domestic product () which refers to the market value of all final goods and services produced within a country in a given period. The term ”gross” refers to particular methodology of enumerating and in economic models we speak simply about product (denoting it ). Hence, the product provides us with the information on how much goods and services is produced (and it is assumed that this amount is also sold) in the economy (of a country/region) in each time period. When trying to express quantitatively, the economists face the problem how to find the representative function. Unlike in mechanics where the time development of the system is uniquely determined by forces and torques acting on the system, no such unequivocal approach exists in economics. Economists state that the characteristics playing a principal role in forming the business cycle are:

investments  the purchased amount of goods which are not consumed but are to be used for future production (e.g. purchases of new machines and buildings intended to be used in a production process),

consumption of households  the amount of money spent by households for goods and services,

total demand for goods and services  the amount of goods and services intended to be purchased in the economy. Let us explain more carefully the concept of total demand using the illustration of money/goods and services flows in an economy (see Fig. 1):
There are three ”players” in an ”economy game”: households that consume goods and services spending the money earned by working in firms, firms who are renting labor from workers to produce goods and services and a government that purchases both goods and services. Hence, the total purchases of goods and services (i.e. demand ) is composed of the household consumption, firm investments and government purchases of goods and services.

autonomous expenditures  expenditures independent from the total income in economy (total income is the amount of money earned by the individuals in economy), i.e. whatever the economic performance of the economy is, these money will be spent (e.g. government expenditures to keep the roads functional). Note that product (which is expressed in monetary units) is identified with total income in the economy^{3}^{3}3 The product is the total amount of goods and services produced in economy. But what was produced is to be sold and all the money are disbursed to the players in economy game (e.g. the employees are paid a wage, the material or semifinished products must be bought from firms  who must paid their employees etc. and the surpluses of corporations and entrepreneurs come as income to individuals in the economy). Hence, what was produced in economy transforms into the total income in economy..
2.1 The Phillips model of business cycle
In the Phillips model of business cycle (developed in 1954), which from the physicist’s point of view is interesting by the methodology used, the product function is deduced from the observed and presumed relationships between the above mentioned characteristics of the economy (see e.g. [11]):

The total demand in economy is composed of the consumption of households , ^{4}^{4}4The individuals divide their income between consumption (spending their money for goods and services) and saving (leaving the money in banks). So for the total income in economy (which can be denoted also since it has the same value as product) we can write , where savings and consumption are given by the income (product) and constant value , which is a given characteristic of an economy. , investments purchases and the autonomous expenditures (we will consider it to be a given constant of government spending):
(1) 
The demand in economy can be satisfied only via the product of the economy (since we do not take into account other economies). Hence, ideally the demand should be equal to supply (product) in any time. But it is assumed (and more realistic) that the product is reacting to the demand with a delay. For example if there is a demand for some fashionable consumer article and there is a lack of it in shops, the producers can react and increase its production, but this will take some time. Thus the supply of the article reflects the past demand. In this way also total product ”persuades” the total demand in time and, particularly in Phillips model it is assumed that the dynamics of product is given by
(2) where is a positive constant.

It is assumed that the potential investment is depending on the product change in time ( being constant), i.e. the investments are needed only if we want to increase production. For example if we want to produce more products in a factory, we will need more inventories, machines, etc. which generates the spendings falling within the ”investments” in the economic language. But the true investment is delayed from the potential one, i.e. in reality the investments do not react on the change in product immediately (in our example we can imagine, that we begin product more products using more extensively the machines we already have and the new machines will be bought later). The change of true investments in time is then supposed to be proportional to the difference , i.e.
(3) where denotes a positive constant.
Now the above mentioned equations describe the model of how business cycle come into existence: the product is delayed from demand, and investment purchases are delayed from their immediate need. These two factors generate in Phillips model the oscillations in time development of product ^{5}^{5}5The sources of the periodic oscillations observed in time development of product still are the matter of investigation. A modern theory of real business cycles suggests that the fluctuations of product can be to a large extent accounted for by real shocks (examples of such shocks include innovations, bad weather, quick oil price increase, stricter environmental regulations, etc.) which appear more or less periodically. The general gist is that something occurs that directly changes the decisions of workers and firms about what they buy and produce and thus eventually affect the product . Note then that Phillips model represents only one of many possible reasoning for the product oscillations.. Combining equations (13) and eliminating and , we obtain a secondorder differential equation for the unknown variable which represents the business cycle:
(4) 
where
are constant values. Equation (4) reminds us of forced oscillations under friction where is the square of the free oscillations frequency and is the external force acting on the system  generally it is a given function of time (this is the case also in Phillips model once the autonomous expenditure is not constant). In physics, the term , , represents the damping arising when the surrounding medium exerts a resistance. In an economic system this term may cause both damped and explosive oscillations depending on the sign of .
Solving the dynamical equation (4) for our economic problem, we can obtain harmonic oscillations, critical damping, damped or even explosive oscillation according to given parameters of the economic system. Given all the other parameters, by changing the constant (which can be found in each of the constants ) we can obtain all types mentioned of time development of the system. The solution of (4) for several values of can be found in Fig. 2^{6}^{6}6 The values of the parameters where chosen according to [11]: and the initial conditions were set as , ..
Note that we do not mention concrete units in the graphs, since the results are derived for artificial parameters and for the purpose of demonstrating the behavior of the product . Visualizing ”nonspecified” graphs is a standard means of qualitative presentation of economic behavior. But if needed, one can scale the product in thousands of USD and the time axis in years.
The manner in which the product is developing in time (harmonic, explosive, damped oscillations or the critical damping) depends on the particular choice of parameters . Given the rest of parameters, the intervals of values corresponding to each possible solution in Fig. 2 can be derived from a general solution of the dynamical equation (4). Analyzing these intervals, the economists conclude that oscillations, will occur if the time lag of product reaction to the demand does not differ significantly from the time lag of the reaction of investment to product change. The damped (explosive) oscillations will occur if the time lag of production to the demand will be significantly higher (lower) than the time lag of the reaction of investments to product change. These explanations following from the mathematical model answer the question of how the business cycle come to existence only to some extent. As we have mentioned there are many other approaches to business cycle modeling and always there may be a debate on the explanations and reasoning of each particular model.
It is worth noting that the final ”equation of motion” (4) for the economic system (we will call it a dynamical equation when talking about economic systems) was obtained based on assumptions (13) about the dynamic relationships among the economic quantities. But there is no mention about the Lagrangian of the economic system^{7}^{7}7Note that in next paragraph we will see that the variational concept has also been adopted by quantitative economists in other models. so as to generate the dynamic equation (4), nor are the terms in (4) interpreted as some economic ”forces”. Thus, although the resulting dynamical equation is well known in physics, the method of obtaining it lacks the systematic spirit of physical methodology (remember Lagrangian approach or Newton equations in mechanics.)
Hence, through certain assumptions Phillips model arrived to the description of oscillating behavior of product by the differential equation of second order. The ideal result of this model would be that once known the constants of the economy () one could predict the future oscillations of . Practically, the political authorities would prefer only steady economical growth (i.e growth of the product), which generates low unemployment, increasing wages and overall increasing life standard of people. Oscillations around this simple growth trend are disturbing and at least knowing how they come into existence or better, what they will look like next period, would bring important support to the policy makers. But as we have mentioned before, such predictions and even the understanding of the oscillation phenomena are not completed and are the subject of ongoing research.
Let us conclude the description of Phillips model of business cycle by the note that using the wellknown and developed physical model (the harmonic oscillations) to describe observed economic behavior (product time development) was a natural initial approach in economics. Although the current philosophy of quantitative economics diverges fundamentally from simply adopting the existing physical models for the description of economic systems, the physical methodology itself is still in the focus of economists. Accordingly, the harmonicoscillationsbased methodology is still a common approach used for the description of economic systems nowadays (for a recent application see e.g.[12]). No matter how historical and simple the approach is, the Phillips model is able to produce at least qualitatively correct solutions (in the sense of correspondence with the observed behavior) as we shall see in Sec. 2.2.
2.2 Two interacting economies
In the preceding model we considered an isolated economy which had no contact with any other economy. In reality, economies are interacting through mutual purchases and money transfers. As typical in economics, a number of models were developed for this more realistic assumption. The most common is the model where, when describing a single economy, we introduce another ”player” (apart from households, firms and government) representing external economies influencing the economy modeled. But let us continue to develop a different, possibly more general, model employing the physical description initialized in the preceding section.
Let us consider economies of two countries which are ”interacting” through purchases and money transfers, i.e. part of the goods and services in an economy is purchased in the another. We can assume that the product in one economy will be affected besides other things also by the product of the other economy (for example if in Germany decreases the total product, i.e. there is less production of goods and services, then the demand of German firms for the subproducts and services imported from Czech Republic decreases and this will naturally affect negatively the product of Czech Republic. And, reciprocally, the less product in Czech Republic means lower income of households and lower consumption of all goods and services, among all those more expensive e.g. those imported from Germany, so the lower demand in Czech Republic will affect to some extent the German product).
Assume that, if there was no such interaction between the economies, each of them could be described in terms of equation (4). Hence, we can describe two economies without interaction as a system with two degrees of freedom using the dynamical equations:
(5)  
(6) 
with constants
where labels the economies. Following the physically motivated design of the Phillips model, we will formulate the interaction between the economies as a constraint binding the two coordinates (). In physics, the constraints are used if the real forces ensuring the observed behavior are not known or are uneasy to quantify. In a constrained system then, the constraint forces arise which ensure the prescribed behavior and have the meaning of real physical forces acting on the system (remember e.g. the rolling of the cylinder without slipping in mechanics where the constraint forces have the meaning of forces and torques stemming from interaction of the cylinder with underlay, see [13] or [14]). Although we have already mentioned that no ”forces” are defined in economics, the constraints seem to be an appropriate methodology since they make it possible to describe the really observed behavior whose ”dynamical generator” is not known. The relationship between the two economies can be stated generally as a nonholonomic constraint. For our artificial example, let us assume this constraint to be
(7) 
where are constants. Since the constraint is linear in velocities and the model itself is physically motivated, we will treat the constrained system in a way common in mechanics  we will use the nonholonomic approach (see e.g. [15] or [16]). Then, the dynamical equations for our constrained system take the form:
(8)  
(9) 
which, together with constraint (7), yields 3 equations for 3 unknown variables , where is the Lagrange multiplier. Eliminating the Lagrange multiplier from (89) and substituting from the constraint (7), we obtain reduced equations of the constrained system:
(10)  
(11) 
where
2.3 Results and discussion
Solving the constrained system (57) for three different sets of parameters in (7): Constraint A: Constraint B: Constraint C:
we obtain the solutions in Fig. 3^{8}^{8}8The parameters used for our example were chosen according to [11] . Hence, the economies are quite similar except the autonomous expenditure . The initial conditions for our problem were the following: , , .. The solutions provide qualitative information about product dynamics for two economies influencing each other. The particular time behavior of products obviously depends on the particular form of the constraint. Note that, for realistic modeling, the particular coefficients should be estimated from observed data. Since the qualitative results of the test case are for guidance only, we can state that, although being simple and based on an (oldfashioned) classical physically motivated model, the solution is in a qualitative correspondence with the observed behavior of the true interacting economies.
At Fig. 4 we can see the GDP time series of Czech Republic and Germany in the period 20072010. As we have already mentioned there is a considerable trading between these two economies which possibly make a connection between the time development of their products as qualitatively visible on Fig. 4. This connection could be qualitatively expressed by the constraint (7) and the model of two intercating economies (1011) could enable to describe the business cycles of the two economies.
Thus, we have presented an example of how an economic model can be designed using the physical motivation. First the classical model of product dynamics based on the mechanical model of an oscillator was employed to describe an isolated economy. Then this model was extended to a model of two interacting economies using an oscillatorbased approach, i.e. we treat each of the economies as a single component of the system with two degrees of freedom. The existing observed relationship between the economies was modeled by a constraint and the qualitative results of the model when compared to the real data (see Fig. 4) are promising. In economic science, such a model design could be successful after careful validation based on a concrete parameter estimation. As in other models, probably only some aspects of the observed behavior could be captured by the newly developed model and, therefore, a debate about the relevance of the assumptions made and the methodology used might be sparked. Nevertheless, such debates take place for any economic model no matter whether with a physical or other background. An important fact is that, so far, hardly any model has been validated in economics in the sense of physics, since there always can be found significant discrepancies between the data predicted by the model and those observed in the real economy world. Indeed, economic science is still trying to find the appropriate theories, methods and even the calculus to achieve the status of an unchallenged source of scholar knowledge.
3 The economic growth model using calculus of variations
While the short run variation in product (measured by gross domestic product GDP per capita ) is usually termed business cycle, the long run increase of per capita GDP is called economic growth. Economic growth is primarily driven by improvements in productivity, which involves producing more goods and services with the same inputs of labor, capital, energy and materials (therefore there is a lasting demand for innovations and technological improvements in developed economies). The longrun path of economic growth is one of the central questions of economics. An increase in GDP of a country greater than population growth is generally taken as an increase in the standard of living of its inhabitants, hence the aim of the policy makers is to achieve steady economic growth. To understand better the phenomena of economic growth and to obtain a quantitative support for making appropriate decisions a number of mathematical models of economic growth was developed during past fifty years. Again as we have already mentioned when generally speaking about economic models all models are representing the economic reality to individual extent and they differ in the assumptions made about how the economic characteristics are mutually interacting and influencing the final growth.
In this section, the CassKoopmansRamsey model of economic growth based on variational approach will be presented and the methodology and motivations behind an economic variational approach in comparison to traditional usage of variational calculus in mechanics will be investigated. Since the aim of this section is not to provide the reader with a thorough economical background of the model, we will describe the philosophy of the model intuitively (for more detailed notes about particular assumptions and economical background see e.g. [18]).
3.1 The CassKoopmansRamsey model
This time, there are two ”players” in the ”economy game”: firms and households (so the flows are similar to the previous model illustrated in Fig. 1 except the missing government purchases). The aim of the model is to answer the question of how much the product of economy at any point of time should be spent for immediate consumption to yield current utility, and how much of it should be saved (and invested)^{2}^{2}2In economics the product (which is considered to be equal to the income) can either be consumed or saved, but what is saved results in investment. For example we can imagine, that the household saves some amount of money in a bank. But the bank uses the money to provide the firms with loans and the firms use the extra money for investment purchases. so as to enhance future production and consumption and, hence, yield future utility. The criterion of optimality is the social welfare which is given by social utility . There are two timedependent variables through which the desired optimal development of the economy can be attained: consumption and capital (measured in units advantageous for the model purposes, for details see [18]). Consumption represents the purchases of goods and services by households while capital is a factor of production, used to produce goods or services, that is not itself significantly consumed (though it may depreciate) in the production process (e.g. machinery, buildings, vehicles). Typically, in economics there is considered a relationship between the investment and the time derivative of the capital. In CassKoopmansRamsey model, the particular relationship takes the form
(12) 
where denotes the product (more precisely it is a known production function prescribing how the total production in economy is dependent on the capital available). Expression has the meaning of investment^{4}^{4}4This expression is connected with relationship used in the previous example (see item (i)). Remember, that it means that the total income (or product, which is equivalent as we have already discussed) in economy is distributed among consumption and savings. Since what is saved is supposed to be invested, the savings define the investments in monetary units. Thus investments savings = product consumption. and is the amount of investment that must be done just to keep at its existing level ( and being constants). The equation (12) states that the change in capital is equal to the investments less the replacement purchases (i.e. if we buy new machinery only to replace the broken ones then, even we made investments, we have not increased capital we still have the same production factors for our production process as before the investment purchases.)
The task in CassKoopmansRamsey growth model is to maximize social utility (which depends only on consumption)^{3}^{3}3It is typically assumed in economics that the utility of the society is given by the consumption of the people only. This assumption can be put into question and actually, several authors did so. It is obvious that if the criterion of optimality in economy is given simply by consumption, then many essential factors considering e.g. natural conditions can be omitted. via time paths of variables and under the constraint (12). For the particular form of the utility function defined in Ramsey model (see [17] or [18]), we can write the overall optimization problem as:
(14)  
where are positive constants with . The form of the consumption function is the matter of intensive investigation and the particular form used in the CassKoopmansRamsey model represents just one approach. The term ensures discounting of economical quantities. It means that the flows of particular quantity (social utility in our problem) to come have smaller weight that its current flow. Note that this approach was doubted by many authors because by this model the priority is given to (our) current consumption at the expense of the consumption of future generations.
From the mathematical point of view this is a constrained variational problem with a velocitydependent constraint (14). Solving (1414), we obtain the optimal growth path, in other words, we realize how and must behave over time to achieve the maximal lifetime utility from the consumption flow given the relation (14) between the product, consumption and investment. Note that in (14) could be substituted from (14) to obtain an unconstrained variational problem with the Lagrangian dependent on both variable and its first derivative . Because the form (1414) of the growth model is typically used in economics (thanks to its better economic interpretation and reasoning), we will not ”simplify” the initial problem (1414) and it will be treated as a constrained variational problem for the solution of which the vakonomic approach will be employed. The Lagrangian for the constrained problem (1414) takes the vakonomic form
(15) 
were is the Lagrange multiplier. Considering three variables instead of the initial pair , we obtain the variational dynamical equations:
(16)  
(17)  
(18) 
Note that the vakonomic approach accounts for the main technique used to solve variational constraint dynamical problems in economics (see [17],[18],[19] for typical examples). For the description and some notes on the comparison of vakonomic and nonholonomic approach see Appendix 1.
3.2 Results and Discussion
Expressing from (17) and substituting it into equation (16), we arrive after some calculation to a system of two equations for two unknown variables :
(19)  
(20) 
Then, it is not difficult to find a solution to our problem, i.e. to find an optimal consumption and capital paths. Nevertheless, in economics the functions are not normally given explicitly, but only assumed to have characteristic qualitative properties (such as that the utility is an increasing non convex function of consumption). Thus the problems have an air of ”theory” rather than ”computation” and reaching a particular solution is neither possible nor necessary. Yet, the content of the problems is meaningful and analyzing the qualitative characteristics of the solutions often generates important insights into economic behavior. In models of economic dynamic optimization, twovariable diagrams are prevalently employed to obtain these qualitative analytical results.
The behavior of an economy in the CassKoopmansRamsey problem is typically described in terms of the evolution of and using a diagram at Fig. 5. This diagram reflects the dynamical equations (1920) where the arrows show the directions of change of both and . In mechanics, we can plot similar diagrams but this is mostly done for a particular solution. In such a case the trajectory of the system in the configuration space is obtained. The diagram in Fig. 5 is used to identify the qualitative optimal behavior of an economy given the initial conditions and additional requirements (such as ). Then, the result consists in a verbal description of the optimal behavior rather then in an analytical solution. As an example, let us discuss the top left ”quadrant” of the graph. The consumption in economy is high and rising and eventually reaches zero. When continues to rise must become negative. But this cannot occur. Since the product is zero when is zero, the consumption must drop to zero. Therefore, such paths can be ruled out from consideration about realistic behavior of an economy. In a similar way, other ”quadrants” in Fig. 5 can be analyzed and we would arrive to a conclusion that the possible time path of the system is the one driving into the point where . Then the economy is said to ”move along the saddle path” to the equilibrium in . The equilibrium is Paretoefficient which means that it is impossible to make anyone better off without making someone else worse off. The explanation lies in the fact that all the households in the model have the same utility. The equilibrium produces the highest possible utility among allocations of capital and consumption that treat all households in the same way. Once the economy arrives to this equilibrium it can not do better and therefore we do not change the level of capital and consumption anymore, i.e. (for details see [18]).
We can see clearly the difference between the application of variational calculus in mechanics where we are searching for a particular solution which prescribes the time evolution of the system precisely, and the economic application, where the behavior of the system is only qualitatively analyzed.
3.3 Lagrange multiplier interpretation
A first look at a nontrivial constrained optimization economic problem may be surprising for the student of physics. Leaving aside the unusual definition of feasible region in economic optimization problem (mostly restricted to nonnegative values of variables), the economic interpretation of the Lagrange multiplier is worth a discussion. Mathematical texts provide no interpretation of the Lagrange multiplier , leaving the student with the impression that has no significance beyond providing an extra variable which transforms a constrained problem into an unconstrained, higher dimensional one. But in economic problems, the Lagrange multiplier can usually be interpreted as the rate of change of optimal value of the criterion function relative to some parameter.
The constraint (14) describes the capital accumulation and thus should correspond to the value of having a tiny bit more capital (see Appendix 2 for simple explanatory example). Hence, we can find the meaning of the Lagrange multiplier by answering the question of how the unit change in capital available in economics will affect the utility in optimum. Let us derive the influence of change in capital onto a ”present value” instantaneous utility (see (14)):
(21) 
From the first variational dynamical equation (16) we obtain
(22) 
and differentiating the constraint (18) we get
(23) 
Substituing into relation (21) from (2223) and having on mind that , we obtain
(24) 
That is, the extra unit of capital will raise the flow of output by an amount each unit of which (without the unit replacement production needed for keeping the capital at existing level and thus not intended for consumption) has a utility value of . Multiplier is referred to as a shadow value of capital (evaluated in present value). New economical interpretation arises then from variational equation (16) and (22), respectively: It states that the marginal utility of consumption^{9}^{9}9utility gained (or lost) from an increase (or decrease) in the consumption is equal to the shadow value of capital, . Thus, at the optimum the consumer is indifferent between consuming an additional unit and investing it. Remember, that the consumer must decide whether to consume or to save (which directly generates the investments) additional unit (of income). If it is consumed it brings direct utility, if it is invested (into capital) then it increases the product and future consumption. If the marginal utility of consumption was larger than the shadow value of capital,
(25) 
then capital would be too high and consuming more and saving less would increase the utility. Similarly, if the marginal utility of consumption was lower that the shadow value of capital
(26) 
then capital is too low and the households should save more to increase utility through higher future consumption. This example of Lagrange multiplier interpretation together with the subsequent discussion of dynamic equations represent the typical approach of optimization problem analysis in economics.
In mechanics, no interpretation is used for multipliers in optimization problems, but remember that the multipliers arise in the constrained forces, which are important from the physical point of view. The concept of constrained forces, on the other hand, has no significance in economic problems although the multipliers do. Hence, although the computational routine remains similar, there are distinct interpretations of the tools used.
4 Conclusion
Classical mechanics has played a significant role in the development of economic thinking influencing it in both principles and the calculus. Currently, different attitudes to the adoption of physical methodology for economic purposes can be distinguished. A strong one supports the philosophy that applying the well developed physical methods could quickly provide the economists with working and applicable models. This strong motion in economic science is reflected by a newly arisen term ”econophysics” which currently designates one of the possible fields where the undergraduates and graduates of physics could be involved.
The examples presented in the paper enable the undergraduates to meet for the first time this highly modern and progressively expanding field, and to get in touch with different philosophy of building scholar knowledge in economics. These examples may serve as a motivation for students to further study the economics  mechanics (or more generally physics) intersections. In addition, studying these problems makes the undergraduate students of physics face the different usage of a known methodology: the constraint variational calculus appears to be a flexible approach which in each discipline provides specific information and the nonholonomic dynamics of systems, based on constraint forces, could be applied even outside mechanics.
In physics education, the economic examples presented can serve not only as a demonstration of interdisciplinary applications of methods typically used in mechanics, but can also provide the teachers with an aid for demonstrating what are general mathematical and what are specifically physical features of the mechanics methodology. In this way, the students can achieve a deeper insight into the physical background of what they have learned in classical mechanics courses.
(27)  
(28) 
where are the Lagrange multipliers. The nonholonomic approach is typically used for solving the mechanical systems with nonholonomic constraints although in mechanics the nonholonomic constraints are frequently affine in velocities. The form of constraint forces for classical ideal constraints is a consequence of d’Alemberts principle. Note, that the nonholonomic approach has been extended also to more general constraints by the work of Chetaev [20] and others (see e.g. [21]; for instructive nonholonomic mechanical problems see e.g. [13], [16], [14]).
The mathematical concept of nonholonomic dynamics had been kept essentially unchanged until 30 years ago when a new dynamics of velocity constrained mechanics system was introduced by Kozlov [22]. This new mechanics was called vakonomic being ”variational axiomatic kind”. If we adopt a variational approach by requiring the motion to be a stationary curve of the action functional among all curves having the same end points and satisfying the nonholonomic constraints, then we get a vakonomic motion, i.e. we search for the solution of (unconstrained) variational problem associated to the Lagrangian function
Namely, the vakonomic motions can be obtained by the Lagrange equations
(29)  
(30) 
which give
(31)  
(32) 
Comparing the nonholonomic (27) and vakonomic equations (31) we arrive to a conclusion that the arising systems of differential equations are not equivalent unless
The author himself in [22] says: ”…vakonomic dynamics, which is an internally consistent model that can be applied to description of the motion of any mechanical systems, is as ”true” as traditional nonholonomic mechanics. The issue of the choice of model for each particular case is ultimately resolved by experiment.” While there are doubts that vakonomic dynamics is a satisfactory model for nonholonomic systems in mechanics (see e.g. [15]), in economics vakonomic approach is the one typically used for the solution of dynamical optimization problems with velocity dependent constraints. Intuitively, the nonholonomic concept based on the definition of constrained forces developed for mechanics could not be valid for the economic problems unless the background of the model is physically motivated. The vakonomic approach, if chosen for solving the constrained systems in economics, can be justified simply by the results it provides (for more examples on vakonomic approach in economics see e.g. [23], [24]).
Appendix 2
In constrained optimization in economics, the value of the Lagrange multiplier at the optimal solution is referred to as a shadow price. Shadow price is the change in the objective value of the optimal solution obtained by relaxing the constraint by one unit.
Consider a simple (static) business optimization problem of profit maximization in a firm. Assume that there are only two products and the firm is deciding about the amount of these two products to be produced in next period. The profit and number of working hours needed per unit of each product are known and denoted by and , respectively. The operating time limit for next period is (for example, we have B working hours available on certain machine for next month). We get the linear programming problem:
(34)  
We introduce a multiplier and form the Lagrangian
Assuming that the firm maximizes profit (given the constraint (34)), the optimal quantities and the multiplier necessarily satisfy the firstorder conditions:
(35)  
(36)  
(37) 
and differentiation of constraint (34) yields
(38) 
Now, using the chain rule and equations (35), (36), (38) we obtain
(39) 
Hence, the Lagrange multiplier measures how the total profit responds to the unit change in total machine operating time available. The shadow price is here the maximum price the manager would be willing to pay for operating the production line for an additional unit of time, based on the benefits he would get from this change. For example, if Euro, then running the machine and producing for an additional hour will gain the profit of Euro and the manager should not pay more for this additional working hour than is the value it produces (for more see e.g. [25]).
References
 [1] E. Ray Canterbery, A Brief History of Economics, World Scientific Publishing (2002).
 [2] GrattanGuinness I, ”How influential was mechanics in the developement of neoclassical economics? A small example of a large question”, Journal of the History of Economic Thought 32, 531–581 (2010).
 [3] H. Stanley et al., ”Anomalous fluctuations in the dynamics of complex systems: From DNA and physiology to econophysics” Physica A 224, 302–321 (1996).
 [4] M. Schulz, Statistical physics and economics, SpringerVerlag New York, Inc., (2003).
 [5] I. Kitov, O. Kitov, Mechanomics: Economics as Classical Mechanics, Lap Lambert Academic Publishing (2010).
 [6] Ch. Schinckus, ”Econophysics and economics: Sister discipline?”, Am. J. Phys.78, 325–327 (2010).
 [7] A. Walstad, ” Comment on Econophysics and economics: Sister discipline?”, Am. J. Phys.78, 789–789 (2010).
 [8] D. Stauffer, ”Teaching computational biosocioeconophysics” Eur. J. Phys. 26, S79–S84 (2005).
 [9] P. Mostardinha, E.J. Durana, F.V. de Abereu, ”The econophysics in the Euromillions lottery”, Eur. J. Phys. 27, 675–684 (2006).
 [10] A. O’Sullivan, S. M. Sheffrin, Economics: Principles in action. Upper Saddle River, Pearson Prentice Hall (2003).
 [11] R.G.D. Allen, Mathematical economics, ST Martin’S Press INC, New York (1963).
 [12] A. Ataullah, M.Tippett, ”Equity prices as a simple harmonic oscillator with noise”, Physica A. 382, 557–564 (2007).
 [13] J. Janová, J. Musilová, J. Bartoš, ”Coupled rolling motion: a student project in nonholonomic mechanics,” Eur. J. Phys. 30, 1257–1269 (2009).
 [14] J.Janová, J. Musilová, ”The streetboard rider: an appealing problem in nonholonomic mechanics”, Eur. J. Phys. 31, 333–345 (2010).
 [15] G. Zampieri, ”Nonholonomic versus Vakonomic Dynamics”, J. Diff. Eq. 163, 335–347 (2000).
 [16] J. Janová, J. Musilová, ”Coupled rolling motion: considering rolling friction in nonholonomic mechanics,” Eur. J. Phys. 32, 245–257 (2011).
 [17] A.C. Chiang, Elements of dynamic optimization, Waveland Press, Inc. (2000).
 [18] D. Romer, Advanced macroeconomics, McGrawHill, (2006).
 [19] P. Aghion, P. Howitt, Endogenous growth theory, The MIT Press, Cambridge (Mass) (1998).
 [20] N.G. Chetaev ”On Gauss Principle” (in russian ), Izv. Kazan. Fiz. Mat. Obs.6, 323–326 (1932).
 [21] M. de Leon, D.M. de Diego, ”On the geometry of nonholonomic Lagrangian systems”, J. Math. Phys. 37, 1–26 (1996).
 [22] V. V. Kozlov, ”Dynamics of systems with nonintegrable constraints”, Part I, Moscow. Univ. Mech. Bull. 37 (1982), 2734; Part II, 7480, Part III 380 , 40–51 (1983).
 [23] A. Grubbstrom, ”Application of Calculus of Variations to Financing Alternatives”, Omega 19, 305–316 (1990).
 [24] A. Molinder, ”Application of calculus of variations to a continuous time aggregate production model”, Int. J. Prod. Econ. 41, 273–280 (1994).
 [25] J.V. Baxley, ”Lagrange Multiplier Problems in Economics”, Am. Math. Mon. 91, 404–412 (1984).