Next Article in Journal
Feature Selection of Non-Dermoscopic Skin Lesion Images for Nevus and Melanoma Classification
Previous Article in Journal
Accurate Sampling with Noisy Forces from Approximate Computing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Theoretical Bases of Heat Engineering, National Research University Moscow Power Engineering Institute, Moscow 111250, Russia
2
Department of Mathematics, Military Educational and Scientific Center of the Air Force “N.E. Zhukovsky and Y.A. Gagarin Air Force Academy”, Voronezh 394064, Russia
3
Viale delle Mimose, 18, 64025 Pineto (TE), Italy
4
MIREA—Russian Technological University, Moscow 119454, Russia
*
Author to whom correspondence should be addressed.
Submission received: 8 April 2020 / Revised: 20 April 2020 / Accepted: 23 April 2020 / Published: 28 April 2020

Abstract

:
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.

1. 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 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.

2. The Method of Symmetric Polynomials for Solving Irrational Equations Using Computer Mathematical Packages

2.1. 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:
a · z + b + c · z + d = f
where z   is the variable, a ,   b , c , d , f are real parameters, with a and c not equal to zero at the same time and f 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 z r 1 ( a , d , c , d , f ) and z r 2 ( a , d , c , d , f ) .
Note that the formulas obtained for the roots z r 1 ( a , d , c , d , f ) and z r 2 ( 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).
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 r 1 ( a , d , c , d , f ) and z r 2 ( 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 r 1 ( a , d , c , d , f ) and z r 2 ( a , d , c , d , f ) , respectively. Obviously, if the roots z r 1 ( a , d , c , d , f ) and z r 2 ( 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 .
We plot 1 ( a , d , c , d , f ) and 2 ( a , d , c , d , f ) as functions of the parameter f for a = 2 ,   b = 1 ,   c = ± 4 ,   d = 1000   (see Figure 3a,b).
Figure 3a,b shows that the roots z r 1 ( a , d , c , d , f ) and z r 2 ( a , d , c , d , f ) do not satisfy Equation (1) for all values of the parameter f . Therefore, for c = 4 (Figure 3a), the root z r 1 ( a , d , c , d , f ) satisfies (1) only for   f 22.4 , and the root z r 2 ( a , d , c , d , f ) satisfies for   f 31.6 . For c = 4   (Figure 3b), the root z r 2 ( a , d , c , d , f ) satisfies (1) only for   f 31.6 , while the root   z r 1 ( a , d , c , d , f ) does not satisfy at all (1).
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).

2.2. 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].

2.2.1. 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 the numerical solution of Equation (1) for the values of the parameters a = 2 ,   b = 1 ,   c = 4 ,   d = 1000 ,   f = 20   the result is “No solution”, which corresponds to the plot in Figure 3a (see the listing of the Mathcad program in Figure 4).
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) is z r 1 ( a , d , c , d , f ) · 2 + 1 = 22.323 , while the value of the second radical is z r 1 ( a , d , c , d , f ) · ( 4 ) + 1000 = 2.323 . Obviously, the correct sum of the radicals f = 20 can be obtained by taking the second radical with the negative sign.

2.2.2. 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:
x = a · z + b ,     y = c · z + d , σ 1 = x + y ,       σ 2 = x · y .  
Obviously the expressions σ 1 = x + y ,     σ 2 = x · y are the simplest symmetric polynomials. Equation (1) becomes:
σ 1 2 = ( a + c ) · z + b + d + 2 · σ 2 , σ 1 = f ,
from this it derives:
σ 2 = f 2   ( a + c ) · z b d 2 .
Because the radicals must be non-negative,   σ 1 0 and σ 2 0 .
Since σ 1 = f , that requirement σ 1 0   is verified for f 0 .
From the condition σ 2 0 , it follows the range of admissible values for the parameter f :
f ( a + c ) · z + b + d .
Finally, we solve the system given by Equation (1) and inequality (5):
a · z + b + c · z + d = f , f ( a + c ) · z + b + d .
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.

