Description:FIELD OF THE INVENTION

The present invention relates to the field of numerical analysis and computational mathematics, with particular emphasis on the development of efficient, stable, and highly accurate numerical methods for solving nonlinear singular boundary value problems such as the Lane–Emden equations. These equations arise prominently in astrophysics, theoretical physics, and engineering, where they are used to model the thermal behavior of polytropic stars, stellar structures, isothermal gas spheres, and other self-gravitating systems. Traditional analytical solutions exist only for specific indices, and conventional numerical approaches often face challenges due to singularities at the origin and nonlinearities inherent in the governing equations. The invention specifically addresses these limitations by introducing a B-spline based collocation method that capitalizes on the smoothness, local support, and approximation versatility of B-spline functions to construct a robust computational framework. By formulating the Lane–Emden problem into a system of algebraic equations through collocation at carefully chosen nodes, the method achieves optimal accuracy and reduces computational complexity compared to standard finite difference, finite element, or spectral approaches. The invention thus belongs to the field of computational techniques designed for scientific computing, numerical simulation, and applied mathematics, with applications extending to astrophysical modeling, fluid mechanics, plasma physics, and other disciplines requiring solutions to nonlinear singular problems. By ensuring high convergence rates, stability, and adaptability to various boundary conditions, the B-spline based collocation framework provides an innovative contribution to computational science, facilitating both theoretical investigations and practical applications in physics, engineering, and applied numerical research.

Background of the proposed invention:

The Lane–Emden problem, a class of singular nonlinear ordinary differential equations, has attracted extensive attention in applied mathematics, physics, and engineering because of its direct relevance to astrophysical modeling and its mathematical complexity that poses significant challenges for accurate numerical approximation. Originating from the study of polytropic processes in gaseous spheres, the Lane–Emden equation describes the dimensionless density distribution of a spherically symmetric self-gravitating polytropic fluid under hydrostatic equilibrium, making it a cornerstone of stellar structure theory and astrophysical simulations. In its general form, the Lane–Emden equation is expressed as
(
1
/
𝑥
2
)
(
𝑑
/
𝑑
𝑥
)
(
𝑥
2
𝑑
𝑦
/
𝑑
𝑥
)
+
𝑦
𝑛
=
0
(1/x
2
)(d/dx)(x
2
dy/dx)+y
n
=0 with specified boundary conditions, where
𝑦
(
𝑥
)
y(x) represents the dimensionless density,
𝑥
x is the dimensionless radial coordinate, and
𝑛
n is the polytropic index. The singularity at the origin (
𝑥
=
0
x=0) complicates the problem because classical solution methods fail to properly handle the division by zero, and nonlinearities exacerbate computational instabilities, making the development of reliable numerical schemes crucial for practical applications. While analytical solutions exist for a few special cases such as
𝑛
=
0
,
1
,
n=0,1, and
5
5, most physically significant cases require robust numerical approximations, motivating researchers to devise efficient computational techniques over the last few decades. Traditional methods like finite difference schemes, shooting techniques, perturbation methods, finite element analysis, and spectral approaches have been widely applied to tackle the Lane–Emden problem, yet each method suffers from inherent limitations. Finite difference schemes, though conceptually simple, struggle with accuracy and stability near singularities and require extremely fine meshes to maintain precision, leading to high computational cost. Shooting methods, which transform the boundary value problem into an initial value problem, often face difficulties with convergence due to sensitivity to initial guesses and instability when dealing with stiff or highly nonlinear systems. Perturbation approaches rely on small parameters, restricting their applicability to cases where perturbation expansions converge reliably, which is not always guaranteed for general Lane–Emden problems. Finite element methods provide flexibility and accuracy but demand significant computational resources and complex implementation, making them less practical for large-scale astrophysical modeling. Spectral methods, leveraging global basis functions like Chebyshev or Legendre polynomials, achieve exponential convergence rates for smooth problems, but their reliance on global basis sets reduces their adaptability to handle local singularities effectively and often results in ill-conditioned systems. These shortcomings underline the necessity of alternative frameworks that combine the accuracy of spectral methods with the flexibility of localized approaches to singular behavior.

B-splines, or basis splines, have emerged as powerful tools in approximation theory and computational mathematics due to their desirable properties of local support, piecewise polynomial structure, high smoothness, and numerical stability. Unlike global polynomials, B-splines avoid oscillations, reduce numerical instabilities, and allow adaptive refinement in regions with singularities or sharp gradients, making them ideal candidates for solving boundary value problems with singularities such as the Lane–Emden equation. Their compact support significantly reduces computational complexity because each basis function affects only a limited portion of the domain, leading to sparse algebraic systems that can be solved efficiently using modern numerical solvers. Furthermore, the continuity and differentiability of B-spline basis functions enable accurate approximation of derivatives, which is essential for differential equation solving. When integrated into a collocation framework, B-splines provide a systematic method for approximating the solution by enforcing the governing equation at discrete collocation points, thereby converting the original singular boundary value problem into a system of algebraic equations. This collocation approach ensures that the nonlinear characteristics of the Lane–Emden equation are captured accurately, while the flexibility of B-splines allows precise handling of boundary conditions and singularities near the origin.

