For a short period of time the missing right hand margin/ word wrap problem was solve but after some "upgrade" weeks ago it reapeared and it is still there for me: Win XP Pro SP1, Firefox 1.01, resolution 1024x768 Regards, Alejandro

I have just updated M10 in a Win XP system at the university, behind a firewall. The first time I did in main menu Tools > Check for Updates, and Standard apparently hanged. Then I killed the maple process, and on a new execution it opened a dialog box offering update. I did, it asked me http proxy and port and then it downloaded/updated smoothly. Apparently it went fine as I get kernelopts(version); Maple 10.01, IBM INTEL NT, Aug 18 2005 Build ID 196612

I have written a new numeric version for chi^2 using piecewise option (chin4, in the new attached worksheet on top). But I need to find yet the trick to use this function once instead of N times, as timings show.

I have written a new numeric version for chi^2 using piecewise option (chin4, in the new attached worksheet on top). But I need to find yet the trick to use this function once instead of N times, as timings show.

Yes, the purpose of this calculation is to minimize this chi^2 function with GOT, as in the worksheet that I have attached in the thread generate code for numerical integral So, the idea is that GOT uses an approximation for chi^2, and its derivatives wrt the parameters, using the interpolating function instead of the integral. Derivatives may be used in the local phase of the search of the minimum, as far as I understand. Hence I guess that some properties of the approximation to the derivatives are also important, eg continuity. For a given dataset, the range of z is fixed by (z_min, z_max), z=2 was just for testing. So, if I understand correctly, I should use the output=Array option, but not the global variable. What do you suggest instead?

Yes, the purpose of this calculation is to minimize this chi^2 function with GOT, as in the worksheet that I have attached in the thread generate code for numerical integral So, the idea is that GOT uses an approximation for chi^2, and its derivatives wrt the parameters, using the interpolating function instead of the integral. Derivatives may be used in the local phase of the search of the minimum, as far as I understand. Hence I guess that some properties of the approximation to the derivatives are also important, eg continuity. For a given dataset, the range of z is fixed by (z_min, z_max), z=2 was just for testing. So, if I understand correctly, I should use the output=Array option, but not the global variable. What do you suggest instead?

Yes, the task is to find the best fit for the parameters beta, lambda, sigma, P, B and M given a data set (in this case test of 54 points) given as lists (zl,ml,sl) for numbers z_i, m_i and errors s_i so that chi^2 is a sum of terms of the form (m_i-m(z_i))^2/s_i^2 and m(z) involves the function DL(z), that is defined in terms of an integral. As said I am using this old 54 data set for testing numerical methods, the actual one being over 150 (and much larger in the future). Yes each run is taking a long time. The function v is used in this syntax of GOT for a constraint among the parameters. For background on the original problem you may wish to give a look at these papers, dealing with simpler models: http://arxiv.org/PS_cache/astro-ph/pdf/0101/0101549.pdf (pp 5-7 and figures at the end) and http://arxiv.org/PS_cache/astro-ph/pdf/0303/0303145.pdf (pp 15-17)

Yes, the task is to find the best fit for the parameters beta, lambda, sigma, P, B and M given a data set (in this case test of 54 points) given as lists (zl,ml,sl) for numbers z_i, m_i and errors s_i so that chi^2 is a sum of terms of the form (m_i-m(z_i))^2/s_i^2 and m(z) involves the function DL(z), that is defined in terms of an integral. As said I am using this old 54 data set for testing numerical methods, the actual one being over 150 (and much larger in the future). Yes each run is taking a long time. The function v is used in this syntax of GOT for a constraint among the parameters. For background on the original problem you may wish to give a look at these papers, dealing with simpler models: http://arxiv.org/PS_cache/astro-ph/pdf/0101/0101549.pdf (pp 5-7 and figures at the end) and http://arxiv.org/PS_cache/astro-ph/pdf/0303/0303145.pdf (pp 15-17)

I am not familiar with Ooura's double exponential integration. Is it better than the NAG method you have suggested before? I have just attached the M9.5 worksheet that generates the stack limit error to the initial message of this thread (it seems that this is the only place where it is possible to attach files). An initial run in an old PC with the addition of method = _d01ajc seemed promising: it did not issued stack error message when it used to. As said before, I am evaluating options whether I can make these calculations from within Maple or eventually export pieces of code and run from outside. And besides fitting I will need to make confidence intervals calculations for the parameters, confidence contour plots, etc. I am not aware of any package for this stuff within Maple, and sounds that I should look for that elsewhere. Any suggestion is also welcome.

I am not familiar with Ooura's double exponential integration. Is it better than the NAG method you have suggested before? I have just attached the M9.5 worksheet that generates the stack limit error to the initial message of this thread (it seems that this is the only place where it is possible to attach files). An initial run in an old PC with the addition of method = _d01ajc seemed promising: it did not issued stack error message when it used to. As said before, I am evaluating options whether I can make these calculations from within Maple or eventually export pieces of code and run from outside. And besides fitting I will need to make confidence intervals calculations for the parameters, confidence contour plots, etc. I am not aware of any package for this stuff within Maple, and sounds that I should look for that elsewhere. Any suggestion is also welcome.

This error message arised using GOT to minimize a chisq object function involving expressions with evalf(Int(...)). I do not know whether NAG or other external libraries are called for the calculation of the integral in this context. The stack limit error appears after some time of calculation (the period varies with the run, apparently some randomness is involved).

This error message arised using GOT to minimize a chisq object function involving expressions with evalf(Int(...)). I do not know whether NAG or other external libraries are called for the calculation of the integral in this context. The stack limit error appears after some time of calculation (the period varies with the run, apparently some randomness is involved).

I find a stack limit problem while using GOT, hence I am exploring options. I would prefer to stay within Maple if possible. I will try again with GOT specifying the NAG method, perhaps it helps to avoid this stack problem. It seems that the integral cannot be expressed in terms of known functions for B not = 0. As a function of z it satisfies a fourth order nonlinear ODE (for B=0 it can be expressed in terms of a hypergeometric function).

I find a stack limit problem while using GOT, hence I am exploring options. I would prefer to stay within Maple if possible. I will try again with GOT specifying the NAG method, perhaps it helps to avoid this stack problem. It seems that the integral cannot be expressed in terms of known functions for B not = 0. As a function of z it satisfies a fourth order nonlinear ODE (for B=0 it can be expressed in terms of a hypergeometric function).

Values of z are in the range 0 up to 1.5 or 2 currrently. For typical values of the parameters I find that a Taylor series departures significantly above z=0.5. An alternative expansion, with uniform convergence in this range would be needed. I will explore that. Any suggestion is welcome.