Evaluating an integral of a numeric function

I want to evaluate an integral over a numerically calculated function. Specifically:

  • = proc (t) options operator, arrow; exp(-rho*t)*c(t)^(1-theta)/(1-theta) end proc
  • = int(Uflow(t), t = 0 .. 400, numeric = true)
  • evalf(Uzero)

where c(t) is calculated as a result of using dsolve/numeric to solve a system of 1st order o.d.e.s (and where the numerical values of rho and theta have been specified). However, Maple's just drawing the integral for me, which is not very helpful. Can anyone suggest what I'm doing wrong?

 

Thanks

David

 

 

 

Robert Israel's picture

Integral involving dsolve solution

You may get better results by solving a larger system of differential equations.  Thus to your system involving c(t), add the differential equation

diff(U(t),t) = exp(-rho*t)*c(t)^(1-theta)/(1-theta)

and boundary condition U(0) = 0.  Your Uzero will then be U(400).

 

Comment viewing options

Select your preferred way to display the comments and click "Save settings" to activate your changes.
}