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

learn more… | top users | synonyms (3)

2
votes
0answers
65 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
52 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
60 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
99 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
93 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
38 views
2
votes
1answer
122 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 \kappa''(s)+\kappa^3(s)...
0
votes
0answers
65 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}}\\ &\dot{z}(t)=-z(t)+\frac{...
3
votes
1answer
68 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
63 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
28 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
84 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
51 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
61 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
70 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 solve....
0
votes
2answers
99 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
88 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
67 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
93 views

How can I solve this PDE?

I have been trying to solve this PDE, but Mathematica didn't give me a solution. Here it is: ...
-1
votes
1answer
57 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
60 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
90 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
102 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
543 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
156 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
98 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
141 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
125 views

Animate spherical pendulum

I want to plot the position of Spherical pendulum. First I tried to plot the simple pendulum: ...
1
vote
1answer
78 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
136 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
73 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
98 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
84 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 ...
3
votes
3answers
167 views

Can I plot the solution to an ODE without using dsolve?

I have the following ODE $\qquad x'= -x^3+\sin(t)$ I want to find numerically the initial condition $x(0)$ corresponding to the $2\pi$ periodic solution. I tried to use ...
1
vote
0answers
54 views

NDSolve cannot consistently handle a simple 2D reaction-diffusion PDE with spiky (and random) initial conditions

I am still fighting with a modified version of the 2D reaction-diffusion PDE that I have already posted in a previous question ( Reaction-diffusion PDE with NDSolve: either very slow or very ...
1
vote
1answer
119 views
0
votes
1answer
44 views

DSolve “Equation or list of equations expected instead of True in the first \ argument ”

I am getting the error: Equation or list of equations expected instead of True in the first argument when I try to evaluate the cell. ...
1
vote
1answer
48 views

Polar and Parametric Plots not showing up

I bet this has asked before, but I can't seem to find an answer anywhere. I have x and y functions which, individually plotted, are showing up fine, but when I try to plot the parametrically, the ...
6
votes
2answers
159 views

NDSolve and matrix formulation do not produce the same result. Why and which is more precise?

Let's assume a vector of variables $\vec a=(a_1(t),a_2(t),a_3(t))$ and let's assume a following DE$$C\ddot{\vec a}=B\dot{\vec a}+A\vec{a}+\vec F.$$ If $A$ and $C$ are both reversible, than the DE can ...
20
votes
2answers
421 views

Numerical optimal control

I was hoping to tackle optimal control using Mathematica in order to learn how I can use Mathematica's built in numerical integration and optimization functions together in order to solve an optimal ...
-2
votes
1answer
57 views

Nonlinear system of ODEs need phase-plane

I have this system of equations : I need to solve this equation for different parameters, for now theta is 1, alpha is 25 and ...
0
votes
1answer
49 views

Failed to extract parameters by fitting data with ODE

I'm trying to extract 4 parameters by fitting one set of data with an ODE via NonlinearModelFit. I originally thought the process should be straightforward, but later on I found I was totally wrong. ...
0
votes
1answer
58 views

NDSolve Unable to continue with complex values or beyond floating-point exceptions error

Initially tried to solve following first degree of differential equations using following method ...
0
votes
0answers
53 views

Solving system of ordinary differential equations where coefficients are matrices

I am attempting to use NDSolve to solve a system of ordinary differential equations defined by matrices. Here's the setup of my problem and what I have tried so far. I'm interested in applying the ...