Exact and inexact differentials in the early development of mechanics and thermodynamics

We give an account and a critical analysis of the use of exact and inexact differentials in the early development of mechanics and thermodynamics, and the emergence of differential calculus and how it was applied to solve some mechanical problems, such as those related to the cycloidal pendulum. The Lagrange equations of motions are presented in the form they were originally obtained in terms of differentials from the principle of virtual work. The derivation of the conservation of energy in differential form as obtained originally by Clausius from the equivalence of heat and work is also examined.


Introduction
It is usual to formulate the basic equations of thermodynamics in terms of differentials. The conservation of energy is written as where dU , dQ, and dW are the differentials of the internal energy, the heat absorbed and the work performed by a system. The conservation of energy is expressed by equation (1) as long as dU is an exact differential, a key feature that is not explicitly stated but is tacitly understood. An exact differential such as dU means that there exists a state function U such that its differential is dU . An inexact differential such as dQ and dW , does not hold this property. It should be pointed out that some authors use a distinction notation to refer do to an inexact differential. One of them is to cross the differential, a procedure that becomes quite confusing and useless. The differential of the work dW is written in terms of exact differential. In the case of mechanical work, for instance, dW = pdV where p is the pressure and dV is the differential of the volume V . In an analogous way the differential of heat can also be written in terms of an exact differential. If a system is in thermodynamic equilibrium than dQ/T is the exact differential dS of the entropy S. This result allows us to write the relation dQ = T dS between heat and entropy, which should be understood as valid as long as the system is in thermodynamic equilibrium.
Differentials appeared with the invention of calculus in the second half of the seventeenth century. The differential dx was understood as a small increment of the * Correspondence email address: oliveira@if.usp.br. variable x. If another variable y depends on the independent variable x, then the resulting increment dy of y is its differential. The quotient of these two differentials, dy/dx, was interpreted geometrically by Leibniz as the ratio of the ordinate y of a point on a curve and the length of the subtangent associated to this point.
The formulation of thermodynamics in terms of differentials seemed to be a natural consequence of the application of differential calculus at the time of the emergence of thermodynamics around the middle of the eighteenth century [1]. Clausius formulated the laws of thermodynamics by the use of differentials [2]. The same can be said about the development of mechanics in the eighteenth century [3], particularly the formulation of analytical mechanics by Lagrange [4].
Here, we give an account and an analysis of the development of the concept of differential and its application during the early development of mechanics and thermodynamics. Our presentation starts with the introduction of the differential and integral calculus by Newton and Leibniz. Then we analyze the developments of the basic rules of differential calculus and the concepts of integrating factor, exact and inexact differentials. Next we show how the application of differential calculus was used to solve mechanical problems such as the brachistochone problem. In the second part, we analyze how the differentials were used by Lagrange in his analytical mechanics and by Clausius in his development of thermodynamics.
As a generic procedure of our exposition we try to use the terminology and notation originally employed by the authors. However, for the sake of a better understanding by the reader, sometimes we use the modern terminology and notation as long as they do not modify the original meaning of the concept being described.

