First of all, let’s give a few words to explain the title of the article, since it contains such concepts as terms, non-classical equations of mathematical physics, hydrodynamics, and equilibrium magneto-hydrodynamics. In fluid dynamics, the Hicks equation, sometimes also called the Bragg-Hawthorne equation (or Squire-Long equation) (see [1] ), is a partial differential equation that describes the distribution of the stream function for an axisymmetric fluid residue. However, the Grad-Shafranov equation [2] , which arises in plasma physics, also has the same form as the Hicks equation (see [3] - [5] ). Therefore, the equilibrium in ideal magneto-hydrodynamics (MHDs) for a two-dimensional plasma, such as the discovered axisymmetric toroidal plasma equation in a tokamak. This equation reflects the form of the Hicks equation (see [3] [4] ) from fluid dynamics and at the same time, it is a two-dimensional nonlinear elliptic partial differential equation obtained by simplifying the ideal MHD equations to two dimensions, usually in the context of toroidal rotational symmetry, as is relevant for tokamaks. Hence, including directly that the Hicks equation from hydrodynamics, also belongs to the category of equations of non-classical mathematical physics. Now, we need to clarify what does it means in equations of non-classical mathematical physics.

Non-classical equations of mathematical physics:

Definition (non-classical equation in mathematical physics) and applicable cases: Non-classical models of mathematical physics are those whose representations in the form of equations or systems of partial differential equations do not fit within the framework of one of the classical types: elliptic, parabolic, or hyperbolic (1) $\Delta U=U_{xx}+U_{yy}=0$ or equation of Poisson: $\Delta U=U_{xx}+U_{yy}=-f\left(x,y\right)$, $\Delta$-operators Laplace, 2) $U_t=a\Delta U$ and 3) $U_{tt}=a{U}_{xx}$). In particular, non-classical models include those described by mixed-type equations (for example, the Tricomi equation [6] $K\left(z\right)u_{xx}\left(x,z\right)+u_{zz}=0$, $K\left(z\right)=z,z^2,z^m$, degenerate equations (for example, the generalized Keldysh equation $K\left(z\right)u_{zz}\left(x,z\right)+u_{xx}\left(x,z\right)+a\left(x,z\right)u_z=0$ [7] ), or Sobolev-type equations (for example, hydrodynamical equations: the Barenblatt-Zheltov-Kochina equation (see [8] [9] )).

The term “non-classical equations of mathematical physics” at first was introduced by V. N. Vragov (see [10] ), and his students (for example, in case of model mixed type equation with changing time direction by M. A. Nurmamedov (see [11] )), and at the same time, applicable cases to astrophysical fluid dynamics are a modern branch of astronomy that includes fluid mechanics, which deals with the movement of fluids such as the gases that make up stars or any fluid that is found in outer space (see based on the point of view of this definition, applicable indicted problems of astronomy and astrophysics which is noted by M. A. Nurmammadov (see [12] [13] )).

Magnetism is found throughout the universe. Magnetic fields are known to exist in planets, stars, accretion disks, the interstellar medium, galaxies, and active galactic nuclei (see the article on galactic magnetic fields). Often, these magnetic fields are generated and maintained by the action of a magneto-hydrodynamic dynamo and in many cases, the magnetic field dynamically dominates, determining the evolution of the object.

Note that the nature of the equilibrium, be it a tokamak, a reversed field pinch, etc., is largely determined by the choice of two functions that depend on the arguments of the flow functions Ψ and how, as well as the boundary conditions. It must be borne in mind that the fundamental condition for equilibrium is that the forces are equal to zero at all points. For the simplest, assuming no plasma resistance (i.e. ideal MHD conditions), we require that the magnetic and pressure forces are balanced at all points, $J\times B=\nabla P$, where, $J$ is the current density, $B$ is the magnetic field and P is the pressure. Frankly, note that the Grad-Shafranov equation is directly derived from the magneto-hydrodynamic (MHD) equilibrium equations, which are given by the following cases:

1) The momentum equation for a plasma in equilibrium:

