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

learn more… | top users | synonyms (3)

0
votes
0answers
44 views

Time dependent/dynamic parameter in the system of ODEs

I am wondering what would be the best way how to code the following situation. I have a solution but I don't think it as elegant as it could be. I have a system with 10 ODEs and lot of parameters. ...
0
votes
1answer
35 views

Define inverse Laplace operator as map from functions to functions

I am trying to implement the inverse Laplace operator through P[g_] := h /.NDSolve[{h''[x] == g[x], h[0] == 0, h[1] == 0}, h, {x, 0, 1}][[1]] If I now define ...
0
votes
2answers
35 views

Beginners problem, Do Loop, Eigenfunction iteration

I am trying to find the first Eigenfunction of the Laplacian (in 1D), i.e. a solution of $$ u''(x)=k u(x)\\ u(0)=u(1)=0 $$ with minimal $k>0$ (in this trivial example, I actually know the analytic ...
1
vote
1answer
143 views

Plotting Phase Diagram

I have the following problem that I'm sure Mathematica can handle, but it's not working for me! In the following code, I'm trying to replicate the Ramsey Model Phase Diagram. In fact, I want now ...
2
votes
1answer
89 views

Dog chases his tail ! - “parametric differential/Integral equation”..?

I have the following situation where I am interested in the function $m(t)$ $$ \frac{dm}{dt}=4T(t)^{3}+T(t)^{2} $$ $$ T(\tau)=T_{0}-\int_{0}^{\tau}(\frac{dm}{dt})dt*Q_{S} $$ Is there a way to solve ...
0
votes
1answer
193 views

Using data from NDSolve into a secondary equation

I have this ODE, f''''[y]+Z1*(6 f''[y]*f'''[y]*f'''[y]+3 f''[y]*f''[y]*f''''[y])- N1*N1*f''[y] == 0 subject to boundary conditions ...
0
votes
0answers
85 views

FourierTransform a differental-equations [duplicate]

I want to use the command FourierTransform to transform a differential equations from time - domain to complex frequency domain but it don't works. My code is ...
0
votes
0answers
26 views

Implementing more accurate boundary conditions in NDSolve

I want to solve a second order differential equation numerically. The boundary conditions are needed to be imposed at z=0 and at ...
1
vote
1answer
130 views

Differential Equations and Unit Notation

What do I get the Syntax::sntxi: Incomplete expression; more input is needed . error when I try to use the math palette for derivatives? What does the dot in ...
5
votes
2answers
272 views

Simpler way to type Derivative[1, 0][u][0, x]

Though not unbearable, I always feel a little nervous when typing partial derivatives at a specified place e.g. $u^{(1,0)}(0,x)$, which happens a lot when setting initial/boundary conditions for PDEs. ...
14
votes
3answers
5k views

Solving a time-dependent Schrödinger equation

I want to solve the time-dependent Schrödinger equation: $$ i\partial_t \psi(t) = H(t)\psi(t)$$ for matrix, time-dependent $H(t)$ and vector $\psi$. What is an efficient way of doing this so that ...
1
vote
1answer
167 views

NDSolve memory usage [closed]

I am trying to solve numerically a system of linear ODEs with some quickly varying driving functions. The basic command is: ...
4
votes
1answer
489 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 ...
2
votes
2answers
124 views

How do I solve Schrödinger equation with a piecewise periodic Hamiltonian? [duplicate]

My Hamiltonian is HI(t) = Piecewise[ {{H1, n τ <= t <= (n + 1/2) τ}, {H2, (n + 1/2) τ <= t <= (n + 1) τ}}] where ...
5
votes
0answers
136 views

Dealing with a nonlinear PDE using FEM

As far as I can see in the Mathematica documentation center, there are no built-in ways for solving nonlinear PDEs using Finite Element Method (FEM). I am dealing with a simple nonlinear ...
4
votes
4answers
254 views

Trying to model a chemical reaction - NDSolve::litarg error

I am trying to model a chemical reaction with NDSolve. I have specified the known differential equations and parameters, ...
2
votes
2answers
66 views

How to simplify a differential equation?

Now I get a differential equation like this: (-1 + x) (1 + x)^2 == (2 (1 + x^2) (1 + x^4))/(-2 x (1 + x + x F[x]) + (-1 + x^3) F'[x]) I want to get the standard ...
1
vote
1answer
103 views

NDSolve with coupled ODE's and unknown singularities

I have two coupled ODEs that I am trying to solve numerically. It appears that there is a singularity in the solution to the equations which I am unsure how to get past. Both functions $\alpha$ and ...
0
votes
0answers
43 views

Check if interpolated function passes through a specified region

Is there an efficient way to check if a curve that was obtained as an Interpolating Function (from solving differential equations numerically) passes through a given region? I am solving for two ...
2
votes
1answer
147 views

solving coupled differential equations with parameters subject to optimization problem

I have an equation, but I don’t know how to program it in Mathematica: x[0]==1; y[0]==1; x'[t]==(a+b) x[t] y[t]; y'[t]==a b y[t]/x[t]; where ...
3
votes
0answers
64 views

ColorBars or PlotLegends for ElementMeshSurfacePlot3D

I want to produce a colorbar legend to the output from ElementMeshSurfacePlot3D (Mathematica 10). The pde was solved using NDSolveValue (FEM pdesolution is the resulting interpolating function) over a ...
0
votes
1answer
37 views

Find a independent variable value from the numerical solutions of a system of differential equations

I have calculated the solutions for the following system of differential equations using NDSolve ...
13
votes
1answer
106 views

Empty WhenEvent action crashes kernel

Bug introduced in 9.0 and persisting through 10.4.1 WhenEvent is new in 9.0. This is an example from the docs, slightly modified (action is wrapped in ...
-1
votes
0answers
68 views

Limitation of MMA in solving complicated functions?

I try to solve nonlinear PDE functions with NDsolve and NDIntegrate in MMA. And it seems that MMA will not work if the functions become too complicated. I do not know why. Below are my codes which ...
1
vote
1answer
167 views

Differential equation in two variables

I have a system of second order differential equations, with two independent variables, which represents a trajectory. I tried solving this system using NDSolve and ...
2
votes
1answer
48 views

Asymptotic Radial Wave Equation

I'm trying to reproduce the following solution to the ODE, $$\frac{d^{2}u}{d\rho^{2}} = \frac{l*(l+1)}{\rho^{2}}u$$ Solution: $$u\left(\rho\right) = C\rho^{l+1}+D\rho^{-l}$$ What I've tried in ...
0
votes
2answers
68 views

Using an interpolating solution of NDSolve into an algebraic expression

I would like to use the interpolating solution of NDSolve into an algebraic expression without creating a system of a differential equation and algebraic equation. The example is the following, the ...
0
votes
0answers
41 views

How to use NDSolve for mixed types of several equations?

I have three equations, one is partial differential equation and another two are not. They are just normal algebra equations but coupled with the differential equation. I wanna use NDSolve to get ...
3
votes
0answers
112 views

Mathematica Newmark Optimization

Here is my Mathematica code which implements the Newmark method to solve a equation of motion. The variable "ag" contains the accelerations from an earthquake record. Is it possible to optimize this ...
6
votes
1answer
148 views

Constrain Locator to specified region

How can I constrain the locator to stay within the region defined by RegionPlot? When the Locator remains within the region, ...
3
votes
2answers
51 views

Numerically solving coupled ODE's with a parameter as initial condition

i'm currently trying to numerically solve a set of coupled ODE's to obtain the functions p(r), h(r) and m(r) in the range of r1 <= r <= r2 with initial conditions m(r1)=a=const and ...
2
votes
1answer
140 views

solving Stokes equations for Force=0 boundary condition

Say you want to solve the movement of a particle in a viscous fluid such that the force on the particle is zero. In general, the solution is trivial since one has no-slip boundary conditions, ie the ...
0
votes
1answer
49 views

Setting Initial Conditions in a System of Differential Equations that don't start at t=0

I have a system of six differential equations (linear for the sake of the discussion). I want to set the initial condition to start at time t=18 and not time t=0. I know the values of the dependent ...
0
votes
0answers
43 views

Extract stiffness matrix from NDSolve

I am solving a set of ODE equations with NDSolve with more than 25 equations. I was wondering is there any way Mathematica gives me the mass and stiffness matrices?
2
votes
1answer
74 views

Does NDSolve iterate faster in a region where the system being solved is in equilibrium?

I was wondering how exactly the time it takes for NDSolve to iterate through a certain region of the domain depends on the local behaviour of the given system being solved. My guess would be the ...
0
votes
0answers
49 views

How to speed up NDSolve?

I have a program where I need to iterate NDSolve thousands of times for different configurations of a diffusion problem. The thing is, I only need about 3 digits of accuracy. I looked thru previous ...
0
votes
0answers
51 views

Lyapunov exponent of non-autonomous differential equation

I am trying to calculate the Lyapunov exponent for a non-autonomous system using M Sandri code, but it always shows me error "infinite expression 1/0 encountered". Does anyone know about this? Mr. ...
3
votes
1answer
201 views

Solve a nonlinear PDE equation with a Neumann boundary condition

I am trying to use Mathematica 10 to solve a PDE $$u_t=u_{xx}+u_{yy}+u(1-u),$$ in the unit disk $(x,y) \in D=\{(x,y):x^2+y^2<1\}$, with the Neumann boundary condtion $$\frac{\partial u}{\partial ...
13
votes
3answers
529 views

Particles bouncing in a 3D box

I made this Manipulate code which shows a pack of particles randomly moving inside a box. It has a huge performance issue. What should be the best way in doing ...
0
votes
1answer
118 views

Trouble solving a system of 2 ODE's with NDSolve

I'm trying to solve a system of two coupled second-order ordinary differential equations using Mathematica's NDSolve. The functions of r that appear in the equations are ...
10
votes
2answers
2k views

Error entering equation in DSolve

I entered a command incorrectly as follows: DSolve[{y'[x]=y[x]},y[x],x] I am now experiencing: ...
0
votes
1answer
34 views

Trying to solve ODE. Need help with error messages [closed]

With this code I get an error-message: Equation or list of equations expected instead of True in the first argument. Apparently y'[0] == v0 gives ...
0
votes
0answers
69 views

Force NDSolve to use finite difference

How to force NDSolve to use finite difference instead of FEM. I have the following code: ...
0
votes
0answers
70 views

Boundary conditions for NDSolve

I get the following message: Boundary values may only be specified for one independent variable. Initial values may only be specified at one value of the other independent variable. >> from ...
0
votes
0answers
48 views

Eigenelements of conductivity equation

I am trying to calculate the eigenvalues and eigenfunctions of the conductivity equation in an annulus. In particular I am looking for $(\lambda, u)$ s.t. $$ \begin{cases} \Delta u = \lambda u & ...
11
votes
1answer
218 views

NDEigensystem returns incorrect eigenvalues for 2D coulomb problem, eigenfunctions contain discontinuity

I posted a similar question a short time ago regarding the 3D Coulomb problem. Jens' excellent answer to this thread allowed me to obtain the correct eigenvalues and eigenenergies for that system. I ...
0
votes
1answer
35 views
0
votes
0answers
32 views

Solving Time-Dependent PDE's Analytically with no non-trivial Boundary Conditions

I'm having some trouble solving a partial differential equation that is supposed to model the propagation of a scalar field in Schwarzschild spacetime using spherical coordinates. I'm using the DSolve ...