
рд╡реИрдЬреНрдЮрд╛рдирд┐рдХ рд╕рдореБрджрд╛рдп рдореЗрдВ, рдбреЗрдЯрд╛ рд╡рд┐рдЬрд╝реБрдЕрд▓рд╛рдЗрдЬрд╝реЗрд╢рди рдФрд░ рд╕рд┐рджреНрдзрд╛рдВрдд рдбрд┐рдЬрд╛рдЗрди рдмрд╣реБрдд рдорд╣рддреНрд╡рдкреВрд░реНрдг рднреВрдорд┐рдХрд╛ рдирд┐рднрд╛рддреЗ рд╣реИрдВред рдЙрдкрдХрд░рдг рдЬреЛ рд▓реИрдЯреЗрдХреНрд╕ рдХрдорд╛рдВрдбреНрд╕ рдХреЛ рд▓рд╛рдЧреВ рдХрд░рддреЗ рд╣реИрдВ, рдЬреИрд╕реЗ рдХрд┐ рдорд╛рд░реНрдХрдбрд╛рдЙрди рдФрд░ рдореИрдердЬреИрдХреНрд╕, рдЕрдХреНрд╕рд░ рдлреЙрд░реНрдореВрд▓реЛрдВ рдХреА рд╕реБрд╡рд┐рдзрд╛рдЬрдирдХ рдФрд░ рд╕реБрдВрджрд░ рдкреНрд░рд╕реНрддреБрддрд┐ рдХреЗ рд▓рд┐рдП рдЙрдкрдпреЛрдЧ рдХрд┐рдП рдЬрд╛рддреЗ рд╣реИрдВред
рдЬреВрд▓рд┐рдпрд╛ рдХреЗ рд▓рд┐рдП, рдкреИрдХреЗрдЬ рдХрд╛ рдПрдХ рд╕реЗрдЯ рднреА рд╣реИ рдЬреЛ рдЖрдкрдХреЛ LaTeX рдХреЗ рд╕рд┐рдВрдЯреИрдХреНрд╕ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рдиреЗ рдХреА рдЕрдиреБрдорддрд┐ рджреЗрддрд╛ рд╣реИ, рдФрд░ рдкреНрд░рддреАрдХрд╛рддреНрдордХ рдмреАрдЬрдЧрдгрд┐рдд рдХреЗ рд╕рд╛рдзрдиреЛрдВ рдХреЗ рд╕рд╛рде рдорд┐рд▓рдХрд░ рд╣рдореЗрдВ рд╕реВрддреНрд░реЛрдВ рдХреЛ рд╕рдВрднрд╛рд▓рдиреЗ рдХреЗ рд▓рд┐рдП рдПрдХ рд╢рдХреНрддрд┐рд╢рд╛рд▓реА рдЙрдкрдХрд░рдг рдорд┐рд▓рддрд╛ рд╣реИред
рдбрд╛рдЙрдирд▓реЛрдб рдХрд░реЗрдВ рдФрд░ рдЖрдЬ рдХреЗ рд▓рд┐рдП рдЖрдкрдХреА рдЬрд╝рд░реВрд░рдд рдХреА рд╣рд░ рдЪреАрдЬрд╝ рдХреЛ рдХрдиреЗрдХреНрдЯ рдХрд░реЗрдВ
using Pkg Pkg.add("Latexify") Pkg.add("LaTeXStrings") Pkg.add("SymEngine") using Latexify, LaTeXStrings, Plots, SymEngine
LaTeXStrings
LaTeXStrings рдПрдХ рдЫреЛрдЯрд╛ рдкреИрдХреЗрдЬ рд╣реИ рдЬреЛ рдЬреВрд▓рд┐рдпрд╛ рд╕реНрдЯреНрд░рд┐рдВрдЧ рд╢рд╛рдмреНрджрд┐рдХ рдореЗрдВ LaTeX рд╕рдореАрдХрд░рдг рджрд░реНрдЬ рдХрд░рдирд╛ рдЖрд╕рд╛рди рдмрдирд╛рддрд╛ рд╣реИред рдЬреВрд▓рд┐рдпрд╛ рдореЗрдВ рдирд┐рдпрдорд┐рдд рд╕реНрдЯреНрд░рд┐рдВрдЧреНрд╕ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддреЗ рд╣реБрдП рдПрдореНрдмреЗрдбреЗрдб рд▓реИрдЯреЗрдХреНрд╕ рд╕рдореАрдХрд░рдгреЛрдВ рдХреЗ рд╕рд╛рде рдПрдХ рд╕реНрдЯреНрд░рд┐рдВрдЧ рд╢рд╛рдмреНрджрд┐рдХ рджрд░реНрдЬ рдХрд░рдиреЗ рдХреЗ рд▓рд┐рдП, рдЖрдкрдХреЛ рдореИрдиреНрдпреБрдЕрд▓ рд░реВрдк рд╕реЗ рд╕рднреА рдмреИрдХрд╕реНрд▓реИрд╢ рдФрд░ рдбреЙрд▓рд░ рдХреЗ рд╕рдВрдХреЗрддреЛрдВ рд╕реЗ рдмрдЪрдирд╛ рдЪрд╛рд╣рд┐рдП: рдЙрджрд╛рд╣рд░рдг рдХреЗ рд▓рд┐рдП, $ \alpha^2 $
рд▓рд┐рдЦрд╛ рд╣реИ \$\\alpha^2\$
ред рдЗрд╕рдХреЗ рдЕрд▓рд╛рд╡рд╛, рд╣рд╛рд▓рд╛рдВрдХрд┐ IJulia LaTeX рд╕реНрд╡рд░реВрдкрд┐рдд рд╕рдореАрдХрд░рдгреЛрдВ ( MathJax рдХреЗ рдорд╛рдзреНрдпрдо рд╕реЗ) рдХреЛ рдкреНрд░рджрд░реНрд╢рд┐рдд рдХрд░рдиреЗ рдореЗрдВ рд╕рдХреНрд╖рдо рд╣реИ, рдпрд╣ рдирд┐рдпрдорд┐рдд рд▓рд╛рдЗрдиреЛрдВ рдХреЗ рд╕рд╛рде рдХрд╛рдо рдирд╣реАрдВ рдХрд░реЗрдЧрд╛ред рдЗрд╕рд▓рд┐рдП, LaTeXStrings рдкреИрдХреЗрдЬ рдкрд░рд┐рднрд╛рд╖рд┐рдд рдХрд░рддрд╛ рд╣реИ:
LaTeXString
рд╡рд░реНрдЧ ( String
рдХрд╛ рдПрдХ рдЙрдкрдкреНрд░рдХрд╛рд░) рдЬреЛ рдПрдХ рд╕реНрдЯреНрд░рд┐рдВрдЧ (рдЕрдиреБрдХреНрд░рдордг, рд░реВрдкрд╛рдВрддрд░рдг, рдЖрджрд┐ рдХреЗ рд▓рд┐рдП) рдХреА рддрд░рд╣ рдХрд╛рдо рдХрд░рддрд╛ рд╣реИ, рд▓реЗрдХрд┐рди рд╕реНрд╡рдЪрд╛рд▓рд┐рдд рд░реВрдк рд╕реЗ IJulia рдореЗрдВ рдкрд╛рда / рд▓реЗрдЯреЗрдХреНрд╕ рдХреЗ рд░реВрдк рдореЗрдВ рдкреНрд░рджрд░реНрд╢рд┐рдд рд╣реЛрддрд╛ рд╣реИред- рд╕реНрдЯреНрд░рд┐рдВрдЧ рдореИрдХреНрд░реЛрдЬрд╝
L"..."
рдФрд░ L"""..."""
рдЬреЛ рдЖрдкрдХреЛ рдмреИрдХрд╕реНрд▓реИрд╢ рдФрд░ рдбреЙрд▓рд░ рдХреЗ рд╕рдВрдХреЗрддреЛрдВ (рдФрд░ рдпрджрд┐ рдЖрдк рдЙрдиреНрд╣реЗрдВ рдЫреЛрдбрд╝ рджреЗрддреЗ рд╣реИрдВ рддреЛ рдЖрдкрдХреЗ рд▓рд┐рдП рдбреЙрд▓рд░ рдХреЗ рд╕рдВрдХреЗрдд рдЬреЛрдбрд╝рддреЗ рд╣реИрдВ) рд╕реЗ рднрд╛рдЧрдиреЗ рдХреЗ рдмрд┐рдирд╛ рд▓рд╛рдЯреЗрдХреНрд╕ рд╕рдореАрдХрд░рдгреЛрдВ рдХреЛ рджрд░реНрдЬ рдХрд░рдиреЗ рдХреА рдЕрдиреБрдорддрд┐ рджреЗрддреЗ рд╣реИрдВред
S = L"1 + \alpha^2"
REPL рдЖрдЙрдЯрдкреБрдЯ рджреЗрдЧрд╛:
"\$1 + \\alpha^2\$"
рдФрд░ рдЬреБрдкрд┐рдЯрд░ рдкреНрд░рджрд░реНрд╢рд┐рдд рдХрд░реЗрдЧрд╛ :
1+ рдЕрд▓реНрдлрд╛2
рдЕрдиреБрдХреНрд░рдордг рд╕рд╛рдорд╛рдиреНрдп рд░реЗрдЦрд╛рдУрдВ рдХреА рддрд░рд╣ рдХрд╛рдо рдХрд░рддрд╛ рд╣реИ:
S[4:7] "+ \\a"
рдЗрд╕ рддрд░рд╣ рдХреА рд▓рд╛рдЗрдиреЗрдВ рд░реЗрдЦрд╛рдВрдХрди рдХреЗ рдбрд┐рдЬрд╛рдЗрди рдореЗрдВ рдЙрдкрдпреЛрдЧреА рд╣реЛ рд╕рдХрддреА рд╣реИрдВред
x = [-3:0.1:3...] y1 = x .^2 ╬▒ = 10 y2 = x .^4 / ╬▒; plot(x,y1, lab = "\$x^2_i\$") plot!(x,y2, lab = L"x^4_i/\alpha")

