Symbolic Computation to Solving an Irrational Equation on Based Symmetric Polynomials Method

In this article, we examine the use of symmetry groups for modeling applied problems through computer symbolic calculus. We consider the problem of solving radical equations symbolically using computer mathematical packages. We propose some methods to obtain a correct analytical solution for this class of equations by means of the Mathcad package. The application of symmetric polynomials is proposed to ensure a correct approach to the solution. Issues on the solvability based on the physical sense of a problem are discussed. Common errors in solving radical equations related to the specificity of the computer usage are analyzed. Provable electrical and geometrical problems are illustrated as example.


Introduction
Currently, an indisputable fact is the widespread and ubiquitous use of computer mathematical packages for the modeling of applied problems in any field [1,2]. As a consequence, the question of the correctness of the solutions obtained using these packages is extremely relevant. This issue is especially critical when performing symbolic calculations, where there is the need to obtain an analytical calculation formula with a large number of variable input parameters. Thus, in this work, some cases are illustrated in which the use of embedded algorithms for computer symbolic calculations leads to the appearance of additional extraneous solutions. The numerical validation of all the obtained solutions described by an analytical formula is not feasible. Therefore, it is necessary to develop a method for obtaining the correct symbolic solution using computer mathematical packages, based on mathematical theory. In this article, it is proposed to use some methods based on the theory of symmetric groups [3].
In this paper, we use such a branch of the theory of symmetric groups as symmetric polynomials [4,5]. Symmetric polynomials are used in a wide class of algebraic problems, in the solution of systems of higher order equations, in the solution of irrational equations, for proving identities and inequalities, for factoring, for simplifying algebraic expressions, for extracting roots and for eliminating irrationalities in the denominator [6,7]. In this article, we consider the use of symmetric polynomials for solving a class of algebraic equations, i.e., the irrational equations, namely, a method is developed for obtaining the correct symbolic solution of this class of equations through the use of mathematical packages. The choice of irrational equations is due to the fact that the specificity of the embedded computer where is the variable, , , , , are real parameters, with and not equal to zero at the same time and ≥ 0.
By solving symbolically Equation (1) using the standard tools of the mathematical package Mathcad, we obtain two roots (see the listing of the Mathcad program in Figure 1). We denote the found roots functions of several parameters as ( , , , , ) and ( , , , , ).
Note that the formulas obtained for the roots ( , , , , ) and ( , , , , ) do not result in root variants for the case = , because for = Equation (1) has a single root due to the reduction of terms with the same coefficients (see the listing of the program part in Mathcad in Figure  2).  We denote the found roots functions of several parameters as z r1 (a, d, c, d, f ) and z r2 (a, d, c, d, f ). Note that the formulas obtained for the roots z r1 (a, d, c, d, f ) and z r2 (a, d, c, d, f ) do not result in root variants for the case a = c, because for a = c Equation (1) has a single root due to the reduction of terms with the same coefficients (see the listing of the program part in Mathcad in Figure 2). embedded computer algorithms for solving this class of equations more often than others leads to the appearance of extraneous solutions. As an example, an electrical problem is examined using the mathematical package Mathcad [8,9]. This package is very often used for solving applied problems by virtue of its intuitive user interface. In addition, the article discusses techniques for analyzing irrational equations, in order to adapt the method based on symmetric polynomials to some specific cases. These techniques are further illustrated by the example of a graphical geometrical problem about the size of a submarine.

