best software package to numerically solve high order nonlinear ODE boundary value problem












1












$begingroup$


I need to numerically solve a 4th and 5th order nonlinear ODE BVP and I was hoping I could get some advice on the best software package to solve these types of problems. I've used MATLABS bvp4c for the 4th order equation and I'm not too confident in the solutions, I am hoping to try another method.
Any advice will be helpful! Thanks!










share|cite|improve this question











$endgroup$








  • 1




    $begingroup$
    I don't know about "best", but you could try Maple, which has a variety of solution methods. Is your problem particularly difficult? For example, is it stiff, or singular, or does it have narrow boundary layers? Perhaps you could give us an example of the type of problem you're looking at.
    $endgroup$
    – Robert Israel
    Dec 19 '18 at 18:17










  • $begingroup$
    Could you also quantify or at least qualify how you determine your non-confidence in the BVP solution? Does it change materially when you change the tolerances? Is the ode45 forward solution far away from the bvp4c solution?
    $endgroup$
    – LutzL
    Dec 19 '18 at 18:44










  • $begingroup$
    What do your mean by "order"? The highest order derivative, the highest non-linear term, or higher order methods like RK5(4) and beyond? That would help.,
    $endgroup$
    – ggcg
    Dec 19 '18 at 21:00










  • $begingroup$
    Also, are you asking about s/w packages or "methods"? This is a big deal as many "packages" will have all the same methods so what's really the point? MAPLE and Mathematica have pretty good ODE solvers and there is the open source ODEINT. If you are not confident in your answer I would recommend posting some code and letting us check it. You can't make a finite solver of a continuous "work" unless you check outputs and use step size control and other features. MATLAB requires you to check, MAPLE for example uses 1000's of intermediate steps for each step to converge.
    $endgroup$
    – ggcg
    Dec 19 '18 at 21:03










  • $begingroup$
    Are you aware of the existence of "stiff solvers" in Matlab like "ode15s" ?
    $endgroup$
    – Jean Marie
    Dec 19 '18 at 21:13


















1












$begingroup$


I need to numerically solve a 4th and 5th order nonlinear ODE BVP and I was hoping I could get some advice on the best software package to solve these types of problems. I've used MATLABS bvp4c for the 4th order equation and I'm not too confident in the solutions, I am hoping to try another method.
Any advice will be helpful! Thanks!










share|cite|improve this question











$endgroup$








  • 1




    $begingroup$
    I don't know about "best", but you could try Maple, which has a variety of solution methods. Is your problem particularly difficult? For example, is it stiff, or singular, or does it have narrow boundary layers? Perhaps you could give us an example of the type of problem you're looking at.
    $endgroup$
    – Robert Israel
    Dec 19 '18 at 18:17










  • $begingroup$
    Could you also quantify or at least qualify how you determine your non-confidence in the BVP solution? Does it change materially when you change the tolerances? Is the ode45 forward solution far away from the bvp4c solution?
    $endgroup$
    – LutzL
    Dec 19 '18 at 18:44










  • $begingroup$
    What do your mean by "order"? The highest order derivative, the highest non-linear term, or higher order methods like RK5(4) and beyond? That would help.,
    $endgroup$
    – ggcg
    Dec 19 '18 at 21:00










  • $begingroup$
    Also, are you asking about s/w packages or "methods"? This is a big deal as many "packages" will have all the same methods so what's really the point? MAPLE and Mathematica have pretty good ODE solvers and there is the open source ODEINT. If you are not confident in your answer I would recommend posting some code and letting us check it. You can't make a finite solver of a continuous "work" unless you check outputs and use step size control and other features. MATLAB requires you to check, MAPLE for example uses 1000's of intermediate steps for each step to converge.
    $endgroup$
    – ggcg
    Dec 19 '18 at 21:03










  • $begingroup$
    Are you aware of the existence of "stiff solvers" in Matlab like "ode15s" ?
    $endgroup$
    – Jean Marie
    Dec 19 '18 at 21:13
















1












1








1





$begingroup$