2.3. 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]:
P 1 = U 1 2 R 1 ,     P 2 = U 2 2 R 2 .
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 :
P 1 · R 1 + ( P P 1 ) · R 2 = U .  
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).
The found roots are functions of several parameters P r 1 ( U , P ,   R 1 , R 2 ) and P r 2 ( 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 r 1 ( U , 8000 ,   1 , 2 ) and P r 2 ( U , 8000 ,   1 , 2 ) are plotted as function of the required total voltage U (Figure 7).
Let us verify whether the functions P r 1 ( U , 8000 ,   1 ,   2 ) and P r 2 ( 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 r 1 ( U , 8000 ,   1 ,   2 ) and P r 2 ( U , 8000 ,   1 ,   2 ) instead of the sought variable P 1 (Figure 8). Obviously the solutions P r 1 ( U , 8000 ,   1 ,   2 ) and P r 2 ( 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.
According to the plot:
  • For U < 89.4 both solutions Pr1 (U, 8000, 1, 2) and Pr2 (U, 8000, 1, 2) are false;
  • For 89.4 ≤ U < 126.5 solution Pr1 (U, 8000, 1, 2) is false, while Pr2 (U, 8000, 1, 2) is true;
  • For U ≥ 126.5 both solutions Pr1 (U, 8000, 1, 2) and Pr2 (U, 8000, 1, 2) are true.
Thus, we have ascertained that the solution of Equation (8) shown in Figure 5 includes extraneous roots.

2.3.1. 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 r 1 ( U , 8000 ,   1 ,   2 ) and P r 2 ( U , 8000 ,   1 ,   2 ) are also not symmetric (Figure 7). We introduce the notation:
x = P 1 · R 1 ,     y = ( P P 1 ) · R 2 , σ 1 = x + y ,       σ 2 = x · y .  
Equation (8) becomes:
σ 1 2 = P · R 2 + ( R 1 R 2 ) · P 1 + 2 · σ 2 , σ 1 = U ,
from this it derives:
σ 2 = U 2   P · R 2 ( R 1 R 2 ) · P 1 2 .
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 :
U P · R 2 + ( R 1 R 2 ) · P 1 .
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):
P 1 · R 1 + ( P P 1 ) · R 2 = U , U P · R 2 + ( R 1 R 2 ) · P 1 .
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 r 1 ( U , 8000 ,   1 ,   2 ) and P r 2 ( U , 8000 ,   1 ,   2 ) .

2.3.2. 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 :
P 1 · R 1 + ( P P 1 ) · R 1 = U .  
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:
P 1 · R 1 + ( P P 1 ) · R 1 = U , U P · R 1 .
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 .
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 r 1 ( U , 8000 ,   1 ,   2 ) and P r 2 ( U , 8000 ,   1 ,   2 ) .
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).

3. 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   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.

3.1. Initial Formulation of the Problem

According to the known geometry formulas, the surface area and the volume are determined by the following relations [16]:
S = 2 · π · R · ( L + H 2 + R 2 ) , V = π · R 2 · ( L + 2 3 H ) .
Expressing the variable L through both Equation (16), we obtain the relations:
L S ( R , S , H ) = S 2 · π · R · ( H 2 + R 2 ) 2 · π · R , L V ( R , V , H ) = V 2 · π · H · R 2 3 π · R 2 .
Equating the ratios
L S ( R , S , H ) = L V ( R , V , H ) ,
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 ) .

3.2. Solution of the Problem in the Initial Formulation

3.2.1. 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 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.

3.2.2. 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.
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 (Figure 16 and Figure 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.

3.3. Problem Solution Based on Symmetric Polynomials

3.3.1. 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:
3 · R · S 6 · V 6 · π · R 2 · H 2 + R 2 + 4 · π · H · R 2 = 0 .
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:
2 · H 2 3 · H 2 + R 2 = 6 · V 3 · R · S 2 · π · R 2 ,
The left side of (20) is not symmetric, however, we select the symmetric part in it and analyze the resulting equation:
( 2 · H 2 2 · H 2 + R 2 ) H 2 + R 2 = 6 · V 3 · R · S 2 · π · R 2 .
Since the values of all radicals in (21) must be non-negative, it follows the inequality:
2 · H 2 2 · H 2 + R 2 < 0 .
Consequently, the entire left side of Equation (21) is also negative, which implies the negativity of the right side (21):
6 · V 3 · R · S 2 · π · R 2 < 0 .
Analogously for H < 0 Equation (19) can be rewritten in the form:
2 · H 2 3 · H 2 + R 2 = 6 · V 3 · R · S 2 · π · R 2 .
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):
R > 2 · V S .
Now, by solving the system given by Equation (18) and inequality (25):
L S ( R , S , H ) = L V ( R , V , H ) , R > 2 · V S ,
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) with the correct functions H 1 ( R , 1500 , 1000 ) and H 2 ( R , 1500 , 1000 ) substituted in place of the sought variable H . Comparing Figure 18a with Figure 16 and Figure 18b with Figure 17, we see that the false solutions were eliminated.

3.3.2. 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):
H > 0 , L V ( R , V , H ) = V 2 · π · H · R 2 3 π · R 2 > 0
or
H > 0 , V > 2 · π · H · R 2 3 .
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 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.

4. 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 , x y ) ,
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.

5. 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].

Author Contributions

