MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Newbie question on FindRoot and NIntergrate

  • To: mathgroup at
  • Subject: [mg80437] Newbie question on FindRoot and NIntergrate
  • From: "Hoa Bui" <hoabui05 at>
  • Date: Wed, 22 Aug 2007 04:49:45 -0400 (EDT)

Dear all,

Please help me with this problem:

I have a series of points:

Define one of my function as the trapezoidal area created by these
points up to some x:
linNx = Interpolation[points, InterpolationOrder -> 1];
ff[x_] := Integrate[linNx[s], {s, 0.0001, x}]

Define the second function as a power law:
gg[x_] := 0.9 x^(1/0.9)

My third function is the root of the equation:
hh[x_] := FindRoot[ff[y] == gg[x], {y, 0.0001, 0.019}];

I can evaluate h at specific values of x, e.g. h[0.001] ({y ->
0.00107654}), h[0.01] ({y -> 0.0101302}), etc...

Now say I want to use hh[x] in an integral,
NIntergrate[hh[x],{x, 0.0001, 0.019}]
obviously it doesn't work, and I have tried using Solve instead of
FindRoot but it also did not output a numerical value for the integral
because the inverse function is not in closed form or something..

Is there a way for me to compute the integral of hh[x] ?

Thank you all so much,
Hoa Bui

  • Prev by Date: Re: But is it worth to use Mathematica ?
  • Next by Date: Re: Re: Cell Bracket Symbols
  • Previous by thread: Re: FWHM, InterpolationFunction & Solve
  • Next by thread: Re: Newbie question on FindRoot and NIntergrate