I need to numerically solve a 4th and 5th order nonlinear ODE BVP and I was hoping I could get some advice on the best software package to solve these types of problems. I've used MATLABS bvp4c for the 4th order equation and I'm not too confident in the solutions, I am hoping to try another method.
Any advice will be helpful! Thanks!










share|cite|improve this question











$endgroup$




I need to numerically solve a 4th and 5th order nonlinear ODE BVP and I was hoping I could get some advice on the best software package to solve these types of problems. I've used MATLABS bvp4c for the 4th order equation and I'm not too confident in the solutions, I am hoping to try another method.
Any advice will be helpful! Thanks!







ordinary-differential-equations numerical-methods boundary-value-problem math-software






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Dec 20 '18 at 5:25







tenichols

















asked Dec 19 '18 at 17:47









tenicholstenichols

83




83








  • 1




    $begingroup$
    I don't know about "best", but you could try Maple, which has a variety of solution methods. Is your problem particularly difficult? For example, is it stiff, or singular, or does it have narrow boundary layers? Perhaps you could give us an example of the type of problem you're looking at.
    $endgroup$
    – Robert Israel
    Dec 19 '18 at 18:17










  • $begingroup$
    Could you also quantify or at least qualify how you determine your non-confidence in the BVP solution? Does it change materially when you change the tolerances? Is the ode45 forward solution far away from the bvp4c solution?
    $endgroup$
    – LutzL
    Dec 19 '18 at 18:44










  • $begingroup$
    What do your mean by "order"? The highest order derivative, the highest non-linear term, or higher order methods like RK5(4) and beyond? That would help.,
    $endgroup$
    – ggcg
    Dec 19 '18 at 21:00










  • $begingroup$
    Also, are you asking about s/w packages or "methods"? This is a big deal as many "packages" will have all the same methods so what's really the point? MAPLE and Mathematica have pretty good ODE solvers and there is the open source ODEINT. If you are not confident in your answer I would recommend posting some code and letting us check it. You can't make a finite solver of a continuous "work" unless you check outputs and use step size control and other features. MATLAB requires you to check, MAPLE for example uses 1000's of intermediate steps for each step to converge.
    $endgroup$
    – ggcg
    Dec 19 '18 at 21:03










  • $begingroup$
    Are you aware of the existence of "stiff solvers" in Matlab like "ode15s" ?
    $endgroup$
    – Jean Marie
    Dec 19 '18 at 21:13
















  • 1




    $begingroup$
    I don't know about "best", but you could try Maple, which has a variety of solution methods. Is your problem particularly difficult? For example, is it stiff, or singular, or does it have narrow boundary layers? Perhaps you could give us an example of the type of problem you're looking at.
    $endgroup$
    – Robert Israel
    Dec 19 '18 at 18:17










  • $begingroup$
    Could you also quantify or at least qualify how you determine your non-confidence in the BVP solution? Does it change materially when you change the tolerances? Is the ode45 forward solution far away from the bvp4c solution?
    $endgroup$
    – LutzL
    Dec 19 '18 at 18:44










  • $begingroup$
    What do your mean by "order"? The highest order derivative, the highest non-linear term, or higher order methods like RK5(4) and beyond? That would help.,
    $endgroup$
    – ggcg
    Dec 19 '18 at 21:00










  • $begingroup$
    Also, are you asking about s/w packages or "methods"? This is a big deal as many "packages" will have all the same methods so what's really the point? MAPLE and Mathematica have pretty good ODE solvers and there is the open source ODEINT. If you are not confident in your answer I would recommend posting some code and letting us check it. You can't make a finite solver of a continuous "work" unless you check outputs and use step size control and other features. MATLAB requires you to check, MAPLE for example uses 1000's of intermediate steps for each step to converge.
    $endgroup$
    – ggcg
    Dec 19 '18 at 21:03










  • $begingroup$
    Are you aware of the existence of "stiff solvers" in Matlab like "ode15s" ?
    $endgroup$
    – Jean Marie
    Dec 19 '18 at 21:13










