# Tag Info

1 vote
Accepted

### Error: Using NeumannValue with NDSolve correctly

The solution isn't unique, you need additional Dirichlet condition. Examplary I take v[0,1]==0. NDSolve command with changed y-range (see my comment) follows to <...

### Shallow-water equations on a rotating sphere by using the NDSolve method with a purpose of earthquake-generated tsunami simulation

Its difficult or hopeless to NDSolve waves on a sphere with boundaries because the governing Laplacian has discrete spectrum. Eigenfunctions of the Laplacian with non-integer eigenvalues belong to ...
• 3,965

### Using Mathematica to solve Newton's Equations

If I understand your question and comment right try NDSolveValue : ...

### How to solve the partial differential equations?

Check the Neumann-conditions! Try ...

### How to solve the partial differential equations?

The inhomogenuous equation $$\Delta_{x,y}\left(-e^x \ \sin (\pi y) +g(x,y)\right) == \ e^x \ \sin (\pi y)$$ results in the homogenous equation. $$\Delta_{x,y} g(x,y)==0$$ Reformulate the other ...
• 3,965

### Speed up finite element method solver with small features

The main problem is that your initial condition for T is not consistent with the PDE. You define: ...
• 40.1k

### Speed up finite element method solver with small features

To solve this problem, we put characteristic scale L=300 and then rescale all geometrical parameters on L, as result we have <...
• 46.4k
1 vote
Accepted

Workaround: ...
• 239k
1 vote

### Unwanted Oscillation in Solution Occurs When Solving Poisson-Nernst-Planck Equation

As one of the possibly solutions we can decrease accuracy and use explicit method as follows ...
• 46.4k
1 vote
Accepted

### Singularity error when solving Poisson-Nernst-Planck equation

We can use FEM to solve this problem as follows ...
• 46.4k
1 vote

### Scaling DiracComb does not work properly in NDSolve

To avoid troubling MMA with complicated functions(DiracComb and DDE), we use WhenEvent instead. ...
• 889

### How to increase the mesh size in NDsolve

You can do that though a Method option but the exact option dependents on whether the DE is stationary or time dependent: ...
• 40.1k
Accepted

### Singularity or stiff of system suspected

There is a typo in predator equation. We can compare Jacobian to that in the paper as follows ...
• 46.4k
1 vote

### Finding and Plotting Minima

You can use ComplexExpand to get rid Abs and so forth, if z, ...
• 3,569

### Solving a stiff non-linear initial value ODE for a set of parameters

There is a typo with WhenEvent usage in the code. In the computation process Z becomes very small and negative that lids to the complex solution. Therefore the ...
• 46.4k
Accepted

### How to Solve a DAE System with Integral Constraint and Moving Boundary in Mathematica?

To compute constraint, we can use code as it is with a small modification as follows ...
• 46.4k

### Boundary conditions issue with Laplace equation

In Version 14.0 this works: ...
• 40.1k
Accepted

### Boundary conditions issue with Laplace equation

We can use OpenCascadeLink to generate mesh for this problem as follows (note, we rescale l, r, rs because the solution does not ...
• 46.4k
Accepted

### Numerically solved PDE of Ornsteinâ€“Uhlenbeck process on 2-Simplex violates conservation of probability

The initial condition was inconsistent and needed to exclude the boundary: ...
1 vote

### Boundary conditions issue with Laplace equation

It works if you remove one of the two balls. For example : : ...
• 18.7k

### Numerically solved PDE of Ornsteinâ€“Uhlenbeck process on 2-Simplex violates conservation of probability

Roughly your equation is $$\ \partial_t P =-\left(\frac{3}{2}\ - \ y \right) \ \partial_y \ P \ -\left(\frac{5}{2}-x \right) \partial_x \ P +\frac{1}{2}\Delta \ P \ + \frac{3}{2} \ P$$ by the fact ...
• 3,965
Accepted

### Numerical solution with both Neuman and Drichlet BC at one point

Using my code from here we have ...
• 46.4k

### NDsolve does not solve equation

Here is an example that actually runs without messages, fixing a number of errors in @zeraoulia rafik's code and adjusting parameter values to make a nice plot: ...
• 239k

### NDsolve does not solve equation

According to the recommendations given by @Nasser in the comment above, I have tried to show your solutions by adding the missing numerical values for your parameters and initial condition values. I ...
• 1,049

### Solving a large number of coupled non-linear equations

For this equation it could be better to use implicit difference scheme in combination with FEM as follows ...
• 46.4k
1 vote

...
• 62.9k
Accepted

### Issue with a Piecewise Diffusion Problem

Change your bc to NeumannValue to force a FiniteElement-solution: ...
1 vote
Accepted

### Young Laplace equation with shooting method

We can use your code as it is with small modification ...
• 46.4k
Accepted

### How to NDSolve a PDE that contains the integral of the solution?

We can solve this problem using method of lines as follows ...
• 46.4k
1 vote
Accepted

### NDSolve::nlnum and NDSolve::bcedge

The first comes from the NDSolve here, where it returns NDSolve::nlnum, To remove NDSolve::nlnum add this method: ...
• 147k

### Are there contributed Mathematica programs which yield the normal form and Lyapunov coefficients of an ODE at a bifurcation point, in dimens 3 or more

This is not an answer to the question or at least only a partial one. It may help to get started though. A Hopf bifurcation at an equilibrium, if I understand correctly, happens when the Jacobian has ...
• 59.6k
Accepted

### Are there contributed Mathematica programs which yield the normal form and Lyapunov coefficients of an ODE at a bifurcation point, in dimens 3 or more

I have no method but Let's refine the script you provided to ensure it handles the computation of the L1 coefficient and the normal form correctly. I'll integrate improvements while maintaining the ...
• 1,049

### Can I solve system of differential equation in a matrix form?

You don't need column vector form for Mma/WL: row will suffice, e.g. v = {x[t], y[t]}; a = {{6, -4}, {1, -2}}; DSolve[D[v, t] == a . v, v, t] -> ...
• 62.9k

### I have a problem with this code because it is time invertible symmetric but it is not working

Although the question is not entirely clear to me, I presume that a solution symmetric about x = 0 is desired. To produce such a solution, both the ODE and its ...
• 61.8k

### Stiff differential equation

Use asymptotic expansion at $x=0$ of linearization of the ODE (the nonlinear term vanishes to high order there) to help get NDSolve started. Then used Brent's ...
• 239k

### Stiff differential equation

You can't solve with BC at infinity numerically. Since this BVP, you can try shooting method. Best I could make run to is up to $x=4$ ...
• 147k

### Integrating over an interval

[H]ow do I integrate v from xmin to x_1, then x_1 to <...
• 239k

### Solving partial differential algebraic equation in Mathematica

The latest packages for PDE use DirchletCondition as a container for the boundary condtions of the unknowns and for Neumann boundary conditions replace the right side $eq==0$ by a source term, that ...
• 3,965

### A differential equation of the form $x''+a^2 x=0$ does not lead to solutions involving sine and cosine

I think the solution you obtained from Mathematica is indeed correct but presented in an exponential form. To obtain the solution in terms of sine and cosine functions, you can use the fact that the ...
• 1,049