# File diffeq.rb, line 575 def compute_order(log10norm, log10abserr, log10relerr) #* # * this is to determine the 'optimal' degree of the taylor expansion. # * # * parameters: # * log10norm: base 10 log of the norm of the initial condition # * log10abserr: base-10 log of the absolute error required # * log10relerr: base-10 log of the relative error required # * # * returned value: 'optimal' degree. #*/ # When I studied it I couldn't really understand Angel Jorba && Maorong Zou 's # method for determining the order. I don't adjust it anymore. # It will be set to the problem file max_terms which defaults to 30 which seems # a good value most of the time. # I've retained $max_taylor in case I find a better way. return($max_taylor) end