In the context of the Lane–Emden problem, the B-spline based collocation method offers several distinct advantages over conventional approaches. First, the method provides high accuracy with relatively coarse discretization, meaning that fewer nodes are required to achieve the same precision compared to finite difference or finite element methods. Second, the locality of B-spline functions minimizes the propagation of errors across the domain, enhancing stability in long-range astrophysical computations. Third, the sparsity of the resulting algebraic system makes the method computationally efficient, even when solving problems with large domains or higher-order polytropic indices. Fourth, the ability to represent smooth curves and functions with adjustable polynomial degrees offers flexibility for refining approximation accuracy without significantly increasing computational cost. Moreover, the B-spline collocation method can be easily extended to multi-dimensional and more complex singular problems, which are often encountered in physical sciences and engineering beyond astrophysics. These include plasma physics, nuclear reaction modeling, stellar oscillations, heat conduction in spherical bodies, isothermal gas spheres, and other nonlinear singular systems.

Historically, numerical approximations of the Lane–Emden problem have sought to balance accuracy, stability, and computational feasibility, but many methods compromise one aspect to achieve another. For example, high-order finite difference methods improve accuracy but at the expense of stability and computational load, while low-order methods are stable but insufficiently precise for astrophysical applications. The B-spline collocation approach overcomes this trade-off by offering both accuracy and stability with reduced computational effort, making it especially suitable for large-scale astrophysical modeling where computational resources are often a limiting factor. Additionally, the modular nature of B-splines enables straightforward implementation of adaptive mesh refinement strategies, allowing the algorithm to concentrate computational effort in critical regions, such as near singularities or steep gradients, while using coarser meshes elsewhere. This adaptivity not only improves accuracy but also enhances efficiency, which is critical when dealing with highly nonlinear systems where iterative solvers must be applied.

The motivation for developing a B-spline based collocation method for solving the Lane–Emden problem therefore stems from both theoretical and practical considerations. From a theoretical standpoint, the method leverages the mathematical properties of B-splines to handle singularities and nonlinearities with optimal approximation accuracy. From a practical standpoint, it provides a computational framework that is simple to implement, efficient in execution, and versatile across different problem types. Numerical experiments conducted on benchmark Lane–Emden problems with various polytropic indices demonstrate the method’s superior performance, as it achieves rapid convergence rates, excellent agreement with known analytical solutions, and high consistency with existing numerical results. Error analysis reveals that the method delivers optimal accuracy in terms of discretization error reduction and maintains stability across a wide range of test cases, thereby validating its reliability for practical applications.

In summary, the background of the proposed invention lies in addressing the inherent challenges posed by the Lane–Emden problem—namely singularities, nonlinearities, and boundary condition enforcement—through the innovative application of B-spline based collocation techniques. Existing numerical methods have proven either insufficiently accurate, computationally demanding, or unstable for certain classes of Lane–Emden problems, which has limited their applicability in astrophysical research and related scientific fields. By exploiting the advantageous properties of B-splines within a collocation framework, the invention overcomes these limitations and provides a robust, efficient, and accurate numerical solution strategy. This not only advances the state of the art in computational methods for singular boundary value problems but also has profound implications for the broader scientific community, particularly in fields such as astrophysics, plasma physics, stellar dynamics, and engineering disciplines where nonlinear singular differential equations frequently arise. The invention thus builds on the foundation of numerical mathematics while extending its application to real-world scientific challenges, contributing significantly to the advancement of computational astrophysics and applied mathematics.

Summary of the proposed invention:

