# Tagged Questions

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

28 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 ...
21 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 ...
40 views

### How to make Mathematica solve an equation in an analytic class?

Here is my small code: ...
144 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 ...
29 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 ...
33 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 ...
68 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 ...
61 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 ...
164 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 ...
69 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 ...
45 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 ...
59 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 ...
23 views

### How to embed NIntegrate inside NDSolve to solve complex functions? [on hold]

I try to solve complex nonlinear PDE functions with intergrations using NDSolve in Mathematica, but failed. To simplify my problems, I wrote the following codes, but it does not work. ...
35 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 ...
89 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 ...
135 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, ...
46 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 ...
125 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 ...
45 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 ...
38 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?
52 views

### Finding roots of Bessel function $y=3J_1(x)+xJ_1'(x)$ is returning inaccurate roots. Not Kernel bug [on hold]

I can't figure out why Mathematica is returning the incorrect roots. The first five should be 2.9496,5.84113,8.87273,11,9561, and 15.0624 according to my textbook. ...
63 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 ...
42 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 ...
45 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. ...
188 views

66 views

### How do I extract the values returned as Root[,] after solving the Sturm Liouville problem to take the square root?

$y''+\lambda y=0, y'(0)=0,y(1)+y'(1)=0$ ...
39 views

### Getting DSolve to return different solutions [on hold]

I am Using 10.0 I am trying to solve the following ODE: $$\frac{d^2 \psi}{d\theta^2}+ \nu \psi =0$$ , $$\psi(0)=0$$,$$\psi(3\pi/4)=0$$ and Mathematica only sppits out the trivial solution PS: I am ...
197 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 ...
79 views

33 views

### Recycle boundary condition

I'm having the following issue that can not find a solution , by using this boundary condition u[0, t] == 2 u[1, t] that Mathematica does not understand, someone ...
76 views

### How to “INJECT” with holding Evaluation, specifically into a WhenEvent

I am importing a set of differential and stochastic equations from an external package, and want to use Mathematica to analyse this. I am struggling building the set of equations that I can use with ...
71 views

### How can I identify which solution triggered the WhenEvent?

I am using NDSolve for problems like this: ...
82 views

### Jump boundary conditions

I'm having the following issue that can not find a solution , by using this boundary condition ...
21 views

### How does d/dx (ln(1/x) + c) = d/dx (c - ln(x))? [migrated]

I have to find the derivative of the following expression ln(1/x) + c it appears the answer is -1/x however i am told it can ...