Basic rules
The differential and integral calculus was invented independently by Newton and Leibniz in the second half of the seventeenth century [5][6][7][8][9][10]. They provided algorithm procedures, including the basic rules of the calculus, which were essential for the development of the concept of derivative and integral [5]. Newton developed the calculus in the years 1665-1666 but his first publication on the subject, which was written in 1676, appeared only in 1704 as a mathematical appendix to his Opticks, although portions of his Principia, published in 1687, contained some of his achievements on calculus [11]. Leibniz developed the subject during the years 1673-1676 and the accounts of his findings were published in two papers [8]. The first, in 1684, on differential calculus [12], the second, in 1686, on integral calculus [13]. Leibniz also created a convenient notation, without which the calculus could not achieve its fundamental role in mathematics [7].
The differential calculus may have its origin in geometric problems such as that of drawing tangents or of locating the maxima and minima of curves [14]. Its origins is found within the domain of analytic geometry where geometry is studied by means of a coordinate system. A curve drawn in a plane is represented in analytic geometry by an equation involving two algebraic variable, the abscissa and the ordinate, which in Cartesian geometry are understood as the distances of a point on the curve from two orthogonal axes. In his first publication on the differential calculus of 1684 [12], Leibniz denoted the differential of a variable x by dx, the same notation that is used today. It was interpreted as a small increment in the abscissa x, which implies an increment, or a differential, dy in the ordinate y, as represented geometrically in figure 1. Leibniz states that the ratio between dy and dx equals the ratio between the ordinate and the subtangent.
The development of Leibniz calculus during the first decades after its discovery was accomplished mainly by Jacob Bernoulli and Johann Bernoulli, and the first book on the subject was published in 1696 by l'Hôpital [15], based on the works of Leibniz, Jacob Bernoulli, and specially on those of Johann Bernoulli. It contained the basic rules of the differential calculus, and included procedures for determining tangents, maxima and minima, inflexion points, involutes and evolutes, and caustics of plane curves.
L'Hôpital starts his exposition on differential calculus by giving the basic rules of the subject, such as the differential of the sum of two variables z = x + y, of the product of two variables z = xy, of the quotient of two variables z = x/y, where m is a rational number, positive or negative. The determination of the tangents of curves, requires the calculation of the subtangent, the length of the segment TP of figure 1. If Pp (dx) is small, then Rm (dy) is also small, the triangles MRm and TPM become similar and as a consequence the length t of the subangent TP is to y as dx is to dy, The right-hand side is understood as the quotient between the two differentials. In fact, this is the expression used by Leibniz to define dy by introducing dx as an arbitrary finite interval [12]. To find the maximum value of the ordinate of a curve it suffices to set the expression of dy obtained for a given curve equal to zero. L'Hôpital gives the following example, which described a certain plane curve. According to the basic rules, Setting this expression equal to zero he finds y = 3x 2 /a which replaced in equation (7), gives x = (a/3) 3 √ 2, the value of x at which the curve has its maximum value.
The determination of inflexion points of curves was also considered by l'Hôpital. In this case it is necessary to find the second differential d 2 y of the ordinate y, understood as the differential of the differential dy and thus holding the same rules of the first differential. As an example, he considers the curve described by From this expression, one finds The second differential is also determined by the basic rules, yielding Setting the numerator equal to zero, he finds the values x = a/ √ 3 for the location of the inflexion point. The corresponding ordinate is y = a/4.

Integrating factor
The differential calculus was by its origin a natural method for solving geometric problems, the results being expressed in terms of algebraic expressions [5]. In this sense it was a method to solve problems in analytic geometry. This changed with Euler, who made the subject a theory of functions without the need of a connection with geometry [5]. Notice however that, for Euler the meaning of function was a formal expression in terms of variables and constants and did not have the present meaning of a conceptual relationship between a quantity and an independent variable [5].
In a paper presented in 1728, published in 1732 [16], concerning a method to reduce second order into first order differential equations, Euler introduces the integrating factor. As an example, he considers the following differential equation To solve equation (12), Euler multiplies the equation by the integrating factor (t − 1) 2 , obtaining The last term is the differential of t − ln t, and the first two terms is identified as the differential of (t − 1) 2 z. Thus the integration gives where a is a constant.
Equation (12) is of the general type where P and Q are functions of t and can be solved in like manner. According to Euler, the integrating factor is Multiplying equation (15) by x, and taking into account that dx = xP dt, one finds which after integration gives This equation gives z as a function of t.

