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

learn more… | top users | synonyms (3)

5
votes
1answer
414 views

How to solve this Integral equation

D[x[t] - x[t - 1]/(2 E), {t, 3}] + Integrate[E^(-δ)*x[t - δ]/5^t, {δ, 2, 2.5}] == 0 I found solve this problem is hard with Mathematica. I also find a article ...
0
votes
1answer
386 views

How to put constraints on NDSolve[] solution

Assumed the following example of double pendulum by Thison (Answered by Jens): Animation of double pendulum Imagine that in the given animation(i.e. http://i.stack.imgur.com/p4TgL.gif ), the second ...
2
votes
3answers
386 views

Mathematica won't plot the a function that is solved from NDSolve

I have used NDSolve to solve for a function (an angle of a triple pendulum), and now I wish to plot the derivative of that function (the angular velocity). ...
0
votes
1answer
106 views

Misunderstanding of Function Notation

I have: sol = DSolve[{y'[t] == y[t]^2 - y[t], y[0] == c}, y[t], t]; y[t_, c_] := y[t] /. First[sol]; Then: y[t,c] Produces ...
0
votes
1answer
83 views

DSolve for first order nonlinear equation is solution numerical or analytical?

Hi I'm new to mathematica so I may have a somewhat trivial question. my code below solves and plots the results of a first order nonlinear ODE. The solutions mathematica gives me contains & and #1 ...
5
votes
3answers
602 views

Numerically Simulating a Simple Differential Equation: RC Circuit

I need to learn to simulate a differential equation in Mathematica. The differential equation that I need to simulate is complicated, without an analytical solution. However, I am brand new to ...
3
votes
1answer
243 views

Breaking out of NDSolve

I am solving a coupled set of differential equations with NDSolve for 6 unknown functions of time. At a certain point in time, the system hits a singular point ...
2
votes
2answers
549 views

NDSolve solution does not satisfy boundary conditions

I am trying to solve a 1st order non-linear ODE ...
1
vote
0answers
242 views

Plotting phase portrait from NDSolve (i.e., must plot x'[x], where $x=x(t)$)

I'm doing a double pendulum system and I wish to plot the phase portraits for each angle. I've successfully used NDSolve to solve for each of the angles. Let's just ...
1
vote
0answers
143 views

How to solve Fuzzy Differential Equations

How to solve Fuzzy Differential Equations using DSolve or NDSolve, Or Differential Equations with fuzzy number as initial conditions
1
vote
1answer
221 views

Plotting ODE solutions with two parameters

I have a system of two first-order ODEs with two parameters, a and b and I'm trying to plot the solution curve in the X-Y plane at a given initial condition for various parameter values. Can someone ...
1
vote
0answers
197 views

Why is my output equal to my input?

I tried to solve a system of differential equations with the code below: ...
3
votes
2answers
208 views

Turning arrowheads off

My thinking was that Arrowheads[0] would turn off the arrowheads. ...
4
votes
1answer
140 views

Can anyone help me to get the solution of this first order ODE?

$\ y'(x)=-\sqrt(y(x)^3)$ with initial condition $\ y(0)=-2$. In fact, this is an example here in page 3-4 ...
1
vote
0answers
139 views

Direction Field for an Ordinary Differential Equation [duplicate]

I have found links for plot the direction field of a differential equation. For example: How can I plot the direction field for a differential equation? And: Plotting the direction field of a ...
8
votes
1answer
418 views

How can I improve speed of code that generates an orbit plot?

I am attempting to generate an orbit for the "time 2 Pi map" of a forced damped pendulum. The example comes from page 57 in the Chaos text by Alligood,Sauer and Yorke, see Figure 2.7. The code ...
2
votes
0answers
52 views

The use of Evaluate in NDSolve and Plotting the answer [duplicate]

I am having a bit of hard time understanding the use of Evaluate when plotting the output of NDSolve. From the documentation, I enter: ...
4
votes
1answer
160 views

System of differential equations on a directed graph

I want to use Mathematica to model blood flow through a capillary network. The approach we are using is going to be to create a directed graph, and have a system of differential equations at each ...
2
votes
0answers
676 views

NDSolve and strange “nonlinear coefficients problem”

I'm stuck solving the following problem. I defined two functions as follows: $$ \varphi(\lambda) = \frac{\left((\lambda-2)^2-1 \right)^2}{4}$$ $$ \gamma(\lambda) = \varepsilon^2 \quad\text{where}\...
2
votes
2answers
210 views

Can't use my own discontinuous coefficient within NDSolve

I want to evaluate the pendulum's motion with changing equilibrium state given by eqilValue[t]: ...
1
vote
2answers
486 views

DSolve bvlim error due to Boundary condition at infinity

I am solving an Ode using Mathematica that has boundary condition at infinity which is why Mathematica is showing error message. Please help to tackle this problem. ...
4
votes
1answer
337 views

Non-separable partial differential equation in polar coordinates

I'm trying to solve the Schroedinger equation in 2D for a system interacting via a dipole potential. This means, in effect, I'm trying to solve the nonlinear PDE $$ -\frac{1}{r} \frac{\partial}{\...
3
votes
1answer
172 views

Troubles with FindFit for ODE

I need to use FindFit for ODE with arbitary number of parameters. This code (from documentation) works: ...
1
vote
2answers
71 views

Strange behavior of the Evaluate function

I am solving a system of differential equation whose unknown are saved under variable vec={C00,C01,C10,C11} and the solutions are saved under variable ...
1
vote
1answer
188 views

using Module in NDSolve

I am trying to solve a set of two wave equations where solution of first equation defines initial conditions for the second one. And I would like to use Module to define those equations as a function ...
1
vote
0answers
301 views

Help deciphering Neumann Value

Version 10 has an option for NDSolve called "NeumannValue" which is supposed to allow you to specify boundary conditions of the form $\vec{n} . (c \nabla u + \alpha u -\gamma) = g-q u$ but there is no ...
0
votes
1answer
299 views

Help solving the following non-linear differential equation

After spending several hours exploring different options, I still couldn't not figure out how to solve the following differential equations. Any suggestions would be highly appreciated. ...
3
votes
1answer
215 views

Hermite method in Mathematica

When I solve the example below, in the results it is mentioned that Mathematica 10 has used the Hermite method. I cannot find anything about this in the Mathematica documentation. Is there nothing ...
1
vote
1answer
120 views

Problem solving a differential equation

I want to solve the following differential equation: y''[x]+a*y[x]=-b*u[x], y[0]==0, y'[0]==0 for y(x). u[x] is a vector, ...
3
votes
2answers
388 views

ParametricPlot3D line not smooth

Multivariate parameters appear to present a jagged appearance of integrands (using default Runge-Kutta ODE integration intervals?) in ParametricPlot3D plotting on a single argument. Higher Mesh.. ...
3
votes
1answer
212 views

Evaluate function defined by DifferentialRoot

I have the following sequence of rationals that I want to find the generating function of: ...
2
votes
2answers
273 views

ParametricNDSolve and plotting problem

I have a problem with ParametricNDSolve: ...
4
votes
1answer
2k views

Combining Gravity Turn and Orbit Models

I have a mathematical model for the motion of an orbiting spacecraft about Earth: ...
0
votes
2answers
107 views

Derivative in function form

I have a function in mathematica defined as b[n_, x_] := x^n/n!/\!\(\*UnderoverscriptBox[\(\[Sum]\), \(k = n\), \(\[Infinity]\)]\*FractionBox[\(x^k\), \(k!\)]\) ...
1
vote
0answers
89 views

Updating procedural fitting algorithm to more efficient style?

For some time now, I have used this procedural programming approach to fit data in my research that is too cumbersome for the built in Mathematica functions. (Although this simple example works fine ...
0
votes
0answers
30 views

Identification of a parameter in a first order differential equation [duplicate]

I wan't to identify a parameter in a first oder differential equation. My goal is proof of the calculation principle i use. First I calculate a response where the paramater is set to 0.3. Later this ...
1
vote
1answer
124 views

NDSolve with WhenEvent, resetting system for a prolonged period

I am currently working on a complex system where I would like to (for the lack of a better description) reset (part) of the system for a certain period. Say I five ODE's: ...
1
vote
0answers
129 views

A Stupid Question Maybe - EigenNDSolve

I would like to find out how accurate and how useful the chebyshev (collocation?) method is in finding many eigenvalues to a second-order ODE in one go. Specifically, I used the Mathematica package "...
5
votes
2answers
183 views

How to add (energy) constraint when using NDSolve to Equation of Motion

To simplify my problem, I will try and solve the Equation of Motion for a particle in a 1D Harmonic Potential. energy[x_, p_, m_, ω_] := p^2/(2 m) + (m ω^2)/2 x^2 ...
1
vote
1answer
133 views

Fit to inverse function

I am trying to fit my data to an inverse function from solving a DGL: ...
0
votes
1answer
84 views

Why does DSolve fail here?

I have the following differential equation: a*L'[t] == 40000*L[t]^(-3/2) - 1.4 - extf I want to fit the solution to data in order to acquire a. I tried DSolve ...
1
vote
1answer
123 views

Tough Calculation, novice mathematica user

I have an equation, that I've been calling $b_N(x)$ that satisfies the following identity: $$-Nb_N(x)^2=(x-N)b_N(x)+xb_N'(x)$$ where $b_N'(x)$ is the first derivative. I take the derivative then of ...
1
vote
0answers
267 views

NDSolve error message: Cannot find starting value for the variable y'

I'm attempting to solve a differential equation using NDSolve. This is the code: ...
2
votes
1answer
193 views

Analytic approximation of NDSOLVE

I have to solve the following differential equation L'[t] == L[t]^(-3/2) - 0.1 I tried DSolve: ...
0
votes
1answer
536 views

StreamPlot in a system x´ = f[x] + g[t] [closed]

StreamPlot can be used in the O.D.E system: x' = 2 x - 3 y + 3 t y' = 5 x + y - t How? ...
7
votes
1answer
566 views

Poisson PDE over a irregular region with FDM

The problem modeled by the Poisson PDE is related to the torsion of prismatic beams and I use the Finite Differences Method (FDM). I've managed to solve the equation over a rectangle region with ...
1
vote
1answer
143 views

Series expansion of InterpolatingFunction obtained from NDSolve

I am trying to obtain a series expansion of the numerical solution of a differential equation. I encounter difficulties going beyond first-order expansions which I believe might be due to my inability ...
20
votes
2answers
910 views

Phase portrait on a cylinder

It is very nice and very easy to make a sketch of a phase portrait with StreamPlot. For example, for the classical pendulum, defined by \begin{eqnarray*} \dot x&...
1
vote
1answer
384 views

Solving system of differential equations using loops

I have $F$ system of differential equations. Out of those $F$ equations except for first and last I have general form for the remaining equations (say $ dP_{i}/dt)$. Let $dP_0/dt,dP_F/dt$ denotes the ...
1
vote
1answer
113 views

NDSolve with a constant

I have a simple differential equation like this s = NDSolve[{y'[x] == A y[x] Cos[x + y[x]], y[0] == 1}, y, {x, 0, 30}] where $A$ is a constant. Mathematica gives ...