Approximating (x^2+a^2)^(1/4) w/ Error < 1e-6

  • Context: Graduate 
  • Thread starter Thread starter dirk_mec1
  • Start date Start date
  • Tags Tags
    Approximate Function
Click For Summary
SUMMARY

The discussion focuses on approximating the function (x^2 + a^2)^(1/4) with an error less than 1e-6 for x in the range [-1, 1] and a in the real numbers. The user considers using Legendre polynomials but finds that convergence requires a high number of polynomials due to the error decreasing at a rate of n^(3/2). They propose using multiple Taylor polynomials over specific intervals to improve accuracy, acknowledging that Taylor polynomials centered at zero yield significant errors away from that point. The conversation highlights the efficiency of rational function approximations, as referenced in Abramowitz and Stegun's work.

PREREQUISITES
  • Understanding of polynomial approximations, specifically Legendre and Taylor polynomials.
  • Familiarity with error analysis in numerical methods.
  • Knowledge of rational function approximations and their applications.
  • Basic concepts of mathematical functions and their properties.
NEXT STEPS
  • Research the use of Legendre polynomials for function approximation in numerical analysis.
  • Explore Taylor series and their convergence properties for functions beyond their center.
  • Study rational function approximations and their implementation in computational programming.
  • Examine the error analysis techniques for polynomial approximations to ensure desired accuracy.
USEFUL FOR

Mathematicians, computer scientists, and engineers involved in numerical analysis, function approximation, and those seeking efficient algorithms for mathematical computations.

dirk_mec1
Messages
755
Reaction score
13
I want to have the following function approximated (in elementary functions) within a error of 1e-6:

[tex](x^2+a^2)^{1/4}\ \forall x\ \in [-1,1],\ a\in \mathbb{R}[/tex]If I use Legendre polynomials I have to use a lot to get convergence (since if I use an estimate of the error the error decreases with n^3/2 with n the number of polynomials). So I suppose that I can use taylor polynomials, right? But taylored around 0 the polynomials tends to have a large error moving away from zero. So my idea is to use several taylor polynomials for certain intervals. Is this an idea or do I have to do something else?
 
Physics news on Phys.org
From the point of view of computer programming, the most efficient approximations of commonly encountered functions are approximations by rational functions ( ratios of polynomials). For example, if you look at the famous book by Abramawitz and Stegun, you find that sort of approximation.
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 11 ·
Replies
11
Views
3K
  • · Replies 19 ·
Replies
19
Views
5K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 14 ·
Replies
14
Views
3K
Replies
6
Views
2K
  • · Replies 8 ·
Replies
8
Views
3K
Replies
4
Views
2K