Complete differential
In another paper, presented in 1734 and published in 1740 [17], Euler considers a function V of two variables x and y. He writes the differential of V as where the first term is obtained by considering y constant and the second by considering x constant. Writing dP = pdx + rdy and dQ = qdx + sdy, Euler argues that q = r [17]. Later, he introduces the notation [18] and writes the equality q = r as Euler warns that not all expressions of the type (19) are differential of a certain function of x and y. If the relation (21) does not hold, no function of x and y exists such that its differential is P dx + Qdy. An example given by Euler is yxdx + x 2 dy.
In a paper published in 1742 [19], Clairaut shows that, if Adx + Bdy represents the differential of a quantity dependent on x and y, then which is the condition (21) found by Euler. According to Clairaut, this condition allows us to known if Adx + Bdy is the differential of a certain quantity in two variables. If the condition is fulfilled, Clairaut calls Adx + Bdy a complete differential. As an example of a complete differential he gives In his book on the figure of the earth published in 1743 [20], Clairaut states that the integral of a complete differential is a function of x and y and he gives the following two examples ydx + xdy, ydx + xdy 2 a 2 + xy , (24) which have as integrals xy and a 2 + xy, respectively.
To integrate a complete differential Adx+Bdy, Clairaut uses the following scheme [19]. One integrates one of the members, say Adx supposing y constant, and to the result one adds a quantity C that depends only on y. Next one calculates the differential of the sum, keeping x constant, and subtract Bdy. The result is either zero or is the differential of C, which can thus be integrated and added to Adx.
If the expression M dx + N dy is not a complete differential, Clairaut explains that it is possible to render a complete differential by multiplying the expression by a factor µ dependent on x and y [19]. Assuming that µM dx + µN dy is a complete differential, it follows from the condition (22) that an equation to be solved for finding µ.

Derivative
The differential and integral calculus as invented by Newton and Leibniz and developed by their immediate successors was very successful, becoming the mathematical basis of physical theories such was the case of Newtonian mechanics. However, criticisms were raised concerning its foundations particularly the employment of infinitely small quantities [5] as occurs in the differential quotient dy/dx. If the differential is a vanishing quantity then the differential quotient would result in an indetermination ratio of the type zero over zero. On the other hand, the rules of the calculus established by Newton and Leibniz give well defined results for the differential quotient. A crucial contribution toward the clarification of this problem was presented by d'Alembert in his article on differential contained in the Encyclopédie, where he explains how to find the differential quotient by using the concept of limit [21]. As we have seen above the differential quotient is equal to the ratio of the ordinate MP of a point M of a curve and the subtangent TP as shown in figure 1. D'Alembert shows this result by considering a straight line QMm, which crosses the curve AMB at the point M. As the point m approaches M, the point Q moves toward the point T and as a consequence the ratio mR/MR, which equals the ratio MP/QT, aproaches the ratio MP/TP of the ordinate and the subtangent. D'Alembert gives the example of a curve described in algebraic terms by ax = y 2 . Denoting mR by z and MR by u, the ratio mR/MR is a/(2y + z) which, as z decreases, approaches the limit a/2y. According to the rules of the calculus, adx = 2ydy, which gives also dy/dx = a/2y.
Up to the last decades of the eighteenth century the differential calculus was centered on differentials. As we have seen above, the differential equations were written in terms of differentials. This state of affairs changed with Lagrange. In his paper on differential and integral calculus [22], published in 1774, and on his theory of analytical functions [23], published in 1797, the prominent role was played by the derivative, or derived function in the terminology introduced by as Lagrange. He also introduced the notation f x for a function of a variable x, without parentheses, but we will write f (x). However, he used parentheses in the case f (x + a).
Wishing to avoid differentials, Lagrange did not use the ratio of differentials as the definition of derivative. Instead he based his definition on the Taylor series. Taylor wrote the series that bears his name by considering two quantities z and x. When z increases to become z + v, the quantity x becomes Taylor published this result in 1715 [24] using the notation and terminology introduced by Newton. A quantity x is called a fluent and its rate of change, denoted byẋ, is called the fluxion of x. Thus, in expression (26), the ratio betweenẋ andż is the derivative of x with respect to z; the ratio betweenẍ andż 2 is the second derivative; and so on. If one wishes to construct a Taylor series for a certain function, one needs to know a priori the derivatives. However, the series can be obtained by other means without reference to derivatives. If this alternative procedure is employed, the derivatives could be obtained by comparing with the Taylor series. This was the reasoning of Lagrange, who considered the Taylor series as a way to define the derivative of a function [5,9]. He writes the fraction and then sets a = 0 after performing the subtraction in the numerator by the use of the power expansion in the increment a [22,23]. He denotes the result f x, but we will write f (x), and calls it the derived function of f (x).