The Core of the Problem of the Extraneous Solutions
Let us demonstrate the essence of the problem of the false solutions to irrational equations when using computer mathematical packages.
Consider an irrational equation with general parameters: where is the variable, , , , , are real parameters, with and not equal to zero at the same time and ≥ 0.
By solving symbolically Equation (1) using the standard tools of the mathematical package Mathcad, we obtain two roots (see the listing of the Mathcad program in Figure 1). We denote the found roots functions of several parameters as ( , , , , ) and ( , , , , ).
Note that the formulas obtained for the roots ( , , , , ) and ( , , , , ) do not result in root variants for the case = , because for = Equation (1) has a single root due to the reduction of terms with the same coefficients (see the listing of the program part in Mathcad in Figure  2).  Similarly, Equation (1) has only one root in the case a = 0 or c = 0. We will consider Equation (1) in the more general case of a c. The analysis of the found roots z r1 (a, d, c, d, f ) and z r2 (a, d, c, d, f ) shows that they do not satisfy Equation (1) for all the parameter values. To demonstrate this, we introduce the quantities ∆ 1 (a, d, c, d, f ) and ∆ 2 (a, d, c, d, f ) which are the differences between the left and right sides of Equation (1) where the variable z is substituted with the values of the found roots z r1 (a, d, c, d, f ) and z r2 (a, d, c, d, f ), respectively. Obviously, if the roots z r1 (a, d, c, d, f ) and z r2 (a, d, c, d, f ) do satisfy Equation (1), then the values ∆ 1 (a, d, c, d, f ) = 0 and ∆ 2 (a, d, c, d, f ) = 0.
(a) (b)  Thus, using mathematical packages it is possible to ascertain the appearance of extraneous roots in the symbolic solution of an irrational equation of the form (1).

Methods of Elimination of the False Solutions
The occurrence of extraneous solutions corresponding to certain values of the parameter , depends on the following fact. The method of solving irrational equations of the type of Equation (1) adopted in the theory of algebra involves the raising of both sides of an equation to the same even degree [10,11]. In this case, the radicals are assumed the principal square root. This occurrence determines the appearance of extraneous roots that correspond to negative values of the radicals, as well as to negative value of the radicands (if we take into account the possibility of complex roots, as implicit in mathematical packages). Two approaches are usually used to eliminate the extraneous roots [10,11].

Validation by the Numerical Solution
The first approach consists in finding all the roots and in their subsequent validation by substitution into the initial equation. This approach is convenient for the numerical solution of equations. In mathematical packages, in particular, in Mathcad, this approach is implemented in the numerical solution of irrational equations, thus ensuring the elimination of extraneous roots. For 31  Thus, using mathematical packages it is possible to ascertain the appearance of extraneous roots in the symbolic solution of an irrational equation of the form (1).

Methods of Elimination of the False Solutions
The occurrence of extraneous solutions corresponding to certain values of the parameter f , depends on the following fact. The method of solving irrational equations of the type of Equation (1) adopted in the theory of algebra involves the raising of both sides of an equation to the same even degree [10,11]. In this case, the radicals are assumed the principal square root. This occurrence determines the appearance of extraneous roots that correspond to negative values of the radicals, as well as to negative value of the radicands (if we take into account the possibility of complex roots, as implicit in mathematical packages). Two approaches are usually used to eliminate the extraneous roots [10,11].

Validation by the Numerical Solution
The first approach consists in finding all the roots and in their subsequent validation by substitution into the initial equation. This approach is convenient for the numerical solution of equations. In mathematical packages, in particular, in Mathcad, this approach is implemented in the numerical solution of irrational equations, thus ensuring the elimination of extraneous roots. For example, in  However, when searching for an analytical solution, it is impossible to perform this test of all the roots for all the conceivable values of the parameters. Therefore, when performing symbolic calculations, mathematical packages do not solve this problem at the beginning and, according to the method of raising the two sides of the equation to the same power degree, they provide a calculation formula for the roots (see Figure 1), which includes both positive and negative values of radicals, as well as of the values of the radicands. Therefore, for example, with parameters values given above , , , , the value of the first radical in Equation (1) is √ 1 ( , , , , ) • 2 + 1 = 22.323, while the value of the second radical is √ 1 ( , , , , ) • (−4) + 1000 = 2.323. Obviously, the correct sum of the radicals = 20 can be obtained by taking the second radical with the negative sign.

