Thank you, and a graph of the truth has become smoother. The question is precisely what can be wrong? I suspect that in the condition I have given the first derivative "Bessel K" and others. I outlined them in the right code? Another question is this. On my system, the variable "R" changes in the gap [0;25] and the variables "Kz" are changed from [0;3.1]. Tried to build a separate S1-S6 and the program does not give me the addiction. I do not understand, it seems the problem in the code. The question is whether I even brought it? For the first time use wolfram mathematica mistakes can be very trivial. If not difficult, review of the at least one "S" may be wrong, I just laid out code. Kind regards, Denis.

Fine. With the default options which you are using the graphics comes out quite irregularly. But adding the options PlotPoints->100 turns it into a system of smooth curves. This, however, does not match your graphics quantitatively. I guess that your definition of S1 ... S6 does not exactly match the form which underlies your graphics. Hope this helps.