Those are some absolutely huge numbers, so it's not surprising that NDSolve can't deal with this over a time range of t=0 to t=17. New in Mathematica 9 Advanced Hybrid and Differential Algebraic Equations Poincaré Sections Compute Poincaré sections using NDSolve with WhenEvent. \c = (\ L^3)/(8 c \^2) (NN (\e11 - \111) V)^2}, Wolfram Community forum discussion about Problem with NDSolve in Mathematica 9.
If you precede your code with alpha11 alpha22 alpha32 alpha42 alpha51, you. Also, NDSolve is strictly a numerical solver, so it doesnt understand or accept unknown symbolic values like alpha1, alpha2, et cetera. Name your constants alpha1 or 1 instead of alpha1. Widely admired for both its technical prowess and elegant ease of. For three decades, Mathematica has defined the state of the art in technical computingand provided the principal computation environment for millions of innovators, educators, students and others around the world.
Na' = -Ip \pa na + (\ea1 c)/\ (Na - na) q + Do not use slashes and underlines in variable names, it does not work. The worlds definitive system for modern technical computing. Stay on top of important topics and build connections by joining Wolfram Community groups relevant to your interests.
The initial conditions for this system are n(0) = NN, na(0) = Na, n'(0) = q(0) = 0, q'(0) = some_constants, na'(0) = Ip(0)×Some_other_constants (shown in my code below). Wolfram Community forum discussion about Extracting data from plot of NDSolve output. I'm trying to solve this system of differential equation (the ones in yellow markers are just constants): 6.7 Some Remarks on NDSolve Since we intend to generally use NDSolve to solve nonlinear ODEs of interest from now on, a few concluding remarks should be.