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

learn more… | top users | synonyms (3)

0
votes
0answers
46 views

Fitting 6 differential equations to 6 sets of data

I am using the newest Mathematica and I am trying to fit a set of 6 linear differential equations to a set of experimental data for each of the six variables. I am a baby beginner at this and could ...
2
votes
3answers
112 views

Strong glitch with Evaluate in Manipulate

I'm experiencing a strong evaluation problem for a Manipulate box and I can't find the best solution yet. Here's a MWE to show the glitch : ...
2
votes
0answers
83 views

NDSolve for unsteady Taylor-Goldstein equation, can't use MethodOfLines

I am trying to solve a relatively simple system of linear PDE in two variables (z,t). The equations are basically the Navier-stokes equations linearized around a constant background velocity profile ...
0
votes
0answers
42 views

Impose conditions on differential equations

This is my first post in Mathematica Stack Exchange, so thanks for your help and sorry for any rookie mistakes. I am solving the Einstein equations for two vectorial pertubations h_Alpha and v_Alpha ...
3
votes
1answer
205 views

NDSolve mixing many scalar and vector equations

I have a set of scalar equations in many unknowns, which I want to combine with a vector equation inside an NDSolve. The equations are a mix of differential and algebraic equations. A set of scalar ...
1
vote
1answer
39 views

“MemoryAllocationFailure” in DSolve