Determining the Domain of the Admissible Solutions through Symmetric Polynomials
The second approach for the elimination of the extraneous roots is to solve simultaneously the system given by the irrational equation and by a number of additional inequality constraints imposed on the roots of the equation and arising from the requirement of using only the principal square root. From the above reasoning, it follows that for the problem in question, these will be inequalities that define the range of admissible values of the parameter .
To determine the range of admissible values of the parameter , it was proposed to use the properties of symmetry, namely the properties of symmetric polynomials. In this case, the requirement of symmetry is not imposed on the equation being solved. We use a method for reducing a polynomial to a symmetric form [6]. We introduce the notation: Obviously the expressions 1 = + , 2 = • are the simplest symmetric polynomials. Equation (1) becomes: from this it derives: Because the radicals must be non-negative, 1 ≥ 0 and 2 ≥ 0. Since 1 = , that requirement 1 ≥ 0 is verified for ≥ 0. From the condition 2 ≥ 0, it follows the range of admissible values for the parameter : However, when searching for an analytical solution, it is impossible to perform this test of all the roots for all the conceivable values of the parameters. Therefore, when performing symbolic calculations, mathematical packages do not solve this problem at the beginning and, according to the method of raising the two sides of the equation to the same power degree, they provide a calculation formula for the roots (see Figure 1), which includes both positive and negative values of radicals, as well as of the values of the radicands. Therefore, for example, with parameters values given above a, d, c, d, f the value of the first radical in Equation (1) Obviously, the correct sum of the radicals f = 20 can be obtained by taking the second radical with the negative sign.

Determining the Domain of the Admissible Solutions through Symmetric Polynomials
The second approach for the elimination of the extraneous roots is to solve simultaneously the system given by the irrational equation and by a number of additional inequality constraints imposed on the roots of the equation and arising from the requirement of using only the principal square root. From the above reasoning, it follows that for the problem in question, these will be inequalities that define the range of admissible values of the parameter f .
To determine the range of admissible values of the parameter f , it was proposed to use the properties of symmetry, namely the properties of symmetric polynomials. In this case, the requirement of symmetry is not imposed on the equation being solved. We use a method for reducing a polynomial to a symmetric form [6]. We introduce the notation: Obviously the expressions σ 1 = x + y, σ 2 = x·y are the simplest symmetric polynomials. Equation (1) becomes: from this it derives: Because the radicals must be non-negative, σ 1 ≥ 0 and σ 2 ≥ 0.
Since From the condition σ 2 ≥ 0, it follows the range of admissible values for the parameter f : Finally, we solve the system given by Equation (1) and inequality (5): The consideration of system (6) instead of Equation (1) makes it possible to obtain a symbolic solution without extraneous roots, as we will show on some applied problems.
Note that system (6) corresponds to a generally accepted scheme for setting the range of admissible values when solving an algebraic equation [4,5]. However, in our case, it is not possible to determine the admissible region of the parameter (5) based only on the general form of the initial Equation (1), since it is required the use of symmetric polynomials.

Application of the Proposed Method for Solving an Electrical Problem
Let us clarify the core of the problem of the extraneous solutions by an example of an applied electrical problem. It is requested to find the optimal distribution of the total power P = P 1 + P 2 to two users of capacities P 1 and P 2 with a required total voltage U = U 1 + U 2 , less than the breakdown voltage. Users resistances are equal respectively to R 1 and R 2 .
According to the theory of electrical circuits, the expressions for the powers are described by the formulas [12]: Then, setting as parameters the resistances R 1 and R 2 , the total voltage U and the total power P, an irrational equation is derived for the sought variable P 1 : By solving symbolically Equation (8) using the standard tools of the mathematical package Mathcad, we obtain two roots (see the listing of the Mathcad program in Figure 5).
Finally, we solve the system given by Equation (1) and inequality (5): The consideration of system (6) instead of Equation (1) makes it possible to obtain a symbolic solution without extraneous roots, as we will show on some applied problems.
Note that system (6) corresponds to a generally accepted scheme for setting the range of admissible values when solving an algebraic equation [4,5]. However, in our case, it is not possible to determine the admissible region of the parameter (5) based only on the general form of the initial Equation (1), since it is required the use of symmetric polynomials.