As an example, Lagrange considers
Setting a = 0, one finds f (x) = 1/2 √ x. The definition of derivative given by Lagrange was based on the assumption that a function f (x) could be expanded in power series in the increment of x. Cauchy noted that this is not always possible. He then advanced a definition based on the concept of limit [9]. D'Alembert had already proposed this procedure but it lacked a proper definition of limit, a concept that was given by Cauchy in the following terms [25]: When the values successively assigned to the same variable approach indefinitely to a fixed value, so as to end by differing from it as little as desired, the latter is called the limit of all the others.
His definition of derivative becomes the limit of the quotient (28) when a approaches zero. This verbal definition, Cauchy translated into the precise language of deltas, epsilons and inequalities [26].
Cauchy also provided a meaning to the differential dy of a function y = f (x) [25]. The differential of the independent variable dx is a finite constant and According to Cauchy, this relation is the reason to call f (x) the differential coefficient [25], a terminology that was introduced by Lacroix [5].

The concept of function
The word function was used by Leibniz as the meaning of any quantity connected to a curve such as the locus, the slope, or the radius of curvature [14]. Johann Bernoulli and Euler regarded a function as an expression or a formula involving variables and constants, which is the concept usually held by the students of elementary mathematics courses [14]. The same concept of function was used by Lagrange. In his treatise on analytical functions, he writes [23]: We call function of one or more quantities, any expression of calculus in which these quantities enter in any manner, mixed or not with other quantities that we consider as having given and invariable values, while the quantities of the function can receive all the possible values.
Fourier, in his investigation on the equation of heat flows, departs from the notion of function as an analytic expression and approaches the modern concept of function. In his treatise on the theory of heat, where he develops the series that bears his name, he writes [27]: In general, the function f x represents a sequence of values, or ordinates, each of which is arbitrary. The abscissa x can receive an infinity of values, there is an equal number of ordinates f x. All have actual numeric values, either positive, or negative, or zero. It is not supposed that these ordinates are subject to a common law; they succeed one another in whatever manner, and each of them is given as it were a single quantity.
Dirichlet improved the definition given by Fourier arriving at a formulation that was based on the idea of function as a relationship between two variables [14,28]. It can be stated as follows: if x and y be two variables such that for each value assigned to x there corresponds, by some rule which need not be an analytical expression, then y is a function of x.