Mathematica 10.4, Windows 10, 6 GB RAM, 64 bits When I trying to solve DSolve[{v'[t] == 100 - 0.15 v[t]^2, v[0] == 30}, v[t], t] I get ...
14
votes
2answers
257 views

What can one do with extremely stiff problem in NDSolve?

Consider the following illustrative problem: $$ \frac {\partial f} {\partial t} = \frac {\partial} {\partial x}(x f) + \frac {\partial} {\partial x}(f \frac {\partial f} {\partial x}) $$ This is ...
2
votes
1answer
45 views

Clearing Problem with DSolve

Consider: This is a classic mistake (using a single equal sign) that has been discussed on Mathematica Stack Exchange before. I know that Clear[Derivative] will cure the problem. However, consider ...
0
votes
1answer
32 views

NDSolve point animation not working [closed]

The below code plots the path of a ball bouncing but loses 90% of its velocity each time it bounces. (adapted from https://reference.wolfram.com/language/ref/WhenEvent.html). ...
0
votes
1answer
59 views

I am having trouble with the DSolve function [closed]

DSolve[{x''[t] + 0.01 x'[t] + x[t] == E^(It), x[0] == 1, x'[0] == 1}, x[t], t] This give the following error :"Equation or list of equations expected instead ...
2
votes
0answers
58 views

How to change the boundary condition near singularity

The following system has a singularity around $x=0.451$. What I want to do is to detect when the solution gets close enough to this singular point (in this example I'm using the threshold of $0.001$), ...
1
vote
0answers
50 views

DSolve does not give an analytic solution for a tutorial input [closed]

I am just starting with Mathematica and I have tried to solve PDE following DSolve, with first three inputs from the Elliptical Partial Differential equations ...
0
votes
0answers
57 views

Nonlinear PDE system in a rectangle

This question is posted by two physics students working together. We are having some trouble while trying to solve Navier-Stokes equations, The main problem is that, while our set of partial ...
0
votes
1answer
95 views

Plotting graphical analysis of chaotic behavior

I am trying to plot some results of chaotic behavior in a system based on a paper for one of my classes. I am a physics undergrad. I found some sample code online springy.nb which I have been trying ...
0
votes
1answer
92 views

How to omit the corrupt value in a program?

I've such errors: Coordinate .. is not a floating-point number and can't correctly fix it. The problem occures when x==Log[2] The program works correctly ...
0
votes
0answers
35 views
2
votes
1answer
119 views

Differential Equation with Numerically Integrated Boundary Conditions

I have a second order differential equation that I'd like to solve numerically, and then integrate its solution twice to get the parametric equations of a curve. The ODE is: $2 ...
0
votes
0answers
63 views

solving system of ODEs in Mathematica

I am using the DSolve to solve the following 2 ODEs \begin{equation} \begin{cases} &\dot{y}(t)=z(t)-y(t)-\frac{c_1c_2y(t)}{-c_1c_2e^{-t}t+c_1(1-c_2)+c_3e^{-t}}\\ ...
3
votes
1answer
66 views

Dsolve svars error: Equations may not give solutions for all “solve” variables

New to Mathematica, trying to solve a set of coupled differential equations related to a geodesics/Calculus of Variations problem. More specifically, I am trying to solve the two Euler-Lagrange ...
0
votes
0answers
55 views

NDSolve for an (slightly) discontinuous Piecewise input

I define a Piecewise discontinuous function V[f[t]], which depends on a function f[t] I want to get numerically using ...
0
votes
0answers
27 views

Solving a set of differential equations with DSolve

I am trying to solve exactly a system of differential equations. The set of differential equations is given by : ...
0
votes
0answers
81 views

How to numerically solve this non-standard Green's function using mathematica?

$$\left(u(1-u^2)\frac{d^2}{du^2}-(u^2+1)\frac{d}{du}-\frac{au}{(1-u^2)}-\frac{bu^3}{(1-u^2)}+c\right)G(u,u')=mu^2\delta(u-u'),$$ where $a$, $b$, $c$, $m$ are constants, $u\in(-1,1)$ and the boundary ...
0
votes
0answers
44 views

How to numerically solve a Green's function using mathematica?

Suppose we have a Green's function $$LG(x)=\delta(x),$$ how to numerically solve it by mathematica? Can the mathematica read the delta function directly? For analytical calculation, we usually ...
0
votes
0answers
56 views

Problem with piecewise function as source function of NDSolve

It seems that something wrong is with my piecewise source function u[t,x], so I get such an error: ...
0
votes
1answer
65 views

Function Evaluation in Table [closed]

Apologies in advance: this is a really basic question, and I am sure it has been answered, but I don't know enough mathematica to even find the answer. I have a differential equation I want to ...
0
votes
2answers
96 views

Solving recursion relation from power series

I am interested in solving differential equations in the form of power series. Let's say we have following equation: $$f^{\prime \prime} (\rho) + \left( \frac{2 e^{-k \rho}}{\rho} - \varepsilon ...
0
votes
2answers
78 views

Need to find value of NDSolve at specific value in range

I am working with NDSolve to plot an ODE as a function of time. I need to find a way to figure out the exact value X of the solution for some specific time, say t*. Even better, a way to figure out ...
1
vote
1answer
62 views

Differential equation [closed]

How does one solve this differential equation? $$y' = \frac{x^3 + 2x^2 y - y^3}{x^3 + x^2 y}$$
2
votes
0answers
58 views
0
votes
0answers
68 views

How to solve this PDE in mathematica?

I have been trying to solve this PDE, but Mathematica didn't give me a solution. Here it is: DSolve[{b*(-(3/2)*(9*Cos[[Theta]] + Cos[3 [Theta]])Csc[[Theta]]^7 D[h[[Theta], t], [Theta]] - (11 ...
-1
votes
1answer
54 views

NDSolve and vectors [closed]

I think I can't find solution of equation below because I need use scalars, not vectors as I have tried. But I don't know how I can do express the problem with scalars. I have read Wolfram Language ...
2
votes
0answers
56 views

Solving a system of first order differential equations [closed]

I've got a set of equations of motion here (describing the interior of a charged black hole, if you're interested). It's just four first order ODEs, so Mathematica should be able to solve this easily ...
0
votes
0answers
34 views

Why differential equation reduces? [duplicate]

I am trying to ND-Solve the following equation: eqns = {x''[t] == -x[t], x'[0] == 0, x[0] == 4} The problem, is when I enter it, Mathematica reduces it to ...
0
votes
1answer
78 views

NDSolve for a system of nonlinear coupled pde

I try to solve a system($n$) of coupled nonlinear 1-d PDE's which is basically coupled sine-Gordon equations (through 2nd spatial derivative) in time and space. I am interested in the dynamics when I ...
4
votes
1answer
89 views

Overriding default ExtrapolationHandler

When modeling systems with WhenEvent and NDSolve I found that the default extrapolation settings used in the returned ...
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 ...
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
57 views

Using DSolve for r(t) y'(t) = y(t) r'(t) [closed]

I am fairly new to Mathematica and I was trying to solve this equation r(t) y'(t) = y(t) r'(t) using DSolve. How can I adapt it to the DSolve function formalism?
0
votes
1answer
86 views

DEigensystem and NDEigensystem simply not working in Mathematica 10.0 [closed]

I am trying to use DEigensystem and NDEigensystem to check some theoretical derivations I have recently done. I am running ...
1
vote
2answers
127 views

How to plot coupled Differential equations with changing parameters?

I'm a new user of mathematica. What I'm trying to do is following: I have two coupled non-linear differential equations. Say $x'(t)$ = $A− (B +1)x + x^2 y$ $y'(t) = Bx − x^2 y$ What I want is to ...
0
votes
0answers
46 views

Overdetermined Differential Equation is not overdetermined

I want to solve a system of two differential equations ...
0
votes
1answer
118 views

Animate spherical pendulum

I want to plot the position of Spherical pendulum. First I tried to plot the simple pendulum: ...
1
vote
1answer
72 views

Numerical integration plot

I am trying to do a numerical integration plot in Mathematica. Specifically, I want to plot $r$ versus $\tau$ where, $$\frac{dr}{d\tau}=\sqrt{2F(r)}$$ $F$ being a quartic polynomial. The limit of $r$ ...
3
votes
1answer
124 views

how to speed up FEM-NDSolve for similar boundary conditions?

I need to solve several stoke's equations for very similar boundary conditions, which represent particles in a fluid with different velocities U[i]. I'm interested in computing the force acting on the ...
2
votes
1answer
72 views

Heat transportation equation 2

I would like to thank you for all the suggestions for the problem with the heat transportation equation I have posted time ago! I rewrited the code along with the implementation of the Neumann ...
3
votes
2answers
89 views

Animating 3D Parametric Plot

I'm trying to create an animation of Lorenz equation as a function of t. Here is the numerical solution of the equation: ...
1
vote
1answer
83 views

Problem with Finding Minimum of a Function Solved With NDSolve

I am trying to find the minimum of a particular function, which is solved using NDSolve. What is making it a little more complicated is that it is not the only function solved within the NDSolve and ...
0
votes
0answers
28 views

Use a Poincare map and Poincare section to show that the equation has a T-periodic solution [duplicate]

The equation is $x'(t)=-x(t)^3+h(t)$, where $h(t)$ is any smooth T-periodic function. Give an argument using a suitably constructed Poincare section/map that shows the equation has a T-periodic ...