Application of the Proposed Method for Solving an Electrical Problem
Let us clarify the core of the problem of the extraneous solutions by an example of an applied electrical problem. It is requested to find the optimal distribution of the total power = + to two users of capacities and with a required total voltage = + , less than the breakdown voltage. Users resistances are equal respectively to and . According to the theory of electrical circuits, the expressions for the powers are described by the formulas [12]: Then, setting as parameters the resistances and , the total voltage and the total power , an irrational equation is derived for the sought variable : By solving symbolically Equation (8) using the standard tools of the mathematical package Mathcad, we obtain two roots (see the listing of the Mathcad program in Figure 5). It could be noted that by executing the calculations using the Mathematica or the Matlab packages, equivalent roots are obtained (see the listings of programs in Wolfram Mathematica [13] and Matlab [14,15] packages in Figure 6). It could be noted that by executing the calculations using the Mathematica or the Matlab packages, equivalent roots are obtained (see the listings of programs in Wolfram Mathematica [13] and Matlab [14,15] packages in Figure 6).    The found roots are functions of several parameters P r1 (U, P, R 1 , R 2 ) and P r2 (U, P, R 1 , R 2 ). After assigning to the parameters the following values P = 8000, R 1 = 1 , R 2 = 2, the roots of Equation (8) P r1 (U, 8000, 1, 2) and P r2 (U, 8000, 1, 2) are plotted as function of the required total voltage U (Figure 7).    Let us verify whether the functions P r1 (U, 8000, 1, 2) and P r2 (U, 8000, 1, 2) are actually solutions of Equation (8) for any value of U. In order to do this, we build a plot of the right side of Equation (8)-simply Y(U) = U, and a plot of the left side of Equation (8), by substituting in it the functions P r1 (U, 8000, 1, 2) and P r2 (U, 8000, 1, 2) instead of the sought variable P 1 (Figure 8). Obviously the solutions P r1 (U, 8000, 1, 2) and P r2 (U, 8000, 1, 2) are true solutions of the problem only for those values of U, where the plots of the left and right sides match.  Thus, we have ascertained that the solution of Equation (8) shown in Figure 5 includes extraneous roots.

Implementation of the Proposed Method
We use the proposed method for the elimination of the extraneous roots. We need to determine the inequalities that define the range of admissible values of the parameter-total voltage .
To determine the range of admissible values of the parameter , it was proposed to use the properties of symmetric polynomials. Note that the Equation (8) here considered is not symmetric with respect to variables , . So the plots of the roots ( , 8000, 1, 2) and ( , 8000, 1, 2) are also not symmetric (Figure 7). We introduce the notation: Equation (8) becomes: from this it derives: Because the radicals must be non-negative, ≥ 0 and ≥ 0. Since = , that requirement ≥ 0 is verified for ≥ 0. From the condition ≥ 0, it follows the range of admissible values for the parameter : Comparing the plots of the roots in Figure 8 and the plot of the function ( ) in Figure 9 for specific values of , , , , we notice that the lower bound for the value of = 89.4 corresponds to the minimum admissible value of the variable = 0, while the point of appearance of the second root = 126.5 corresponds to the maximum admissible value = = 8000. According to the plot: • For U < 89.4 both solutions P r1 (U, 8000, 1, 2) and P r2 (U, 8000, 1, 2) are false; • For 89.4 ≤ U < 126.5 solution P r1 (U, 8000, 1, 2) is false, while P r2 (U, 8000, 1, 2) is true; • For U ≥ 126.5 both solutions P r1 (U, 8000, 1, 2) and P r2 (U, 8000, 1, 2) are true.
Thus, we have ascertained that the solution of Equation (8) shown in Figure 5 includes extraneous roots.

Implementation of the Proposed Method
We use the proposed method for the elimination of the extraneous roots. We need to determine the inequalities that define the range of admissible values of the parameter-total voltage U.
To determine the range of admissible values of the parameter U, it was proposed to use the properties of symmetric polynomials. Note that the Equation (8) here considered is not symmetric with respect to variables P 1 , P 2 . So the plots of the roots P r1 (U, 8000, 1, 2) and P r2 (U, 8000, 1, 2) are also not symmetric (Figure 7). We introduce the notation: Equation (8) becomes: from this it derives: Because the radicals must be non-negative, σ 1 ≥ 0 and σ 2 ≥ 0. Since σ 1 = U, that requirement σ 1 ≥ 0 is verified for U ≥ 0. From the condition σ 2 ≥ 0, it follows the range of admissible values for the parameter U: Comparing the plots of the roots in Figure 8 and the plot of the function σ 2 (U) in Figure 9 for specific values of P, R 1 , R 2 , P 1 , we notice that the lower bound for the value of U = 89.4 corresponds to the minimum admissible value of the variable P 1 = 0, while the point of appearance of the second root U = 126.5 corresponds to the maximum admissible value P 1 = P = 8000.  Finally, we solve the system given by Equation (8) and inequality (12): The plot of the correct analytical solution obtained using Mathcad is shown in Figure 10. Comparing it with the plot in Figure 7, we see that the extraneous roots are eliminated, and, as noted earlier, in the interval 89.4 ≤ < 126.5 Equation (8) has only one root ( , 8000, 1, 2), while in the interval ≥ 126.5-two roots ( , 8000, 1, 2) and ( , 8000, 1, 2).