Investigation, V.O.; conceptualization, V.O. and K.O.; methodology, V.O., I.V. and E.N.; formal analysis, I.V.; visualization, V.O. and M.N.; writing—original draft preparation, I.V., M.N., and E.N.; writing—review and editing, V.O., K.O., M.N., and E.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NRU "MPEI" and the Ministry of Science and Higher Education of the Russian Federation (unique identifier RFMEFI60719X0323).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fre, P.G.; Fedotov, A. Groups and Manifolds: Lectures for Physicists with Examples in Mathematica; Walter de Gruyter: Berlin, Germany, 2018; 475p. [Google Scholar]
  2. Brenner, A.; Shacham, M.; Cutlip, B. Applications of mathematical software packages for modeling and simulations in environmental engineering education. Environ. Model. Softw. 2005, 20, 1307–1313. [Google Scholar] [CrossRef]
  3. Steeb, W.-H.; Tanski, I.; Hardy, Y. Problems and Solutions for Groups, Lie Groups, Lie Algebras with Applications; World Scientific Publishing Company: Singapore, 2012. [Google Scholar] [CrossRef]
  4. Meier, J.; Smith, D. Algebra and symmetry. In Exploring Mathematics: An Engaging Introduction to Proof (Cambridge Mathematical Textbooks); Cambridge University Press: Cambridge, MA, USA, 2017; pp. 229–249. [Google Scholar] [CrossRef]
  5. Goodman, F.M. Algebra: Abstract and Concrete; SemiSimple Press: Iowa City, IA, USA, 2015. [Google Scholar]
  6. Kuang, Y.; Zheng, Y.; Åström, K. Partial Symmetry in Polynomial Systems and Its Applications in Computer Vision. In Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition, Columbus, OH, USA, 23–28 June 2014; pp. 438–445. [Google Scholar] [CrossRef] [Green Version]
  7. Larsson, V.; Åström, K. Uncovering Symmetries in Polynomial Systems. In Proceedings of the 14th European Conference Computer Vision—ECCV 2016, Amsterdam, The Netherlands, 11–14 October 2016; Part III. pp. 252–267. [Google Scholar] [CrossRef]
  8. Maxfield, B. Essential Mathcad for Engineering, Science, and Math ISE, 2nd ed.; Academic Press: Cambridge, MA, USA, 2009. [Google Scholar]
  9. Maxfield, B. Essential PTC Mathcad Prime 3.0. A Guide for New and Current Users; Academic Press: Cambridge, MA, USA, 2013; 584p. [Google Scholar] [CrossRef]
  10. Dugopolski, M. Elementary and Intermediate Algebra, 4th ed.; McGraw-Hill: New York, NY, USA, 2012; 1104p. [Google Scholar]
  11. Aufmann, R.N.; Barker, V.N.; Lockwood, J. Algebra: Introductory and Intermediate, 4th ed.; Houghton Mifflin: Boston, MA, USA, 2007. [Google Scholar]
  12. Bird, J. Electrical Circuit Theory and Technology, 6th ed.; Routledge: London, UK, 2017. [Google Scholar]
  13. Napolitano, J.A. Mathematica Primer for Physicists; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  14. Gilat, A. Matlab: An Introduction with Applications, 6th ed.; Wiley: Hoboken, NJ, USA, 2017. [Google Scholar]
  15. Michalowski, T. Applications of Matlab in Science and Engineering; InTech: London, UK, 2011. [Google Scholar] [CrossRef]
  16. Swokowski, E.; Cole, J.A. Algebra and Trigonometry with Analytic Geometry, 12th ed.; Brooks Cole: Belmont, CA, USA, 2010. [Google Scholar]
  17. Bodrova, E.V.; Golovanova, N.B. Modernization of the higher technical school: Historical experience and prospects. Russ. Technol. J. 2017, 5, 73–97. (In Russian) [Google Scholar] [CrossRef]
  18. Ochkov, V. 25 Problems for STEM Education; Chapman and Hall/CRC: Boca Raton, FL, USA, 2020. [Google Scholar]