Latexify
рдПрдХ рдЕрдзрд┐рдХ рдХрд╛рд░реНрдпрд╛рддреНрдордХ рдкреИрдХреЗрдЬ рд▓реЗрдЯреЗрдХреНрд╕рд┐рдлрд╝ ( рдЧрд╛рдЗрдб ) рд╣реИред рдЗрд╕реЗ рдЬреВрд▓рд┐рдпрд╛ рдСрдмреНрдЬреЗрдХреНрдЯреНрд╕ рд╕реЗ рд▓рд╛рдЯреЗрдХреНрд╕ рдЧрдгрд┐рдд рдЙрддреНрдкрдиреНрди рдХрд░рдиреЗ рдХреЗ рд▓рд┐рдП рдбрд┐рдЬрд╝рд╛рдЗрди рдХрд┐рдпрд╛ рдЧрдпрд╛ рд╣реИред рдпрд╣ рдкреИрдХреЗрдЬ LaTeX рдкреНрд░рд╛рд░реВрдк рдореЗрдВ рднрд╛рд╡реЛрдВ рдХреЛ рд╕реНрдЯреНрд░рд┐рдВрдЧ рдореЗрдВ рдмрджрд▓рдиреЗ рдХреЗ рд▓рд┐рдП рдЬреВрд▓рд┐рдпрд╛ рдХреА рд╣реЛрдореЛрд▓реЙрдЬреА рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реИред рд▓реЗрдЯреЗрдХреНрд╕рд┐рдлрд╛рдИ.рдЬреНрд▓ рд╡рд┐рднрд┐рдиреНрди рдЬреВрд▓рд┐рдпрд╛ рд╡рд╕реНрддреБрдУрдВ рдХреЛ рдкрд░рд┐рд╡рд░реНрддрд┐рдд рдХрд░рдиреЗ рдХреЗ рд▓рд┐рдП рдХрд╛рд░реНрдп рдкреНрд░рджрд╛рди рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕рдореЗрдВ рд╢рд╛рдорд┐рд▓ рд╣реИрдВ:
- рднрд╛рд╡
- рд▓рд╛рдЗрди,
- рд╕рдВрдЦреНрдпрд╛ (рддрд░реНрдХрд╕рдВрдЧрдд рдФрд░ рдЬрдЯрд┐рд▓ рд╕рд╣рд┐рдд),
- рдкреНрд░рддреАрдХрд╛рддреНрдордХ рднрд╛рд╡ SymEngine.jl рд╕реЗ,
- ParameterizedFunctions рдФрд░ ReactionNetworks рдбрд┐рдлрд╝рд░реЗрдВрд╢рд┐рдпрд▓ рдИрдХреНрд╢рдВрд╕ред Jl рд╕реЗ
рд╕рд╛рде рд╣реА рд╕рднреА рд╕рдорд░реНрдерд┐рдд рдкреНрд░рдХрд╛рд░реЛрдВ рдХреЗ рд╕рд░рдгрд┐рдпрд╛рдБред
ex = :(x/(y+x)^2)
fracx left(y+x right)2
str = "x/(2*k_1+x^2)"
fracx2 cdotk1+x2
рд╡рд┐рд╖рдо рддрддреНрд╡реЛрдВ рдХреА рд╕рд░рдгреА:
m = [2//3 "e^(-c*t)" 1+3im; :(x/(x+k_1)) "gamma(n)" :(log10(x))] latexify(m)
\ left [\ start {array} {ccc} \ frac {2} {3} & e ^ {- c \ cdot t} & 1 + 3 \ textit {i} \\ \ frac {x} {x + k_ {1}} рдФрд░ \ Gamma \ left (n \ рджрд╛рдПрдБ) рдФрд░ \ log_ {10} \ рдмрд╛рдПрдБ (x \ рджрд╛рдПрдБ) \\ \ end {рд╕рд░рдгреА} \ right]
\ left [\ start {array} {ccc} \ frac {2} {3} & e ^ {- c \ cdot t} & 1 + 3 \ textit {i} \\ \ frac {x} {x + k_ {1}} рдФрд░ \ Gamma \ left (n \ рджрд╛рдПрдБ) рдФрд░ \ log_ {10} \ рдмрд╛рдПрдБ (x \ рджрд╛рдПрдБ) \\ \ end {рд╕рд░рдгреА} \ right]
рдЖрдк рдПрдХ рдРрд╕рд╛ рдлрд╝рдВрдХреНрд╢рди рдирд┐рд░реНрджрд┐рд╖реНрдЯ рдХрд░ рд╕рдХрддреЗ рд╣реИрдВ рдЬреЛ рдлрд╝рд╛рд░реНрдореБрд▓реЛрдВ рдХреЛ рдкреНрд░рджрд░реНрд╢рд┐рдд рдХрд░рддрд╛ рд╣реИ рдФрд░ рдЙрдиреНрд╣реЗрдВ рдЙрд╕ рд░реВрдк рдореЗрдВ рдмрдлрд░ рдореЗрдВ рдХреЙрдкреА рдХрд░рддрд╛ рд╣реИ рдЬреЛ рдХрд┐ рд╣рдмрд░ рдХреЗ рд▓рд┐рдП рд╕рдордЭ рдореЗрдВ рдЖрддрд╛ рд╣реИ:
function habr(formula) l = latexify(formula) res = "\$\$display\$$l\$display\$\$\n" clipboard(res) return l end habr(ex)
fracx left(y+x right)2
<p>$<!-- math>$$display$$\frac{x}{\left( y + x \right)^{2}}$$display$$</math -->$</p>
рдзреНрдпрд╛рди рд░рдЦреЗрдВ
latexify("x/y") |> display
fracxy
latexify("x/y") |> print
$\frac{x}{y}$
SymEngine
SymEngine рдПрдХ рдкреИрдХреЗрдЬ рд╣реИ рдЬреЛ рдкреНрд░рддреАрдХрд╛рддреНрдордХ рдЧрдгрдирд╛ рдкреНрд░рджрд╛рди рдХрд░рддрд╛ рд╣реИ рдЬрд┐рд╕реЗ рдЖрдк рдЕрдкрдиреЗ рдмреГрд╣рд╕реНрдкрддрд┐ рдореЗрдВ рд▓реЗрдЯреЗрдХреНрд╕рд┐рдлрд╝ рдХреЗ рд╕рд╛рде рдХрд▓реНрдкрдирд╛ рдХрд░ рд╕рдХрддреЗ рд╣реИрдВред
рдЖрдк рд▓рд╛рдЗрдиреЛрдВ рдФрд░ рдЙрджреНрдзрд░рдг ( quote
) рдореЗрдВ рд╡рд░реНрдг рдирд┐рд░реНрджрд┐рд╖реНрдЯ рдХрд░ рд╕рдХрддреЗ рд╣реИрдВ:
julia> a=symbols(:a); b=symbols(:b) b julia> a,b = symbols("ab") (a, b) julia> @vars ab (a, b)
рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХреЛ рдкрд░рд┐рднрд╛рд╖рд┐рдд рдХрд░реЗрдВ рдФрд░ рдЦреВрдмрд╕реВрд░рддреА рд╕реЗ рдЗрд╕реЗ рдкреНрд░рджрд░реНрд╢рд┐рдд рдХрд░реЗрдВ
u = [symbols("u_$i$j") for i in 1:3, j in 1:3] 3├Ч3 Array{Basic,2}: u_11 u_12 u_13 u_21 u_22 u_23 u_31 u_32 u_33 u |> habr
\ рдЫреЛрдбрд╝ рджрд┐рдпрд╛ [\ start {array} {ccc} u_ {11} & u_ {12} & u_ {13} \\ u_ {21} & u_ {22} & u_ {23} \\ u_ {31} рдФрд░ u_ {32} & u_ {33} \\ \ end {рд╕рд░рдгреА} \ right]
\ рдЫреЛрдбрд╝ рджрд┐рдпрд╛ [\ start {array} {ccc} u_ {11} & u_ {12} & u_ {13} \\ u_ {21} & u_ {22} & u_ {23} \\ u_ {31} рдФрд░ u_ {32} & u_ {33} \\ \ end {рд╕рд░рдгреА} \ right]
рдорд╛рди рд▓реАрдЬрд┐рдП рд╣рдорд╛рд░реЗ рдкрд╛рд╕ рд╡реИрдХреНрдЯрд░ рд╣реИрдВ
C = symbols("╬й_b/╬й_l") J = [symbols("J_$i") for i in ['x','y','z'] ] h = [0, 0, symbols("h_z")] 3-element Array{Basic,1}: 0 0 h_z
рдЬреЛ рд╡реЗрдХреНрдЯрд░-рдЧреБрдгрд╛ рд╣реЛрдирд╛ рдЪрд╛рд╣рд┐рдП
using LinearAlgebra ├Ч = cross latexify(J├Чh, transpose = true)
$ $ рдкреНрд░рджрд░реНрд╢рди $ $ \ _ рд╢реБрд░реВ {рд╕рдореАрдХрд░рдг} \ рдЫреЛрдбрд╝ [\ рд╢реБрд░реВ {рд╕рд░рдгреА} {c} J_ {y} \ cdot h_ {z} \\ - J_ {x} \ cdot h_ {z} \\ 0 \\ \ _ рдЕрдВрдд {array} \ right] \ end {рд╕рдореАрдХрд░рдг} $ $ рдбрд┐рд╕реНрдкреНрд▓реЗ $$
рдкреВрд░реНрдг рдореИрдЯреНрд░рд┐рдХреНрд╕ рдЧрдгрдирд╛:
dJ = C*(u*J.^3)├Чh latexify( dJ, transpose = true) habr(ans)
$ $ рдкреНрд░рджрд░реНрд╢рди $ $ \ _ рд╢реБрд░реВ {рд╕рдореАрдХрд░рдг} \ рдЫреЛрдбрд╝ [\ рд╢реБрд░реВ {рд╕рд░рдгреА} {c} \ frac {h_ {z} \ cdot \ left (u_ {21} \ cdot J_ {x} ^ {3} + u_ { 22} \ cdot J_ {y} ^ {3} + u_ {23} \ cdot J_ {z} ^ {3} \ right) \ cdot \ Omega_ {b}} {\ Omega_ {l}} \\ \ frac { - h_ {z} \ cdot \ left (u_ {11} \ cdot J_ {x} ^ {3} + u_ {12} \ cdot J_ {y} ^ {3} + u_ {13} \ cdot \ _ {z} ^ {3} \ right) \ cdot \ Omega_ {b}} {\ Omega_ {l}} \\ 0 \\ \ end {рд╕рд░рдгреА} \ right] \ end {рд╕рдореАрдХрд░рдг} $ $ $ $ рдкреНрд░рджрд░реНрд╢рди
рд▓реЗрдХрд┐рди рдЗрд╕ рддрд░рд╣ рдХреА рдПрдХ рд╕рд░рд▓ рд╢реНрд░реГрдВрдЦрд▓рд╛ рд╕реЗ рдЖрдк рдирд┐рд░реНрдзрд╛рд░рдХ рдХреЛ рдЦреЛрдЬ рд╕рдХрддреЗ рд╣реИрдВ рдФрд░ рдЗрд╕реЗ рд╣реИрдмрд░ рдХреЛ рднреЗрдЬ рд╕рдХрддреЗ рд╣реИрдВ
u |> det |> habr
u_ {11} \ cdot \ left (u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ рдмрд╛рдПрдБ (u_ {33}} \ _ frac) {u_ {13} \ cdot u_ {31}} {u_ {11}} - \ frac {\ _ рдмрд╛рдИрдВ (u_ {23} - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {32} - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}} \ right)
u_ {11} \ cdot \ left (u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ рдмрд╛рдПрдБ (u_ {33}} \ _ frac) {u_ {13} \ cdot u_ {31}} {u_ {11}} - \ frac {\ _ рдмрд╛рдИрдВ (u_ {23} - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {32} - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}} \ right)
Rekursivnenko! рд╡реНрдпреБрддреНрдХреНрд░рдо рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХреА рдЧрдгрдирд╛ рд╕рдВрднрд╡рддрдГ рдЗрд╕реА рддрд░рд╣ рд╕реЗ рдХреА рдЬрд╛рддреА рд╣реИ:
u^-1 |> habr
рд╕реНрдкреЙрдЗрд▓рд░$ $ $ $ $ \ рдмрд╛рдпрд╛рдБ [\ "{рд╕рд░рдгреА} {ccc} \ frac {1 - \ frac {u_ {12} \ cdot \ left (\ frac {- u_ {21}}} {u_ {11}}} \ _ рдкреНрд░рджрд░реНрд╢рд┐рдд рдХрд░реЗрдВ frac {\ left (\ frac {- u_ {31}} {u_ {11}} + \ frac {u_ {21} \ cdot \ left (u_ {32} - \ frac {u_ {12} \ _ \ _ u_ {31) }} {u_ {11}} \ right)} {u_ {11} \ cdot \ left (u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}} рд╕рд╣реА) } \ right) \ cdot \ left (u_ {23} - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right)} {u_ {33} - \ frac [u_ {13] } \ cdot u_ {31}} {u_ {11}} - \ frac {\ _ рдмрд╛рдИрдВ (u_ {23} - \ frac {u_ {13} \ _ cdot u_ {21}} {u_ {11}} рд╕рд╣реА) \ _ cdot \ left (u_ {32} - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} / cdot u_ { 21}} {u_ {11}}}} рд╕рд╣реА)} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}} - \ frac {u_ [13]} \ cdot \ left (\ frac {- u_ {31}} {u_ {11}} + \ frac {u_ {21} \ cdot \ left (u_ {32} - \ frac {u_ {12} \ _ \ _ u_ {31) }} {u_ {11}} \ right)} {u_ {11} \ cdot \ left (u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}} рд╕рд╣реА) } \ right)} {u_ {33} - \ frac {u_ {13} \ cdot u_ {31}} {u_ {11}} - \ frac {\ _ рдмрд╛рдПрдВ (23_) - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {32} - \ frac {u _ {12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}}}} {u_ {11}} & \ frac {\ frac {- u_ {12} \ cdot \ left (1 + \ _ frac {\ _ рдмрд╛рдИрдВ рдУрд░ (u_ {23} - \ frac {u_ {13} \ _ cades u-{21}) {u_ {11}} \ рд░рд╛рдЗрдЯ) \ cdot \ left (u_ {32} - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} \ рд░рд╛рдЗрдЯ)} {\ _ рдмрд╛рдПрдВ (u) {22 } - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {33} - \ frac {u_ {13} \ cdot u_ [31}} {{ u_ {11}} - \ frac {\ left (u_ {23} - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ / 32) - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}} } (рджрд╛рдпрд╛рдБ)} \ рд╕рд╣реА)} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}} + \ frac {u_ {13} \ cdot рдХрд╛ рдмрд╛рдпрд╛рдБ (u_) {32} - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} \ right)} {\ _ рдмрд╛рдПрдВ (u_ {22} - \ frac {u_ {12} \ _ cades u_ {21} } {u_ {11}} \ right) \ cdot \ left (u_ {33} - \ frac {u_ {13} \ cdot u_ {31}} {u_ {11}} - \ frac {\ _ рдмрд╛рдПрдВ (23) } - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {32} - \ frac {u_ {12} \ cdot u_ [31}} {{ u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}} \ rig ht)}} {u_ {11}} & \ frac {\ _ frac {- u_ {13}} {u_ {33} - \ frac {u_ {13} \ cdot u_ {31}} {u_ [11}} - \ frac {\ _ рдмрд╛рдИрдВ (u_ {23} - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {32) - \ frac [u_ {12] } \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}}} + frac { u_ {12} \ cdot \ left (u_ {23} - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right)} {\ _ рдмрд╛рдПрдВ (u_ {22} - frac { u_ {12} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {33} - \ frac {u_ {13} \ cdot u_ {31}) {u_ {11}} - \ _ рдлрд╝реНрд░рд╛рдХ {\ _ рдмрд╛рдПрдВ (u_ {23} - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {32) - \ f \ _ \ _ { 12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}} {рджрд╛рдИрдВ рдУрд░}} } {u_ {11}} \\ \ frac {\ _ frac {- u_ {21}} {u_ {11}} - \ frac {\ _ рдмрд╛рдИрдВ (\ frac {- u_ {31}} {u_ {11}} + \ frac {u_ {21} \ cdot \ left (u_ {32} - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {11} / cdot \ left (u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}} \ рд╕рд╣реА)} \ рджрд╛рдПрдБ) \ cdot \ рдмрд╛рдПрдБ (u_ {23} - \ frac {u_ {13}) \ cdot u_ {21}} {u_ {11}} \ right)} {u_ {33} - \ frac {u_ {13} \ cdot u_ {31}} {u_ {11}} - \ _ рдлрд╝реНрд░рд╛рдХ {\ _ рдмрд╛рдПрдВ (u_ {23} - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {32) - \ f \ _ \ _ { 12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}}} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}} & \ frac {1 + \ _ frac {\ _ рдмрд╛рдИрдВ рдУрд░ (u_ {23} - \ frac {u_ / 13}} cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {32} - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} {рджрд╛рдИрдВ рдУрд░}} {} \ рдмрд╛рдПрдВ (u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {33} - \ frac {u_] 13 # \ _ cdot u_ {31}} {u_ {11}} - \ frac {\ left (u_ {23} - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ _ \ _ \ _ \ _ \ _ (u_ {32} - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} \ cdot u {21}} {u_ {11}}} \ right)}} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}} & \ frac {- \ left (u) {23 } - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right)} {\ _ рдмрд╛рдПрдВ (u_ {22} - \ frac {u_ {12} \ _ cdot u_ [21}}} { u_ {11}} \ right) \ cdot \ left (u_ {33} - \ frac {u_ {13} \ cdot u_ {31}} {u_ {11}} - \ frac {\ _ рдмрд╛рдПрдВ (23) -) \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {32} - \ frac { u_ {12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}} \ right )} \\ \ frac {\ frac {- u_ {31}} {u_ {11}} + \ frac {u_ {21} \ cdot \ left (u_ {32} - \ frac {u_ {12] \ _ cdot u_ {рейрез}} {u_ {резрез}} \ рджрд╛)} {u_ {резрез} \ cdot \ left (u_ {реиреи} - \ frac {u_ {резреи} \ _ cdot u_ {реирез}}} {u_ {резрез}}} рджрд╛рдПрдБ)}} {u_ {33} - \ frac {u_ {13} \ cdot u_ {31}} {u_ {11}} - \ frac {\ рдмрд╛рдИрдВ (u_ {23} - \ frac {u_ / 13} \ _) cdot u_ {21}} {u_ {11}} \ right) \ cdot \ left (u_ {32} - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} {рджрд╛рдИрдВ рдУрд░}} {} u_ {22} - \ frac {u_ {12} \ cdot u_ {21}} {u_ {11}}}} & \ frac {- \ left (u_ {32} - \ frac {u_ {12} \ _ cdot u_ {рейрез}} {u_ {резрез}} \ рджрд╛)} {\ рд╡рд╛рдо (u_ {реиреи}} - \ frac {u_ {резреи} \ _ рдХреВрдЯрдирд╛ u_ {реирез}} {u_ {резрез}} {рджрд╛рдПрдВ) \ _ \ _ \ _ \ _} (u_ {33} - \ frac {u_ {13} \ cdot u_ {31}} {u_ {11}} - \ frac {\ left (u_ {23} - \ frac {u_ {13} \ _dades u_ {21] }} {u_ {11}} \ right) \ cdot \ left (u_ {32} - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ [22} { - {frac {u_ {12} \ cdot u_ {21}} {u_ {11}}} \ right)} рдФрд░ \ рдмрд╛рдПрдВ (u_ {33} - \ frac {u_ {13} \ cdot u_ [31}} {{ u_ {11}} - \ frac {\ left (u_ {23} - \ frac {u_ {13} \ cdot u_ {21}} {u_ {11}} \ right) \ c рдбреЙрдЯ \ рд▓реЗрдлреНрдЯ (u_ {32} - \ frac {u_ {12} \ cdot u_ {31}} {u_ {11}} \ right)} {u_ {22} - \ frac {u_ {12} \ _ \ _ udot { 21}} {u_ {11}}} \ right) ^ {- 1} \\ \ end {рд╕рд░рдгреА} \ right] $ $ рдкреНрд░рджрд░реНрд╢рди $ $
рдпрджрд┐ рдЖрдк рдЕрдкрдиреЗ рдмреНрд░рд╛рдЙрдЬрд░ рдХреЗ рдордЯрдЬреИрдХреНрд╕ рдХреЛ рдЪреЛрдЯ рдкрд╣реБрдВрдЪрд╛рдирд╛ рдЪрд╛рд╣рддреЗ рд╣реИрдВ, рддреЛ рдорд╛рдЗрдирд╕ рджреВрд╕рд░реА рдбрд┐рдЧреНрд░реА (рдЙрд▓рдЯрд╛ рдореИрдЯреНрд░рд┐рдХреНрд╕ рдХрд╛ рд╡рд░реНрдЧ) рдбрд╛рд▓реЗрдВ
рд╡реИрд╕реЗ, SymEngine рдбреЗрд░рд┐рд╡реЗрдЯрд┐рд╡реНрд╕ рдкрд░ рд╡рд┐рдЪрд╛рд░ рдХрд░рддрд╛ рд╣реИ:
dJ[1] |> habr
\ frac {h_ {z} \ cdot \ left (u_ {21} \ cdot J_ {x} ^ {3} + u_ {22} \ _ cdot J_ {y} ^ {3} + u_ {23] \ cdot J_ {z} ^ {3} \ right) \ cdot \ Omega_ {b}} {\ Omega_ {l}}
\ frac {h_ {z} \ cdot \ left (u_ {21} \ cdot J_ {x} ^ {3} + u_ {22} \ _ cdot J_ {y} ^ {3} + u_ {23] \ cdot J_ {z} ^ {3} \ right) \ cdot \ Omega_ {b}} {\ Omega_ {l}}
diff(dJ[1], J[1]) |> habr
frac3 cdotJ2x cdothz cdotu21 cdot Omegab Omegal
рд╡реИрд╕реЗ, рдЬреВрд▓рд┐рдпрд╛ рди рдХреЗрд╡рд▓ рдлреЙрд░реНрдореВрд▓реЗ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░ рд╕рдХрддреА рд╣реИрдВ, рдмрд▓реНрдХрд┐ рд▓рд╛рдЯреЗрдХреНрд╕ рдХреЗ рдП рд╕реЗ рдЧреНрд░рд╛рдлрд┐рдХреНрд╕ рднреАред рдФрд░ рдпрджрд┐ рдЖрдкрдиреЗ рдорд┐рдХрдЯреЗрдХреНрд╕ рд╕реНрдерд╛рдкрд┐рдд рдХрд┐рдпрд╛ рдФрд░ рдкрд╣рд▓реЗ рд╕реЗ рд╣реА рдбрд╛рдЙрдирд▓реЛрдб рдХрд┐рдП рдЧрдП pgfplots , рддреЛ рдЙрдкрдпреБрдХреНрдд рд╡рд╛рддрд╛рд╡рд░рдг рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рдХреЗ рдЖрдк рдЬреВрд▓рд┐рдпрд╛ рдХреЗ рд╕рд╛рде рджреЛрд╕реНрдд рдмрдирд╛ рд╕рдХрддреЗ рд╣реИрдВ, рдЬреЛ рд╣рд┐рд╕реНрдЯреЛрдЧреНрд░рд╛рдо, рддреНрд░рд┐-рдЖрдпрд╛рдореА рд░реЗрдЦрд╛рдВрдХрди, рдЖрдХреГрддрд┐ рдХреЗ рд╕рд╛рде рддреНрд░реБрдЯрд┐рдпреЛрдВ рдФрд░ рд░рд╛рд╣рддреЗрдВ рдмрдирд╛рдиреЗ рдХрд╛ рдЕрд╡рд╕рд░ рдкреНрд░рджрд╛рди рдХрд░реЗрдЧрд╛, рдФрд░ рдлрд┐рд░ рдЗрд╕реЗ рд▓рд╛рдЯреЗрдХ рджрд╕реНрддрд╛рд╡реЗрдЬрд╝ рдореЗрдВ рдПрдХреАрдХреГрдд рдХрд░реЗрдЧрд╛ред
рдпрд╣ рд╕рдм рдлрд╝рд╛рд░реНрдореБрд▓реЛрдВ рдХреЗ рд╕рд╛рде рд╣реИ, рд▓реЗрдХрд┐рди рдкреНрд░рддреАрдХрд╛рддреНрдордХ рдЧрдгрдирд╛рдУрдВ рдХреЗ рд╕рд╛рде рдирд╣реАрдВ: рдЬреВрд▓рд┐рдпрд╛ рдореЗрдВ рдЕрднреА рднреА рдкреНрд░рддреАрдХрд╛рддреНрдордХ рдмреАрдЬрдЧрдгрд┐рдд рдХреЗ рд▓рд┐рдП рдЕрдзрд┐рдХ рдЬрдЯрд┐рд▓ рдФрд░ рджрд┐рд▓рдЪрд╕реНрдк рд╕рдорд╛рдзрд╛рди рд╣реИрдВ, рдЬреЛ рд╣рдо рдирд┐рд╢реНрдЪрд┐рдд рд░реВрдк рд╕реЗ рдЕрдЧрд▓реА рдмрд╛рд░ рдХреБрдЫ рд╕рдордп рдХреЗ рд╕рд╛рде рдирд┐рдкрдЯреЗрдВрдЧреЗред