A Special Case
Let us consider a special case of Equation (8), when the user resistances are equal = . In this case, the left side of Equation (8) is already initially a symmetric polynomial with respect to the variables , : The plots of the solution (14), shown in Figure 11 and similar to the plots in Figure 7, have a symmetric appearance. Finally, we solve the system given by Equation (8) and inequality (12): The plot of the correct analytical solution obtained using Mathcad is shown in Figure 10. Comparing it with the plot in Figure 7, we see that the extraneous roots are eliminated, and, as noted earlier, in the interval 89.4 ≤ U < 126.5 Equation (8) has only one root P r (U, 8000, 1, 2), while in the interval U ≥ 126.5-two roots P r1 (U, 8000, 1, 2) and P r2 (U, 8000, 1, 2).  Finally, we solve the system given by Equation (8) and inequality (12): The plot of the correct analytical solution obtained using Mathcad is shown in Figure 10. Comparing it with the plot in Figure 7, we see that the extraneous roots are eliminated, and, as noted earlier, in the interval 89.4 ≤ < 126.5 Equation (8) has only one root ( , 8000, 1, 2), while in the interval ≥ 126.5-two roots ( , 8000, 1, 2) and ( , 8000, 1, 2).

A Special Case
Let us consider a special case of Equation (8), when the user resistances are equal = . In this case, the left side of Equation (8) is already initially a symmetric polynomial with respect to the variables , : The plots of the solution (14), shown in Figure 11 and similar to the plots in Figure 7, have a symmetric appearance.

A Special Case
Let us consider a special case of Equation (8), when the user resistances are equal R 1 = R 2 . In this case, the left side of Equation (8) is already initially a symmetric polynomial with respect to the variables P 1 , P 2 : The plots of the solution (14), shown in Figure 11 and similar to the plots in Figure 7, have a symmetric appearance.  System (13) in this case reduces to the form: The overall form of the analytical solution of system (15), obtained using Mathcad, is shown in Figure 12. The analytical solution is written by using the conditional if operator, which allows to determine a different type of solution for different ranges of values . The plot of the correct analytical solution of the system (15), shown in Figure 13, has obviously a symmetric appearance, i.e., on the entire allowable interval ≥ 126.5 there are two symmetric roots ( , 8000, 1, 2) and ( , 8000, 1, 2).  System (13) in this case reduces to the form: The overall form of the analytical solution of system (15), obtained using Mathcad, is shown in Figure 12. The analytical solution is written by using the conditional if operator, which allows to determine a different type of solution for different ranges of values U.  Figure 11. The initial set of the solutions of Equation (14).
System (13) in this case reduces to the form: The overall form of the analytical solution of system (15), obtained using Mathcad, is shown in Figure 12. The analytical solution is written by using the conditional if operator, which allows to determine a different type of solution for different ranges of values .   The plot of the correct analytical solution of the system (15), shown in Figure 13, has obviously a symmetric appearance, i.e., on the entire allowable interval U ≥ 126.5 there are two symmetric roots P r1 (U, 8000, 1, 2) and P r2 (U, 8000, 1, 2).  System (13) in this case reduces to the form: The overall form of the analytical solution of system (15), obtained using Mathcad, is shown in Figure 12. The analytical solution is written by using the conditional if operator, which allows to determine a different type of solution for different ranges of values . The plot of the correct analytical solution of the system (15), shown in Figure 13, has obviously a symmetric appearance, i.e., on the entire allowable interval ≥ 126.5 there are two symmetric roots ( , 8000, 1, 2) and ( , 8000, 1, 2). Figure 13. Set of the correct solutions of system (15). Figure 13. Set of the correct solutions of system (15).
Thus, in this section, we have described a method for obtaining the correct analytical solution of an irrational equation using the properties of symmetry, demonstrating it using as example the solution of Equation (8).