Figure 1. Listing of the Mathcad program for the symbolic solution of Equation (1).
Figure 1. Listing of the Mathcad program for the symbolic solution of Equation (1).
Computation 08 00040 g001
Figure 2. Listing of the Mathcad program for the symbolic solution of Equation (1) if a = c .
Figure 2. Listing of the Mathcad program for the symbolic solution of Equation (1) if a = c .
Computation 08 00040 g002
Figure 3. Determining the range of permissible values of the parameter f for Equation (1): (a) c = −4; (b) c = 4.
Figure 3. Determining the range of permissible values of the parameter f for Equation (1): (a) c = −4; (b) c = 4.
Computation 08 00040 g003
Figure 4. Listing of the program for the numerical solution of Equation (1).
Figure 4. Listing of the program for the numerical solution of Equation (1).
Computation 08 00040 g004
Figure 5. Listing of the Mathcad program for the symbolic solution of Equation (8).
Figure 5. Listing of the Mathcad program for the symbolic solution of Equation (8).
Computation 08 00040 g005
Figure 6. Listings of the programs for the symbolic solution of Equation (2) with the packages: (a) Wolfram Mathematica; (b) Matlab.
Figure 6. Listings of the programs for the symbolic solution of Equation (2) with the packages: (a) Wolfram Mathematica; (b) Matlab.
Computation 08 00040 g006
Figure 7. The initial set of the solutions of Equation (8).
Figure 7. The initial set of the solutions of Equation (8).
Computation 08 00040 g007
Figure 8. Determining the range of admissible values of the parameter for Equation (8).
Figure 8. Determining the range of admissible values of the parameter for Equation (8).
Computation 08 00040 g008
Figure 9. Plot of the function σ 2 ( U ) .
Figure 9. Plot of the function σ 2 ( U ) .
Computation 08 00040 g009
Figure 10. The set of the correct solutions of Equation (8).
Figure 10. The set of the correct solutions of Equation (8).
Computation 08 00040 g010
Figure 11. The initial set of the solutions of Equation (14).
Figure 11. The initial set of the solutions of Equation (14).
Computation 08 00040 g011
Figure 12. Listing of the program for the symbolic solution of system (15).
Figure 12. Listing of the program for the symbolic solution of system (15).
Computation 08 00040 g012
Figure 13. Set of the correct solutions of system (15).
Figure 13. Set of the correct solutions of system (15).
Computation 08 00040 g013
Figure 14. Geometrical dimensions of a submarine.
Figure 14. Geometrical dimensions of a submarine.
Computation 08 00040 g014
Figure 15. Listing for the solution of Equation (18).
Figure 15. Listing for the solution of Equation (18).
Computation 08 00040 g015
Figure 16. Initial set of the solutions of Equation (18).
Figure 16. Initial set of the solutions of Equation (18).
Computation 08 00040 g016
Figure 17. Determination of the range of the permissible values of the parameter for Equation (18).
Figure 17. Determination of the range of the permissible values of the parameter for Equation (18).
Computation 08 00040 g017
Figure 18. The correct set of the solutions of Equation (18): (a) solutions H ( R , V , S ) ; (b) values of L S ( R , S , H ) = L V ( R , V , H ) corresponding to the solutions.
Figure 18. The correct set of the solutions of Equation (18): (a) solutions H ( R , V , S ) ; (b) values of L S ( R , S , H ) = L V ( R , V , H ) corresponding to the solutions.
Computation 08 00040 g018
Figure 19. The set of the correct solutions of Equation (18), taking into account the physical meaning: (a) first subdomain of solutions   H ( R , V , S ) ; (b) values of L S ( R , S , H ) = L V ( R , V , H ) , corresponding to the first subdomain; (c) the second subdomain of the solutions H ( R , V , S ) ; (d) values of L S ( R , S , H ) = L V ( R , V , H ) corresponding to the second subdomain.
Figure 19. The set of the correct solutions of Equation (18), taking into account the physical meaning: (a) first subdomain of solutions   H ( R , V , S ) ; (b) values of L S ( R , S , H ) = L V ( R , V , H ) , corresponding to the first subdomain; (c) the second subdomain of the solutions H ( R , V , S ) ; (d) values of L S ( R , S , H ) = L V ( R , V , H ) corresponding to the second subdomain.
Computation 08 00040 g019
Figure 20. A 3d-plot of the solution area ( H , L , R ) : (a) region of the correct solutions; (b) region of the solutions with physical meaning.
Figure 20. A 3d-plot of the solution area ( H , L , R ) : (a) region of the correct solutions; (b) region of the solutions with physical meaning.
Computation 08 00040 g020

Share and Cite

MDPI and ACS Style

Ochkov, V.; Vasileva, I.; Nori, M.; Orlov, K.; Nikulchev, E. Symbolic Computation to Solving an Irrational Equation on Based Symmetric Polynomials Method. Computation 2020, 8, 40. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8020040

AMA Style

Ochkov V, Vasileva I, Nori M, Orlov K, Nikulchev E. Symbolic Computation to Solving an Irrational Equation on Based Symmetric Polynomials Method. Computation. 2020; 8(2):40. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8020040

Chicago/Turabian Style

Ochkov, Valery, Inna Vasileva, Massimiliano Nori, Konstantin Orlov, and Evgeny Nikulchev. 2020. "Symbolic Computation to Solving an Irrational Equation on Based Symmetric Polynomials Method" Computation 8, no. 2: 40. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8020040

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop