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

learn more… | top users | synonyms

0
votes
0answers
5 views

How To Use D[ ] to Define a System of PDEs

We have two variables $u_{1}$ and $u_{2}$ which are functions of $x,t$. That is, $u_{1}=u_{1}(x,t)$ and $u_{2}=u_{2}(x,t)$. Consider the the system of partial differential equations $$\frac{\partial ...
5
votes
1answer
87 views

How to plot the solution of a Partial Differential Equation?

My attempt. I need to solve numerically the Complex Ginzburg-Laudau Equation (CGLE): $$ \frac{\partial A}{\partial t}=\epsilon A-(1+i\beta)|A|^2A+(1+i\alpha)\nabla^2A $$ I'm using a uniform initial ...
1
vote
1answer
105 views

Solving coupled differential equations with an eigenvalue

I have 2 coupled differential equations with an eigenvalue Ei and want to solve them ...
2
votes
2answers
89 views

Numerical solution of IVP for linear ODE with variable coefficient runs wild soon

Cross posted in scicomp.SE. A friend of mine showed me this initial value problem (IVP) for a linear ordinary differential equation (ODE) with variable coefficient: $$y''(x)=\left(x^2-1\right) ...
1
vote
1answer
64 views

Constants from differential equations

I have a problem, I've solved 2 diff equations and I have: ...
2
votes
3answers
540 views

DSolve will not apply assumption m ∈ Integers

I am trying to solve a linear second order ODE using DSolve which involves an arbitrary integer m. ...
-1
votes
0answers
23 views

Difference between DSolve and NDSolve [on hold]

What's the difference between DSolve and NDSolve? I know both can be used but I'm not sure of the exact difference.
1
vote
3answers
73 views

DSolve with assumptions

I'm new to mathematica and was wondering if someone could help me. I have a differential equation which is: $\frac{dy}{dx} = 1/2 y(1 - y)(a - y + a y)$ It has initial condition $y(0) = y0$ and ...
4
votes
3answers
104 views

How to fit parameters in a system of coupled differential equations

Problem description: I'm studying a chemical reaction that can be modeled by a system of coupled differential equations, and I want to use measured data to determine the parameters that appear in the ...
16
votes
1answer
417 views

Analogue for Maple's dchange - change of variables in differential expressions

Maple owns an interesting function called dchange which can change the variables of differential equations, but there seems to be no such function in Mathematica. ...
3
votes
0answers
32 views

Strange behaviour of suggestion bar to solve an ODE, potential bug

InputE^(-(1/2) x[t]^2) (Derivative[1][x][t]^2 - (x^\[Prime]\[Prime])[t]) == 0and evaluate, then press the solve ode button on the suggestion bar. It will return ...
3
votes
0answers
47 views

Diverging solution to coupled second order ODEs from NDSolve

For a physics application I am considering the radial Bogoliubov-de Gennes equations in two dimensions. In dimensionless form they read as follows $$ \begin{cases} ...
1
vote
1answer
92 views

How to combine VectorPlot and Plot3D

I have a problem with VectorPlot and Plot3D, and I have a code for both but I dont know how to plot the same figure. Here is the code for Plot3D: ...
2
votes
0answers
74 views

Creating a nonlinear phase portrait

I am putting all of my code in one block for easy copy and paste. ...
0
votes
0answers
30 views

How to solve a second order non-linear differential equation analytically?

I have the thermoelectric second order non-linear differential equation: Can some one please help me solve this equation. I need to solve for $T$. The two boundary conditions for $T$ are: $\quad ...
0
votes
1answer
35 views
0
votes
0answers
33 views

Strange difficulty with DSolveValue

Consider: sys = {u (1 - u + a v), r v (1 - v + b u)} (* {u (1 - u + a v), r (1 + b u - v) v} *) And: ...
0
votes
0answers
58 views

ODE using Mathematica [closed]

A man throws a pebble straight up. Its height $y(t)$ satisfies the differential equation $y"+\gamma y'=-g$, where $g$ is the acceleration of gravity and $\gamma$ is the damping rate due to frictional ...
1
vote
1answer
45 views

A With and ContourPlot difficulty

Consider: sys = {u (1 - u + a v), r v (1 - v + b u)}; With[{a = 2, b = 3, r = 1}, Evaluate[Thread[sys == {0, 0}]]] Which gives the output: ...
1
vote
0answers
34 views
5
votes
1answer
94 views

Detecting nearly simultaneous WhenEvents in NDSolve

I am trying to solve a system of (many) coupled nonlinear ODEs. I need to decouple some of the equations (i.e. set the time derivatives of some of the dependent variables to zero) at various points in ...
0
votes
1answer
38 views

How I can make the StreamPlot of this differential equation?

I need the StreamPlot of this differential equation but I don't know how. dp/dt=0.4p(1 - p/30), 0 <= t < 5 and 0.4p(1 - p/30) - 0.25p, t > 5
1
vote
2answers
134 views

How can I create function which contains an ODE system and its numeric solution?

I have created a function like the following, which contains an ODE system. And I want to know the effect of different parameters. ...
5
votes
2answers
117 views

Kernel quits without error in NDSolveValue

I am using the latest 10.1. The following single line command makes the kernel quit without error message (just a beep): ...
3
votes
2answers
443 views

BC for transport equation using NDSolve

First I can solve a transport equation with a source (Is it still called transport equation?) using DSolve. The form of the source serves only as an example. It can ...
9
votes
1answer
295 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 ...
6
votes
2answers
165 views

Plotting separatrices for nonlinear system

Consider the system: \begin{align*} x'&=(1-x-y)x\\ y'&=(4-7x-3y)y \end{align*} The system has a saddle point at (1/4,3/4). How can I plot the separatrices on the phase portrait having domain ...
1
vote
1answer
40 views

Plotting 2nd order ODE solution for specific values of constants

Suppose that I have a 2nd order, linear ODE and I want to plot the solution I found using DSolve. How can I plot giving specific values of ...
0
votes
1answer
59 views

Plot an ODE solution for a list of specific constants [closed]

Let's say I have a 1st order ODE and want to plot the solution. What I have been doing until now is: ...
0
votes
1answer
72 views

Plotting a NDSolve solution

I'm working on a neutrino oscillations problem, and am running into issues ploting the solution to my differential equations. The equations have non-analytic solutions so I'm using the ...
1
vote
1answer
58 views

Error with using NDSolve [duplicate]

I am currently trying to numerically solve a system of differential equations. Here is my code: ...
0
votes
1answer
105 views

Time evolution of a wave packet from the time-independent Schroedinger equation [duplicate]

Starting off with the time-independent Schroedinger equation (TISE) $\quad \quad -\frac{\hbar^2}{2m} \nabla^2 \psi + V(r, \theta) \psi = E\, \psi,$ I would like to study the time evolution of an ...
0
votes
2answers
36 views

Create a random 2x2 matrix with a repeated eigenvalue and single eigenvecor

Does anyone have a good technique for creating a random $2\times 2$ matrix that has one eigenvalue of multiplicity two, but only a single eigenvector? And more generally, has anyone written a ...
0
votes
0answers
105 views

Issue in NDSolve and unknown error

I try to solve a simple BVP: $A {v}''(t)+B v(t)=f(t)$ where $A, B$ — $5\times 5$ matrices, $f(t)$ — 5-dimentional vector-function with the use of NDSolve. Each component of $f(t)$ is a combination ...
4
votes
1answer
113 views

NDSolve to solve a PDF

I would like to solve numerically the following PDE: $$ \frac{\partial}{\partial t} p(x,t)=-\frac{\partial}{\partial x}p(x,t)+\frac{1}{2}\frac{\partial^2}{\partial x^2}[x^2\,p(x,t)],\;x\ge0,\;t\ge0, ...
0
votes
1answer
43 views
0
votes
0answers
39 views

Integro differential eq boundary difficulties

I'm trying numericaly solve simple integro-differential equation, but have some problems with boundary conditions. System: ...
2
votes
1answer
110 views

How can I invoke the solution of NDSolve to determine a parameter in my equation just inside NDSlove?

I am trying to solve a differential equation by NDSlove for $h(x,t)$. It reads $$h_t=h_{xx}-V_h-\lambda(t)$$ where $V_h$ is a given function of $h(x,t)$ denoted by ...
0
votes
0answers
106 views

What is the largest system of PDE you have solved with Mathematica?

I need to solve a system of partial differential equations. In the simplest case it consists of 2000 of PDE equations. In the most complex situation it consists up to 1 million of such equations. ...
-3
votes
1answer
61 views
0
votes
2answers
79 views

Using parts of piecewise function

I define: f[x_] := Piecewise[{{5, 0 <= x < 10}, {g[x], 10 <= x < 15}, {h[x], x > 15}}] Then I try to solve ODE: ...
0
votes
1answer
86 views

How to adjust excel data for NonlinearModelFit?

I'm really new to Mathematica and I'm facing a problem and would really appreciate some help! Maybe my question is similar to one described in How to fit 3 data sets to a model of 4 differential ...
0
votes
2answers
47 views

Problems solving/plotting a differential equation

I am trying the following: Plot[DSolve[{1220*x''[t] + 1000*x'[t] + 35600 x[t] + 4500*x[t]^3 + 2135 == 0, x[0] == 0, x'[0] == -5}, x, t], {t, 0, 10}] Could ...
2
votes
2answers
127 views

Fighting the Current

I am working a problem I am working on. A boat start across a river at the point (c,0) and points its bow directly across the river to the point (0,0). The parameter a is the river current velocity, b ...
1
vote
3answers
155 views

Solving differential equation

I want to solve the following differential equation DSolve[y''[x] == (λ y[x])/Sqrt[-1 + x], y[x],x] But do not know how to actually solve it. Any suggestion?
2
votes
5answers
178 views

Simplify an integral involving trigonometric functions

We were doing variation of parameters in differential equations tonight and had to do the following integral, the result given by hand calculations. $$\int \cos t\tan^2 t\,dt=\ln|\sec t+\tan t|-\sin ...
-1
votes
1answer
73 views

Using InterpolatingFunction in equations

I tried to use InterpolatingFunction in another ODE, but it doesn't work, because it seems for Mathematica that 2 functions are unknown instead of one. ...
0
votes
1answer
71 views

NDSolve::bcedge: Boundary condition not specified on a single edge

NDSolve::bcedge: Boundary condition c[t,5]==Cout is not specified on a single edge of the boundary of the computational domain. >> I'd like to plot $\frac{\partial}{\partial ...
1
vote
1answer
86 views

Solving ODES in order to find unknown constants

I am trying to solve a set of ODES which represent a second order consecutive chemical reaction. ...