## 3 Numeric Experiments and Results Discussion
In the case where heat flux influenced the left, bottom and top of the plate (boundary condition (12), (Table 7c)), almost all inclusions were identified (Table 7f). In Problems 2, 3 and 4, the best identification occurred for the boundary condition (12). It should be noted that changing the inclusions' geometry does not have an effect on the identification process under the boundary condition (12). Problem 6: Identification of the 18 Inclusions in the Form of a Square under the Influence of Temperature and Heat Fluxes  
Let us consider the case when the plate contains 18 inclusions in the form of squares, each with a side size of 0.08 m. The coordinates of these inclusions are the same as the coordinates of the rectangles (Table 1). The boundary conditions and the form of temperature and heat fluxes are the same as in Problem 2. Investigation of the convergence of FEM to solve the direct problem of identifying inclusions in the form of squares were presented in Table 8. The boundary conditions and  
Analysis of the results obtained (Table 8) showed that to identify an area with 18 square inclusions, 1330 finite elements are sufficient, with an error of \(=0.5112\%\). With the influence of the heat flux only on the left border of the plate (boundary condition ((10), (Table 9a)), not all inclusions are fully identified (Table 9d). This is due to the action of heat flux. It should be noted that the top and bottom inclusions are not fully identified. In the case when heat flux is applied to the left and upper boundaries of the plate (boundary condition (11), (Table 9b)), then the identification results (Table 9e) are better in contrast to Problem 2, but the upper and lower inclusions are also not completely identified. In the case when the heat flux is applied to the top, right and left boundaries of the plate (boundary condition (12), (Table 9c)), all inclusions are identified with good accuracy (Table 9f). An analysis of the results, which are presented in Tables 3-5, 7 and 9, leads to the following conclusions. Under the action of the heat flux only on the left boundary of the plate (boundary condition (10), (Tables 3a, 4a, 5a, 7a and 9a)) only the left column and the middle of the middle column were defined (Tables 3d, 4d, 5d, 7d and 9d). It should be noted that the best identification result was obtained for small square inclusions; for rhombuses, it was slightly worse and for large square inclusions, identification was minimal. Thus, the type and size of the hole affects the process of determining the location and type of inclusions. In the case of action of the heat fluxes on the left and upper border of the plate (boundary condition (11), (Tables 3b, 4b, 5b, 7b and 9b) parts of the inclusions in the lower rows were not defined (Tables 3e, 4e, 5e, 7e and 9e). In the case when heat flux acted on the left, upper and lower boundaries of the plate (boundary condition (12), (Tables 3f, 4f, 5f, 7f and 9f), then all inclusions were well-defined. As discussed in the above study, we can conclude that for all studies considered, the identification of all inclusions was fully achieved only when the heat flux was applied to three boundaries and boundary condition (12) was considered.