Application of Additional Methods of Analysis for Solving Irrational Equations Using Computer Mathematical Packages
In many applications, it is possible to simplify the method of obtaining the correct analytical solution of an irrational equation, if in the process of reduction to a symmetric form, the initial irrational equation is analyzed.
Let us consider the problem of finding the geometrical dimensions of a submarine, whose shape consists of a straight circular cylinder and two identical straight circular cones ( Figure 14). Thus, in this section, we have described a method for obtaining the correct analytical solution of an irrational equation using the properties of symmetry, demonstrating it using as example the solution of Equation (8).

Application of Additional Methods of Analysis for Solving Irrational Equations Using Computer Mathematical Packages
In many applications, it is possible to simplify the method of obtaining the correct analytical solution of an irrational equation, if in the process of reduction to a symmetric form, the initial irrational equation is analyzed.
Let us consider the problem of finding the geometrical dimensions of a submarine, whose shape consists of a straight circular cylinder and two identical straight circular cones ( Figure 14). The surface area and the volume of the vessel are given. It is necessary to find the length of the cylinder, the radius of the base of the two cones and of the cylinder and the height of the two cones.

Initial Formulation of the Problem
According to the known geometry formulas, the surface area and the volume are determined by the following relations [16]: Expressing the variable through both Equation (16), we obtain the relations: Equating the ratios we obtain an equation in two variables and . Considering as parameters the values of , , , we could solve Equation (18) with respect to the variable ( , , ).

Analytical Solution
As in the previous example ( Figure 5), the analytical solution of the irrational Equation (18) in the Mathcad package consists of two roots (see listing in Figure 15). The surface area S and the volume V of the vessel are given. It is necessary to find the length L of the cylinder, the radius R of the base of the two cones and of the cylinder and the height H of the two cones.

Initial Formulation of the Problem
According to the known geometry formulas, the surface area and the volume are determined by the following relations [16]: Expressing the variable L through both Equation (16), we obtain the relations: Equating the ratios we obtain an equation in two variables H and R. Considering as parameters the values of S, V, R, we could solve Equation (18) with respect to the variable H(R, V, S).

Analytical Solution
As in the previous example ( Figure 5), the analytical solution of the irrational Equation (18) in the Mathcad package consists of two roots (see listing in Figure 15).  The plots of the roots ( , , ) and ( , , ) with some specified values for the parameters = 1000 м and = 1500 м are shown in Figure 16. For the generality of the considered example, we will not yet take into account the physical meaning of the variable ( , , ), allowing for negative and even complex values.   The plots of the roots H 1 (R, V, S) and H 2 (R, V, S) with some specified values for the parameters S = 1000 м 2 and V = 1500 м 2 are shown in Figure 16. For the generality of the considered example, we will not yet take into account the physical meaning of the variable H(R, V, S), allowing for negative and even complex values.  The plots of the roots ( , , ) and ( , , ) with some specified values for the parameters = 1000 м and = 1500 м are shown in Figure 16. For the generality of the considered example, we will not yet take into account the physical meaning of the variable ( , , ), allowing for negative and even complex values.