1




1




$begingroup$
I don't know about "best", but you could try Maple, which has a variety of solution methods. Is your problem particularly difficult? For example, is it stiff, or singular, or does it have narrow boundary layers? Perhaps you could give us an example of the type of problem you're looking at.
$endgroup$
– Robert Israel
Dec 19 '18 at 18:17




$begingroup$
I don't know about "best", but you could try Maple, which has a variety of solution methods. Is your problem particularly difficult? For example, is it stiff, or singular, or does it have narrow boundary layers? Perhaps you could give us an example of the type of problem you're looking at.
$endgroup$
– Robert Israel
Dec 19 '18 at 18:17












$begingroup$
Could you also quantify or at least qualify how you determine your non-confidence in the BVP solution? Does it change materially when you change the tolerances? Is the ode45 forward solution far away from the bvp4c solution?
$endgroup$
– LutzL
Dec 19 '18 at 18:44




$begingroup$
Could you also quantify or at least qualify how you determine your non-confidence in the BVP solution? Does it change materially when you change the tolerances? Is the ode45 forward solution far away from the bvp4c solution?
$endgroup$
– LutzL
Dec 19 '18 at 18:44












$begingroup$
What do your mean by "order"? The highest order derivative, the highest non-linear term, or higher order methods like RK5(4) and beyond? That would help.,
$endgroup$
– ggcg
Dec 19 '18 at 21:00




$begingroup$
What do your mean by "order"? The highest order derivative, the highest non-linear term, or higher order methods like RK5(4) and beyond? That would help.,
$endgroup$
– ggcg
Dec 19 '18 at 21:00












$begingroup$
Also, are you asking about s/w packages or "methods"? This is a big deal as many "packages" will have all the same methods so what's really the point? MAPLE and Mathematica have pretty good ODE solvers and there is the open source ODEINT. If you are not confident in your answer I would recommend posting some code and letting us check it. You can't make a finite solver of a continuous "work" unless you check outputs and use step size control and other features. MATLAB requires you to check, MAPLE for example uses 1000's of intermediate steps for each step to converge.
$endgroup$
– ggcg
Dec 19 '18 at 21:03




$begingroup$
Also, are you asking about s/w packages or "methods"? This is a big deal as many "packages" will have all the same methods so what's really the point? MAPLE and Mathematica have pretty good ODE solvers and there is the open source ODEINT. If you are not confident in your answer I would recommend posting some code and letting us check it. You can't make a finite solver of a continuous "work" unless you check outputs and use step size control and other features. MATLAB requires you to check, MAPLE for example uses 1000's of intermediate steps for each step to converge.
$endgroup$
– ggcg
Dec 19 '18 at 21:03












$begingroup$
Are you aware of the existence of "stiff solvers" in Matlab like "ode15s" ?
$endgroup$
– Jean Marie
Dec 19 '18 at 21:13






$begingroup$
Are you aware of the existence of "stiff solvers" in Matlab like "ode15s" ?
$endgroup$
– Jean Marie
Dec 19 '18 at 21:13












0






active

oldest

votes











Your Answer





StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});

function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});














draft saved

draft discarded


















StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3046671%2fbest-software-package-to-numerically-solve-high-order-nonlinear-ode-boundary-val%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown

























0






active

oldest

votes








0






active

oldest

votes









active

oldest

votes






active

oldest

votes
















draft saved

draft discarded




















































Thanks for contributing an answer to Mathematics Stack Exchange!


  • Please be sure to answer the question. Provide details and share your research!

But avoid



  • Asking for help, clarification, or responding to other answers.

  • Making statements based on opinion; back them up with references or personal experience.


Use MathJax to format equations. MathJax reference.


To learn more, see our tips on writing great answers.




draft saved


draft discarded














StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3046671%2fbest-software-package-to-numerically-solve-high-order-nonlinear-ode-boundary-val%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

Quarter-circle Tiles

build a pushdown automaton that recognizes the reverse language of a given pushdown automaton?

Mont Emei