Help with Newton-Raphson Method

  • Thread starter GSXtuner21
  • Start date
  • Tags
    Method
In summary, the problem is to find the lower limit of a function h(x) integrated from two given bounds, a and b. The Newton-Raphson method is used to find a.
  • #1
GSXtuner21
3
0
Hello All,
I am new to this forum, and was wondering if anyone could help me with a homework problem I am having.

1. Consider the Integral, int(h(x)dx)=c, integrated from lower limit "a" and upper limit "b".
where "a" is unknown and 0<a<b, c>0, h(x)>0, a<=x<=b.

Task: Develop a general expression for determining "a" by a Newton-Raphson procedure, assuming that b,c,h(x) are known.


2. Homework Equations

I know I need to derive the NR equation below in order to solve for the unknown lower limit.

A(p+1)=A(p) -[f(Ap)/f ' (Ap)]

I am just not sure how to apply the integral of a function to the NR method. Any help will be greatly appreciated!
 
Physics news on Phys.org
  • #2
GSXtuner21 said:
Hello All,
I am new to this forum, and was wondering if anyone could help me with a homework problem I am having.

1. Consider the Integral, int(h(x)dx)=c, integrated from lower limit "a" and upper limit "b".
where "a" is unknown and 0<a<b, c>0, h(x)>0, a<=x<=b.

Task: Develop a general expression for determining "a" by a Newton-Raphson procedure, assuming that b,c,h(x) are known.


2. Homework Equations

I know I need to derive the NR equation below in order to solve for the unknown lower limit.

A(p+1)=A(p) -[f(Ap)/f ' (Ap)]

I am just not sure how to apply the integral of a function to the NR method. Any help will be greatly appreciated!

You want to set it up so you're finding the root of an equation. So you want to find [tex]a[/tex] to solve
[tex] \int_a^b h(x)dx - c = 0 [/tex]
So if you know h(x), you could solve
[tex] F(b) - F(a) - c = 0 [/tex]
where F is the anitderivative of h. Here you get the derivative for free because the derivative of F is h, which you know. Or if you don't know F, you could do the integral numerically. For example, using the trapezoidal rule you get something like
[tex] \frac{h(b)}{2} + \sum_j h(x_j) + \frac{h(a)}{2} - c = 0 [/tex]
In this case you would be computing derivatives of h(a) with respect to a.
I don't know if any of those are correct, but they're some suggestions to start with.
 
  • #3
Okay, thanks for the input.

I am not given b,c, or h(x), I am just supposed to assume that they are provided, and then use that to develop an expression for determining "a" using the Newton-Raphson method.

So if I say that F(b)-F(a) -c = 0, can I also say that

[tex]\int[/tex] h*(b - a) - c = 0 ? (integral of h multiplying b subtracted by integral of h multiplying a)

I think I need to use Leibnitz's Rule, but I am a little confused how to apply this equation to the rule. Thanks!
 
Last edited:
  • #4
My take on it

Let

[tex] F(x) = \int_a^b {h(x) dx} = c [/tex]

which is also

[tex] F(b) - F(a) = c [/tex]

Then, for our zero function, let

[tex]G(a) = F(b) - F(a) - C = 0 [/tex]

Note that here, our variable is "a". Thus, our Newton search will be expressed as:

[tex] \[
a_{\nu+1} = a_{\nu} - \frac{G(a_{\nu})}{G'(a_{\nu})}
\]
[/tex]

Here we note that the prime on G indicates a derivative wrt a

For the derivative, we have (assuming b and c are constants)

[tex]\[
G'(a) = \frac{d G(a)}{da} = -h(a)
\]
[/tex]

Thus, the iteration scheme is

[tex] \[
a_{\nu+1} = a_{\nu} + \frac{\int^b_{a_{\nu}} {h(x) dx} - c}{h(a_{\nu})}
\]
[/tex]
 
  • #5
Looks like you want to write your own code. Try looking up a book called Numerical Recipies in C. It's free online.
 
  • #6
hi all, I've had this same question. and when i go to apply the method it seems to make sense that the f(a)/f'(a) term should be smaller than the a term in the equation such that there is actual convergence.

well, no matter what i do i can't get f(a)/f'(a) to be smaller than my inital guess. I'm ripping my hair out. anyone have any suggestions?
 
  • #7
I've had it work for some contrived examples. But you have to be careful as, just like all Newton-Raphson implementations, the starting guess must be sufficiently close to the root.

Also, if a numerical procedure is being used to evaluate the integral, care must be exercised such that the evaluation is performed with sufficient accuracy to allow convergence. Of course, for the iteration process, you don't want to go too crazy here - pick an efficient solver. A delicate balance here.

Can you provide some more detail on the actual problem you are trying to solve and methods you are using?

P.S. I don't think the Numerical Recipies one line books are free anymore. However, efficient solvers are readily available - may need to select based on problem, but Simp's Rule may be adequate for most depending on the precision of "a" one is trying to solve for.
 
  • #8
the function is x^2*exp(x)

integrated from a to b=3 set =c which is 97.704

then we have to find a convergence to a

so i integrated, and solved for the subsequent function of a=0

I plotted the function and saw a root near 1. so it seems I'm on the right track. Now i just need to figure out how to make a loop in matlab...
 
  • #9
Yep, sounds like you're on the right track.

GSXtuner21 - did you have any other questions? (you're OP after all).
 
  • #10
I think you guys have sufficiently answered my questions and then some, so thank you very much for the help!

Cheers
 

What is the Newton-Raphson Method?

The Newton-Raphson Method is an algorithm used to find the roots of a mathematical equation. It is a numerical method that uses an iterative process to approximate the root of the equation.

How does the Newton-Raphson Method work?

The method starts with an initial guess for the root of the equation. Then, using the derivative of the equation, it calculates the slope of the curve at that point and finds the point where the tangent line intersects with the x-axis. This new point becomes the next guess for the root, and the process continues until the desired level of accuracy is achieved.

What are the advantages of using the Newton-Raphson Method?

The Newton-Raphson Method is a fast and efficient way to find the roots of an equation. It converges quickly and can handle complex equations with multiple roots. It also has a high accuracy compared to other numerical methods.

What are the limitations of the Newton-Raphson Method?

The method may not converge if the initial guess is not close enough to the actual root or if the equation has multiple roots that are close together. It also requires knowledge of the derivative of the equation, which may not always be readily available.

How can I use the Newton-Raphson Method to solve a problem?

To use the method, you need to have an equation and an initial guess for the root. You also need to know the derivative of the equation. Then, you can plug in the values and use the iterative process to approximate the root. There are also many calculators and software programs available that can perform the calculations for you.

Similar threads

  • Engineering and Comp Sci Homework Help
Replies
2
Views
1K
  • Engineering and Comp Sci Homework Help
Replies
2
Views
2K
  • Engineering and Comp Sci Homework Help
Replies
4
Views
1K
  • Engineering and Comp Sci Homework Help
Replies
9
Views
1K
  • Engineering and Comp Sci Homework Help
Replies
3
Views
1K
  • General Math
Replies
7
Views
1K
  • Engineering and Comp Sci Homework Help
Replies
6
Views
798
  • Engineering and Comp Sci Homework Help
Replies
2
Views
2K
  • Engineering and Comp Sci Homework Help
Replies
7
Views
1K
  • Engineering and Comp Sci Homework Help
Replies
1
Views
4K
Back
Top