Analysis of the Solution of the Problem in the Initial Formulation
Let us check whether the functions H 1 (R, 1500, 1000) and H 2 (R, 1500, 1000) are solutions of Equation (18) for all the values of the parameter R. In order to do this check, we construct the plots of the left L S and right L V sides of Equation (18) with the sought variable H substituted by the functions H 1 (R, 1500, 1000) and H 2 (R, 1500, 1000) ( Figure 17). Obviously, the solutions are true only for those values of R, where the plots of the left and right sides coincide.  The plots of the roots ( , , ) and ( , , ) with some specified values for the parameters = 1000 м and = 1500 м are shown in Figure 16. For the generality of the considered example, we will not yet take into account the physical meaning of the variable ( , , ), allowing for negative and even complex values.   As can be seen from the plots in Figure 17, the solutions H 1 (R, 1500, 1000) and H 2 (R, 1500, 1000) are not valid for some values of R, that is, the resulting analytical solution contains extraneous roots (these are the left branches in Figure 16). It is difficult to determine the lower limit of the values of the R parameter so that the solutions are true, because the plots (Figures 16 and 17) show only real values of the roots H 1 (R, 1500, 1000) and H 2 (R, 1500, 1000). In this case, generally speaking, the analytical solution presented in Figure 16 includes complex values. This is due to the fact that if the initially real numerical parameters of the equation are not specified in the mathematical package, the symbolic solution is performed over the entire field of the complex numbers. Therefore, for example, when R = 3, we get the imaginary root H 1 (R, V, S) = 4.025 i, which satisfies Equation (18), that is, the values of L V and L S are complex and coincide.

Determining the Range of Admissible Solutions
In order to determine the range of admissible values for the parameter R of Equation (18), we use also the reduction method to symmetric polynomials. We note that after introducing a restriction for the variable R in the region of the real numbers, we thus proceed to solve the equation in the real domain, and the sought analytical solution will no longer contain complex values.
Transform Equation (18) to the form: In order to bring the left side of Equation (19) to a symmetric form, it is necessary to introduce a second irrationality of the variable H.
For the generality of the considered case we allow negative values of the variable H , neglecting for the time the physical sense. Therefore, we will consider the two options: H ≥ 0 and H < 0.
For H ≥ 0 Equation (19) can be rewritten in the form: The left side of (20) is not symmetric, however, we select the symmetric part in it and analyze the resulting equation: Since the values of all radicals in (21) must be non-negative, it follows the inequality: Consequently, the entire left side of Equation (21) is also negative, which implies the negativity of the right side (21): Analogously for H < 0 Equation (19) can be rewritten in the form: In this case, even without isolating the symmetric part, it is obvious that the left side of Equation (24) is negative due to the use of the principal square root. Therefore, for H < 0 also (17) is true.
Thus, by analyzing Equation (19), we obtain an inequality that determines the range of admissible values of the parameter R, which follows directly from (23): Now, by solving the system given by Equation (18) and inequality (25): we get the correct analytical solution of Equation (18) without extraneous roots. The general view of the correct analytical solution obtained with Mathcad is exactly identical to the solution presented in Figure 18. Figure 18a shows the plots of the true analytical solution. Figure 18b shows the plots of the left L S and right L V sides of Equation (18) we get the correct analytical solution of Equation (18) without extraneous roots. The general view of the correct analytical solution obtained with Mathcad is exactly identical to the solution presented in Figure 18. Figure 18a shows the plots of the true analytical solution. Figure 18b shows the plots of the left and right sides of Equation (18) with the correct functions ( , 1500,1000) and ( , 1500,1000) substituted in place of the sought variable . Comparing Figure 18a with Figure  16 and Figure 18b with Figure 17, we see that the false solutions were eliminated.

Eliminating Solutions without Physical Meaning
Now let's return to the physical meaning of the geometrical problem under consideration. It is necessary to take into account the positivity of the sought values of , , . The positivity of follows from inequality (25). Therefore, it is sufficient to add to system (26) the corresponding inequalities for the quantities , using (17): or > 0, Then, solving together (26) and (28), we obtain the correct analytical solution that takes into account the physical meaning of the geometrical quantities. The solutions ( , , ) together with the corresponding values ( , , ) = ( , , ) are shown in Figure 19. Obviously, they are nonnegative sub-regions of the solutions shown in Figure 18.