The proposed invention introduces a novel B-spline based collocation method specifically designed to solve the Lane–Emden problem with optimal accuracy, stability, and computational efficiency, addressing long-standing challenges associated with the singular and nonlinear nature of this class of equations. The Lane–Emden equation, which arises prominently in astrophysics, fluid mechanics, and plasma physics, is used to model the structure of polytropic stars, self-gravitating gaseous spheres, and similar systems, but its singularity at the origin and nonlinear behavior complicate numerical treatment. Traditional numerical methods such as finite difference, shooting, perturbation, finite element, and spectral techniques, while widely used, have inherent limitations such as slow convergence, instability near singular points, difficulty handling nonlinearities, or excessive computational demands. The invention overcomes these drawbacks by employing B-splines—piecewise polynomial functions with local support and high smoothness—within a collocation framework to approximate the solution. The essence of the method is to represent the unknown solution of the Lane–Emden equation as a linear combination of B-spline basis functions defined over a partitioned domain, and then enforce the governing differential equation at carefully chosen collocation points, transforming the problem into a sparse algebraic system that can be solved efficiently using modern solvers. This design ensures that the singularity at the origin is managed smoothly, as B-splines’ local flexibility allows for accurate representation of the solution behavior near critical points, while their global smoothness ensures stability and precision throughout the computational domain. The method is further adaptable in terms of polynomial degree and mesh refinement, allowing users to achieve desired levels of accuracy without incurring significant computational overhead, and the sparse system structure drastically reduces storage requirements and accelerates solution times compared to dense matrix formulations. Unlike global polynomial or spectral approaches that are highly sensitive to oscillations and ill-conditioning, the B-spline collocation method is numerically stable and maintains error localization, preventing propagation of approximation errors across the domain. Extensive numerical experiments conducted on benchmark Lane–Emden problems with various polytropic indices confirm the superiority of the proposed method, as results exhibit high agreement with exact analytical solutions available for specific cases, and demonstrate improved precision compared to conventional numerical techniques for general cases. Convergence analysis reveals that the error decreases rapidly with refinement, confirming the optimal accuracy of the framework. Moreover, the invention’s efficiency makes it suitable not only for single-equation astrophysical models but also for larger, more complex systems of nonlinear singular equations encountered in applied physics and engineering contexts. Another strength of the invention lies in its modularity: the B-spline based collocation framework can be extended seamlessly to higher-dimensional problems, systems of equations, or cases involving additional physical constraints, thus broadening its applicability beyond the classical Lane–Emden formulation. The invention provides a balance of accuracy, stability, and efficiency that is rarely achieved simultaneously in singular boundary value problem solving, enabling researchers and practitioners to model astrophysical and physical systems with greater confidence and reduced computational cost. From a practical standpoint, the method can be implemented with straightforward coding efforts, as B-spline libraries are widely available, and the collocation approach translates directly into algebraic formulations that are compatible with standard numerical solvers. The invention therefore not only contributes to the theoretical advancement of numerical analysis but also offers a ready-to-use computational tool for scientists, engineers, and researchers working in astrophysics, plasma dynamics, and related fields. By ensuring high convergence rates, low computational complexity, robustness against singularities, and adaptability to various problem classes, the proposed B-spline based collocation method establishes itself as a transformative technique in computational mathematics. In conclusion, the summary of the proposed invention highlights its central contribution: the design of a computationally efficient, accurate, and stable framework that exploits the advantageous properties of B-splines within a collocation scheme to solve the Lane–Emden problem and related singular nonlinear boundary value problems with optimal accuracy. This innovation bridges the gap left by existing methods, paving the way for more reliable simulations in astrophysical modeling and extending to diverse applications across physical sciences and engineering, thereby setting a new benchmark in the field of numerical computation for singular differential equations.
Brief description of the proposed invention:

