The possibility of using Neumann's method to solve the boundary problems for thin elastic shells is studied. The variational statement of the static problems for the shells allows for a problem examination within the distribution space. The convergence of Neumann's method is proven for the shells with holes when the boundary of the domain is not completely fixed. The numerical implementation of Neumann's method normally requires significant time before any reliable results can be achieved. This paper suggests a way to improve the convergence of the process, and allows for parallel computing and evaluation during the calculations.