Questions on the symbolic (DSolve, DifferentialRoot) and numerical (NDSolve) solutions of differential equations in Mathematica.

learn more… | top users | synonyms (3)

2
votes
0answers
438 views

NDSolve Issue with initial conditions

I am trying to implement a Poincare section for a gravitational movement on the plane $(x,y)$. Here is the code I wrote ...
2
votes
1answer
1k views

Nonlinear least squares

I have a typical nonlinear least square problem and it seems very similar to the example provided in the section on parameter fitting at the Wolfram website. Unfortunately I always get some error ...
4
votes
2answers
793 views

FindFit for a differential equation model with multiple variables?

I'm trying to figure out how to use FindFit with a multivariable differential equation model and data. I've successfully made it work for the one-variable version ...
3
votes
1answer
301 views

NDSolve much slower when using InterpolatingFunctions inside

I am solving an ODE system with Mathematica. Let's say the system is of the form (to simplify). sol = NDSolve[{x'[t] == f[t] x[t] + g[t], x[0] == 0}, x, {t, 0, 1}] ...
3
votes
2answers
295 views

DSolve vs a system of differential equations

I'm trying to solve a simple system of differential equations. ...
-1
votes
1answer
200 views

Linking two differential equations

The solution and plot of the Differential equation is $\frac{d^2S(r)}{dr^2}+\frac{1}{r}\frac{dS(r)}{dr}-S(r)+3S^3(r)=0$ ...
-1
votes
1answer
894 views

plotting the solutions of system of differential equations [duplicate]

The following code include the system of differential equations: ...
1
vote
1answer
306 views

Dsolve too slow — is there anyway around?

I am trying to solve: ...
10
votes
1answer
502 views

Modelling the effect of a structure on a “tsunami” (hyperbolic wave equation)

So, the hyperbolic wave equation can be quite easily solved in Mathematica like this: ...
2
votes
2answers
221 views

How to pick out specific parts of the solution of a differential equation efficiently?

I'm trying to solve a set of differential equations (in parallel) and then select specific parts of the solution. What I have at the moment (see below-I've simplified the differential equations, in ...
2
votes
1answer
83 views

Incorrect information displayed when slider is moved (intermittent)

All, the following code graphs 10 different solutions for a first order ODE and allows the user to modify certain equation parameters to see how the graph is affected: ...
4
votes
2answers
810 views

DSOLVE for tenth order boundary value problem

I am currently using Mathematica 9 to find solutions for higher order ordinary differential equations. The following was my input for a tenth order boundary value problem (BVP): ...
2
votes
3answers
954 views

How to generate 2D FEM mesh in Mathematica?

To generate 3D mesh TetGen can be easily used. Are there similar functions (or a way to use TetGen) to generate 2d mesh? I know that such functionality can be easily implemented but I would like to ...
16
votes
3answers
887 views

Only final result from NDSolve

Finally, I started to play with differential equations in Mathematica. And I have faced the problem, which seems to me so basic that I'm afraid this question is going to be closed soon. However, I'...
0
votes
0answers
620 views

Plot basins of attraction using separatrices for particular pendulum [duplicate]

How can I plot the stable and unstable separatrices to make a graphic of the basins of attraction of this differential equation? ...
3
votes
1answer
565 views

Use output of NDSolve as an integration limit of NIntegrate

I've been working on a ODE, for a function y[t], that, due to it's complexity, has to be solved numerically. I then need to use the results of NDSolve, which I ...
3
votes
1answer
403 views

Transform recursion for coefficients into differential equation for generating function

Assume, one is given a linear recursion with polynomial coefficients for a sequence $(a_i)_i$, such as a[i] == i a[i-1] I would like to convert this recursion ...
7
votes
1answer
844 views

Modelling a Rocket Launch using NDSolve

I'm trying to model a rocket launch with Mathematica but I've run into a little problem since I don't know how to turn the thrust off. I'm using Newton's Law of Universal Gravitation plus an added ...
4
votes
1answer
295 views

Create ParametricFunction object without ParametricNDSolve

Documentation of v9 states that ParametricFunction is generated by ParametricNDSolve and ...
0
votes
1answer
958 views

Bifurcation diagrams for system of equations

I am interested in constructing a bifurcation diagram for parameter a in the dynamical system given in the code below. I want to see how parameter changes affect the stability of the system. The ...
0
votes
1answer
1k views

Can graph interpolating function, but not use it. Any thoughts?

When solving a system of differential equations numerically, I get two interpolating functions. I can graph both of these, but when I try to evaluate either of the functions at a point, it just spits ...
0
votes
2answers
242 views

Interpolating a ParametricFunction

Original domain of my independent variable t is 0 to 10 nanometers, but in Mathematica (as in most numerical software), it's good idea to solve differential ...
1
vote
0answers
590 views

NDSolve for a large system of coupled (complex valued) ODEs

The system under consideration here is a generalization of the one I discussed in this problem. Qualitatively, I am trying to solve a system of coupled nonlinear (up to order 4 in the dep variables) ...
6
votes
1answer
437 views

Changing a wave simulation

I want to adjust some code that I found on this Wolfram Demonstration: http://demonstrations.wolfram.com/WaterRipplesPassingThroughASlit/ Basically it shows a wave passing through a single slit. I ...
3
votes
2answers
671 views

Ploting a 3D vector field using Manipulate, VectorPlot3D and NDSolve