Laws of motion
The fundamental laws of motion were established by Newton in his treatise on mechanics, known as the Principia, published in 1687 [29]. A translation from the original in Latin was published later in 1729 [30]. Newton was not the first to formulate laws of motion. Prior to Newton, Galileo had formulated the law of inertia, which is one the fundamental laws of Newtonian mechanics, and also studied the motion of falling bodies, the trajectory of projectiles and the oscillations of a pendulum.
In his Dialogues Concerning two New Sciences, Galileo describes his observations of oscillations in the following terms [31,32]: Thousands of times I have observed vibrations especially in churches where lamps, suspended by long cords, had been inadvertently set into motion.
These occurrences might have stimulated him to study the oscillations of a pendulum. Galileo found that the period of oscillations of a simple pendulum is proportional to the length of the pendulum. He might have considered that the period is independent of the amplitude of oscillations but there is no mention concerning this issue on his writings [33].
In his treatise on the pendulum clock, published in 1673, Huygens showed that a cycloidal pendulum is isochronous, that is, the period is independent of its amplitude [34]. In contrast to a simple pendulum, which follows a circular path, a cycloidal pendulum follows a cycloid, a curve traced by a point on the edge of a circle rolling without sliding along a straight line, as shown in figure 2. When the point K goes into the point M, the disk rotates by an angle LOK so that MN equals the displacement of the center of the disk, which is the arc LK, plus LN, which gives ML equals the arc LK. Considering that CN is equal to the semi-circumference, it follows that CM is equal to the arc LG minus LN, which defines geometrically the cycloid. In modern notation and in parametric form the cycloid is described by where x and y denotes the distances BM and CM, θ is the angle LOG and r is the radius of the disk. Huygens demonstrated the isochronous properties of a cycloid by showing that it is a tautochrone, a curve holding the property that the time it takes for a body to descend to the lowest point is the same regardless of the starting point.
The cycloid is also the solution of the brachistochrone problem, posed by Johann Bernoulli in 1696 [35]. The problem is to find the curve which gives the fastest descent of a sliding body moving between two points that are not in the same vertical line. The solution given by Johann Bernoulli himself was based on the Fermat principle of the least time, which he adapted to his mechanical problem [36]. Accordingly, he used the sine law which is obtained from this principle. This law states that the sine of the angle of inclination of the curve at a certain point with respect to the vertical is proportional to the velocity. Referring to the figure 2, it is the ratio of the differentials dy =nm and ds =Mm, which is thus proportional to the velocity v, that is, or ady = vds, which after squaring gives a 2 dy 2 = v 2 ds 2 . Taking into account that ds 2 = dx 2 +dy 2 , where dx =Mn, one finds (a 2 − v 2 )dy 2 = v 2 dx 2 from which follows the equation To relate velocity and position, Johann Bernoulli uses the relation found by Galileo in his investigation on the fall of bodies according to which the velocity is proportional to the square root of the altitude traversed from the rest, v = √ ax. Employing this relation, Johann Bernoulli reaches the differential equation To shown that this equation describes a cycloid, he writes this equation in the equivalent form The integral of the second term is √ ax − x 2 , which is LN of figure 2, and the integral of the first term is the length of the arc GL. The total integral y = GL -LN, which is the statement that AMK is cycloid as we have seen above.
In 1687, Leibniz posed the following problem. Find a line of descent, in which a heavy body descends uniformly, and also approaches the horizontal line in equal time, that is, with a constant vertical velocity [37]. The solution is a semi-cubic parabola and a solution was given by Jacob Bernoulli in a paper published in 1690 [38]. A differential equation describing the curve was set up by Jacob Bernoulli as follows. The straight lines AL and MG are tangents to the points L and G of the curve NLG shown in figure 3. The ratio of LH= dy and FL= dx 2 + dy 2 is equal to the ratio of the vertical component of the velocity v y and the velocity v at the point L, dy A similar relation can be written for the point G, Considering that the vertical velocity is the same, v y = V y , we find the following relation b a Jacob Bernoulli uses the Galileo result according to which the velocity is proportional to the square root of the vertical altitude, to reach the result a b . After squaring, the following differential equation is found that describes the desired curve. The integral found by Bernoulli is Replacing y − a 3 /b 2 = z, which describes a semi-cubic parabola.

