I have a problem: I have an ODE, dT/dy, that depends on y and k, but k itself depends on rho and T. The only problem is rho cannot be obtained analytically--it's a complicated function. However, I can use a rootfinding algorithm,(fsolve), to determine rho at given values of y, thereby making my ODE solvable. I want to use dsolve to numerically solve this ode. The problem is I can't seem to be able to input the fsolve equation as an argument to the dsolve command. It keeps telling me : Error in fsolve, y is not solved for. But my intention is to use y to solve for T and rho in dsolve. If anyone has any idea I would greatly appreciate it.

Hey guys Does anyone know how to include an fsolve command incorporated into the dsolve command? In other words, as dsolve iterates the runge kutta procedure, it uses the inital values to use fsolve which would then be used over again in the rk procedure.. Thanks B-rad

Does it surprise anyone that Maple can't run a loop from i=10^8 to i=10^9 ? Or is my computer just too slow? Thanks Brandon

Does anyone know how to incorporate a rootfinding algorithm with a numerical dsolve? What I mean is to solve an ODE numerically, but having to use a rootfinding technique at the same time... Thanks for your input guys, Brado

Hi guys, I just have a quick question: Does anyone know if there is a way to access to actual code Maple uses for solving ode's numerically? Like, for instance, if you use the Runge-Kutta method and you use dsolve(numeric), is there anyone to actually see the code that Maple uses to calculate the ODE? Because I need to add something inside the loop used to solve the ODE. If anyone knows an easier way to add something to the method that Maple uses to solve ODE's numerically, I'd love to know!!! Thanks guys, Brandon