My aim is to plot a vector field of the following system with a few trajectories: $$r'(t)=i-l.r(t)-\text{ux}. r(t). x(t)-\text{uy}. r(t). y(t) \\ x'(t)=\text{ex}. \text{ux}. r(t). x(t)-\text{mx}. x(...
6
votes
2answers
582 views

NDSolve: ProcessEquations and Reinitialize with Piecewise functions

I am having trouble with using NDSolve`Reinitialize when the system consists of a piecewise function. If we define the ODE system ...
1
vote
1answer
2k views

3D Phase Portrait with direction fields and sample solution curves [duplicate]

I have seen NDSolve and ParamtericPlot3D used for things like the Lorenz Attractor, but was wondering if there is a way to draw a 3D phase portrait for a system. Something akin to this: http://www....
1
vote
0answers
199 views

Having Mathematica solve or help solve a set of 8 equations

There is a set of DDE's that arise in leukemia dynamics and I'd like to find the equilibrium points of the system by setting all of the derivatives equal to zero. The system can be described in ...
1
vote
1answer
1k views

Problem in solving PDE

I have a problem in solving a type of PDE. The actual problem involves diffusion and adsorption in a trench, but here is a simpler example. I have two independent variables $t$ and $x$. There are two ...
3
votes
1answer
615 views

Limitations of ParametricNDSolve family w.r.t objective functions

Observation: I can see even for very simple modification in case of an scalar objective involving an definite integral in time ParametricNDSolve fails. Here is an ...
4
votes
1answer
298 views

How to demonstrate lack of stability with advection equation

I am trying to that using a coarse grid with an explicit method for, say, the advection equation leads to an unstable solution. The trouble is Mathematica avoids unstable solutions with good ...
3
votes
2answers
358 views

How do I find the best fit?

I am trying to solve following system of three differential equations ...
2
votes
1answer
639 views

Crank-Nicolson with NDSolve?

As far as I understand, the Crank-Nicolson method (a.k.a. trapezoidal method) can be expressed as a second order implicit Runge-Kutta method. It's Butcher tableau is: ...
0
votes
2answers
802 views

Manipulating differential equations

I would be very grateful if anyone could help with manipulating the solution to the system of differential equations below. I've tried just about everything I can think of all to no avail. The code ...
3
votes
1answer
252 views

Using WhenEvent for derivative of discontinuous function

I have a discontinuous function ($u(t)$, a square wave) and I would like WhenEvent to trigger when the signal goes high/low, i.e. when the value of $u(t)$ changes. ...
3
votes
1answer
283 views

Utterly confounding behaviour in Mathematica 9 DSolve

I've been baffled by a lot of the behaviours of Mathematica 9 — so much that I regret upgrading my copy of Mathematica 8 (I would downgrade but I'm currently on vacation). Here is an example: ...
1
vote
1answer
123 views

Incorporating a Manipulation in an Animation

I am having trouble incorporating a manipulation into my animation. I want to be able to maniuplate the 5 constants at the beginning. Am I going to need to combine all of the code into one string? ...
1
vote
1answer
237 views

Using new symbol in WhenEvent in NDSolve

This is an example in Mathematica's help document. ...
1
vote
1answer
183 views

How can I eliminate $x_1$ from my two ODEs?

How can I get Mathematica to eliminate $x_1$ in the equations in the first box to achieve the equation in the second box? My code: ...
2
votes
1answer
2k views

DSolve doesn't find a solution

I have this partial differential equation (G and M are constants): and a solution is: Mathematica doesn't find any solutions. My code for DSolve is: ...
9
votes
2answers
430 views

Adding a constant vector to a vector differential equation seems to break NDSolve. Why?

I'm trying to solve a differential equation that's phrased in terms of matrices and vectors. My minimum working example is this: ...
1
vote
1answer
315 views
1
vote
2answers
784 views

How to define a function based on the output of DSolve? [duplicate]

I would like to define a function based on the output of the following DSolve call: ...
-3
votes
1answer
257 views

NDSolve doesn't solve my system of ODEs

What is causing the error in this code? ...
0
votes
2answers
523 views

Autonomous equations' critical points and eigenvalues

I have three-dimensional autonomous equations and want to find their critical points and corresponding eigenvalues. They have 6 critical points each. x, ...
1
vote
1answer
264 views

Changing ParametricNDSolveValue options for only part of the run

Here is an example of the system that I am running, ...
4
votes
1answer
492 views

How can I deal with a non-numerical value for a derivative at $t = 0$ when using NDSolve?

I want to solve two coupled equations with NDSolve, ...
20
votes
2answers
4k views

How to fit 3 data sets to a model of 4 differential equations?

I'm a biologist and a newbie in Mathematica. I want to fit three data sets to a model consisting of four differential equations and 10 parameters. I want to find the parameters best fitting to my ...
9
votes
1answer
671 views

Complex valued 2+1D nonlinear PDE using NDSolve

I am trying to follow the main ideas presented in this question, applying it to my own problem, which is a complex, time-dependent, nonlinear PDE: $$i \frac{\partial \psi}{\partial t} = \left[ -\...
12
votes
1answer
388 views

solve ODE with divergencies

The solutions of the second order differential equation $$\frac{1}{\eta}\frac{d}{d\eta}\left(\eta \frac{df}{d\eta}\right)+\left(1-\frac{s^2}{\eta^2}\right)f-f^3=0$$ is shown in Fig. 5.2 below, for ...