Principle of d'Alembert
Analytical mechanics was the name that Lagrange chose for his book on mechanics to emphasize the analytical framework on which rests his treatise on the subject [4]. He stressed his point of view by remarking that [4] No figures will be found in this work. The methods that expose in it require neither constructions nor geometrical or mechanical reasonings, but only the algebraic operations inherent to a regular and uniform process. Those who love analysis will, with joy, see mechanics become a new branch, and I am grateful that I have extended its domain.
The analytical mechanics of Lagrange is founded on the principle of virtual work, which was developed by Johann Bernoulli, in the static version, and by d'Alembert, in the dynamic version. The principle of virtual work stated by Johann Bernoulli generalized the previous formulations given by his predecessor [3]. The principle gives the condition for the static equilibrium of a mechanical system acted by several forces. In a letter written in 1717 to Varignon, which was reproduced in the treatise of Varignon on statics [39], Johann Bernoulli states the principle of virtual work, which he calls principle of virtual velocities, as follows: In all equilibrium of any forces, in whatever way they are applied, and in whatever direction that they act on each other, indirectly or directly, the sum of the positive energies will be equal to the sum of the negative energies taken positively.
In this statement, energy means the product of the magnitude of a force by a small displacement of the point acted by the force. The small displacement is parallel to the force, being positive if in the direction of the force and negative in the opposite direction. It should be noted that the small displacements Johann Bernoulli called virtual velocities.
In his treatise on analytical mechanics, Lagrange adopts the principle of virtual work, as a fundamental principle of statics. He writes the principle in the analytical form as where Q i are forces acting along the lines that define center of foces, and the differentials q i are small displacements along the lines of forces, described by the variables q i . Lagrange states that the equation (43) is the condition of equilibrium of a system of forces. It should be remarked that forces of reaction, which occurs for instance when a body rests on a surface, do not enter this equation because their works perform no work. Lagrange supposes next that the expression in the lefthand side of (43) is an exact differential dΦ of a function Φ, which means that Φ is a function of the variables q i . In modern terms, Lagrange assumes that the forces are conservative, Φ being understood as the potential energy. The condition for equilibrium becomes dΦ = 0, which means that the systems is configured in such as way that Φ is a maximum or a minimum. Lagrange shows furthermore that if the function is a minimum then the equilibrium is stable and the system will display small oscillations. If on the contrary the function is a maximum, the equilibrium is unstable and, being once disturbed, the system will depart from equilibrium.
The static principle of virtual work of Johann Bernoulli was extended to dynamic problems by d'Alembert [40]. According to d'Alembert, in dynamic problems one should take into account the inertial forces in addition to the actual forces. An analytical expression of the d'Alembert principle was given by Lagrange, and was obtained as follows. Lagrange argues that, the laws of motion of a body will be reduced to the laws of equilibrium if, following d'Alembert, one includes the inertial force, which is the mass of the body multiplied by the acceleration [4].
If we denote by x i a Cartesian coordinate of a body of mass m, the work of the inertial force along this coordinate is m(d 2 x i /dt 2 )δx i . Adding the work of all bodies and all Cartesian coordinates to the left hand side of (43), which is the analytical expression of the d'Alembert principle according to Lagrange [4]. Lagrange remarks that the differential signalized by δ, called variation, describes arbitrary increments or decrements and should be distinguished from the usual differential, signalized by d, which describes increments or decrements caused by the actual motion of the body.

Equations of Lagrange
The equations of motion were obtained by Lagrange from the principle expressed by equation (44) as follows. One starts by performing a change of variables from x i to new variables ξ j , not necessarily Cartesian. The relation between the differentials is where A ij = ∂x i /∂ξ j are functions of the new variables ξ j , and an analogous relation holds between the variations The first relation gives From the last two relations one obtains where which depends on the new variables ξ j only. Defining the quantity the following relations are obtained Comparing (48) with (51) and (52) one reaches the result The crucial step in the derivation of equation (53) is found in the second equality of equation (52). To show its validity it suffices to determine δ(dx i ) from (45) and d(δx i ) from (46). Since these two quantities are equal, the following relation is obtained where the equality δ(dξ j ) = d(δξ j ) has taken into account. Using this result and the relation (49) between B ij and A ij , it is straightforward to reach the second equality of equation (52).
Lagrange considers next the case where the forces Q i in equation (44) are such that is a complete differential and V will be a function of q i . If V is expressed in terms of the variables ξ i , we may write The replacement of this expression and the expression (53) in equation (44), yields the result where Lagrange concludes that if the variables ξ i are independent, each coefficient of δξ i will vanish, that is, which are the Lagrange equations of motion.