$J\times B=\nabla P$ (1)

No time variation ( $\upsilon$ is velocity), therefore must satisfy equality: $\frac{\text{d}\upsilon }{\text{d}t}=0$.

2) Amperes law relating the current density $J$ to the curl of the magnetic flux density $B$.

${\mu }_{0}J=B=\nabla \times B$ (2)

where ${\mu }_0$ is the magnetic permeability, $P=P\left(\psi \right)$ is the pressure.

3) The divergence-free postulate, stating that there are no sources of magnetic flux, that is, no magnetic monopoles:

$\nabla \times B=0$ (3)

Since Equation (1) establishes the force balance needed for equilibrium, the pressure gradient (expansion force) needs to be equal to the magnetic force (confinement force). In this way, the plasma is in equilibrium. It is important to note that the plane is defined by $J$ and everywhere tangent to the surfaces of P.

Note that the above-mentioned two physical effects are individually expressed as follows:

1) The first physical effect occurs when a good conductor is placed in a magnetic field and an electric current is induced in the conductor. Then, according to Lenz’s law, this current creates its own magnetic field. Consequently, this induced magnetic field tends to neutralize the original, and externally supported, field. Thus, effectively eliminating the magnetic field line from the conductor. Then, the result is that the lines of force seem to stretch along with the conductor and they “float with the flow”. Moreover, if the conductor is like a liquid with complex movements, the resulting distribution of the magnetic field can become very complex. In this case, the current will increase until its growth is balanced by Ohmic dissipation.

2) The second key physical effect is dynamic. That is, the Lorentz force (or $J\times B$) will act on the fluid only when the currents are induced by the movement of a conducting fluid through a magnetic field. This will allow it to change the movement of the fluid through the magnetic field. In order to derive the basic MHD equations, first of all, we should always consider the choice of the main variables.

Therefore, in the theory of electromagnetism, we define spatial and temporal changes in either the electromagnetic field or its source, electric charge density and current density (for example, in Figure 1 ). However, it can be calculated from another, namely using Maxwell’s equations, supplemented by suitable boundary conditions.

Finally, in terminology MHD and plasma note that, the fluid is electrically conductive and moves in a magnetic field, which is a type of flow known as magneto-hydrodynamics or MHD for short. However, the main application of MHD is in plasma physics (plasma is a hot ionized gas containing free electrons and ions). It is not at all obvious that plasma can be considered a liquid, since the mean free paths for Coulomb collisions between electrons and ions are macroscopically long. Now we need to give magneto-hydrodynamics concept and for this equilibrium:

Definition (about of magneto-hydrodynamic): Magneto-hydrodynamics (MHDs) (magneto-hydrodynamics or hydro magnetics) is the study of the dynamics of electrically conductive fluids. Examples of such fluids include plasma, liquid metals, and salt water or electrolytes. The word magneto-hydrodynamics comes from the words: “magneto”—magnetic field, “hydro”—fluid and “dynamic”—motion. The field of MHD was initiated by Hannes Alfven, for which he received the Nobel Prize in Physics in 1970. The fundamental concept of MHD is that magnetic fields can induce currents in a moving conducting fluid, which in turn creates forces on the fluid and changes the magnetic field itself. Currently, there is no other global solution to the system of equations, which includes a complex of equations: describing the steady flow of an inviscid magnetic fluid with high electrical conductivity, the equation Faraday’s law of electromagnetic induction, the equation of momentum, the equation of conservation of magnetic flux, with the addition of an equation of state, a detailed energy balance equation, and Poisson’s law of gravity. An important class of problems studied so far is the magnetostatic equilibrium of plasma in magnetic and gravitational fields.

The sphere in the center of the image is a neutron star, curved lines indicate the magnetic field lines of the pulsar (see Figure 2 and Figure 3 ) decrease by many orders of magnitude, and since the magnetic moment of the star is completely preserved, the strength of the magnetic field of the neutron star increases proportionally by many orders of magnitude. The rapid rotation of neutron stars turns them into a pulsar, which emits a narrow beam of energy. It is important to express magneto-hydrodynamic conditions and mixed magnetic field equations in the following.

