{"id":395958,"date":"2024-06-29T12:32:51","date_gmt":"2024-06-29T12:32:51","guid":{"rendered":"http:\/\/savepearlharbor.com\/?p=395958"},"modified":"-0001-11-30T00:00:00","modified_gmt":"-0001-11-29T21:00:00","slug":"","status":"publish","type":"post","link":"https:\/\/savepearlharbor.com\/?p=395958","title":{"rendered":"<span>\u0420\u0438\u0441\u043e\u0432\u0430\u043b\u043a\u0430 \u0434\u043b\u044f \u0430\u0442\u043e\u043c\u043d\u044b\u0445 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439 \u043d\u0430 Python<\/span>"},"content":{"rendered":"<div><!--[--><!--]--><\/div>\n<div id=\"post-content-body\">\n<div>\n<div class=\"article-formatted-body article-formatted-body article-formatted-body_version-2\">\n<div xmlns=\"http:\/\/www.w3.org\/1999\/xhtml\">\n<p>\u041d\u0430\u0447\u0430\u043b\u0441\u044f \u043d\u043e\u0432\u044b\u0439 \u0443\u0447\u0435\u0431\u043d\u044b\u0439 \u0433\u043e\u0434, \u0438 \u043f\u0440\u0435\u043f\u043e\u0434\u0430\u0432\u0430\u0442\u0435\u043b\u044f\u043c, \u0441\u0442\u0443\u0434\u0435\u043d\u0442\u0430\u043c \u0438 \u0448\u043a\u043e\u043b\u044c\u043d\u0438\u043a\u0430\u043c, \u0432\u043e\u0437\u043c\u043e\u0436\u043d\u043e, \u0442\u0440\u0435\u0431\u0443\u0435\u0442\u0441\u044f (\u0438\u043b\u0438 \u043f\u0440\u043e\u0441\u0442\u043e \u0445\u043e\u0447\u0435\u0442\u0441\u044f) \u043f\u043e\u0441\u043c\u043e\u0442\u0440\u0435\u0442\u044c \u043d\u0430 \u0442\u043e, \u043a\u0430\u043a \u0432\u044b\u0433\u043b\u044f\u0434\u044f\u0442 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0438, \u043d\u0430 \u043a\u043e\u0442\u043e\u0440\u044b\u0445 \u0441\u0438\u0434\u044f\u0442 \u044d\u043b\u0435\u043a\u0442\u0440\u043e\u043d\u044b \u0432 \u0430\u0442\u043e\u043c\u0430\u0445: \u0432\u0441\u0435 \u044d\u0442\u0438 \u0437\u0430\u0432\u043e\u0440\u0430\u0436\u0438\u0432\u0430\u044e\u0449\u0438\u0435 \u0431\u0443\u043a\u043e\u0432\u043a\u0438 s, p, d, f, \u0438 \u0442.\u0434. \u0414\u0430, \u043a\u0430\u0440\u0442\u0438\u043d\u043e\u043a \u043f\u043e\u043b\u043d\u043e \u043a\u0430\u043a \u0432 \u0443\u0447\u0435\u0431\u043d\u0438\u043a\u0430\u0445, \u0442\u0430\u043a \u0438 \u0432 \u0418\u043d\u0442\u0435\u0440\u043d\u0435\u0442\u0435, \u043d\u043e \u043f\u043e\u043a\u0440\u0443\u0442\u0438\u0442\u044c \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0438 \u043d\u0430 \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0435 \u043d\u0435 \u043f\u043e\u043b\u0443\u0447\u0438\u0442\u0441\u044f, \u0430 \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0443 \u0438\u0437 \u0443\u0447\u0435\u0431\u043d\u0438\u043a\u0430\/\u0441 \u043b\u0435\u0432\u043e\u0433\u043e \u0441\u0430\u0439\u0442\u0430 \u0432 \u043f\u0440\u0435\u0437\u0435\u043d\u0442\u0430\u0446\u0438\u044e\/\u0440\u0435\u0444\u0435\u0440\u0430\u0442 \u0431\u0435\u0437 \u043c\u043e\u0440\u043e\u043a\u0438 \u0441 \u043b\u0438\u0446\u0435\u043d\u0437\u0438\u0435\u0439 \u043f\u0438\u0445\u0430\u0442\u044c (\u043f\u043e-\u0445\u043e\u0440\u043e\u0448\u0435\u043c\u0443) \u043d\u0435 \u0441\u0442\u043e\u0438\u0442. \u041f\u043e\u044d\u0442\u043e\u043c\u0443 \u0432 \u044d\u0442\u043e\u043c \u043f\u043e\u0441\u0442\u0435 \u043c\u044b \u0440\u0430\u0437\u0431\u0435\u0440\u0451\u043c \u043e\u0434\u043d\u0443 \u0438\u0437 \u0432\u043e\u0437\u043c\u043e\u0436\u043d\u044b\u0445 \u0440\u0435\u0430\u043b\u0438\u0437\u0430\u0446\u0438\u0439 \u0440\u0438\u0441\u043e\u0432\u0430\u043b\u043a\u0438 \u0434\u043b\u044f \u044d\u0442\u0438\u0445 \u0441\u0430\u043c\u044b\u0445 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u043a. <\/p>\n<hr\/>\n<p>\u0421\u0430\u043c\u0443 \u043a\u043e\u043d\u0446\u0435\u043f\u0446\u0438\u044e \u0442\u043e\u0433\u043e, \u0447\u0442\u043e \u0442\u0430\u043a\u043e\u0435 <a href=\"https:\/\/en.wikipedia.org\/wiki\/Atomic_orbital\" rel=\"noopener noreferrer nofollow\">\u0430\u0442\u043e\u043c\u043d\u044b\u0435 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438<\/a> \u044f \u0443\u0436\u0435 \u0440\u0430\u0437\u0431\u0438\u0440\u0430\u043b (<a href=\"https:\/\/habr.com\/ru\/post\/443232\/\" rel=\"noopener noreferrer nofollow\">https:\/\/habr.com\/ru\/post\/443232\/<\/a>), \u0442\u043e \u043a\u0430\u043a \u0440\u0438\u0441\u043e\u0432\u0430\u0442\u044c \u044d\u0442\u0438 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 \u043f\u0440\u0438 \u043f\u043e\u043c\u043e\u0449\u0438 \u043c\u0435\u0442\u043e\u0434\u0430 \u041c\u043e\u043d\u0442\u0435-\u041a\u0430\u0440\u043b\u043e \u0432 \u0432\u0430\u0440\u0438\u0430\u043d\u0442\u0435 <a href=\"https:\/\/ru.wikipedia.org\/wiki\/%D0%90%D0%BB%D0%B3%D0%BE%D1%80%D0%B8%D1%82%D0%BC_%D0%9C%D0%B5%D1%82%D1%80%D0%BE%D0%BF%D0%BE%D0%BB%D0%B8%D1%81%D0%B0_%E2%80%94_%D0%93%D0%B0%D1%81%D1%82%D0%B8%D0%BD%D0%B3%D1%81%D0%B0\" rel=\"noopener noreferrer nofollow\">\u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u0430 \u041c\u0435\u0442\u0440\u043e\u043f\u043e\u043b\u0438\u0441\u0430 <\/a>\u0442\u043e\u0436\u0435 (<a href=\"https:\/\/habr.com\/ru\/post\/548704\/\" rel=\"noopener noreferrer nofollow\">https:\/\/habr.com\/ru\/post\/548704\/<\/a>),<\/p>\n<details class=\"spoiler\">\n<summary>c\u043e\u0431\u0441\u0442\u0432\u0435\u043d\u043d\u043e, \u0432\u043e\u0442 \u0438 \u0441\u0430\u043c \u043a\u043e\u0434 \u0440\u0438\u0441\u043e\u0432\u0430\u043b\u043a\u0438 \u043d\u0430 Python (\u0441 \u043a\u043e\u043c\u043c\u0435\u043d\u0442\u0430\u0440\u0438\u044f\u043c\u0438),<\/summary>\n<div class=\"spoiler__content\">\n<pre><code class=\"python\">#! \/usr\/bin\/python3  import numpy as np  # NumPy \u043d\u0443\u0436\u0435\u043d \u0434\u043b\u044f \u0432\u0441\u0435\u0433\u043e import argparse     # \u042d\u0442\u043e \u043c\u043e\u0434\u0443\u043b\u044c, \u0434\u043b\u044f \"\u043e\u0431\u0449\u0435\u043d\u0438\u044f\" \u0441 \u044e\u0437\u0435\u0440\u043e\u043c    # \u042d\u0442\u0430 \u0444\u0443\u043d\u043a\u0446\u0438\u044f \u0434\u0430\u0451\u0442 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \"\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438\" \u0432 \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u043e\u0439 \u0442\u043e\u0447\u043a\u0435 \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0430 # \u041e\u0440\u0431\u0438\u0442\u0430\u043b\u044c \u044f\u0432\u043b\u044f\u0435\u0442\u0441\u044f \u0441\u0443\u043c\u043c\u043e\u0439 \u0442\u0440\u0451\u0445 \u0443\u0433\u043b\u043e\u0432\u044b\u0445 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442 \u0441 \u0442\u0440\u0435\u043c\u044f \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442\u0430\u043c\u0438, \u0434\u043e\u043c\u043d\u043e\u0436\u0435\u043d\u043d\u043e\u0439 \u043d\u0430 \u0435\u0434\u0438\u043d\u0443\u044e  # \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u0443\u044e \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u0443. # r  -- \u044d\u0442\u043e \u0442\u0440\u0451\u0445\u043c\u0435\u0440\u043d\u044b\u0439 \u0432\u0435\u043a\u0442\u043e\u0440, \u0440\u0430\u0441\u043f\u043e\u043b\u043e\u0436\u0435\u043d\u0438\u0435 \u044d\u043b\u0435\u043a\u0442\u0440\u043e\u043d\u0430 \u0432 \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0435 # n1,n2,n3 -- \u044d\u0442\u043e \u0442\u0440\u0451\u0445\u043c\u0435\u0440\u043d\u044b\u0439 \u0432\u0435\u043a\u0442\u043e\u0440 \u0441\u0442\u0435\u043f\u0435\u043d\u0435\u0439 x,y,z \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442 \u0432\u0435\u043a\u0442\u043e\u0440\u0430 \u0432 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0447\u0430\u0441\u0442\u0438 \u043a\u0430\u0436\u0434\u043e\u0439 \u0438\u0437 \u0442\u0440\u0451\u0445 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442 # coeff1,coeff2,coeff3 -- \u044d\u0442\u043e \u0441\u043e\u0431\u0441\u0442\u0432\u0435\u043d\u043d\u043e \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442\u044b \u043f\u0435\u0440\u0435\u0434 \u0443\u0433\u043b\u043e\u0432\u044b\u043c\u0438 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u0430\u043c\u0438 # r0 -- \u044d\u0442\u043e \"\u0440\u0430\u0437\u043c\u0435\u0440\" \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u043e\u0439 \u0447\u0430\u0441\u0442\u0438, \u0434\u0430\u0432\u0430\u0435\u043c\u043e\u0439 \u0432\u044b\u0440\u0430\u0436\u0435\u043d\u0438\u0435\u043c exp(-|r|\/r0) def orb(r, n1, n2=None, n3=None, coeff1=1., coeff2=1., coeff3=1., r0 = 1.0):      # \u0441\u0447\u0438\u0442\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b Ang     Ang = coeff1 * np.prod(np.power(r,n1))     if not n2 is None:          Ang += coeff2 * np.prod(np.power(r,n2))     if not n3 is None:          Ang += coeff3 * np.prod(np.power(r,n3))      Arg = np.sqrt(np.dot(r,r))\/r0 # \u044d\u0442\u043e \u0430\u0440\u0433\u0443\u043c\u0435\u043d\u0442 \u0434\u043b\u044f \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b      Rad = np.exp(-Arg)   # \u0441\u0447\u0438\u0442\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438          return  Rad * Ang # \u043e\u0442\u0432\u0435\u0442 -- \u044d\u0442\u043e \u043f\u0440\u043e\u0438\u0437\u0432\u0435\u0434\u0435\u043d\u0438\u0435 \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u043e\u0439 \u0438 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442      # \u0434\u0430\u043d\u043d\u0430\u044f \u0444\u0443\u043d\u043a\u0446\u0438\u044f \u043d\u0443\u0436\u043d\u0430, \u0447\u0442\u043e\u0431\u044b \u0441\u0434\u0435\u043b\u0430\u0442\u044c \u043d\u0430\u0437\u0432\u0430\u043d\u0438\u0435 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 # \u043d\u0430\u043c \u043d\u0443\u0436\u043d\u044b \u0442\u043e\u043b\u044c\u043a\u043e \u0441\u0442\u0435\u043f\u0435\u043d\u0438 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b (n) \u0438 \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043f\u0435\u0440\u0435\u0434 \u043d\u0435\u0439 def PowToName(n, coeff=1.):     FunctionName = \"\"      # \u043d\u0435\u0442 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b -- \u043d\u0435\u0442 \u0438 \u0444\u0443\u043d\u043a\u0446\u0438\u0438     if n is None or coeff is None:         return \"\"              # \u0441\u043c\u043e\u0442\u0440\u0438\u043c \u043d\u0430 \u0437\u043d\u0430\u043a \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442\u0430 \u0434\u043b\u044f \u0443\u0433\u043b\u043e\u0432\u043e\u0433\u043e \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u0430         if not coeff is None:         if coeff > 0:             FunctionName += \"+%.2f\"  % coeff         elif coeff &lt;0:             FunctionName += \"%.2f\" % coeff         else:             return \"\" # \u0435\u0441\u043b\u0438 \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043d\u0443\u043b\u0435\u0432\u043e\u0439, \u044d\u0442\u0430 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u0430 \u043d\u0438\u0447\u0435\u0433\u043e \u0438 \u043d\u0435 \u0434\u0430\u0451\u0442          # \u0421\u0442\u0435\u043f\u0435\u043d\u044c \u043f\u043e\u043b\u0438\u043d\u043e\u043c\u0430 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b -- \u044d\u0442\u043e \u0442\u0438\u043f \u0444\u0443\u043d\u043a\u0446\u0438\u0438 (s\/p\/d\/f...), \u0435\u0433\u043e \u043c\u044b \u0438 \u043f\u0440\u043e\u0432\u0435\u0440\u044f\u0435\u043c     if np.sum(n) == 0:         FunctionName += \"s\"     elif np.sum(n) == 1:         FunctionName += \"p\"     elif np.sum(n) == 2:         FunctionName += \"d\"     elif np.sum(n) == 3:         FunctionName += \"f\"     elif np.sum(n) == 4:         FunctionName += \"g\"     elif np.sum(n) == 5:         FunctionName += \"h\"     else:         FunctionName += \"X3\" # \u043c\u043e\u0436\u043d\u043e \u043f\u043e\u0439\u0442\u0438 \u0438 \u0434\u0430\u043b\u044c\u0448\u0435, \u043d\u043e \u043a\u043e\u043c\u0443 \u044d\u0442\u043e \u043d\u0443\u0436\u043d\u043e?               # \u0438 \u043f\u043e\u043a\u0430\u0437\u044b\u0432\u0430\u0435\u043c \u043a\u0430\u043a\u0430\u044f \u043a\u043e\u043d\u043a\u0440\u0435\u0442\u043d\u043e \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u0430\u044f \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u0430 \u0443 \u043d\u0430\u0441 \u0442\u0443\u0442 \u0435\u0441\u0442\u044c (\u0435\u0441\u0442\u0435\u0441\u0442\u0432\u0435\u043d\u043d\u043e, \u043d\u0443\u043b\u0435\u0432\u044b\u0435 \u0441\u0442\u0435\u043f\u0435\u043d\u0438 \u043c\u044b \u043d\u0435 \u043f\u0435\u0447\u0430\u0442\u0430\u0435\u043c)     if n[0]>0:         FunctionName += \"x\"         if n[0] > 1:             FunctionName += \"%i\" %  (n[0])     if n[1]>0:         FunctionName += \"y\"         if n[1] > 1:             FunctionName += \"%i\" %  (n[1])     if n[2]>0:         FunctionName += \"z\"         if n[2] > 1:             FunctionName += \"%i\" %  (n[2])                     return FunctionName       # argument parser to get options from command line parser = argparse.ArgumentParser(description='\u0414\u0430\u0432\u0430\u0439\u0442\u0435-\u043a\u0430 \u043d\u0430\u0440\u0438\u0441\u0443\u0435\u043c \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0438 \u043f\u0440\u0438 \u043f\u043e\u043c\u043e\u0449\u0438 \u043c\u0435\u0442\u043e\u0434\u0430 \u041c\u043e\u043d\u0442\u0435-\u041a\u0430\u043b\u043e (\u041c\u041a)!')    parser.add_argument('--c1', type=float, default=1.0,                      help=\"\u041a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043f\u0435\u0440\u0435\u0434 \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0435\u0439\")     parser.add_argument('--c2', type=float, default=0.0,                      help=\"\u041a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043f\u0435\u0440\u0435\u0434 \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0435\u0439\")     parser.add_argument('--c3', type=float, default=0.0,                      help=\"\u041a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043f\u0435\u0440\u0435\u0434 \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0435\u0439\")    parser.add_argument('--nx1', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 x**n \u0434\u043b\u044f \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--nx2', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 x**n \u0434\u043b\u044f \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--nx3', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 x**n \u0434\u043b\u044f \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--ny1', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 y**n \u0434\u043b\u044f \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--ny2', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 y**n \u0434\u043b\u044f \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--ny3', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 y**n \u0434\u043b\u044f \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--nz1', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 z**n \u0434\u043b\u044f \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--nz2', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 z**n \u0434\u043b\u044f \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--nz3', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 z**n \u0434\u043b\u044f \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--NumSteps', type=int, default=50000,                      help=\"\u0427\u0438\u0441\u043b\u043e \u0448\u0430\u0433\u043e\u0432 \u0432 \u041c\u041a \u0441\u0438\u043c\u0443\u043b\u044f\u0446\u0438\u0438\")   parser.add_argument('--R0', type=float, default=1.0,                      help=\"\u0420\u0430\u0437\u043c\u0435\u0440 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 R0\")   parser.add_argument('--PartToIgnore', type=float, default=0.1,                      help=\"\u0414\u043e\u043b\u044f \u041c\u041a \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438 \u0432 \u043d\u0430\u0447\u0430\u043b\u0435, \u043a\u043e\u0442\u043e\u0440\u0430\u044f \u0438\u0433\u043d\u043e\u0440\u0438\u0440\u0443\u0435\u0442\u0441\u044f \u0432 \u0440\u0430\u0441\u0447\u0451\u0442\u0430\u0445, \u043a\u0430\u043a \u0441\u0442\u0430\u0434\u0438\u044f \u0443\u0440\u0430\u0432\u043d\u043e\u0432\u0435\u0448\u0438\u0432\u0430\u043d\u0438\u044f\")    parser.add_argument('--MaxBox', type=float, default=0.1,                      help=\"\u0420\u0430\u0437\u043c\u0435\u0440 \u0441\u0442\u0430\u0440\u0442\u043e\u0432\u043e\u0439 \u043e\u0431\u043b\u0430\u0441\u0442\u0438, \u0438\u0437 \u043a\u043e\u0442\u043e\u0440\u043e\u0439 \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u043e \u0432\u044b\u0431\u0438\u0440\u0430\u0435\u0442\u0441\u044f \u0441\u0442\u0430\u0440\u0442\u043e\u0432\u0430\u044f \u0442\u043e\u0447\u043a\u0430 \u041c\u041a-\u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438\")    parser.add_argument('--MaxStep', type=float, default=2.0,                      help=\"\u041c\u0430\u043a\u0441\u0438\u043c\u0430\u043b\u044c\u043d\u0430\u044f \u0434\u043b\u0438\u043d\u0430 \u0448\u0430\u0433\u0430 \u0432 \u041c\u041a-\u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438\")    parser.add_argument('--MaxSmallShift', type=float, default=0.1,                      help=\"\u041c\u0430\u043b\u0435\u043d\u044c\u043a\u0430\u044f \u0434\u043e\u0431\u0430\u0432\u043e\u0447\u043a\u0430 \u043a \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438, \u043a\u043e\u0433\u0434\u0430 \u043c\u044b \u043e\u0441\u0442\u0430\u0451\u043c\u0441\u044f \u0432 \u0442\u043e\u0439 \u0436\u0435 \u0442\u043e\u0447\u043a\u0435 (\u0434\u043b\u044f \u043a\u0440\u0430\u0441\u043e\u0442\u044b)\")    parser.add_argument('--plot',action='store_true',                      help=\"\u0424\u043b\u0430\u0433, \u0447\u0442\u043e\u0431\u044b \u043d\u0430\u0440\u0438\u0441\u043e\u0432\u0430\u0442\u044c \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0443 \u043f\u0440\u0438 \u043f\u043e\u043c\u043e\u0449\u0438 matplotlib\")     # \u0430 \u0442\u0435\u043f\u0435\u0440\u044c \u0440\u0430\u0441\u043f\u0430\u0440\u0441\u0438\u043c \u0430\u0440\u0433\u0443\u043c\u0435\u043d\u0442\u044b args,unknown = parser.parse_known_args()   # \u0441\u043e\u0441\u0442\u0430\u0432\u043b\u044f\u0435\u043c \u0432\u0435\u043a\u0442\u043e\u0440\u0430 \u0438\u0437 \u0441\u0442\u0435\u043f\u0435\u043d\u0435\u0439 \u0434\u043b\u044f \u043a\u0430\u0436\u0434\u043e\u0439 \u0438\u0437 \u0442\u0440\u0451\u0445 \u0443\u0433\u043b\u043e\u0432\u044b\u0445 \u0444\u0443\u043d\u043a\u0446\u0438\u0439 \u0438\u0437 \u0442\u043e\u0433\u043e, \u0447\u0442\u043e \u0432\u0432\u0451\u043b \u044e\u0437\u0435\u0440 NPow1 = np.array([args.nx1,args.ny1,args.nz1] ,dtype=int) NPow2 = np.array([args.nx2,args.ny2,args.nz2] ,dtype=int) NPow3 = np.array([args.nx3,args.ny3,args.nz3] ,dtype=int)               # \u0432\u044b\u0431\u0438\u0440\u0430\u0435\u043c \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u0443\u044e \u0441\u0442\u0430\u0440\u0442\u043e\u0432\u0443\u044e \u0442\u043e\u0447\u043a\u0443 \u0432 \u0437\u0430\u0434\u0430\u043d\u043d\u043e\u0439 \u044e\u0437\u0435\u0440\u043e\u043c \u043e\u0431\u043b\u0430\u0441\u0442\u0438 \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0430 (\u043a\u0443\u0431\u0438\u043a\u0435 \u0432\u043e\u043a\u0440\u0443\u0433 r=(0,0,0) \u0441\u043e \u0441\u0442\u043e\u0440\u043e\u043d\u043e\u0439 2*MaxBox) currentR = np.random.uniform(low=-args.MaxBox, high=args.MaxBox, size=3)  # \u0441\u0447\u0438\u0442\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043e\u0440\u0438\u0431\u0438\u0442\u0430\u043b\u0438 \u0432 \u044d\u0442\u043e\u0439 \u0442\u043e\u0447\u043a\u0435 currentOrbP = orb(r=currentR, n1=NPow1, n2=NPow2, n3=NPow3, coeff1=args.c1, coeff2=args.c2, coeff3=args.c3, r0 = args.R0)  # \u0432 \u044d\u0442\u043e\u0442 \u0444\u0430\u0439\u043b \u043c\u044b \u0431\u0443\u0434\u0435\u043c \u0441\u043e\u0445\u0440\u0430\u043d\u044f\u0442\u044c \u0441\u0430\u043c\u0443 \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u044e FileName=\"orbital_%s%s%s.dat\" % (PowToName(NPow1, args.c1),PowToName(NPow2, args.c2),PowToName(NPow3, args.c3)) outf = open(FileName, \"w\")  AccRate = 0 # \u042d\u0442\u043e \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u043f\u0440\u0438\u043d\u0438\u043c\u0430\u0435\u043c\u044b\u0445 \u0448\u0430\u0433\u043e\u0432 \u0432 \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438, \u0447\u0442\u043e\u0431\u044b \u0432\u044b\u0447\u0438\u0441\u043b\u0438\u0442\u044c \u0434\u043e\u043b\u044e \u043f\u0440\u0438\u043d\u044f\u0442\u0438\u044f (\u0434\u043e\u043b\u0436\u043d\u0430 \u0431\u044b\u0442\u044c \u0432 \u0434\u0438\u0430\u043f\u0430\u0437\u043e\u043d\u0435 20-60 %)  # \u043d\u0430\u0447\u0438\u043d\u0430\u0435\u043c \u0438\u0442\u0435\u0440\u0430\u0446\u0438\u0438 \u043c\u0435\u0442\u043e\u0434\u0430 \u041c\u0435\u0442\u0440\u043e\u043f\u043e\u043b\u0438\u0441\u0430  for nstep in range(0,args.NumSteps):     # \u0433\u0435\u043d\u0435\u0440\u0438\u0440\u0443\u0435\u043c \u043d\u043e\u0432\u0443\u044e \u0442\u043e\u0447\u043a\u0443, \u043f\u0440\u0438\u043f\u0438\u0441\u044b\u0432\u0430\u044f \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u043e\u0435 \u0441\u043c\u0435\u0449\u0435\u043d\u0438\u0435 \u043a \u043d\u044b\u043d\u0435\u0448\u043d\u0435\u043c\u0443 \u043f\u043e\u043b\u043e\u0436\u0435\u043d\u0438\u044e \u0432 \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0435     trialR = currentR + np.random.uniform(low=-args.MaxStep, high=args.MaxStep, size=3)      # \u0441\u0447\u0438\u0442\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 \u0432 \u044d\u0442\u043e\u0439 \u0442\u043e\u0447\u043a\u0435     trialOrbP = orb(r=trialR, n1=NPow1, n2=NPow2, n3=NPow3, coeff1=args.c1, coeff2=args.c2, coeff3=args.c3, r0 = args.R0)             # \u0432\u044b\u0431\u0440\u0430\u0441\u044b\u0432\u0430\u0435\u043c \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u044c \u043f\u0440\u0438\u043d\u044f\u0442\u0438\u044f \u043d\u043e\u0432\u043e\u0439 \u0442\u043e\u0447\u043a\u0438     Ptrial = np.random.uniform(low=0, high=1)          # \u0430 \u0437\u0434\u0435\u0441\u044c \u043c\u044b \u0441\u043e\u0445\u0440\u0430\u043d\u0438\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0442\u043e\u0433\u043e, \u043f\u0440\u0438\u043c\u0435\u043c \u043b\u0438 \u043c\u044b \u043d\u043e\u0432\u0443\u044e \u0442\u043e\u0447\u043a\u0443, \u0438\u043b\u0438 \u043d\u0435\u0442     DoWeAccept = False          # \u0441\u043d\u0430\u0447\u0430\u043b\u0430 \u043f\u0440\u043e\u0441\u0442\u043e \u043f\u0440\u043e\u0432\u0435\u0440\u0438\u043c, \u0447\u0442\u043e \u043d\u043e\u0432\u0430\u044f \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u044c (\u043a\u0432\u0430\u0434\u0440\u0430\u0442 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438) \u0432\u044b\u0448\u0435, \u0447\u0435\u043c \u0431\u044b\u043b\u0430, \u0438 \u0435\u0441\u043b\u0438 \u0434\u0430, \u0442\u043e \u043f\u0440\u0438\u043d\u0438\u043c\u0430\u0435\u043c \u0442\u043e\u0447\u043a\u0443     if (trialOrbP**2 > currentOrbP**2):         DoWeAccept = True     else:         # \u0435\u0441\u043b\u0438 \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u044c \u043d\u043e\u0432\u043e\u0439 \u0442\u043e\u0447\u043a\u0438 \u043d\u0438\u0436\u0435, \u0441\u0440\u0430\u0432\u043d\u0438\u0432\u0430\u0435\u043c \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u044c \u043f\u0435\u0440\u0435\u0445\u043e\u0434\u0430 \u0432 \u044d\u0442\u0443 \u0442\u043e\u0447\u043a\u0443 \u0441 \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u044c\u044e \u043f\u0440\u0438\u043d\u044f\u0442\u0438\u044f         if trialOrbP**2\/currentOrbP**2 > Ptrial:             DoWeAccept = True          # \u0435\u0441\u043b\u0438 \u043d\u043e\u0432\u0443\u044e \u0442\u043e\u0447\u043a\u0443 \u043f\u0440\u0438\u043d\u0438\u043c\u0430\u0435\u043c, \u0442\u043e \u043e\u043d\u0430 \u0441\u0442\u0430\u043d\u043e\u0432\u0438\u0442\u0441\u044f \u043d\u043e\u0432\u044b\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0435\u043c \u0442\u043e\u0447\u043a\u0438 \u0432 \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438, \u0430 \u0435\u0449\u0451 \u0443\u0432\u0435\u043b\u0438\u0447\u0438\u0432\u0430\u0435\u043c \u0447\u0438\u0441\u043b\u043e \u043f\u0440\u0438\u043d\u044f\u0442\u044b\u0445 \u0442\u043e\u0447\u0435\u043a     if DoWeAccept:         currentR = trialR         currentOrbP = trialOrbP         AccRate += 1         shift = np.zeros(currentR.shape) # \u0432 \u044d\u0442\u043e\u043c \u0441\u043b\u0443\u0447\u0430\u0435 \u0441\u0434\u0432\u0438\u0433 \u043d\u0438\u043a\u0430\u043a\u043e\u0439 \u043d\u0430\u043c \u043d\u0435 \u043d\u0443\u0436\u0435\u043d     else:         # \u0435\u0441\u043b\u0438 \u0442\u043e\u0447\u043a\u0443 \u043d\u0435 \u043f\u0440\u0438\u043d\u044f\u043b\u0438, \u0447\u0442\u043e\u0431\u044b \u0431\u044b\u043b\u043e \u043d\u0435 \u043e\u0447\u0435\u043d\u044c \u0443\u043d\u044b\u043b\u043e, \u0441\u043e\u0445\u0440\u0430\u043d\u0438\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0441\u0442\u0430\u0440\u043e\u0439 \u0442\u043e\u0447\u043a\u0438 \u0441 \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u044b\u043c \u043c\u0430\u043b\u044b\u043c \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u044b\u043c \u0441\u0434\u0432\u0438\u0433\u043e\u043c         shift = np.random.uniform(low=-args.MaxSmallShift, high=args.MaxSmallShift, size=3)          # \u0438 \u0435\u0441\u043b\u0438 \u0447\u0430\u0441\u0442\u044c \u0438\u0433\u043d\u043e\u0440\u0438\u0440\u043e\u0432\u0430\u043d\u0438\u044f \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438 \u0437\u0430\u043a\u043e\u043d\u0447\u0438\u043b\u0430\u0441\u044c, \u0442\u043e \u0441\u043e\u0445\u0440\u0430\u043d\u044f\u0435\u043c \u043d\u043e\u0432\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435         if nstep > int(args.NumSteps * args.PartToIgnore):         outf.write(\" %15.10f %15.10f %15.10f \" % tuple(currentR + shift))         outf.write(\"  %15.5e   %10i\\n\" % (currentOrbP, nstep))             outf.flush()           # \u0432 \u043a\u043e\u043d\u0446\u0435 \u0432\u044b\u0434\u0430\u0451\u043c \u043f\u0440\u043e\u0446\u0435\u043d\u0442 \u043f\u0440\u0438\u043d\u044f\u0442\u044b\u0445 \u0442\u043e\u0447\u0435\u043a \u0432 \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u044e (\u0445\u043e\u0442\u0435\u043b\u043e\u0441\u044c \u0431\u044b, \u0447\u0442\u043e\u0431\u044b \u043e\u043d\u043e \u0431\u044b\u043b\u043e \u0432 \u0440\u0430\u0439\u043e\u043d\u0435 20-60 %) print(\"Acc rate = %15.10f %%\" % (100.*float(AccRate)\/float(args.NumSteps)) ) outf.close()       ############## \u0410 \u0442\u0435\u043f\u0435\u0440\u044c \u043f\u043e\u0441\u0442\u0440\u043e\u0438\u043c \u0432\u0441\u044e \u044d\u0442\u0443 \u043a\u0440\u0430\u0441\u043e\u0442\u0443, \u0435\u0441\u043b\u0438 \u044e\u0437\u0435\u0440 \u0437\u0430\u0445\u043e\u0442\u0435\u043b if args.plot:     import matplotlib.pyplot as plt     from mpl_toolkits.mplot3d import Axes3D          fig = plt.figure()     ax = fig.add_subplot(projection='3d')          data = np.loadtxt(FileName)     ax.scatter(data[:,0], data[:,1], data[:,2], c=data[:,3],cmap = 'coolwarm',s=args.MaxSmallShift*0.1)     ax.set_xlabel('X')     ax.set_ylabel('Y')     ax.set_zlabel('Z')          plt.show()  <\/code><\/pre>\n<\/div>\n<\/details>\n<p>\u0441\u043f\u0430\u0441\u0438\u0431\u043e \u0437\u0430 \u0432\u043d\u0438\u043c\u0430\u043d\u0438\u0435. <\/p>\n<p>\u041b\u0430\u0434\u043d\u043e, \u0448\u0443\u0447\u0443, \u0437\u0434\u0435\u0441\u044c \u043c\u044b \u0435\u0449\u0451 \u0440\u0430\u0437 \u043d\u0435\u043c\u043d\u043e\u0433\u043e \u043f\u0440\u043e\u0439\u0434\u0451\u043c\u0441\u044f \u043f\u043e \u0442\u043e\u043c\u0443, \u0447\u0442\u043e \u0442\u0430\u043a\u043e\u0435 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438, \u043a\u0430\u043a\u0438\u0435 \u043e\u043d\u0438 \u0431\u044b\u0432\u0430\u044e\u0442, \u0438 \u043a\u0430\u043a \u0438\u0445 \u0441\u0442\u0440\u043e\u0438\u0442\u044c \u0441 \u043f\u043e\u043c\u043e\u0449\u044c\u044e \u0432\u044b\u0448\u0435\u043f\u0440\u0438\u0432\u0435\u0434\u0451\u043d\u043d\u043e\u0433\u043e \u0441\u043a\u0440\u0438\u043f\u0442\u0430. <\/p>\n<h2>\u041a\u0430\u043a \u0432\u044b\u0433\u043b\u044f\u0434\u044f\u0442 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438?<\/h2>\n<p>\u0410\u0442\u043e\u043c\u043d\u0430\u044f \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c &#8212; \u044d\u0442\u043e \u0432\u043e\u043b\u043d\u043e\u0432\u0430\u044f \u0444\u0443\u043d\u043a\u0446\u0438\u044f \u0434\u043b\u044f \u0432\u044b\u0431\u0440\u0430\u043d\u043d\u043e\u0433\u043e \u044d\u043b\u0435\u043a\u0442\u0440\u043e\u043d\u0430 \u0432 \u0430\u0442\u043e\u043c\u0435, \u043f\u043e\u044d\u0442\u043e\u043c\u0443 \u0431\u0443\u0434\u0435\u043c \u043e\u0431\u043e\u0437\u043d\u0430\u0447\u0430\u0442\u044c \u0438\u0445 \u0431\u0443\u043a\u043e\u0432\u043a\u043e\u0439 <em>\u03c8<\/em>. \u0421\u043e\u0433\u043b\u0430\u0441\u043d\u043e \u043a\u0432\u0430\u043d\u0442\u043e\u0432\u043e\u0439 \u043c\u0435\u0445\u0430\u043d\u0438\u043a\u0435, \u043c\u044b <a href=\"https:\/\/ru.wikipedia.org\/wiki\/%D0%9F%D1%80%D0%B8%D0%BD%D1%86%D0%B8%D0%BF_%D0%BD%D0%B5%D0%BE%D0%BF%D1%80%D0%B5%D0%B4%D0%B5%D0%BB%D1%91%D0%BD%D0%BD%D0%BE%D1%81%D1%82%D0%B8\" rel=\"noopener noreferrer nofollow\">\u043d\u0435 \u043c\u043e\u0436\u0435\u043c \u0437\u043d\u0430\u0442\u044c \u0433\u0434\u0435 \u044d\u043b\u0435\u043a\u0442\u0440\u043e\u043d \u043d\u0430\u0445\u043e\u0434\u0438\u0442\u0441\u044f \u0438 \u043a\u0430\u043a \u0431\u044b\u0441\u0442\u0440\u043e \u043e\u043d \u043a\u0443\u0434\u0430-\u0442\u043e \u043b\u0435\u0442\u0438\u0442,<\/a> \u043d\u043e \u043d\u0430\u043c \u0434\u043e\u0441\u0442\u0443\u043f\u043d\u043e \u0437\u043d\u0430\u043d\u0438\u0435 (\u043d\u0430\u043f\u0440\u0438\u043c\u0435\u0440) \u043e \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u0438 \u043d\u0430\u0445\u043e\u0436\u0434\u0435\u043d\u0438\u044f \u044d\u0442\u043e\u0433\u043e \u044d\u043b\u0435\u043a\u0442\u0440\u043e\u043d\u0430 \u0432 \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0435. \u0418 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 \u044d\u0442\u043e \u0437\u043d\u0430\u043d\u0438\u0435 \u043d\u0430\u043c \u0438 \u043f\u0440\u0435\u0434\u043e\u0441\u0442\u0430\u0432\u043b\u044f\u044e\u0442: \u043d\u0430\u0439\u0442\u0438 \u044d\u043b\u0435\u043a\u0442\u0440\u043e\u043d \u0432 \u043c\u0435\u0441\u0442\u0435 \u0441 \u043a\u043e\u043e\u0440\u0434\u0438\u043d\u0430\u0442\u043e\u0439 <em>x<\/em> (<em>p<\/em>(<em>x<\/em>)) \u0434\u0430\u0451\u0442\u0441\u044f \u043a\u0432\u0430\u0434\u0440\u0430\u0442\u043e\u043c \u0432\u043e\u043b\u043d\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438: <em>p<\/em>(<em>x<\/em>)=|<em>\u03c8<\/em>(<em>x<\/em>)|<sup>2<\/sup>. \u0421\u0430\u043c\u0438 \u0430\u0442\u043e\u043c\u043d\u044b\u0435 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 \u043c\u043e\u0436\u043d\u043e \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u0438\u0442\u044c \u0432 \u0432\u0438\u0434\u0435 \u0434\u0432\u0443\u0445 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442: \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u043e\u0439 (<em>R<\/em>, \u0437\u0430\u0432\u0438\u0441\u044f\u0449\u0435\u0439 \u0442\u043e\u043b\u044c\u043a\u043e \u043e\u0442 \u0442\u043e\u0433\u043e, \u043a\u0430\u043a \u0434\u0430\u043b\u0435\u043a\u043e \u043e\u0442 \u044f\u0434\u0440\u0430 \u0443\u0431\u0435\u0436\u0430\u043b \u044d\u043b\u0435\u043a\u0442\u0440\u043e\u043d) \u0438 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 (<em>A<\/em>, \u0437\u0430\u0432\u0438\u0441\u044f\u0449\u0435\u0439 \u0442\u043e\u043b\u044c\u043a\u043e \u043e\u0442 \u0442\u043e\u0433\u043e, \u0441 \u043a\u0430\u043a\u043e\u0439 \u0438\u043c\u0435\u043d\u043d\u043e \u0441\u0442\u043e\u0440\u043e\u043d\u044b \u043e\u0442 \u044f\u0434\u0440\u0430 \u043d\u0430\u0445\u043e\u0434\u0438\u0442\u0441\u044f \u043d\u0430\u0448\u0430 \u0447\u0430\u0441\u0442\u0438\u0446\u0430), \u0442\u0430\u043a \u0447\u0442\u043e <em>\u03c8<\/em>=<em>R<\/em>\u00b7<em>A<\/em>.<\/p>\n<p>\u0414\u043b\u044f \u043d\u0430\u0448\u0438\u0445 \u0446\u0435\u043b\u0435\u0439 \u0432\u0438\u0437\u0443\u0430\u043b\u0438\u0437\u0430\u0446\u0438\u0438 \u043c\u044b \u0432\u043f\u043e\u043b\u043d\u0435 \u043c\u043e\u0436\u0435\u043c \u0441\u0435\u0431\u0435 \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u0438\u0442\u044c \u043a\u0430\u0436\u0434\u0443\u044e \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c \u0432 \u0432\u0438\u0434\u0435:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"formula\" source=\"\\psi(x,y,z) = \\underbrace{\\exp\\left(-\\frac{\\sqrt{x^2 + y^2 +z^2}}{R_0}\\right)}_{R} \\cdot \\underbrace{x^{n_x} \\cdot y^{n_y} \\cdot z^{n_z}}_{A}\" alt=\"\\psi(x,y,z) = \\underbrace{\\exp\\left(-\\frac{\\sqrt{x^2 + y^2 +z^2}}{R_0}\\right)}_{R} \\cdot \\underbrace{x^{n_x} \\cdot y^{n_y} \\cdot z^{n_z}}_{A}\" src=\"https:\/\/habrastorage.org\/getpro\/habr\/upload_files\/fae\/2ea\/e79\/fae2eae79675669cc5869c097856f178.svg\" width=\"433\" height=\"91\"\/><\/p>\n<p>\u0437\u0434\u0435\u0441\u044c \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u0430\u044f \u0447\u0430\u0441\u0442\u044c \u0434\u0430\u043d\u0430 \u043a\u0430\u043a \u044d\u043a\u0441\u043f\u043e\u043d\u0435\u043d\u0446\u0438\u0430\u043b\u044c\u043d\u043e \u0443\u0431\u044b\u0432\u0430\u044e\u0449\u0430\u044f \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u044c \u043d\u0430\u0439\u0442\u0438 \u044d\u043b\u0435\u043a\u0442\u0440\u043e\u043d \u0432\u0434\u0430\u043b\u0438 \u043e\u0442 \u044f\u0434\u0440\u0430, \u0433\u0434\u0435 \u043f\u0430\u0440\u0430\u043c\u0435\u0442\u0440 <em>R<\/em><sub>0<\/sub> \u043a\u043e\u043d\u0442\u0440\u043e\u043b\u0438\u0440\u0443\u0435\u0442 &#171;\u0440\u0430\u0437\u043c\u0435\u0440&#187; \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438, \u0430 \u0443\u0433\u043b\u043e\u0432\u0430\u044f \u0447\u0430\u0441\u0442\u044c \u0434\u0430\u043d\u0430 \u043a\u0430\u043a \u043f\u0440\u043e\u0438\u0437\u0432\u0435\u0434\u0435\u043d\u0438\u0435 \u043a\u043e\u043e\u0440\u0434\u0438\u043d\u0430\u0442 \u0432 \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u044b\u0445 \u0441\u0442\u0435\u043f\u0435\u043d\u044f\u0445 <em>n<\/em>\u22650. \u0421\u043e\u0431\u0441\u0442\u0432\u0435\u043d\u043d\u043e, \u0442\u0438\u043f \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438, \u043a\u043e\u0442\u043e\u0440\u044b\u0439 \u0434\u0430\u0451\u0442\u0441\u044f \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043d\u044b\u043c \u043a\u0432\u0430\u043d\u0442\u043e\u0432\u044b\u043c \u0447\u0438\u0441\u043b\u043e\u043c <em>l<\/em> \u0432 \u0434\u0430\u043d\u043d\u043e\u043c \u043f\u0440\u0435\u0434\u0441\u0442\u0430\u0432\u043b\u0435\u043d\u0438\u0438 &#8212; \u044d\u0442\u043e \u043f\u0440\u043e\u0441\u0442\u043e \u0441\u0443\u043c\u043c\u0430 \u0432\u0441\u0435\u0445 \u0441\u0442\u0435\u043f\u0435\u043d\u0435\u0439 \u0434\u043b\u044f x,y,z \u043a\u043e\u043e\u0440\u0434\u0438\u043d\u0430\u0442:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"formula\" source=\"l = n_x + n_y + n_z\" alt=\"l = n_x + n_y + n_z\" src=\"https:\/\/habrastorage.org\/getpro\/habr\/upload_files\/885\/6ef\/72d\/8856ef72d8ae645bd8a99f05bf1cd967.svg\" width=\"143\" height=\"22\"\/><\/p>\n<p>\u0422.\u0435. \u0434\u043b\u044f s-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439 (<em>l<\/em>=0) \u043d\u0430\u043c \u043d\u0443\u0436\u043d\u044b <em>n<sub>x<\/sub><\/em>=<em>n<sub>y<\/sub><\/em>=<em>n<sub>z<\/sub><\/em>=0, \u0434\u043b\u044f p-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439, \u043a\u043e\u0442\u043e\u0440\u044b\u0445 \u0443 \u043d\u0430\u0441 \u0442\u0440\u0438 \u0448\u0442\u0443\u043a\u0438, \u043d\u0430\u043c \u043d\u0443\u0436\u043d\u044b \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0438\u0435 \u043a\u043e\u043c\u0431\u0438\u043d\u0430\u0446\u0438\u0438:<\/p>\n<ul>\n<li>\n<p><em>n<sub>x<\/sub><\/em>=1, <em>n<sub>y<\/sub><\/em>=<em>n<sub>z<\/sub><\/em>=0 (p<sub>x<\/sub>-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c),<\/p>\n<\/li>\n<li>\n<p><em>n<sub>y<\/sub><\/em>=1, <em>n<sub>x<\/sub><\/em>=<em>n<sub>z<\/sub><\/em>=0 (p<sub>y<\/sub>-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c),<\/p>\n<\/li>\n<li>\n<p><em>n<sub>z<\/sub><\/em>=1, <em>n<sub>x<\/sub><\/em>=<em>n<sub>y<\/sub><\/em>=0 (p<sub>z<\/sub>-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c).<\/p>\n<\/li>\n<\/ul>\n<p>\u041d\u043e \u0434\u0430\u043b\u044c\u0448\u0435 \u0432\u0441\u0451 \u0431\u0443\u0434\u0435\u0442 \u043f\u043b\u043e\u0445\u043e. \u0414\u043b\u044f \u043a\u0430\u0436\u0434\u043e\u0433\u043e \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043d\u043e\u0433\u043e \u043a\u0432\u0430\u043d\u0442\u043e\u0432\u043e\u0433\u043e \u0447\u0438\u0441\u043b\u0430 <em>l<\/em> \u0443 \u043d\u0430\u0441 \u0434\u043e\u043b\u0436\u043d\u0430 \u0431\u044b\u0442\u044c <em>N<\/em>=2<em>l<\/em>+1 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c, \u0438 \u0434\u043b\u044f s \u0438 p \u0443 \u043d\u0430\u0441 \u0432\u0441\u0451 \u0443\u0434\u043e\u0432\u043b\u0435\u0442\u0432\u043e\u0440\u044f\u043b\u043e \u044d\u0442\u043e\u043c\u0443 \u0441\u043e\u043e\u0442\u043d\u043e\u0448\u0435\u043d\u0438\u044e. \u0410 \u0432\u043e\u0442 \u0434\u0430\u043b\u044c\u0448\u0435, \u0438\u0437 \u0443\u0441\u043b\u043e\u0432\u0438\u044f <em>l<\/em>=<em>n<sub>x<\/sub><\/em>+<em>n<sub>y<\/sub><\/em>+<em>n<sub>z<\/sub><\/em> \u0443 \u043d\u0430\u0441 \u0431\u0443\u0434\u0435\u0442 \u043f\u043e\u043b\u0443\u0447\u0430\u0442\u044c\u0441\u044f \u0441\u043b\u0438\u0448\u043a\u043e\u043c \u043c\u043d\u043e\u0433\u043e \u0434\u043e\u043f\u0443\u0441\u0442\u0438\u043c\u044b\u0445 \u043a\u043e\u043c\u0431\u0438\u043d\u0430\u0446\u0438\u0439. \u041d\u0430\u043f\u0440\u0438\u043c\u0435\u0440, \u0434\u043b\u044f d-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439 (<em>l<\/em>=2) \u043f\u043e \u0444\u043e\u0440\u043c\u0443\u043b\u0435 \u0443 \u043d\u0430\u0441 \u0434\u043e\u043b\u0436\u043d\u043e \u0431\u044b\u0442\u044c 5 \u0444\u0443\u043d\u043a\u0446\u0438\u0439, \u043d\u043e \u0432\u043e\u0437\u043c\u043e\u0436\u043d\u044b\u0445 \u043a\u043e\u043c\u0431\u0438\u043d\u0430\u0446\u0438\u0439 \u0443 \u043d\u0430\u0441 \u0431\u0443\u0434\u0435\u0442 \u0443\u0436\u0435 6 \u0448\u0442\u0443\u043a. \u0427\u0442\u043e \u0434\u0435\u043b\u0430\u0442\u044c \u0441 \u043b\u0438\u0448\u043d\u0438\u043c\u0438? <\/p>\n<p>\u0412\u0435\u0441\u044c \u0441\u044b\u0440-\u0431\u043e\u0440 \u0431\u0443\u0434\u0435\u0442 \u0437\u0430\u043a\u043b\u044e\u0447\u0430\u0442\u044c\u0441\u044f \u0432 \u0442\u043e\u043c, \u0447\u0442\u043e \u043b\u0438\u0448\u043d\u0438\u0435 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 &#8212; \u044d\u0442\u043e \u0431\u0443\u0434\u0443\u0442 \u043a\u0430\u043a\u0438\u0435-\u0442\u043e \u043a\u043e\u043c\u0431\u0438\u043d\u0430\u0446\u0438\u0438 \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u044b\u0445 <em>N<\/em>=2<em>l<\/em>+1-\u0438\u0437\u0431\u0440\u0430\u043d\u043d\u044b\u0445, \u043d\u0430\u0441\u0442\u043e\u044f\u0449\u0438\u0445 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439, \u043a\u043e\u0442\u043e\u0440\u044b\u0435 \u0442\u043e \u043d\u0430\u0441 \u0438 \u0438\u043d\u0442\u0435\u0440\u0435\u0441\u0443\u044e\u0442. \u041f\u043e\u044d\u0442\u043e\u043c\u0443, \u0447\u0442\u043e\u0431\u044b \u0438\u0445 \u043f\u043e\u0441\u0442\u0440\u043e\u0438\u0442\u044c, \u043d\u0443\u0436\u043d\u043e \u0431\u0443\u0434\u0435\u0442 \u043d\u0435\u0441\u043a\u043e\u043b\u044c\u043a\u043e \u043e\u0431\u043e\u0431\u0449\u0438\u0442\u044c \u0443\u0433\u043b\u043e\u0432\u0443\u044e \u0447\u0430\u0441\u0442\u044c \u0432\u043e\u043b\u043d\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438, \u0432\u0437\u044f\u0432 \u0432\u0437\u0432\u0435\u0448\u0435\u043d\u043d\u0443\u044e \u0441\u0443\u043c\u043c\u0443 \u0440\u0430\u0437\u043b\u0438\u0447\u043d\u044b\u0445 \u0444\u043e\u0440\u043c \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"formula\" source=\"A_\\text{\u043a\u043e\u043c\u0431} = c_1 x^{n_{x1}} y^{n_{y1}} z^{n_{z1}} + c_2 x^{n_{x2}} y^{n_{y2}} z^{n_{z}} + c_3 x^{n_{x3}} y^{n_{y3}} z^{n_{z3}} + \\ldots\" alt=\"A_\\text{\u043a\u043e\u043c\u0431} = c_1 x^{n_{x1}} y^{n_{y1}} z^{n_{z1}} + c_2 x^{n_{x2}} y^{n_{y2}} z^{n_{z}} + c_3 x^{n_{x3}} y^{n_{y3}} z^{n_{z3}} + \\ldots\" src=\"https:\/\/habrastorage.org\/getpro\/habr\/upload_files\/71f\/010\/7fa\/71f0107fa87699011726653af7dcf239.svg\" width=\"495\" height=\"28\"\/><\/p>\n<p>\u0433\u0434\u0435 <em>c<\/em> &#8212; \u044d\u0442\u043e \u043d\u0435\u043a\u0438\u0435 \u043f\u043e\u0441\u0442\u043e\u044f\u043d\u043d\u044b\u0435 \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442\u044b. \u0418 \u0432\u043e\u0442 \u0434\u043b\u044f \u043e\u043a\u043e\u043d\u0447\u0430\u0442\u0435\u043b\u044c\u043d\u043e\u0433\u043e \u0438\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u044f \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u043a \u043c\u044b \u0438 \u0431\u0443\u0434\u0435\u043c \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u0442\u044c \u0432\u044b\u0440\u0430\u0436\u0435\u043d\u0438\u0435 \u0434\u043b\u044f \u0432\u043e\u043b\u043d\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"formula\" source=\"\\psi(x,y,z) = A_\\text{\u043a\u043e\u043c\u0431} \\cdot \\overbrace{ \\exp\\left( -\\frac{\\sqrt{x^2 + y^2 +z^2}}{R_0} \\right) }^{R}\" alt=\"\\psi(x,y,z) = A_\\text{\u043a\u043e\u043c\u0431} \\cdot \\overbrace{ \\exp\\left( -\\frac{\\sqrt{x^2 + y^2 +z^2}}{R_0} \\right) }^{R}\" src=\"https:\/\/habrastorage.org\/getpro\/habr\/upload_files\/e77\/909\/924\/e77909924de5efe6ccc511e20678bac5.svg\" width=\"365\" height=\"89\"\/><\/p>\n<p>\u0412 \u0441\u043a\u0440\u0438\u043f\u0442\u0435 \u0432\u0441\u0435 \u044d\u0442\u0438 \u0441\u0442\u0435\u043f\u0435\u043d\u0438 <em>n<\/em> \u0438 \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442\u044b c \u0437\u0430\u0434\u0430\u044e\u0442\u0441\u044f \u0440\u0430\u0437\u043b\u0438\u0447\u043d\u044b\u043c\u0438 \u0430\u0440\u0433\u0443\u043c\u0435\u043d\u0442\u0430\u043c\u0438 (\u043e \u0447\u0451\u043c \u043c\u043e\u0436\u043d\u043e \u0443\u0437\u043d\u0430\u0442\u044c \u0437\u0430\u043f\u0443\u0441\u0442\u0438\u0432 \u0441\u043a\u0440\u0438\u043f\u0442 \u0441 \u043e\u043f\u0446\u0438\u044f\u043c\u0438 -h \u0438\u043b\u0438 &#8212;help):<\/p>\n<pre><code class=\"bash\">  --c1 C1               \u041a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043f\u0435\u0440\u0435\u0434 \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0435\u0439   --c2 C2               \u041a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043f\u0435\u0440\u0435\u0434 \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0435\u0439   --c3 C3               \u041a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043f\u0435\u0440\u0435\u0434 \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0435\u0439   --nx1 NX1             \u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 x**n \u0434\u043b\u044f \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438   --nx2 NX2             \u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 x**n \u0434\u043b\u044f \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438   --nx3 NX3             \u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 x**n \u0434\u043b\u044f \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438   --ny1 NY1             \u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 y**n \u0434\u043b\u044f \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438   --ny2 NY2             \u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 y**n \u0434\u043b\u044f \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438   --ny3 NY3             \u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 y**n \u0434\u043b\u044f \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438   --nz1 NZ1             \u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 z**n \u0434\u043b\u044f \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438   --nz2 NZ2             \u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 z**n \u0434\u043b\u044f \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438   --nz3 NZ3             \u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 z**n \u0434\u043b\u044f \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438<\/code><\/pre>\n<p>\u0421 \u044d\u0442\u0438\u043c\u0438 \u0437\u043d\u0430\u043d\u0438\u044f\u043c\u0438 \u043f\u043e\u043f\u0440\u043e\u0431\u0443\u0435\u043c \u043d\u0430\u0440\u0438\u0441\u043e\u0432\u0430\u0442\u044c \u043f\u0440\u043e\u0441\u0442\u0435\u0439\u0448\u0438\u0435 \u0442\u0438\u043f\u044b \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u043a.<\/p>\n<h2>S-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438<\/h2>\n<p>\u041a\u0430\u043a \u043c\u043d\u043e\u0433\u0438\u0435 \u043c\u043e\u0433\u0443\u0442 \u043f\u043e\u043c\u043d\u0438\u0442\u044c (\u0438\u0437 \u043a\u0443\u0440\u0441\u0430 \u0448\u043a\u043e\u043b\u044c\u043d\u043e\u0439 \u0445\u0438\u043c\u0438\u0438) &#8212; \u044d\u0442\u043e \u0442\u0430\u043a\u0438\u0435 &#171;\u0448\u0430\u0440\u0438\u043a\u0438&#187;, \u043a\u043e\u0433\u0434\u0430 \u044d\u043b\u0435\u043a\u0442\u0440\u043e\u043d \u043b\u0435\u0442\u0430\u0435\u0442 \u0432\u043e\u043a\u0440\u0443\u0433 \u0430\u0442\u043e\u043c\u0430 \u0432 \u0432\u0438\u0434\u0435 \u0448\u0430\u0440\u043e\u043e\u0431\u0440\u0430\u0437\u043d\u043e\u0433\u043e \u043e\u0431\u043b\u0430\u043a\u0430. \u041d\u0430\u0448\u0430 \u0432\u043e\u043b\u043d\u043e\u0432\u0430\u044f \u0444\u0443\u043d\u043a\u0446\u0438\u044f \u0438\u043c\u0435\u0435\u0442 \u0432\u0438\u0434 <em>c<\/em><sub>1<\/sub>=1, <em>n<sub>x<\/sub><\/em><sub>1<\/sub>=<em>n<sub>y<\/sub><\/em><sub>1<\/sub>=<em>n<sub>z<\/sub><\/em><sub>1<\/sub>=0. \u0427\u0442\u043e\u0431\u044b \u043f\u043e\u0441\u0442\u0440\u043e\u0438\u0442\u044c \u044d\u0442\u0443 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0443 \u0441\u043a\u0440\u0438\u043f\u0442\u043e\u043c, \u0434\u0430\u043d\u043d\u044b\u043c \u0432 \u0441\u0430\u043c\u043e\u043c \u043d\u0430\u0447\u0430\u043b\u0435 \u043f\u043e\u0441\u0442\u0430, \u043d\u0430\u0434\u043e \u0431\u0443\u0434\u0435\u0442 \u0432\u044b\u043f\u043e\u043b\u043d\u0438\u0442\u044c \u043a\u043e\u043c\u0430\u043d\u0434\u0443,<\/p>\n<pre><code class=\"bash\">python3 gen-plot_orb.py --plot<\/code><\/pre>\n<p>\u0433\u0434\u0435 <code>gen-plot_orb.py<\/code> &#8212; \u044d\u0442\u043e \u043d\u0430\u0437\u0432\u0430\u043d\u0438\u0435 \u0441\u043a\u0440\u0438\u043f\u0442\u0430, \u0430 &#171;<code>--plot<\/code>&#187;  &#8212; \u044d\u0442\u043e \u043e\u043f\u0446\u0438\u044f, \u0447\u0442\u043e\u0431\u044b \u0432\u043a\u043b\u044e\u0447\u0438\u0442\u044c \u043e\u0442\u0440\u0438\u0441\u043e\u0432\u043a\u0443 \u043f\u043e\u043b\u0443\u0447\u0435\u043d\u043d\u043e\u0439 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438. \u0412 \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442\u0435 \u043c\u044b \u043f\u043e\u043b\u0443\u0447\u0438\u043c \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0443, \u0433\u0434\u0435 \u0447\u0438\u0441\u043b\u043e \u0442\u043e\u0447\u0435\u043a \u0438 \u0438\u0445 \u0446\u0432\u0435\u0442 \u0431\u0443\u0434\u0435\u0442 \u043e\u0442\u043e\u0431\u0440\u0430\u0436\u0430\u0442\u044c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0432\u043e\u043b\u043d\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438:<\/p>\n<figure class=\"bordered full-width\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/habrastorage.org\/r\/w1560\/getpro\/habr\/upload_files\/243\/d9c\/0d5\/243d9c0d5dd5d5ac1e6d56654cab6793.png\" alt=\"\u0418\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u0435 s-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438, \u043f\u043e\u043b\u0443\u0447\u0435\u043d\u043d\u043e\u0435 \u0441\u043a\u0440\u0438\u043f\u0442\u043e\u043c \u0432\u044b\u0448\u0435\" title=\"\u0418\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u0435 s-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438, \u043f\u043e\u043b\u0443\u0447\u0435\u043d\u043d\u043e\u0435 \u0441\u043a\u0440\u0438\u043f\u0442\u043e\u043c \u0432\u044b\u0448\u0435\" width=\"640\" height=\"480\" data-src=\"https:\/\/habrastorage.org\/getpro\/habr\/upload_files\/243\/d9c\/0d5\/243d9c0d5dd5d5ac1e6d56654cab6793.png\"\/><figcaption>\u0418\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u0435 s-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438, \u043f\u043e\u043b\u0443\u0447\u0435\u043d\u043d\u043e\u0435 \u0441\u043a\u0440\u0438\u043f\u0442\u043e\u043c \u0432\u044b\u0448\u0435<\/figcaption><\/figure>\n<p>\u042f \u043d\u0438\u043a\u0430\u043a\u0443\u0449\u0438\u0439 \u0441\u043f\u0435\u0446\u0438\u0430\u043b\u0438\u0441\u0442 \u043f\u043e Matplotlib, \u043f\u043e\u044d\u0442\u043e\u043c\u0443 \u043e\u0442\u0441\u044e\u0434\u0430 \u0438 \u0434\u0430\u043b\u0435\u0435 \u044f \u0431\u0443\u0434\u0443 \u043f\u0440\u0438\u0432\u043e\u0434\u0438\u0442\u044c \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0438, \u043f\u043e\u0441\u0442\u0440\u043e\u0435\u043d\u043d\u044b\u0435 \u043f\u0440\u0438 \u043f\u043e\u043c\u043e\u0449\u0438 Gnuplot \u043a\u043e\u0441\u044b\u043c\u0438 \u0441\u043a\u0440\u0438\u043f\u0442\u0430\u043c\u0438, \u0442\u0438\u043f\u0430 \u0442\u043e\u0433\u043e, \u0447\u0442\u043e \u0443\u0431\u0440\u0430\u043d \u043f\u043e\u0434 \u0441\u043f\u043e\u0439\u043b\u0435\u0440.  <\/p>\n<details class=\"spoiler\">\n<summary>\u0421\u043a\u0440\u0438\u043f\u0442 \u0434\u043b\u044f Gnuplot<\/summary>\n<div class=\"spoiler__content\">\n<pre><code>set view 82.,32.  set xlabel \"x\" set xtics 0.,1e5,1e5   set ylabel \"y\" set ytics 0.,1e5,1e5  set zlabel \"z\" set ztics 0.,1e5,1e5  set cblabel \"\u03c8\" set cbtics 0.,1e5,1e5  set xyplane at 0 set view equal xyz  unset key  set terminal pngcairo size 1000,1000 font \"FreeMono,24\"   ##### s ################################# set title \"s\" FNAME=\"orbital_+1.00s.dat\" set output \"s.png\" maxWFN=0.9089320000 ###########################################  set cbrange [-maxWFN:maxWFN] splot FNAME u 1:2:3:4 w p pt 7 ps 0.1 lc palette<\/code><\/pre>\n<\/div>\n<\/details>\n<p>\u0421\u043e\u043e\u0442\u0432\u0435\u0442\u0441\u0442\u0432\u0435\u043d\u043d\u043e, \u043f\u0435\u0440\u0435\u0441\u0442\u0440\u043e\u0435\u043d\u043d\u0430\u044f Gnuplot-\u043e\u043c s-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c \u0431\u0443\u0434\u0435\u0442 \u0432\u044b\u0433\u043b\u044f\u0434\u0435\u0442\u044c \u0442\u0430\u043a:<\/p>\n<figure class=\"bordered full-width\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/habrastorage.org\/r\/w1560\/getpro\/habr\/upload_files\/dda\/480\/2b8\/dda4802b8ee6cdd094e364fecc31a3dd.png\" alt=\"s-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c, \u043f\u043e\u0441\u0442\u0440\u043e\u0435\u043d\u043d\u0430\u044f Gnuplot-\u043e\u043c.\" title=\"s-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c, \u043f\u043e\u0441\u0442\u0440\u043e\u0435\u043d\u043d\u0430\u044f Gnuplot-\u043e\u043c.\" width=\"1000\" height=\"1000\" data-src=\"https:\/\/habrastorage.org\/getpro\/habr\/upload_files\/dda\/480\/2b8\/dda4802b8ee6cdd094e364fecc31a3dd.png\"\/><figcaption>s-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c, \u043f\u043e\u0441\u0442\u0440\u043e\u0435\u043d\u043d\u0430\u044f Gnuplot-\u043e\u043c.<\/figcaption><\/figure>\n<p>\u0412\u0441\u0451, \u043a\u0430\u043a \u043c\u044b \u0438 \u043f\u0440\u0438\u0432\u044b\u043a\u043b\u0438, \u0448\u0430\u0440\u0438\u043a.<\/p>\n<h2>P-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438<\/h2>\n<p>\u0417\u0434\u0435\u0441\u044c <em>l<\/em>=1, \u0438, \u0441\u043e\u043e\u0442\u0432\u0435\u0442\u0441\u0442\u0432\u0435\u043d\u043d\u043e, \u0442\u0440\u0438 p-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438, \u043c\u044b \u043c\u043e\u0436\u0435\u043c \u043f\u043e\u0441\u0442\u0440\u043e\u0438\u0442\u044c \u043a\u043e\u043c\u0430\u043d\u0434\u0430\u043c\u0438:<\/p>\n<\/p>\n<pre><code class=\"bash\">python3 gen-plot_orb.py --nx1 1 --plot python3 gen-plot_orb.py --ny1 1 --plot python3 gen-plot_orb.py --nz1 1 --plot<\/code><\/pre>\n<p>\u0412 \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442\u0435 \u0447\u0435\u0433\u043e \u043c\u044b \u043f\u043e\u043b\u0443\u0447\u0438\u043c \u0441\u043b\u0435\u0434\u0443\u044e\u0449\u0438\u0435 \u0438\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u044f \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439:<\/p>\n<figure class=\"bordered full-width\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/habrastorage.org\/r\/w1560\/getpro\/habr\/upload_files\/b71\/8ca\/0cb\/b718ca0cbedf21b14d52b750d53dbb51.png\" alt=\"\u0422\u0440\u0438 p-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438.\" title=\"\u0422\u0440\u0438 p-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438.\" width=\"1325\" height=\"442\" data-src=\"https:\/\/habrastorage.org\/getpro\/habr\/upload_files\/b71\/8ca\/0cb\/b718ca0cbedf21b14d52b750d53dbb51.png\"\/><figcaption>\u0422\u0440\u0438 p-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438.<\/figcaption><\/figure>\n<p>\u0412\u043d\u043e\u0432\u044c, \u0432\u0441\u0451 \u0432\u044b\u0433\u043b\u044f\u0434\u0438\u0442 \u043a\u0430\u043a \u0438 \u043e\u0436\u0438\u0434\u0430\u043b\u043e\u0441\u044c \u043f\u043e\u0441\u043b\u0435 \u0438\u0437\u0443\u0447\u0435\u043d\u0438\u044f \u0448\u043a\u043e\u043b\u044c\u043d\u043e\u0439 \u043f\u0440\u043e\u0433\u0440\u0430\u043c\u043c\u044b: \u0442\u0440\u0438 &#171;\u0433\u0430\u043d\u0442\u0435\u043b\u044c\u043a\u0438&#187;, \u043d\u0430\u043f\u0440\u0430\u0432\u043b\u0435\u043d\u043d\u044b\u0435 \u0432 \u0442\u0440\u0438 \u0440\u0430\u0437\u043d\u044b\u0435 \u0441\u0442\u043e\u0440\u043e\u043d\u044b, \u0438 \u043a\u0430\u0436\u0434\u0430\u044f \u043f\u043e\u043b\u043e\u0432\u0438\u043d\u043a\u0430 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 \u0438\u043c\u0435\u0435\u0442 \u0437\u043d\u0430\u043a, \u043f\u0440\u043e\u0442\u0438\u0432\u043e\u043f\u043e\u043b\u043e\u0436\u043d\u044b\u0439 \u0434\u0440\u0443\u0433\u043e\u0439 \u043f\u043e\u043b\u043e\u0432\u0438\u043d\u043a\u0435 \u0442\u043e\u0439 \u0436\u0435 \u0433\u0430\u043d\u0442\u0435\u043b\u0438. <\/p>\n<h2>D-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438<\/h2>\n<p><em>l<\/em>=3, \u0438 \u0442\u0443\u0442 \u0434\u043e\u043b\u0436\u043d\u043e \u0431\u044b\u0442\u044c \u043f\u044f\u0442\u044c \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439. \u0422\u0440\u0438 \u0438\u0437 \u043d\u0438\u0445 \u0441\u0442\u0440\u043e\u044f\u0442\u0441\u044f \u043f\u0440\u043e\u0441\u0442\u043e (\u0442.\u043d. d<sub>xy<\/sub>, d<sub>xz<\/sub>, d<sub>yz<\/sub>):<\/p>\n<pre><code class=\"bash\">python3 gen-plot_orb.py --nx1 1 --ny1 1 --plot python3 gen-plot_orb.py --nx1 1 --nz1 1 --plot python3 gen-plot_orb.py --ny1 1 --nz1 1 --plot<\/code><\/pre>\n<p>\u0412\u044b\u0433\u043b\u044f\u0434\u044f\u0442 \u043e\u043d\u0438 \u043a\u0430\u043a \u0434\u0432\u0435 \u0441\u043f\u0430\u0440\u0435\u043d\u043d\u044b\u0435 \u0433\u0430\u043d\u0442\u0435\u043b\u0438 p-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439 \u043d\u043e \u0441 \u0438\u0437\u043c\u0435\u043d\u0451\u043d\u043d\u044b\u043c\u0438 \u0437\u043d\u0430\u043a\u0430\u043c\u0438 \u043f\u043e\u043b\u043e\u0432\u0438\u043d\u043e\u043a, \u043d\u0430\u043f\u0440\u0438\u043c\u0435\u0440:<\/p>\n<figure class=\"bordered full-width\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/habrastorage.org\/r\/w1560\/getpro\/habr\/upload_files\/519\/c28\/1cb\/519c281cb1a3aac216bf58aa03c8f286.png\" alt=\"dxz-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0430\" title=\"dxz-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0430\" width=\"1000\" height=\"1000\" data-src=\"https:\/\/habrastorage.org\/getpro\/habr\/upload_files\/519\/c28\/1cb\/519c281cb1a3aac216bf58aa03c8f286.png\"\/><figcaption>dxz-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0430<\/figcaption><\/figure>\n<p>\u0415\u0449\u0451 \u043e\u0434\u043d\u0430 \u0438\u0437 \u0442\u0430\u043a\u0438\u0445 \u0441\u043f\u0430\u0440\u0435\u043d\u043d\u044b\u0445 \u0433\u0430\u043d\u0442\u0435\u043b\u0435\u0439 \u043f\u043e\u043b\u0443\u0447\u0430\u0435\u0442\u0441\u044f \u043a\u0430\u043a \u0440\u0430\u0437\u043d\u043e\u0441\u0442\u044c <em>x<\/em><sup>2<\/sup> \u0438 <em>y<\/em><sup>2<\/sup> \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439 \u043a\u043e\u043c\u0430\u043d\u0434\u043e\u0439 <\/p>\n<pre><code class=\"bash\">python3 gen-plot_orb.py --nx1 2 --ny2 2 --c1 1 --c2 -1 --plot<\/code><\/pre>\n<p>\u0410 \u0432\u043e\u0442 \u043f\u043e\u0441\u043b\u0435\u0434\u043d\u044f\u044f \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c, \u043d\u0430\u0437\u044b\u0432\u0430\u0435\u043c\u0430\u044f d(<em>z<\/em><sup>2<\/sup>) &#8212; \u0441\u0430\u043c\u0430\u044f \u0441\u043b\u043e\u0436\u043d\u0430\u044f, \u043e\u043d\u0430 \u043d\u0430 \u0441\u0430\u043c\u043e\u043c \u0434\u0435\u043b\u0435 \u0441\u043e\u0441\u0442\u043e\u0438\u0442 \u0438\u0437 \u0442\u0440\u0451\u0445 \u043d\u0430\u0448\u0438\u0445 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439-\u043a\u0438\u0440\u043f\u0438\u0447\u0438\u043a\u043e\u0432. \u041a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442\u044b \u043c\u0435\u0436\u0434\u0443 \u043f\u043e\u043b\u043e\u0436\u0438\u0442\u0435\u043b\u044c\u043d\u044b\u043c\u0438 \u0438 \u043e\u0442\u0440\u0438\u0446\u0430\u0442\u0435\u043b\u044c\u043d\u044b\u043c\u0438 \u043a\u0443\u0441\u043a\u0430\u043c\u0438 \u043d\u0430 \u0441\u0430\u043c\u043e\u043c \u0434\u0435\u043b\u0435 \u043c\u043e\u0436\u043d\u043e \u043d\u0430\u0439\u0442\u0438, \u043d\u043e \u044d\u0442\u043e \u043d\u0430\u0434\u043e \u0438\u043b\u0438 \u043a\u043e\u043f\u0430\u0442\u044c \u043b\u0438\u0442\u0435\u0440\u0430\u0442\u0443\u0440\u0443, \u0438\u043b\u0438 \u043f\u043e\u0434\u043d\u0438\u043c\u0430\u0442\u044c \u0437\u0443\u0431\u043e\u0434\u0440\u043e\u0431\u0438\u0442\u0435\u043b\u044c\u043d\u044b\u0439 \u043c\u0430\u0442\u0430\u043d, \u0442\u0430\u043a \u0447\u0442\u043e \u044f \u0438\u0445 \u043f\u0440\u043e\u0441\u0442\u043e \u043f\u043e\u0434\u043e\u0431\u0440\u0430\u043b, \u0447\u0442\u043e\u0431\u044b \u0431\u044b\u043b\u043e \u043f\u043e\u0445\u043e\u0436\u0435 \u043d\u0430 \u0442\u043e, \u043a\u0430\u043a \u044d\u0442\u0430 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c \u0434\u043e\u043b\u0436\u043d\u0430 \u0432\u044b\u0433\u043b\u044f\u0434\u0435\u0442\u044c:<\/p>\n<pre><code class=\"bash\">python3 gen-plot_orb.py --nz1 2 --nx2 2 --ny3 2 --c1 3 --c2 -1 --c3 -1 --plot<\/code><\/pre>\n<p>\u0447\u0442\u043e \u0434\u0430\u0451\u0442 \u0432 \u0440\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442\u0435 \u0438\u0441\u043a\u043e\u043c\u0443\u044e \u0444\u043e\u0440\u043c\u0443 \u043e\u0434\u043d\u043e\u0437\u043d\u0430\u0447\u043d\u043e\u0439 \u0433\u0430\u043d\u0442\u0435\u043b\u044c\u043a\u0438 \u0441 \u043d\u0430\u0434\u0435\u0442\u044b\u043c \u043d\u0430 \u043d\u0435\u0451 <s>\u0431\u0443\u0431\u043b\u0438\u043a\u043e\u043c<\/s> \u0442\u043e\u0440\u043e\u043c:<\/p>\n<figure class=\"bordered full-width\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/habrastorage.org\/r\/w1560\/getpro\/habr\/upload_files\/649\/c16\/f65\/649c16f658971fdd6868080b151fbc4f.png\" alt=\"dz2-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0430\" title=\"dz2-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0430\" width=\"1000\" height=\"1000\" data-src=\"https:\/\/habrastorage.org\/getpro\/habr\/upload_files\/649\/c16\/f65\/649c16f658971fdd6868080b151fbc4f.png\"\/><figcaption>dz2-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0430<\/figcaption><\/figure>\n<p>\u0410 \u0431\u043e\u043b\u0435\u0435 \u0441\u043b\u043e\u0436\u043d\u044b\u0435 \u0432\u0438\u0434\u044b \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439 \u043e\u0441\u0442\u0430\u0432\u043b\u044f\u044e \u043d\u0430 \u0441\u0430\u043c\u043e\u0441\u0442\u043e\u044f\u0442\u0435\u043b\u044c\u043d\u044b\u0439 \u0440\u0430\u0437\u0431\u043e\u0440, \u0438\u0431\u043e f-\u0444\u0443\u043d\u043a\u0446\u0438\u0439 \u0431\u0443\u0434\u0435\u0442 \u0443\u0436\u0435 7 \u0448\u0442\u0443\u043a, \u043a\u043e\u043d\u0441\u0442\u0440\u0443\u0438\u0440\u0443\u0435\u043c\u044b\u0445 \u0438\u0437 10 \u043f\u0440\u043e\u0441\u0442\u0435\u0439\u0448\u0438\u0445 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u043a. <\/p>\n<h2>\u041f\u0430\u0440\u043e\u0447\u043a\u0430 \u0434\u0438\u0441\u043a\u043b\u0435\u0439\u043c\u0435\u0440\u043e\u0432 \u0438 \u0430\u043b\u044c\u0442\u0435\u0440\u043d\u0430\u0442\u0438\u0432<\/h2>\n<p>\u041f\u0440\u043e\u0441\u0442\u043e \u043f\u0430\u0440\u043e\u0447\u043a\u0430 \u043a\u043e\u043c\u043c\u0435\u043d\u0442\u0430\u0440\u0438\u0435\u0432, \u043a\u043e\u0442\u043e\u0440\u044b\u0435 \u044f \u043d\u0435 \u0437\u043d\u0430\u044e \u043a\u0443\u0434\u0430 \u0437\u0430\u043f\u0438\u0445\u0430\u0442\u044c \u0435\u0449\u0451.<\/p>\n<ul>\n<li>\n<p>\u041d\u0430 \u0441\u0430\u043c\u043e\u043c \u0434\u0435\u043b\u0435, \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u0430\u044f \u0447\u0430\u0441\u0442\u044c \u0430\u0442\u043e\u043c\u043d\u044b\u0445 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439 \u0438\u043c\u0435\u0435\u0442 \u0441\u0443\u0449\u0435\u0441\u0442\u0432\u0435\u043d\u043d\u043e \u0431\u043e\u043b\u0435\u0435 \u0441\u043b\u043e\u0436\u043d\u044b\u0439 \u0432\u0438\u0434, \u043d\u043e \u0434\u043b\u044f \u0446\u0435\u043b\u0435\u0439 \u0432\u0438\u0437\u0443\u0430\u043b\u0438\u0437\u0430\u0446\u0438\u0438 \u044d\u0442\u043e \u043d\u0435 \u043e\u0447\u0435\u043d\u044c \u0432\u0430\u0436\u043d\u043e, \u0438 \u0432 \u0442\u0435\u0445 \u0436\u0435 \u043a\u0432\u0430\u043d\u0442\u043e\u0432\u043e-\u0445\u0438\u043c\u0438\u0447\u0435\u0441\u043a\u0438\u0445 \u0440\u0430\u0441\u0447\u0451\u0442\u0430\u0445 \u043f\u043e \u043c\u0435\u0442\u043e\u0434\u0443 \u041c\u041e \u041b\u041a\u0410\u041e (\u043c\u043e\u043b\u0435\u043a\u0443\u043b\u044f\u0440\u043d\u044b\u0435 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 \u043a\u0430\u043a \u043b\u0438\u043d\u0435\u0439\u043d\u044b\u0435 \u043a\u043e\u043c\u0431\u0438\u043d\u0430\u0446\u0438\u0438 \u0430\u0442\u043e\u043c\u043d\u044b\u0445 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439) \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u044e\u0442\u0441\u044f \u0438\u043c\u0435\u043d\u043d\u043e \u0432\u044b\u0440\u0430\u0436\u0435\u043d\u0438\u044f, \u0430\u043d\u0430\u043b\u043e\u0433\u0438\u0447\u043d\u044b\u0435 \u0442\u043e\u043c\u0443, \u0447\u0442\u043e \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u043e\u0432\u0430\u043b\u0438 \u043c\u044b, \u0442\u0430\u043a \u0447\u0442\u043e \u0441\u0438\u043b\u044c\u043d\u043e \u043c\u044b \u043e\u0442 \u0442\u0430\u043a\u043e\u0433\u043e \u0443\u043f\u0440\u043e\u0449\u0435\u043d\u0438\u044f \u043d\u0435 \u0441\u043e\u0432\u0440\u0430\u043b\u0438.<\/p>\n<\/li>\n<li>\n<p>\u042f \u043f\u043e\u043d\u0438\u043c\u0430\u044e, \u0447\u0442\u043e \u043c\u043e\u0433 \u0443\u0436\u0435 \u0434\u043e\u0441\u0442\u0430\u0442\u044c \u0432\u0441\u044f\u043a\u0438\u043c\u0438 \u0448\u0442\u0443\u043a\u0430\u043c\u0438 \u0442\u0438\u043f\u0430 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u043a \u0438 \u041c\u043e\u043d\u0442\u0435-\u041a\u0430\u0440\u043b\u043e. \u0415\u0441\u043b\u0438 \u044d\u0442\u043e \u0442\u0430\u043a, \u0434\u0430\u0439\u0442\u0435 \u0437\u043d\u0430\u0442\u044c \u0432 \u043a\u043e\u043c\u043c\u0435\u043d\u0442\u0430\u0445 (\u0442\u0430\u043c \u043c\u043e\u0436\u043d\u043e \u043a\u0430\u043a \u0432\u044b\u0440\u0430\u0437\u0438\u0442\u044c \u0441\u0432\u043e\u0451 &#171;\u0444\u0438&#187;, \u0442\u0430\u043a \u0438 \u0434\u043e\u0431\u0430\u0432\u0438\u0442\u044c \u0431\u043e\u043b\u0435\u0435 \u043a\u043e\u043d\u0441\u0442\u0440\u0443\u043a\u0442\u0438\u0432\u043d\u044b\u0445 \u043a\u043e\u043c\u043c\u0435\u043d\u0442\u043e\u0432, \u043d\u0430\u043f\u0440\u0438\u043c\u0435\u0440, \u043e \u0447\u0451\u043c \u0431\u044b \u0445\u043e\u0442\u0435\u043b\u043e\u0441\u044c \u0432\u0438\u0434\u0435\u0442\u044c \u0441\u0442\u0430\u0442\u044c\u0438).<\/p>\n<\/li>\n<li>\n<p>\u041c\u043d\u043e\u0433\u043e \u043a\u043e\u043c\u0443 \u043d\u0435 \u043d\u0440\u0430\u0432\u0438\u0442\u0441\u044f \u044d\u0442\u043e\u0442 \u0441\u0442\u0438\u043b\u044c \u0438\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u044f \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439 (\u0434\u043e\u0441\u0442\u0430\u0442\u043e\u0447\u043d\u043e \u043f\u043e\u0441\u043c\u043e\u0442\u0440\u0435\u0442\u044c \u043d\u0430 \u043e\u0442\u043a\u043b\u0438\u043a\u0438 \u043a \u043f\u043e\u0441\u0442\u0443 <a href=\"https:\/\/habr.com\/ru\/post\/548704\/\" rel=\"noopener noreferrer nofollow\">https:\/\/habr.com\/ru\/post\/548704\/<\/a>). \u041c\u043d\u0435 \u044d\u0442\u043e\u0442 \u0441\u0442\u0438\u043b\u044c, \u043e\u0447\u0435\u0432\u0438\u0434\u043d\u043e, \u043d\u0440\u0430\u0432\u0438\u0442\u0441\u044f, \u043f\u043e\u0441\u043a\u043e\u043b\u044c\u043a\u0443 \u043e\u043d, \u043a\u0430\u043a \u043c\u043d\u0435 \u043a\u0430\u0436\u0435\u0442\u0441\u044f, \u043e\u0442\u0440\u0430\u0436\u0430\u0435\u0442 \u0441\u0442\u0430\u0442\u0438\u0441\u0442\u0438\u0447\u0435\u0441\u043a\u0443\u044e \u043f\u0440\u0438\u0440\u043e\u0434\u0443 \u0432\u043e\u043b\u043d\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438, \u0430 \u0435\u0449\u0451 \u0432 \u043c\u043e\u0451\u043c \u0431\u043e\u043b\u044c\u043d\u043e\u043c \u043c\u043e\u0437\u0433\u0443, \u0432\u0438\u0434\u0438\u043c\u043e, \u043e\u043d \u0430\u0441\u0441\u0438\u0446\u0438\u0438\u0440\u0443\u0435\u0442\u0441\u044f \u0441 \u0447\u0435\u043c-\u0442\u043e \u0442\u0438\u043f\u0430 \u0438\u043c\u043f\u0440\u0435\u0441\u0441\u0438\u043e\u043d\u0438\u0437\u043c\u0430. \u042d\u0442\u043e\u0442 \u0441\u043a\u0440\u0438\u043f\u0442 \u0438 \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0438 \u044f \u0434\u0435\u043b\u0430\u043b \u0434\u043b\u044f \u0441\u0432\u043e\u0435\u0439 \u043e\u0431\u0440\u0430\u0437\u043e\u0432\u0430\u0442\u0435\u043b\u044c\u043d\u043e\u0439 \u0434\u0435\u044f\u0442\u0435\u043b\u044c\u043d\u043e\u0441\u0442\u0438, \u0442\u0430\u043a \u0447\u0442\u043e \u043c\u043e\u0433\u0443 \u0442\u043e\u043b\u044c\u043a\u043e \u0441\u043a\u0430\u0437\u0430\u0442\u044c &#171;\u043d\u0435 \u043d\u0440\u0430\u0432\u0438\u0442\u0441\u044f &#8212; \u043d\u0435 \u0435\u0448\u044c\u0442\u0435&#187;. <\/p>\n<\/li>\n<li>\n<p>\u0412 \u043a\u0430\u0447\u0435\u0441\u0442\u0432\u0435 \u043f\u0440\u043e\u0441\u0442\u043e\u0439 (\u0438 \u0431\u0435\u0441\u043f\u043b\u0430\u0442\u043d\u043e\u0439!) \u0430\u043b\u044c\u0442\u0435\u0440\u043d\u0430\u0442\u0438\u0432\u044b, \u043c\u043e\u0433\u0443 \u043f\u043e\u0441\u043e\u0432\u0435\u0442\u043e\u0432\u0430\u0442\u044c \u043c\u043e\u043b\u0435\u043a\u0443\u043b\u044f\u0440\u043d\u044b\u0439 \u0432\u0438\u0437\u0443\u0430\u043b\u0438\u0437\u0430\u0442\u043e\u0440 Jmol (<a href=\"http:\/\/jmol.sourceforge.net\/\" rel=\"noopener noreferrer nofollow\">http:\/\/jmol.sourceforge.net\/<\/a>). \u041f\u043e\u043c\u0438\u043c\u043e \u043a\u0443\u0447\u0438 \u0440\u0435\u0430\u043b\u044c\u043d\u043e \u043f\u043e\u043b\u0435\u0437\u043d\u044b\u0445 \u0432 \u0440\u0430\u0431\u043e\u0442\u0435 \u0444\u0443\u043d\u043a\u0446\u0438\u0439, \u043e\u043d \u0443\u043c\u0435\u0435\u0442 \u0438 \u043f\u0440\u043e\u0441\u0442\u043e \u0441\u0442\u0440\u043e\u0438\u0442\u044c \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 \u0432\u043e\u0434\u043e\u0440\u043e\u0434\u043e\u043f\u043e\u0434\u043e\u0431\u043d\u043e\u0433\u043e \u0430\u0442\u043e\u043c\u0430 \u0432 \u0441\u0442\u0430\u043d\u0434\u0430\u0440\u0442\u043d\u043e\u043c \u0432\u0438\u0434\u0435. \u041d\u0430\u043f\u0440\u0438\u043c\u0435\u0440, \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 \u043c\u043e\u0436\u043d\u043e \u043f\u043e\u0441\u0442\u0440\u043e\u0438\u0442\u044c \u0438\u0437 \u043a\u043e\u043d\u0441\u043e\u043b\u0438 \u043f\u0430\u0440\u043e\u0447\u043a\u043e\u0439 \u043a\u043e\u043c\u0430\u043d\u0434 (\u0441\u043c. \u043d\u0430\u043f\u0440\u0438\u043c\u0435\u0440 <a href=\"https:\/\/chemapps.stolaf.edu\/jmol\/docs\/examples-11\/surfacedemos.htm\" rel=\"noopener noreferrer nofollow\">\u0437\u0434\u0435\u0441\u044c<\/a>).<\/p>\n<\/li>\n<li>\n<p>\u0410 \u0435\u0449\u0451 \u043d\u0430 YouTube \u0435\u0441\u0442\u044c \u043a\u043b\u0430\u0441\u0441\u043d\u044b\u0439 \u0440\u043e\u043b\u0438\u043a \u043e\u0431 \u0438\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u0438 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u0439 \u0432 \u0432\u0438\u0434\u0435 \u0411\u043e\u043c\u043e\u0432\u0441\u043a\u0438\u0445 \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0439: <a href=\"https:\/\/youtu.be\/W2Xb2GFK2yc\" rel=\"noopener noreferrer nofollow\">https:\/\/youtu.be\/W2Xb2GFK2yc<\/a> . <\/p>\n<\/li>\n<\/ul>\n<h2>P.S.<\/h2>\n<p>\u0427\u0442\u043e\u0431\u044b \u043d\u0435 \u0431\u044b\u0442\u044c \u0433\u043e\u043b\u043e\u0441\u043b\u043e\u0432\u043d\u044b\u043c \u043f\u043e \u043f\u043e\u0432\u043e\u0434\u0443 Jmol, \u0432\u043e\u0442 \u0438\u0437\u043e\u0431\u0440\u0430\u0436\u0435\u043d\u0438\u0435 d-\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0438, \u043f\u043e\u043b\u0443\u0447\u0435\u043d\u043d\u043e\u0435 \u043a\u043e\u043c\u0430\u043d\u0434\u0430\u043c\u0438 \u043d\u0438\u0436\u0435:<\/p>\n<figure class=\"bordered\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/habrastorage.org\/r\/w780q1\/getpro\/habr\/upload_files\/5b3\/cc5\/326\/5b3cc5326a118d58350d31dee7227507.jpg\" width=\"508\" height=\"492\" data-src=\"https:\/\/habrastorage.org\/getpro\/habr\/upload_files\/5b3\/cc5\/326\/5b3cc5326a118d58350d31dee7227507.jpg\" data-blurred=\"true\"\/><figcaption><\/figcaption><\/figure>\n<pre><code class=\"java\">isosurface phase atomicOrbital 3 2 0 set axesMolecular;set axesScale 0.5;axes on moveto 1.0 { 462 -868 -180 47.18} 141<\/code><\/pre>\n<p>\u0421\u043f\u0430\u0441\u0438\u0431\u043e \u0437\u0430 \u0432\u043d\u0438\u043c\u0430\u043d\u0438\u0435. <\/p>\n<\/div>\n<\/div>\n<\/div>\n<p><!----><!----><\/div>\n<p><!----><!----><br \/> \u0441\u0441\u044b\u043b\u043a\u0430 \u043d\u0430 \u043e\u0440\u0438\u0433\u0438\u043d\u0430\u043b \u0441\u0442\u0430\u0442\u044c\u0438 <a href=\"https:\/\/habr.com\/ru\/articles\/577450\/\"> https:\/\/habr.com\/ru\/articles\/577450\/<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"<div><!--[--><!--]--><\/div>\n<div id=\"post-content-body\">\n<div>\n<div class=\"article-formatted-body article-formatted-body article-formatted-body_version-2\">\n<div xmlns=\"http:\/\/www.w3.org\/1999\/xhtml\">\n<p>\u041d\u0430\u0447\u0430\u043b\u0441\u044f \u043d\u043e\u0432\u044b\u0439 \u0443\u0447\u0435\u0431\u043d\u044b\u0439 \u0433\u043e\u0434, \u0438 \u043f\u0440\u0435\u043f\u043e\u0434\u0430\u0432\u0430\u0442\u0435\u043b\u044f\u043c, \u0441\u0442\u0443\u0434\u0435\u043d\u0442\u0430\u043c \u0438 \u0448\u043a\u043e\u043b\u044c\u043d\u0438\u043a\u0430\u043c, \u0432\u043e\u0437\u043c\u043e\u0436\u043d\u043e, \u0442\u0440\u0435\u0431\u0443\u0435\u0442\u0441\u044f (\u0438\u043b\u0438 \u043f\u0440\u043e\u0441\u0442\u043e \u0445\u043e\u0447\u0435\u0442\u0441\u044f) \u043f\u043e\u0441\u043c\u043e\u0442\u0440\u0435\u0442\u044c \u043d\u0430 \u0442\u043e, \u043a\u0430\u043a \u0432\u044b\u0433\u043b\u044f\u0434\u044f\u0442 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0438, \u043d\u0430 \u043a\u043e\u0442\u043e\u0440\u044b\u0445 \u0441\u0438\u0434\u044f\u0442 \u044d\u043b\u0435\u043a\u0442\u0440\u043e\u043d\u044b \u0432 \u0430\u0442\u043e\u043c\u0430\u0445: \u0432\u0441\u0435 \u044d\u0442\u0438 \u0437\u0430\u0432\u043e\u0440\u0430\u0436\u0438\u0432\u0430\u044e\u0449\u0438\u0435 \u0431\u0443\u043a\u043e\u0432\u043a\u0438 s, p, d, f, \u0438 \u0442.\u0434. \u0414\u0430, \u043a\u0430\u0440\u0442\u0438\u043d\u043e\u043a \u043f\u043e\u043b\u043d\u043e \u043a\u0430\u043a \u0432 \u0443\u0447\u0435\u0431\u043d\u0438\u043a\u0430\u0445, \u0442\u0430\u043a \u0438 \u0432 \u0418\u043d\u0442\u0435\u0440\u043d\u0435\u0442\u0435, \u043d\u043e \u043f\u043e\u043a\u0440\u0443\u0442\u0438\u0442\u044c \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0438 \u043d\u0430 \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0435 \u043d\u0435 \u043f\u043e\u043b\u0443\u0447\u0438\u0442\u0441\u044f, \u0430 \u043a\u0430\u0440\u0442\u0438\u043d\u043a\u0443 \u0438\u0437 \u0443\u0447\u0435\u0431\u043d\u0438\u043a\u0430\/\u0441 \u043b\u0435\u0432\u043e\u0433\u043e \u0441\u0430\u0439\u0442\u0430 \u0432 \u043f\u0440\u0435\u0437\u0435\u043d\u0442\u0430\u0446\u0438\u044e\/\u0440\u0435\u0444\u0435\u0440\u0430\u0442 \u0431\u0435\u0437 \u043c\u043e\u0440\u043e\u043a\u0438 \u0441 \u043b\u0438\u0446\u0435\u043d\u0437\u0438\u0435\u0439 \u043f\u0438\u0445\u0430\u0442\u044c (\u043f\u043e-\u0445\u043e\u0440\u043e\u0448\u0435\u043c\u0443) \u043d\u0435 \u0441\u0442\u043e\u0438\u0442. \u041f\u043e\u044d\u0442\u043e\u043c\u0443 \u0432 \u044d\u0442\u043e\u043c \u043f\u043e\u0441\u0442\u0435 \u043c\u044b \u0440\u0430\u0437\u0431\u0435\u0440\u0451\u043c \u043e\u0434\u043d\u0443 \u0438\u0437 \u0432\u043e\u0437\u043c\u043e\u0436\u043d\u044b\u0445 \u0440\u0435\u0430\u043b\u0438\u0437\u0430\u0446\u0438\u0439 \u0440\u0438\u0441\u043e\u0432\u0430\u043b\u043a\u0438 \u0434\u043b\u044f \u044d\u0442\u0438\u0445 \u0441\u0430\u043c\u044b\u0445 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0435\u043a. <\/p>\n<hr\/>\n<p>\u0421\u0430\u043c\u0443 \u043a\u043e\u043d\u0446\u0435\u043f\u0446\u0438\u044e \u0442\u043e\u0433\u043e, \u0447\u0442\u043e \u0442\u0430\u043a\u043e\u0435 <a href=\"https:\/\/en.wikipedia.org\/wiki\/Atomic_orbital\" rel=\"noopener noreferrer nofollow\">\u0430\u0442\u043e\u043c\u043d\u044b\u0435 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438<\/a> \u044f \u0443\u0436\u0435 \u0440\u0430\u0437\u0431\u0438\u0440\u0430\u043b (<a href=\"https:\/\/habr.com\/ru\/post\/443232\/\" rel=\"noopener noreferrer nofollow\">https:\/\/habr.com\/ru\/post\/443232\/<\/a>), \u0442\u043e \u043a\u0430\u043a \u0440\u0438\u0441\u043e\u0432\u0430\u0442\u044c \u044d\u0442\u0438 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 \u043f\u0440\u0438 \u043f\u043e\u043c\u043e\u0449\u0438 \u043c\u0435\u0442\u043e\u0434\u0430 \u041c\u043e\u043d\u0442\u0435-\u041a\u0430\u0440\u043b\u043e \u0432 \u0432\u0430\u0440\u0438\u0430\u043d\u0442\u0435 <a href=\"https:\/\/ru.wikipedia.org\/wiki\/%D0%90%D0%BB%D0%B3%D0%BE%D1%80%D0%B8%D1%82%D0%BC_%D0%9C%D0%B5%D1%82%D1%80%D0%BE%D0%BF%D0%BE%D0%BB%D0%B8%D1%81%D0%B0_%E2%80%94_%D0%93%D0%B0%D1%81%D1%82%D0%B8%D0%BD%D0%B3%D1%81%D0%B0\" rel=\"noopener noreferrer nofollow\">\u0430\u043b\u0433\u043e\u0440\u0438\u0442\u043c\u0430 \u041c\u0435\u0442\u0440\u043e\u043f\u043e\u043b\u0438\u0441\u0430 <\/a>\u0442\u043e\u0436\u0435 (<a href=\"https:\/\/habr.com\/ru\/post\/548704\/\" rel=\"noopener noreferrer nofollow\">https:\/\/habr.com\/ru\/post\/548704\/<\/a>),<\/p>\n<details class=\"spoiler\">\n<summary>c\u043e\u0431\u0441\u0442\u0432\u0435\u043d\u043d\u043e, \u0432\u043e\u0442 \u0438 \u0441\u0430\u043c \u043a\u043e\u0434 \u0440\u0438\u0441\u043e\u0432\u0430\u043b\u043a\u0438 \u043d\u0430 Python (\u0441 \u043a\u043e\u043c\u043c\u0435\u043d\u0442\u0430\u0440\u0438\u044f\u043c\u0438),<\/summary>\n<div class=\"spoiler__content\">\n<pre><code class=\"python\">#! \/usr\/bin\/python3  import numpy as np  # NumPy \u043d\u0443\u0436\u0435\u043d \u0434\u043b\u044f \u0432\u0441\u0435\u0433\u043e import argparse     # \u042d\u0442\u043e \u043c\u043e\u0434\u0443\u043b\u044c, \u0434\u043b\u044f \"\u043e\u0431\u0449\u0435\u043d\u0438\u044f\" \u0441 \u044e\u0437\u0435\u0440\u043e\u043c    # \u042d\u0442\u0430 \u0444\u0443\u043d\u043a\u0446\u0438\u044f \u0434\u0430\u0451\u0442 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \"\u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438\" \u0432 \u043d\u0435\u043a\u043e\u0442\u043e\u0440\u043e\u0439 \u0442\u043e\u0447\u043a\u0435 \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0430 # \u041e\u0440\u0431\u0438\u0442\u0430\u043b\u044c \u044f\u0432\u043b\u044f\u0435\u0442\u0441\u044f \u0441\u0443\u043c\u043c\u043e\u0439 \u0442\u0440\u0451\u0445 \u0443\u0433\u043b\u043e\u0432\u044b\u0445 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442 \u0441 \u0442\u0440\u0435\u043c\u044f \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442\u0430\u043c\u0438, \u0434\u043e\u043c\u043d\u043e\u0436\u0435\u043d\u043d\u043e\u0439 \u043d\u0430 \u0435\u0434\u0438\u043d\u0443\u044e  # \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u0443\u044e \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u0443. # r  -- \u044d\u0442\u043e \u0442\u0440\u0451\u0445\u043c\u0435\u0440\u043d\u044b\u0439 \u0432\u0435\u043a\u0442\u043e\u0440, \u0440\u0430\u0441\u043f\u043e\u043b\u043e\u0436\u0435\u043d\u0438\u0435 \u044d\u043b\u0435\u043a\u0442\u0440\u043e\u043d\u0430 \u0432 \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0435 # n1,n2,n3 -- \u044d\u0442\u043e \u0442\u0440\u0451\u0445\u043c\u0435\u0440\u043d\u044b\u0439 \u0432\u0435\u043a\u0442\u043e\u0440 \u0441\u0442\u0435\u043f\u0435\u043d\u0435\u0439 x,y,z \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442 \u0432\u0435\u043a\u0442\u043e\u0440\u0430 \u0432 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0447\u0430\u0441\u0442\u0438 \u043a\u0430\u0436\u0434\u043e\u0439 \u0438\u0437 \u0442\u0440\u0451\u0445 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442 # coeff1,coeff2,coeff3 -- \u044d\u0442\u043e \u0441\u043e\u0431\u0441\u0442\u0432\u0435\u043d\u043d\u043e \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442\u044b \u043f\u0435\u0440\u0435\u0434 \u0443\u0433\u043b\u043e\u0432\u044b\u043c\u0438 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u0430\u043c\u0438 # r0 -- \u044d\u0442\u043e \"\u0440\u0430\u0437\u043c\u0435\u0440\" \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u043e\u0439 \u0447\u0430\u0441\u0442\u0438, \u0434\u0430\u0432\u0430\u0435\u043c\u043e\u0439 \u0432\u044b\u0440\u0430\u0436\u0435\u043d\u0438\u0435\u043c exp(-|r|\/r0) def orb(r, n1, n2=None, n3=None, coeff1=1., coeff2=1., coeff3=1., r0 = 1.0):      # \u0441\u0447\u0438\u0442\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b Ang     Ang = coeff1 * np.prod(np.power(r,n1))     if not n2 is None:          Ang += coeff2 * np.prod(np.power(r,n2))     if not n3 is None:          Ang += coeff3 * np.prod(np.power(r,n3))      Arg = np.sqrt(np.dot(r,r))\/r0 # \u044d\u0442\u043e \u0430\u0440\u0433\u0443\u043c\u0435\u043d\u0442 \u0434\u043b\u044f \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b      Rad = np.exp(-Arg)   # \u0441\u0447\u0438\u0442\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438          return  Rad * Ang # \u043e\u0442\u0432\u0435\u0442 -- \u044d\u0442\u043e \u043f\u0440\u043e\u0438\u0437\u0432\u0435\u0434\u0435\u043d\u0438\u0435 \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u043e\u0439 \u0438 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442      # \u0434\u0430\u043d\u043d\u0430\u044f \u0444\u0443\u043d\u043a\u0446\u0438\u044f \u043d\u0443\u0436\u043d\u0430, \u0447\u0442\u043e\u0431\u044b \u0441\u0434\u0435\u043b\u0430\u0442\u044c \u043d\u0430\u0437\u0432\u0430\u043d\u0438\u0435 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 # \u043d\u0430\u043c \u043d\u0443\u0436\u043d\u044b \u0442\u043e\u043b\u044c\u043a\u043e \u0441\u0442\u0435\u043f\u0435\u043d\u0438 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b (n) \u0438 \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043f\u0435\u0440\u0435\u0434 \u043d\u0435\u0439 def PowToName(n, coeff=1.):     FunctionName = \"\"      # \u043d\u0435\u0442 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b -- \u043d\u0435\u0442 \u0438 \u0444\u0443\u043d\u043a\u0446\u0438\u0438     if n is None or coeff is None:         return \"\"              # \u0441\u043c\u043e\u0442\u0440\u0438\u043c \u043d\u0430 \u0437\u043d\u0430\u043a \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442\u0430 \u0434\u043b\u044f \u0443\u0433\u043b\u043e\u0432\u043e\u0433\u043e \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u0430         if not coeff is None:         if coeff > 0:             FunctionName += \"+%.2f\"  % coeff         elif coeff &lt;0:             FunctionName += \"%.2f\" % coeff         else:             return \"\" # \u0435\u0441\u043b\u0438 \u043a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043d\u0443\u043b\u0435\u0432\u043e\u0439, \u044d\u0442\u0430 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u0430 \u043d\u0438\u0447\u0435\u0433\u043e \u0438 \u043d\u0435 \u0434\u0430\u0451\u0442          # \u0421\u0442\u0435\u043f\u0435\u043d\u044c \u043f\u043e\u043b\u0438\u043d\u043e\u043c\u0430 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u044b -- \u044d\u0442\u043e \u0442\u0438\u043f \u0444\u0443\u043d\u043a\u0446\u0438\u0438 (s\/p\/d\/f...), \u0435\u0433\u043e \u043c\u044b \u0438 \u043f\u0440\u043e\u0432\u0435\u0440\u044f\u0435\u043c     if np.sum(n) == 0:         FunctionName += \"s\"     elif np.sum(n) == 1:         FunctionName += \"p\"     elif np.sum(n) == 2:         FunctionName += \"d\"     elif np.sum(n) == 3:         FunctionName += \"f\"     elif np.sum(n) == 4:         FunctionName += \"g\"     elif np.sum(n) == 5:         FunctionName += \"h\"     else:         FunctionName += \"X3\" # \u043c\u043e\u0436\u043d\u043e \u043f\u043e\u0439\u0442\u0438 \u0438 \u0434\u0430\u043b\u044c\u0448\u0435, \u043d\u043e \u043a\u043e\u043c\u0443 \u044d\u0442\u043e \u043d\u0443\u0436\u043d\u043e?               # \u0438 \u043f\u043e\u043a\u0430\u0437\u044b\u0432\u0430\u0435\u043c \u043a\u0430\u043a\u0430\u044f \u043a\u043e\u043d\u043a\u0440\u0435\u0442\u043d\u043e \u0440\u0430\u0434\u0438\u0430\u043b\u044c\u043d\u0430\u044f \u043a\u043e\u043c\u043f\u043e\u043d\u0435\u043d\u0442\u0430 \u0443 \u043d\u0430\u0441 \u0442\u0443\u0442 \u0435\u0441\u0442\u044c (\u0435\u0441\u0442\u0435\u0441\u0442\u0432\u0435\u043d\u043d\u043e, \u043d\u0443\u043b\u0435\u0432\u044b\u0435 \u0441\u0442\u0435\u043f\u0435\u043d\u0438 \u043c\u044b \u043d\u0435 \u043f\u0435\u0447\u0430\u0442\u0430\u0435\u043c)     if n[0]>0:         FunctionName += \"x\"         if n[0] > 1:             FunctionName += \"%i\" %  (n[0])     if n[1]>0:         FunctionName += \"y\"         if n[1] > 1:             FunctionName += \"%i\" %  (n[1])     if n[2]>0:         FunctionName += \"z\"         if n[2] > 1:             FunctionName += \"%i\" %  (n[2])                     return FunctionName       # argument parser to get options from command line parser = argparse.ArgumentParser(description='\u0414\u0430\u0432\u0430\u0439\u0442\u0435-\u043a\u0430 \u043d\u0430\u0440\u0438\u0441\u0443\u0435\u043c \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0438 \u043f\u0440\u0438 \u043f\u043e\u043c\u043e\u0449\u0438 \u043c\u0435\u0442\u043e\u0434\u0430 \u041c\u043e\u043d\u0442\u0435-\u041a\u0430\u043b\u043e (\u041c\u041a)!')    parser.add_argument('--c1', type=float, default=1.0,                      help=\"\u041a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043f\u0435\u0440\u0435\u0434 \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0435\u0439\")     parser.add_argument('--c2', type=float, default=0.0,                      help=\"\u041a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043f\u0435\u0440\u0435\u0434 \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0435\u0439\")     parser.add_argument('--c3', type=float, default=0.0,                      help=\"\u041a\u043e\u044d\u0444\u0444\u0438\u0446\u0438\u0435\u043d\u0442 \u043f\u0435\u0440\u0435\u0434 \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0435\u0439\")    parser.add_argument('--nx1', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 x**n \u0434\u043b\u044f \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--nx2', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 x**n \u0434\u043b\u044f \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--nx3', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 x**n \u0434\u043b\u044f \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--ny1', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 y**n \u0434\u043b\u044f \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--ny2', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 y**n \u0434\u043b\u044f \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--ny3', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 y**n \u0434\u043b\u044f \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--nz1', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 z**n \u0434\u043b\u044f \u043f\u0435\u0440\u0432\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--nz2', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 z**n \u0434\u043b\u044f \u0432\u0442\u043e\u0440\u043e\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--nz3', type=int, default=0,                      help=\"\u0421\u0442\u0435\u043f\u0435\u043d\u044c n \u0432 \u0444\u0443\u043d\u043a\u0446\u0438\u0438 z**n \u0434\u043b\u044f \u0442\u0440\u0435\u0442\u044c\u0435\u0439 \u0443\u0433\u043b\u043e\u0432\u043e\u0439 \u0444\u0443\u043d\u043a\u0446\u0438\u0438\")   parser.add_argument('--NumSteps', type=int, default=50000,                      help=\"\u0427\u0438\u0441\u043b\u043e \u0448\u0430\u0433\u043e\u0432 \u0432 \u041c\u041a \u0441\u0438\u043c\u0443\u043b\u044f\u0446\u0438\u0438\")   parser.add_argument('--R0', type=float, default=1.0,                      help=\"\u0420\u0430\u0437\u043c\u0435\u0440 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 R0\")   parser.add_argument('--PartToIgnore', type=float, default=0.1,                      help=\"\u0414\u043e\u043b\u044f \u041c\u041a \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438 \u0432 \u043d\u0430\u0447\u0430\u043b\u0435, \u043a\u043e\u0442\u043e\u0440\u0430\u044f \u0438\u0433\u043d\u043e\u0440\u0438\u0440\u0443\u0435\u0442\u0441\u044f \u0432 \u0440\u0430\u0441\u0447\u0451\u0442\u0430\u0445, \u043a\u0430\u043a \u0441\u0442\u0430\u0434\u0438\u044f \u0443\u0440\u0430\u0432\u043d\u043e\u0432\u0435\u0448\u0438\u0432\u0430\u043d\u0438\u044f\")    parser.add_argument('--MaxBox', type=float, default=0.1,                      help=\"\u0420\u0430\u0437\u043c\u0435\u0440 \u0441\u0442\u0430\u0440\u0442\u043e\u0432\u043e\u0439 \u043e\u0431\u043b\u0430\u0441\u0442\u0438, \u0438\u0437 \u043a\u043e\u0442\u043e\u0440\u043e\u0439 \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u043e \u0432\u044b\u0431\u0438\u0440\u0430\u0435\u0442\u0441\u044f \u0441\u0442\u0430\u0440\u0442\u043e\u0432\u0430\u044f \u0442\u043e\u0447\u043a\u0430 \u041c\u041a-\u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438\")    parser.add_argument('--MaxStep', type=float, default=2.0,                      help=\"\u041c\u0430\u043a\u0441\u0438\u043c\u0430\u043b\u044c\u043d\u0430\u044f \u0434\u043b\u0438\u043d\u0430 \u0448\u0430\u0433\u0430 \u0432 \u041c\u041a-\u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438\")    parser.add_argument('--MaxSmallShift', type=float, default=0.1,                      help=\"\u041c\u0430\u043b\u0435\u043d\u044c\u043a\u0430\u044f \u0434\u043e\u0431\u0430\u0432\u043e\u0447\u043a\u0430 \u043a \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438, \u043a\u043e\u0433\u0434\u0430 \u043c\u044b \u043e\u0441\u0442\u0430\u0451\u043c\u0441\u044f \u0432 \u0442\u043e\u0439 \u0436\u0435 \u0442\u043e\u0447\u043a\u0435 (\u0434\u043b\u044f \u043a\u0440\u0430\u0441\u043e\u0442\u044b)\")    parser.add_argument('--plot',action='store_true',                      help=\"\u0424\u043b\u0430\u0433, \u0447\u0442\u043e\u0431\u044b \u043d\u0430\u0440\u0438\u0441\u043e\u0432\u0430\u0442\u044c \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u044c\u043a\u0443 \u043f\u0440\u0438 \u043f\u043e\u043c\u043e\u0449\u0438 matplotlib\")     # \u0430 \u0442\u0435\u043f\u0435\u0440\u044c \u0440\u0430\u0441\u043f\u0430\u0440\u0441\u0438\u043c \u0430\u0440\u0433\u0443\u043c\u0435\u043d\u0442\u044b args,unknown = parser.parse_known_args()   # \u0441\u043e\u0441\u0442\u0430\u0432\u043b\u044f\u0435\u043c \u0432\u0435\u043a\u0442\u043e\u0440\u0430 \u0438\u0437 \u0441\u0442\u0435\u043f\u0435\u043d\u0435\u0439 \u0434\u043b\u044f \u043a\u0430\u0436\u0434\u043e\u0439 \u0438\u0437 \u0442\u0440\u0451\u0445 \u0443\u0433\u043b\u043e\u0432\u044b\u0445 \u0444\u0443\u043d\u043a\u0446\u0438\u0439 \u0438\u0437 \u0442\u043e\u0433\u043e, \u0447\u0442\u043e \u0432\u0432\u0451\u043b \u044e\u0437\u0435\u0440 NPow1 = np.array([args.nx1,args.ny1,args.nz1] ,dtype=int) NPow2 = np.array([args.nx2,args.ny2,args.nz2] ,dtype=int) NPow3 = np.array([args.nx3,args.ny3,args.nz3] ,dtype=int)               # \u0432\u044b\u0431\u0438\u0440\u0430\u0435\u043c \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u0443\u044e \u0441\u0442\u0430\u0440\u0442\u043e\u0432\u0443\u044e \u0442\u043e\u0447\u043a\u0443 \u0432 \u0437\u0430\u0434\u0430\u043d\u043d\u043e\u0439 \u044e\u0437\u0435\u0440\u043e\u043c \u043e\u0431\u043b\u0430\u0441\u0442\u0438 \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0430 (\u043a\u0443\u0431\u0438\u043a\u0435 \u0432\u043e\u043a\u0440\u0443\u0433 r=(0,0,0) \u0441\u043e \u0441\u0442\u043e\u0440\u043e\u043d\u043e\u0439 2*MaxBox) currentR = np.random.uniform(low=-args.MaxBox, high=args.MaxBox, size=3)  # \u0441\u0447\u0438\u0442\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043e\u0440\u0438\u0431\u0438\u0442\u0430\u043b\u0438 \u0432 \u044d\u0442\u043e\u0439 \u0442\u043e\u0447\u043a\u0435 currentOrbP = orb(r=currentR, n1=NPow1, n2=NPow2, n3=NPow3, coeff1=args.c1, coeff2=args.c2, coeff3=args.c3, r0 = args.R0)  # \u0432 \u044d\u0442\u043e\u0442 \u0444\u0430\u0439\u043b \u043c\u044b \u0431\u0443\u0434\u0435\u043c \u0441\u043e\u0445\u0440\u0430\u043d\u044f\u0442\u044c \u0441\u0430\u043c\u0443 \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u044e FileName=\"orbital_%s%s%s.dat\" % (PowToName(NPow1, args.c1),PowToName(NPow2, args.c2),PowToName(NPow3, args.c3)) outf = open(FileName, \"w\")  AccRate = 0 # \u042d\u0442\u043e \u043a\u043e\u043b\u0438\u0447\u0435\u0441\u0442\u0432\u043e \u043f\u0440\u0438\u043d\u0438\u043c\u0430\u0435\u043c\u044b\u0445 \u0448\u0430\u0433\u043e\u0432 \u0432 \u0442\u0440\u0430\u0435\u043a\u0442\u043e\u0440\u0438\u0438, \u0447\u0442\u043e\u0431\u044b \u0432\u044b\u0447\u0438\u0441\u043b\u0438\u0442\u044c \u0434\u043e\u043b\u044e \u043f\u0440\u0438\u043d\u044f\u0442\u0438\u044f (\u0434\u043e\u043b\u0436\u043d\u0430 \u0431\u044b\u0442\u044c \u0432 \u0434\u0438\u0430\u043f\u0430\u0437\u043e\u043d\u0435 20-60 %)  # \u043d\u0430\u0447\u0438\u043d\u0430\u0435\u043c \u0438\u0442\u0435\u0440\u0430\u0446\u0438\u0438 \u043c\u0435\u0442\u043e\u0434\u0430 \u041c\u0435\u0442\u0440\u043e\u043f\u043e\u043b\u0438\u0441\u0430  for nstep in range(0,args.NumSteps):     # \u0433\u0435\u043d\u0435\u0440\u0438\u0440\u0443\u0435\u043c \u043d\u043e\u0432\u0443\u044e \u0442\u043e\u0447\u043a\u0443, \u043f\u0440\u0438\u043f\u0438\u0441\u044b\u0432\u0430\u044f \u0441\u043b\u0443\u0447\u0430\u0439\u043d\u043e\u0435 \u0441\u043c\u0435\u0449\u0435\u043d\u0438\u0435 \u043a \u043d\u044b\u043d\u0435\u0448\u043d\u0435\u043c\u0443 \u043f\u043e\u043b\u043e\u0436\u0435\u043d\u0438\u044e \u0432 \u043f\u0440\u043e\u0441\u0442\u0440\u0430\u043d\u0441\u0442\u0432\u0435     trialR = currentR + np.random.uniform(low=-args.MaxStep, high=args.MaxStep, size=3)      # \u0441\u0447\u0438\u0442\u0430\u0435\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438 \u0432 \u044d\u0442\u043e\u0439 \u0442\u043e\u0447\u043a\u0435     trialOrbP = orb(r=trialR, n1=NPow1, n2=NPow2, n3=NPow3, coeff1=args.c1, coeff2=args.c2, coeff3=args.c3, r0 = args.R0)             # \u0432\u044b\u0431\u0440\u0430\u0441\u044b\u0432\u0430\u0435\u043c \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u044c \u043f\u0440\u0438\u043d\u044f\u0442\u0438\u044f \u043d\u043e\u0432\u043e\u0439 \u0442\u043e\u0447\u043a\u0438     Ptrial = np.random.uniform(low=0, high=1)          # \u0430 \u0437\u0434\u0435\u0441\u044c \u043c\u044b \u0441\u043e\u0445\u0440\u0430\u043d\u0438\u043c \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435 \u0442\u043e\u0433\u043e, \u043f\u0440\u0438\u043c\u0435\u043c \u043b\u0438 \u043c\u044b \u043d\u043e\u0432\u0443\u044e \u0442\u043e\u0447\u043a\u0443, \u0438\u043b\u0438 \u043d\u0435\u0442     DoWeAccept = False          # \u0441\u043d\u0430\u0447\u0430\u043b\u0430 \u043f\u0440\u043e\u0441\u0442\u043e \u043f\u0440\u043e\u0432\u0435\u0440\u0438\u043c, \u0447\u0442\u043e \u043d\u043e\u0432\u0430\u044f \u0432\u0435\u0440\u043e\u044f\u0442\u043d\u043e\u0441\u0442\u044c (\u043a\u0432\u0430\u0434\u0440\u0430\u0442 \u043e\u0440\u0431\u0438\u0442\u0430\u043b\u0438) \u0432\u044b\u0448\u0435, \u0447\u0435\u043c \u0431\u044b\u043b\u0430, \u0438 \u0435\u0441\u043b\u0438 \u0434\u0430, \u0442\u043e \u043f\u0440\u0438\u043d\u0438\u043c\u0430\u0435\u043c \u0442\u043e\u0447\u043a\u0443     <\/code><\/pre>\n<\/div>\n<\/details>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[],"tags":[],"class_list":["post-395958","post","type-post","status-publish","format-standard","hentry"],"_links":{"self":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/395958","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=395958"}],"version-history":[{"count":0,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=\/wp\/v2\/posts\/395958\/revisions"}],"wp:attachment":[{"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=395958"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=395958"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/savepearlharbor.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=395958"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}