Thermodynamics
Thermodynamics emerged around the middle of the nineteenth century and presented two changes in the way the heat was conceived [1]. The first was the recognition that heat should be understood as a form of work, which lead to the law of conservation of energy. The second was the way in which heat was transformed in mechanical work, which allowed Clausius to define entropy and lead to the law of the increase in entropy. The heat absorbed minus the work performed by a system along a thermodynamic process is equal to the increase in the internal energy, and is independent of the trajectory connecting the the final and initial states. The energy is a state function which means to say that differential dQ − dW is a exact differential. In his first paper on thermodynamics [41], Clausius showed that dQ − dW is an exact differential by the use of the equivalence between heat and mechanical work.
Clausius starts his reasoning by assuming that a small quantity dQ of heat exchanged when the volume V and temperature T of a gas changes by dV and dT is given by where M and N are functions of V and T . Then he considers the total heat exchanged in a clockwise cycle. dQ = dW = pdV or dU = 0 along an isotherm. In an equivalent form (∂U/∂V ) T = 0, that is, the function U is independent of V , depending only on T .
Recalling that dQ = M dV + N dT , equation (74) is written in the form dQ = dU + pdV and we see that the heat capacity at constant volume is where C = dU/dT , and depends only on T , and that the heat capacity at constant pressure is from which follows C p = C v + R, that is, the difference in the heat capacities of ideal gas is a constant. The way in which heat is transformed in mechanical work was the main subject of research on heat carried out by Carnot. His investigations lead him to the following fundamental principle. When a system undergoes a cyclic process composed by two isothermal and two adiabatic processes the ratio of the work produced and the heat depends only on the two temperatures. In this principle of Carnot, heat was understood as a conserved quantity which in this case descends from a high temperature to a low temperature. Clausius modified this principle by replacing heat by heat absorbed. In addition he uses the law of Mayer and Joule to state that the work W is the heat absorbed Q 1 by the system at high temperature minus the heat Q 2 released by the system at low temperature, or more precisely, Q = Q 1 − Q 2 .
The modified principle of Carnot was written by Clausius in the form Q 1 /T 1 = Q 2 /T 2 where T 1 and T 2 are the absolute temperatures corresponding to the two isotherms. The generalization of this expression for any cycle lead Clausius to the result that dQ/T is a exact differential [43]. Later on, he wrote dS = dQ/T and called S the entropy [44]. The infinitesimal heat absorbed by a system in equilibrium becomes related to the differential of entropy by dQ = T dS and which is the conservation of energy in differential form, where all differentials involved are exact differentials.

Conclusion
We have analyzed the role of exact and inexact differentials in the early developments of mechanics and thermodynamics. We have also examined the evolution of differential calculus in relation to the concepts related to exact and inexact differentials. Euler introduced the concept of integrating factor, which he used to solve an ordinary linear differential equation of the first-order. Euler also found the condition for a complete differential by examining the differential of a function of two variables. In an independent way, Clairaut also reached the same condition. When this condition is not fulfilled the differential is incomplete or inexact. In the analytical treatment of mechanics, Lagrange considered forces whose differential work is an exact differential, in which case it is possible to define a work function and reach the conservation of energy. When the differential is inexact, it is possible to transform it into and an exact differential as long as an integrating factor can found. The fundamental law of equilibrium thermodynamics introduced by Clausius stating that dQ/T = dS can equally be formulated by declaring that 1/T is an integrating factor related to the inexact differential dQ. The law of conservation of energy was also written by Clausius in terms of exact differentials. The process of finding exact differentials was a fundamental procedure which allowed the formulation of thermodynamics in terms of state functions such as energy and entropy.