In this model, the hydrostatic equilibrium equation is modified to include a magnetic term due to the Lorentz force as:

$\nabla P\rho +\nabla \Phi =\left(\nabla \times B\right)\times B4\pi \rho =4\pi \rho F_L,$ (4)

$\nabla \Phi =4\pi \rho G$ (5)

$\Phi =-\frac{GM}{r}$ (6)

where M is the mass of the star, r is the distance from the center, P is pressure, ρ is density, G is constant of gravity, Φ is gravitational potential, $B$ is the vector magnetic field, and is the Lorentz force. We refer to Equation (4) as the hydromagnetic equilibrium equation. In this case, magnetic field configurations satisfying Equation (4) must satisfy an additional constraint determined by taking the rotor of the hydromagnetic equilibrium equation. Since the gradient of the vortex is zero, i.e. $\nabla \times \nabla f=0$, where f is any scalar function, the left side of Equation (4) and the $\nabla P\rho +\nabla {\Phi }_g=f$ disappears with such an operation, and we arrive at the following restriction:

$\nabla \times B\times \left(\nabla \times B\right)\times \rho =0,$ (7)

where ρ is the barotropic equation of state. Since the density ρ is present in this constraint, Equation (6) must be solved together with Equation (4), and mixed magnetic fields, including both poloidal and toroidal components, have stable configurations (see Figure 4 ).

Magnetic field lines lie on nested surfaces of constant magnetic flux. Stellar winds and jets from accretion discs are examples of outflows in which rotation and magnetic fields have important or essential roles [16] - [18] . Using cylindrical polar coordinates (r, φ, z), we examine steady (∂/∂t = 0), axisymmetric (∂/∂φ = 0) models based on the equations of ideal MHD.

$B_r=-\frac{1}{r}\frac{\partial \psi }{\partial z}$, $B_z=\frac{1}{r}\frac{\partial \psi }{\partial r}$, where $\psi \left(r,z\right)$ is the magnetic flux function (as corresponding similar to Figure 5 ). This is related to the magnetic vector potential by $\psi \left(r,z\right)=r{A}_{\varphi }$. The magnetic flux contained inside the circle (r = constant, z = constant). The poloidal (meridional) and toroidal (azimuthal) parts of the magnetic field: $B=\nabla \psi \times \nabla \Phi +B_{\varphi }{e}_{\varphi }=B_p+B_{\varphi }{e}_{\varphi }=-\frac{1}{r}{e}_{\varphi }\times \nabla \psi +B_{\varphi }{e}_{\varphi }$. Note that, in this equation, first addend $-\frac{1}{r}{e}_{\varphi }\times \nabla \psi$ means poloidal (see Figure 6 ), but the second addend $B_{\varphi }{e}_{\varphi }$ means toroidal (see Figure 7 ) magnetic field of $B$, which is mixed magnetic field.

In this model, the magnetic field is divided into poloidal ( $B_p$) and ( $B_t$) toroidal components, where $B_p=\left(B_r,B_{\theta },0\right)$ and $B_t=\left(0,0,B_{\Phi }\right)$. is the axisymmetric mixed model of the poloidal-toroidal magnetic field in Tokmak’s terminology describes mixed elliptic-hyperbolical, when, for example, the following conditions are met for a star (for example):

