# Questions tagged [differential-equations]

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

698 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
0answers
802 views

### DSolve returns an answer that is not a solution

Investigating DSolve misses a solution of a differential equation, I came across this odd behavior of DSolve. The following DSolve command returns an answer to the ...
0answers
161 views

### Workaround for DSolve in V 12 when it gives undefined as solution to 1D heat PDE?

Reported to WRI CASE: 4278450 Comparing the following, all done from clean kernel The strange thing is that V 12 can solve this same PDE without the assumptions ...
0answers
502 views

### From where to learn finite element method?

Can anyone reccomend me a book or site for learning finite element method with mathematica, besides wolfram language official site, https://reference.wolfram.com/language/FEMDocumentation/tutorial/...
0answers
104 views

### What does DSolveDsolveLaplaceTransformDump mean in solution of DSolve?

Reported to Wolfram with CASE number 4028951 on 11.3, ...
0answers
301 views

### Spurious Error Messages from NDSolve when Using WhenEvent with Time Delays

Bug introduced between versions 10.1 and 10.4, and resolved in 11.3. Using either 11.2 or 10.4 on Windows 10 (64 bit), I am unable to reproduce the answer by March to question 99576. Specifically, ...
0answers
487 views

### benchmark of Mathematica's FEM?

Does anyone have any numbers on how mathematica compares to other commercial (eg ANSYS) or free FEM softwares (eg. FreeFEM++, FEniCS, elmer)? If this is too vague say solving the diffusion equation in ...
0answers
3k views

### Precision and accuracy in NDSolve and NMinimize

I use NDSolve quite a lot and have noticed that setting values for AccuracyGoal and ...
0answers
92 views

### How to modify NDSolveStateData without crashing the kernel?

Probably a hard question, but it's better to cry out loud. Reminded by Chris K, I noticed my fix function has been broken since v11.3. After some checking, I ...
0answers
176 views

### DEigenvalues and NDEigenvalues return different values

In the following example, DEigenvalues and NDEigenvalues return different results despite having identical arguments. Does anyone know why? (I use Mathematica 11.3) ...
0answers
476 views

### How to put 2 or 3 NeumannValue conditions?

I want to solve numerically the heat exchange of a fin in transient state and two-dimentional, here is my attempt but instead of put DirichletCondition on ...
0answers
86 views

### DSolve fails with four variables

How is that DSolve has no trouble solving this: DSolve[{ D[G[x1,x2,x3],x1]==0, D[G[x1,x2,x3],x2]==0, D[G[x1,x2,x3],x3]==0 },G[x1,x2,x3],{x1,x2,x3}] but it fails ...
0answers
86 views

### Higher order Laplacian flows

Given disjoint surfaces $q_i$ in 3D and their 1D boundary curves $\partial q_i = \gamma_i$, I seek a surface $p$ that joins the $q_i$, where $p \cup q_i$ forms a (piecewise) $C^k$ surface that ...
0answers
75 views

### why DSolve gives different looking solution when using u vs. u[x,t] in the call?

Should the answer by DSolve be different when calling it as DSolve[ode,u,{x,t}] vs. ...
0answers
284 views

### How to use periodic boundary conditions on the derivative of u[x,y,z]?

I have a structural mechanical problem, where i try to calculate the strain field in the modelled material. I simplified the problem in the following example: There is a material with two phases, ...
0answers
374 views

### Mathematica packages for numerically solving the Klein-Gordon and Dirac equations?

The book "Advanced Visual Quantum Mechanics" by Thaller includes Mathematica software packages for the numerical solution of the Klein-Gordon equation and the Dirac equation (subject to user-defined ...
0answers
97 views

### Order of boundary condition inside DSolve matters?

Bug introduced in 11.1.1 or earlier. Consider the following example: ...
0answers
404 views

### Symbolic Weak Form

Usually I write the weak form by hand for my FEM code, but it's a little annoying and mechanic sometimes. So, I wonder, is there any way to generate the symbolic weak form in Mathematica? For ...
0answers
74 views

### EventLocator in NDSolve seemingly hogs memory in Windows, but not Mac

I have written a function (ftest in my code) which does the following: Takes a vector {x,vx,vy} as input Numerically solves (with high precision) a specific system ...
0answers
267 views

### Applying method of dominant balance directly to an equation

The non-linear equation $y''(x) = \sin k x y$ for $k\in \mathbb R$ cannot be solved analytically to yield a closed form solution. However, by applying the method of dominant balance one may obtain ...
0answers
227 views

### DSolve breaks when the ordering of independent variables aren't proper?

Bug introduced in 5.2 or earlier and persisting through 12.0. I encountered this when trying to solve this problem with DSolve: ...
0answers
1k views

### Modelling Hysteresis with a Differential Equation

I want to implement the bulk ferromagnetic hysteresis model (mostly the Jiles-Atherton Model), see http://drum.lib.umd.edu/bitstream/1903/6043/1/PhD_99-1.pdf page 44 equation (30). The needed ...
0answers
2k views

0answers
720 views

### Alternatives to FiniteElement as Spatial Discretization Method for NDSolve

Finite Element Programming: [...] It is possible to skip this section and continue with the discretization stage and make use of the initialized data structures ...
0answers
58 views

### ParametricNDSolveValue with SparseArrays

I want to solve large, sparse matrix systems with ParametricNDSolveValue and was hoping to use SparseArray to speed up the computation. However I find that ParametricNDSolveValue doesn't seem to work ...
0answers
50 views

### SystemException[“MemoryAllocationFailure”] with DSolve and using TimeConstrained

Report to WRI CASE:4238154 There is a similar question where SystemException["MemoryAllocationFailure"] is generated by DSolve ...
1answer
193 views

### Neumann boundary condition is not satisfied

I want to solve the diffusion equation on a disk centered at (0,0) with a radius of 1. I also want the flux at a radius of 0.8 to be zero. I have this initial condition at time zero: ...
0answers
40 views

### DSolveValue gives a confusing result for a simple vector-valued equation

If we use DSolveValue on just about the simplest vector ODE I can think of, it gives the following result. ...
0answers
111 views

### DSolve, NDSolve with WhenEvent Give Incorrect Solution for Simple ODE

NDSolve Results On the course of addressing question 181974, I encountered the following problem. ...
0answers
267 views

### Spectral problem for differential vector operator (calculation of EM field in a cavity)

I know that mathematica has a DEigensystem and NDEigensystem which allow one to find eigenfunction and ...
0answers
168 views

### How to find an exact solution to this ODE?

I am trying to find an exact solution to this differential equation, ode = f'''[x] + f[x]*f''[x] - f'[x]^2 == 0 For which, I already know the exact solution ...
0answers
171 views

### Why is this non-vectorized NDSolve faster than vectorized?

I'm trying to get a feeling for stochastic population models. A simple starting place is a birth-death process whose deterministic limit is the logistic equation. The idea is to model the probability <...
0answers
175 views

### Laplace equation in square region

Suppose we have square region with side $l$ and rectangular with sides $l$ and $l+e$, $e=0.001$. ...