The proposed invention provides a comprehensive and efficient computational framework that employs a B-spline based collocation method to solve the Lane–Emden problem, which is a singular, nonlinear boundary value problem of fundamental importance in astrophysics and applied mathematics, and the brief description of the invention centers around how this framework is designed, how it operates, and why it achieves superior performance compared to existing numerical techniques. At its core, the Lane–Emden problem models the dimensionless density distribution in polytropic stars or self-gravitating gaseous spheres under hydrostatic equilibrium, and it is mathematically expressed as
(
1
/
𝑥
2
)
(
𝑑
/
𝑑
𝑥
)
(
𝑥
2
𝑑
𝑦
/
𝑑
𝑥
)
+
𝑦
𝑛
=
0
(1/x
2
)(d/dx)(x
2
dy/dx)+y
n
=0 with initial boundary conditions at the origin
𝑦
(
0
)
=
1
,
𝑦
′
(
0
)
=
0
y(0)=1,y
′
(0)=0, where
𝑛
n is the polytropic index, but the singularity at
𝑥
=
0
x=0 and the nonlinear term
𝑦
𝑛
y
n
make direct analytical or standard numerical treatments challenging. The invention addresses this challenge by constructing the approximate solution
𝑦
(
𝑥
)
y(x) as a linear combination of B-spline basis functions defined on a partitioned interval, thereby ensuring piecewise polynomial representation with desirable smoothness properties, while leveraging the local support of B-splines to minimize computational complexity. The collocation method is then applied, where the Lane–Emden differential equation is enforced at strategically selected collocation points within the domain, transforming the continuous singular differential equation into a finite system of sparse nonlinear algebraic equations. These algebraic equations are solved using robust iterative numerical solvers, with the sparsity structure resulting from the local nature of B-splines allowing for efficient computation even on large domains. The design is inherently modular, meaning that the degree of the B-spline basis can be adjusted to balance smoothness and computational load, while mesh refinement can be selectively applied in regions of interest such as near the singular origin or in areas with steep solution gradients, ensuring accuracy without excessive computational burden. Unlike finite difference methods that require very fine discretization or spectral methods that struggle with local singularities, the B-spline collocation technique maintains high accuracy with relatively coarse meshes, as the localized polynomial approximation captures essential features of the solution without oscillatory artifacts. This not only reduces memory and time requirements but also enhances numerical stability, as local errors remain contained instead of propagating across the entire domain. The invention also provides flexibility for generalization beyond the classical Lane–Emden equation, as the same framework can be extended to higher-order singular equations, multi-dimensional problems, or coupled nonlinear systems frequently encountered in astrophysics, plasma dynamics, and thermal modeling of spherical bodies. Extensive testing of the invention demonstrates its practical value: for benchmark Lane–Emden problems with indices where analytical solutions exist (such as 𝑛 = 0,1,5), the results produced by the invention match exactly with theoretical values, validating the method’s precision, and for more general cases where analytical solutions are unavailable, the results exhibit close agreement with high-accuracy numerical solutions available in the literature, confirming its reliability. Furthermore, error and convergence analysis shows that the proposed invention achieves optimal convergence rates, with errors decreasing systematically as the mesh is refined or the polynomial degree of the B-splines is increased, establishing its theoretical robustness. Another important feature of the invention is its adaptability to computational environments, as it can be implemented with existing mathematical software and optimized linear algebra solvers, requiring no special-purpose hardware or overly complex programming infrastructure, making it accessible to a wide community of researchers. By addressing the deficiencies of traditional approaches—such as instability in shooting methods, slow convergence in finite differences, high computational cost in finite elements, and ill-conditioning in spectral methods—the invention provides a balanced solution that delivers accuracy, efficiency, and stability in one unified framework. Its novelty lies not only in the application of B-splines to this singular problem but also in the careful integration of collocation techniques that maximize the strengths of B-splines while directly tackling the singularity and nonlinearity. The broader impact of the invention extends to numerous scientific and engineering applications: in astrophysics, it can aid in the modeling of stellar interiors, polytropic gas spheres, and white dwarf stars; in physics, it can be applied to plasma confinement, thermodynamic equilibrium problems, and spherical diffusion models; and in engineering, it can be used in heat transfer problems involving spherical geometries and nonlinear reaction-diffusion systems. The invention thus not only advances the mathematical toolbox available for solving nonlinear singular boundary value problems but also provides a practical computational methodology that can be readily employed across disciplines. In essence, the brief description of the proposed invention emphasizes that it is an innovative, adaptable, and computationally efficient method that exploits the powerful approximation properties of B-splines within a collocation framework to deliver optimal accuracy and stability for the Lane–Emden problem, overcoming the drawbacks of earlier methods and opening pathways for broader applications in science and engineering.
, Claims:We Claim:
1. A numerical method for solving the Lane–Emden problem, comprising the steps of:
(a) approximating the unknown solution of the Lane–Emden equation using a linear combination of B-spline basis functions;
(b) selecting collocation points across the computational domain; and
(c) enforcing the Lane–Emden differential equation at said collocation points to generate a system of algebraic equations that represents the solution with optimal accuracy.

2. The method of claim 1, wherein the B-spline basis functions are chosen for their local support, smoothness, and piecewise polynomial structure, enabling accurate handling of singularities near the origin.

3. The method of claim 1, wherein the system of algebraic equations generated is sparse due to the compact support of B-spline functions, thereby reducing memory requirements and computational cost.

4. The method of claim 1, wherein adaptive mesh refinement is applied, providing finer partitions near singularities or steep gradients and coarser partitions elsewhere to balance accuracy and efficiency.

5. The method of claim 1, wherein the polynomial degree of the B-splines is adjustable, allowing optimization of convergence rate and computational performance for different polytropic indices.

6. The method of claim 1, wherein the resulting nonlinear algebraic system is solved using iterative solvers, including Newton–Raphson or quasi-Newton methods, ensuring stable and efficient convergence.

7. The method of claim 1, wherein the framework is extended to multi-dimensional singular boundary value problems, enabling applications beyond the classical one-dimensional Lane–Emden formulation.

8. The method of claim 1, wherein error analysis confirms that the numerical solution converges optimally, with error reduction directly proportional to mesh refinement and spline order.

9. The method of claim 1, wherein the computational framework is implementable in standard scientific programming environments using available B-spline libraries and linear algebra packages, ensuring broad accessibility.

10. The method of claim 1, wherein the application extends to astrophysical modeling, plasma physics, thermodynamics of self-gravitating systems, and engineering problems involving nonlinear singular differential equations in spherical geometries.