$\left|B\right|=B<\{\begin{array}{l}5\times {10}^{15}G,\text{purely-poloidal}\\ 2\times {10}^{16}G,\text{purely-toriodal}\end{array}$ (8)

For example, let us take the model that we have considered for the case of rotating pulsars (purpose as strong magnetic field and satisfying conditions of mixed field), where conditions (8) are satisfied. In the physical and mechanical sense, the situation is not so pessimistic, if we consider more slowly rotating pulsars. Consider, for example, the current LIGO upper limit of for the Crab pulsar [19] . For the models we have considered, we must then have. At first glance, these limits may seem less interesting than the limits for a more rapidly rotating system. However, we must remember that young pulsars, such as the Crab Nebula, are believed to have much stronger magnetic fields. Maybe to take data observation of another class of magnetic stars. Along with these additional examples, we will simultaneously solve a specific problem together with a general class of problems based on the studies of the model presented in this paper.

The ideal MHD equilibrium equation for an axisymmetric toroidal plasma can be written as a differential equation for the poloidal magnetic flux function Ψ and is known as the Grad-Shafranov equation. It contains two arbitrary poloidal magnetic flux functions, $P\left(\psi \right)$ and $g\left(\psi \right)$ respectively associated with the pressure and current profiles (see Appendix A),

$\begin{array}{l}\left(\frac{1}{r}\frac{\partial }{\partial r}r\frac{\partial }{\partial r}+\frac{1}{r^2}\frac{{\partial }^2}{\partial {\theta }^2}\right)\psi -\frac{1}{R_0+r\cos \theta }\left(\frac{\partial \psi }{\partial r}\cos \theta -\frac{\partial \psi }{\partial \theta }\frac{1}{r}\sin \theta \right)\\ =-{\mu }_0{\left(R_0+r\cos \theta \right)}^2\frac{\partial P}{\partial \psi }-\frac{\partial g}{\partial \psi }g\left(\psi \right)\end{array}$ (9)

Equilibrium is partly determined by external conditions such as the total current, the applied toroidal magnetic field and by notation, we have:

$K_1\left(r\right)=r^3,K_2\left(\theta \right)={\sin }^2\theta ,a\left(r,\theta \right)=r\left(r-\frac{1}{R_0+r\cos \theta }\cos \theta \right),b\left(r,\theta \right)=-\sin \theta ,$

$K_1\left(r\right)\frac{{\partial }^2\psi }{\partial r^2}+K_2\left(\theta \right)\frac{{\partial }^2\psi }{\partial {\theta }^2}+a\left(r,\theta \right)\frac{\partial \psi }{\partial r}+b\left(r,\theta \right)\frac{\partial \psi }{\partial \theta }=rf\left(r,\theta ,\psi \right)=f_1\left(r,\theta ,\psi \right),$ (10)

$K_1\left(r\right){\psi }_{rr}+K_2\left(\theta \right){\psi }_{\theta \theta }+a\left(r,\theta \right){\psi }_{\theta }+b{\left(r,\theta \right)}_{\theta }=rf\left(r,\theta ,\psi \right)=f_1\left(r,\theta ,\psi \right),$ (11)

$K_1\left(x\right)u_{xx}\left(x,z\right)+K_2\left(z\right)u_{zz}+a\left(x,z\right)u_x+b\left(x,z\right)u_z=0,$ (12)

$K_1\left(x\right)u_{xx}\left(x,z\right)+K_2\left(z\right)u_{zz}+a\left(x,z\right)u_x+b\left(x,z\right)u_z+c\left(x,z\right)=f\left(x,z\right),$ (13)

When $x{K}_1\left(x\right)>0$, for $x\ne 0$, $z{K}_1\left(z\right)<0$, for $z\ne 0$. Equation (13) is called mixed elliptic-hyperbolic type of M. V. Keldysh in two-dimensional problems (see [11] - [13] ). Two-dimensional inhomogeneities of the kind represented by equations of (11) and (12), (13) can be expected to arise in toroidal fields, such as those created in tokamaks. Finally, frankly speaking, the cold plasma and toroidal-poloidal plasma, which in tokamaks, describe mathematical models corresponding to elliptical-hyperbolic equation of mixed type, both Keldysh and Tricomi or Keldysh-Trikomi types (see Figure 8 , in case of sign of coefficients).

In Figure 8 , using sign of coefficients for Equations (12) and (13), the physical, hydrodynamical, aerodynamical and gas dynamic means in terminology number Mache is described (see Figure 9 ).

In this sense, as a study from a mathematical point of view, the presented model, described by the use of non-classical equations of mathematical physics (physical mean see Figure 9 ), which included the Tricomi and Keldysh equations, was studied for the first time. In order to have motivation for non-classical models, we first illustrate the Tricomi, Keldysh, and Tricomi-Keldysh type model equations:

$z{u}_{xx}\left(x,z\right)+u_{zz}=0,K\left(z\right)u_{xx}\left(x,z\right)+u_{zz}=0,K\left(z\right)=z,z^2,z^m$

Equations (10)-(12) are the particular case of Equation (13) by M. A. Nurmammadov (see [11] - [13] ). Therefore, in any case, pure poloidal or pure toroidal magnetic field similar to the so-called tokamak plasma model are also described by equations in frequent derivatives. Introducing stream function as $S\left(r,\theta \right)$ in the following form:

$B_r=\frac{1}{r^2\sin \theta }\frac{\partial S}{\partial \theta },$ $B_{\theta }=-\frac{1}{r\sin \theta }\frac{\partial S}{\partial r},$ (14)

and accordance to the Hopf’s theorem, from Maxwell’s equation and magnetic force balance equation, we get $B_P\cdot \nabla \left(r\sin \left(\theta \right)B_t\right)=0$, and it follows that $B_{\Phi }=\frac{\beta \left(S\right)}{r\sin \left(\theta \right)}$. Where β is some function of the stream function $S\left(r,\theta \right)$. This means that the toroidal part of the field is a function of the poloidal part. In this case for steam function, $S\left(r,\theta \right)$ elliptical equation can be written as:

$\Delta S\left(r,\theta \right)=\frac{B}{R^2}{r}^2{\sin }^2\left(\theta \right)$ (15)

where $\Delta$-is Laplace operator respect to independent variables $\left(r,\theta \right)$. Hence, by means of boundary conditions, the stream function $S\left(r,\theta \right)$ can be found and finally, the $B_r,B_{\theta }$ solution may be expressed by analytical formulas. In this case, if consider tokamak poloidal-toroidal instead of $\left|B\right|=B$ replacing in (8) we obtain mixed poliadal-toriadal potential equation of magnetic field (which is physical corresponding as mathematical model to the poloidal-toroidal field equation):

$\Delta S\left(r,\theta \right)=\{\begin{array}{l}\frac{4.5\times {10}^{15}G}{R^2}{r}^2{\sin }^2\left(\theta \right)\\ \frac{1.9\times {10}^{16}G}{R^2}{r}^2{\sin }^2\left(\theta \right)\end{array}$

If take $B\propto r^{-3}$ as chosen ordinary case (see [20] ) and by astronomical coordinate with $\mu =\cos \left(\theta \right)$, then equation of (15) can be written:

$\Delta S\left(r,\theta \right)=\frac{1}{r{R}^2}\left(1-{\mu }^2\right)=f\left(r,\mu \right)$ or $\frac{{\partial }^{2}S}{\partial r^2}+\frac{1-{\mu }^2}{r^2}\frac{{\partial }^{2}S}{\partial {\mu }^2}=\frac{1}{r{R}^2}\left(1-{\mu }^2\right)=f\left(r,\mu \right)$

$\frac{{\partial }^{2}S}{\partial r^2}+\frac{1-{\mu }^2}{r^2}\frac{{\partial }^{2}S}{\partial {\mu }^2}=\{\begin{array}{l}\frac{4.5\times {10}^{15}G}{R^2}{r}^2\left(1-{\mu }^2\right)\\ \frac{1.9\times {10}^{16}G}{R^2}{r}^2\left(1-{\mu }^2\right)\end{array}$ (16)

In this case when $\theta =0^{\circ }$ $\mu =1$, the equation is degenerating and belongs to the class equation Keldysh-Tricomi type:

$\frac{{\partial }^{2}S}{\partial r^2}+\frac{1-{\mu }^2}{r^2}\frac{{\partial }^{2}S}{\partial {\mu }^2}=f\left(r,\mu \right)$ (17)

If $r\to \infty$ $\theta ={90}^{\circ }$ this equation belongs to the class of Tricomi type and at the same time for $0<\theta \le {90}^{\circ }$ having singularity:

$\frac{{\partial }^{2}S}{\partial r^2}=f\left(r,\mu \right),$ (18)

when $\theta =0^{\circ }$ $\mu =1$, then we have the following class equations:

$\frac{{\partial }^{2}S}{\partial r^2}+\frac{1}{r^2}\frac{{\partial }^{2}S}{\partial {\mu }^2}=f\left(r,\mu \right),$ (19)

Note that the Equation (18) has a general solution:

$S\left(r,\mu \right)={\displaystyle \underset{r_0}{\overset{r}{\int }}{\displaystyle \underset{r_0}{\overset{r}{\int }}f\left(\xi ,\mu \right)\text{d}\xi \text{d}\xi }}+{\displaystyle \underset{r_0}{\overset{r}{\int }}{C}_1\left(\mu \right)\text{d}\xi }+C_2\left(\mu \right)$

where $C{}_1,C_2$ are arbitrary constants. Physical means the corresponding singular point of the sonic line, cutoff, and resonance phenomena having notations—parallel propagation: the principal resonance to be these which occurs at $\theta =0^{\circ }$ $\theta =\frac{\pi }{2}$. As an initial condition, it is assumed that the magnetic field is a dipole, described using the notation as:

$U_{initial}=U_{dipol}=\frac{1-{\cos }^2\left(\theta \right)}{r},$ (20)

If $\theta =0^{\circ }$, then $U_{initial}=U_{dipol}=0$, this means that beyond the critical radius, the magnetic field may be too weak to support plasma coronation and coronation violations can be calculated for magnetic star from the height of the integral conductivity and external plasma flux.

Finally, common case, Equations (17) and (19) can be solved by means of the following recurrent formulas (see Section 6), which is at first applied in the theory of astrophysics.

$\frac{{\partial }^{2}U}{\partial r^2}+\frac{1-{\cos }^2\left(\theta \right)}{r^2}\frac{{\partial }^{2}U}{\partial {\theta }^2}=f\left(r,\mu ,U\right),K\left(r,\theta \right)=\frac{1-{\cos }^2\left(\theta \right)}{r^2}$ (21)

Then, it is easy to see that Equations (18) and (19) are Tricomi type:

$\frac{{\partial }^{2}U}{\partial r^2}+K_1\left(r\right)\frac{{\partial }^{2}U}{\partial {\theta }^2}=f\left(r,0,U\right),\frac{{\partial }^{2}U}{\partial r^2}+\frac{1}{r^2}\frac{{\partial }^{2}U}{\partial {\theta }^2}=f\left(r,0,U\right),K_1\left(r\right)=\frac{1}{r^2}.$ (22)

The same analogical as above from equation of (14), we have equation of Keldysh type:

$r^2\frac{{\partial }^2\psi }{\partial r^2}+\frac{{\partial }^2\psi }{\partial {\mu }^2}=f\left(r,\mu ,\psi \right),$ (23)

Note. The method of straight lines can be based on two concepts. One of them considers this method as a development of the finite difference method. The other concept is based on the ideas of projecting infinite-dimensional spaces of solutions to problems for partial differential equations onto spaces with a finite functional basis.

Problem 1. Find solution of Equation (23) such that satisfies the following conditions:

$\psi \left(r,{\mu }^0\right)={\varphi }^0\left(r\right),\left(r,{\mu }^0+l\right)={\varphi }^1\left(r\right),\left({\alpha }_0\le r\le {\alpha }_1\right),$

$\psi \left({\alpha }_0,\mu \right)={\varphi }^2\left(\mu \right),\psi \left({\alpha }_1,\mu \right)={\varphi }^3\left(r\right),\left({\mu }^0\le \mu \le {\mu }^0+l\right).$

After application of method of straight lines (see note), we obtain recurrent formulas of Equation (22):

$r^2{U^{″}}_k\left(r\right)+\frac{1}{h^2}\left[U_{k+1}\left(r\right)-2U_k\left(r\right)+U_{k-1}\left(r\right)\right]=0,U_0\left(r\right)=U_{n+1}\left({\alpha }_1\right)=0,$ (24)

Main Result: Finally, we found general solution with control tests:

$\begin{array}{l}{U}_k\left(r\right)={\displaystyle \underset{s=1}{\overset{n}{\sum }}\sin \frac{s\pi \left({\mu }^s-{\mu }^0\right)}{l}}\left[A_s{\text{e}}^{{\delta }_{s}r}+B_s{\text{e}}^{-{\delta }_{s}r}\right],\\ {{\delta }^{\prime }}_s^2=\frac{24{\sin }^2\frac{\pi \left({\mu }^s-{\mu }^0\right)}{2l},}{h^2\left(5+\cos \frac{\pi \left({\mu }^s-{\mu }^0\right)}{l}\right)}\left(s=0,1,2,\cdots ,n\right).\end{array}$

where $A_i,B_i$ are arbitrarily constants. Application in case of $n=3,h=0.25$ and $f_n$ right hand of Equation (23) in case for Equation (23) can be written in the following:

$r^2{U^{″}}_1\left(r\right)+16\left[U_3\left(r\right)-2U_2\left(r\right)\right]=f_1,U_0\left(r\right)=U_4\left(r\right)=0,$

$r^2{U^{″}}_2\left(r\right)+16\left[U_3\left(r\right)-2U_2\left(r\right)+U_1\left(r\right)\right]=f_2,$

$r^2{U^{″}}_3\left(r\right)+16\left[U_2\left(r\right)-2U_3\left(r\right)\right]=f_3,U_i\left(0\right)=U_i\left(0,5\right)=0.$

Let’s the particular solution in the form $U_i\left(r\right)=A_i=const,i=1,2,3$. For given $f_k\left(r\right)=r^2{f}_{1,k}\left(r\right)$-right hand of Equation (14) or Equation (23). Using the value of $A_1,A_2,A_3$ and obtained formulas of general solution, we have:

$\begin{array}{c}{U}_1\left(r\right)=\sin 0.25\pi \left(C_1{\text{e}}^{{\delta }_{1}r}+D_1{\text{e}}^{-{\delta }_{1}r}\right)+\sin 0.5\pi \left(C_2{\text{e}}^{{\delta }_{2}r}+D_2{\text{e}}^{-{\delta }_{2}r}\right)\\ +\sin 0.75\pi \left(C_3{\text{e}}^{{\delta }_{3}r}+D_3{\text{e}}^{-{\delta }_{3}r}\right)+\frac{f_2-f_1}{48},\end{array}$

$\begin{array}{c}{U}_2\left(r\right)=\sin 0.5\pi \left(C_1{\text{e}}^{{\delta }_{1}r}+D_1{\text{e}}^{-{\delta }_{1}r}\right)+\sin \pi \left(C_2{\text{e}}^{{\delta }_{2}r}+D_2{\text{e}}^{-{\delta }_{2}r}\right)\\ +\sin 1.5\pi \left(C_3{\text{e}}^{{\delta }_{3}r}+D_3{\text{e}}^{-{\delta }_{3}r}\right)+\frac{2f_2+f_3}{48},\end{array}$

$\begin{array}{c}{U}_3\left(r\right)=\sin 0.75\pi \left(C_1{\text{e}}^{{\delta }_{1}r}+D_1{\text{e}}^{-{\delta }_{1}r}\right)+\sin 1.5\pi \left(C_2{\text{e}}^{{\delta }_{2}r}+D_2{\text{e}}^{-{\delta }_{2}r}\right)\\ +\sin 2.25\pi \left(C_3{\text{e}}^{{\delta }_{3}r}+D_3{\text{e}}^{-{\delta }_{3}r}\right)+\frac{7f_1+2f_3}{48},\end{array}$

${\delta }_1^2=4\ast 16{\sin }^2\frac{\pi }{8}=9.3726;{\delta }_1=3.0611;$

${\delta }_2^2=4\ast 16{\sin }^2\frac{\pi }{4}=32;{\delta }_2=5.6568,$

${\delta }_3^2=4\ast 16{\sin }^2\frac{3\pi }{8}=54.6274;{\delta }_2=\mathrm{7.3910.}$

Hence, account into boundary value conditions and symmetrically, take $C_i=D_i$ obtain inhomogeneous algebraically system equations respect to $C_1,C_2,C_3$ and after substituting its value, which is found above as the general solution $U_i\left(r\right),i=1,2,3$. Accordance to testing values of $f_1=8$; $f_2=3.5$; $f_3=-10$ finding constants are $C_1=-0.0266$; $C_2=0$; $C_3=-0.00009$.

Thus, we have $U_1\left(r\right)=U_3\left(r\right)=0.0938-0.0376Ch\left(3.0611r\right)-0.00013Ch\left(7.3910r\right)$; $U_2\left(r\right)=0.01250-0.0532Ch\left(3.0611r\right)+0.00018Ch\left(7.3910r\right)$; hence, approximately solution in the center of square will be $U_2\left(0\right)=0.0720$; but at the same time, the exact solution is $\psi \left(r,\mu \right)=u\left(0,0\right)=0.07360$; $error=0.0016$. It means that this method is best and at first applied for the stars model.

After substituting these flux band stream functions, it can be found axisymmetric mixed model of the poloidal-toroidal magnetic field. Hence, the poloidal (meridional) and toroidal (azimuthal) parts of the magnetic field may be expressed by means above general solution and using boundary conditions special solution which choosing near to exact solution. That is, at first, we applied new methods and determined initial boundary conditions in the following: As an initial condition, it is assumed that the magnetic field is a dipole, described using the previous notation as: $U_{initial}=0$, $U_{initial}=\frac{1-{\mu }^2}{r}$. Finally, these obtained solution provide the conditions of magneto-hydrodynamic equilibrium on the magnetosphere of magnetic stars.

Accordance equality of force balance in the magnetosphere, standard form defines the plasma and the Alfven Mach number for magneto-hydrodynamic equilibrium can used (Appendix B).

Remark (about open question): In almost all cases of numerical modeling, the recreation of real physical conditions in the depths of stars will remain unattainable in the foreseeable future, regardless of the expected increase in computing power and therefore need to create retrospective physical and mathematical models.

Note. However, as already seen by the applications of mathematical models, that in the theory of magnetic field still has many shortcomings, so it is necessary to make more and more new theoretical assumptions and then compare the reliability of the practical applied. Therefore, our duty is to seek and try to find new theoretical assumptions on the basis of the existing justified theory and the facts confirming these theories of observational results. Although, in this work, more or less we tried to observe these necessary only that we noted approaches, which confirms the above all assumptions.

In this work, first of all, we give the definition of non-classical equations of mathematical physics and their applicable cases presented, which are used by the author’s works. Considering the concept of MHD equilibrium and the field of basic applications of MHD, the model of the MHD equilibrium field of a magnetic star is constructed. Using the main physical effects in MHD to model MHD equilibrium and the field of a magnetic star, we chose a mathematical model. The main results consist of remonstration magneto-hydrodynamic conditions of mixed magnetic field equations. A new equation of model and axisymmetric mixed magnetic field model is obtained into the Keldysh-Tricomi type equations for magnetic stars. Finally, the application of some new methods to the solution of the equation for the hydromagnetic-dynamic equilibrium of magnetic stars is introduced. For the linear and quasi-linear (in the case of the right hand of the equation) cases, we apply the straight-line method to obtain a recurrent formula for the calculation of the stream and magnetic flux functions. For this reason, the corresponding problems are solved and given numerical calculations, which are tested.