Eliminating Solutions without Physical Meaning
Now let's return to the physical meaning of the geometrical problem under consideration. It is necessary to take into account the positivity of the sought values of H, L, R. The positivity of R follows from inequality (25). Therefore, it is sufficient to add to system (26) the corresponding inequalities for the quantities H, L using (17): or Then, solving together (26) and (28), we obtain the correct analytical solution that takes into account the physical meaning of the geometrical quantities. The solutions H(R, V, S) together with the corresponding values L S (R, S, H) = L V (R, V, H) are shown in Figure 19. Obviously, they are non-negative sub-regions of the solutions shown in Figure 18.  It is possible to define up to third digit the boundaries of the first sub-region of solutions as 3.146 ≤ ≤ 4.578 , and with 3.146 ≤ < 3.208 there are two roots ( , 1500,1000) and ( , 1500,1000), while with ≥ 3.208 only one ( , 1500,1000). The second subdomain of the solutions is defined for 10.708 ≤ ≤ 12.157, and here there is only one root ( , 1500,1000).
The corresponding 3d-plot of the solution area ( , , ), calculated analytically taking into account the physical meaning of the variables, is shown in Figure 20. It is possible to define up to third digit the boundaries of the first sub-region of solutions as 3.146 ≤ R ≤ 4.578, and with 3.146 ≤ R < 3.208 there are two roots H 1 (R, 1500, 1000) and H 2 (R, 1500, 1000), while with R ≥ 3.208 only one H 1 (R, 1500, 1000). The second subdomain of the solutions is defined for 10.708 ≤ R ≤ 12.157, and here there is only one root H(R, 1500, 1000).
The corresponding 3d-plot of the solution area (H, L, R), calculated analytically taking into account the physical meaning of the variables, is shown in Figure 20.

Discussions of the Results
The results obtained in the article give the opportunity to indicate the directions of further research in the field of identified problems.
Primarily, it is necessary to determine the usage field of the proposed method, i.e., the determination of the irrational equations that, can be reduced to a form that allows the use of such a method. The starting point of the study can be the statement [4,5], that there is an asymmetric (globally) polynomial ( , ), such as: for any symmetric polynomial ( , ). This problem contains the need of research of method usage possibility for irrational equations with high degree radicals using symmetric polynomials of the form = + . This can be realized due to the statement [4,5], that any sum of power is representable as a polynomial in variables = + , = • . Certainly, this will require additional analysis of the inequalities for the determination of the admissible region of parameter values. Since one of the main applications of symmetric polynomials in algebra is their use for solving systems of equations [6,7], an extension of the proposed method for solving systems of irrational equations seems to be a promising direction.
A separate study requires the application of the method for system of irrational equations, nonlinear in the variables. In Section 3 of the article, we considered a special case of applying the method for nonlinear irrational equation. However, a generalization of the issue requires additional research, which we plan to conduct in future research.

Discussions of the Results
The results obtained in the article give the opportunity to indicate the directions of further research in the field of identified problems.
Primarily, it is necessary to determine the usage field of the proposed method, i.e., the determination of the irrational equations that, can be reduced to a form that allows the use of such a method. The starting point of the study can be the statement [4,5], that there is an asymmetric (globally) polynomial F(σ 1 , σ 2 ), such as: F(x, y) = F(x + y, xy), for any symmetric polynomial F(x, y). This problem contains the need of research of method usage possibility for irrational equations with high degree radicals using symmetric polynomials of the form S n = x n + y n . This can be realized due to the statement [4,5], that any sum of power S n is representable as a polynomial in variables σ 1 = x + y, σ 2 = x·y. Certainly, this will require additional analysis of the inequalities for the determination of the admissible region of parameter values.
Since one of the main applications of symmetric polynomials in algebra is their use for solving systems of equations [6,7], an extension of the proposed method for solving systems of irrational equations seems to be a promising direction.
A separate study requires the application of the method for system of irrational equations, nonlinear in the variables. In Section 3 of the article, we considered a special case of applying the method for nonlinear irrational equation. However, a generalization of the issue requires additional research, which we plan to conduct in future research.

Conclusions
The article presented critical errors in the symbolic solution of irrational equations using mathematical packages. A method of obtaining a correct analytical solution using the properties of symmetry, namely symmetric polynomials, is proposed. Some methods of further analysis of irrational equations with the aim of obtaining a correct analytical solution are presented. The solution methods are established by taking into account the physical meaning of the problems. This algorithm can also be used in innovative engineering education [17], in STEM (Science, Technology, Engineering and Math including Computer Science) education [18].