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

learn more… | top users | synonyms (3)

1
vote
1answer
151 views

Problem with convergence of limiting value in Manipulate

I created a nice Manipulate box of the magnetic field outside and inside a spherical shell placed in an homogeneous magnetic field. The code is working great, but it has a problem at a limiting ...
1
vote
0answers
62 views

DSolve throws an error about Solve and returns unevaluated [duplicate]

I evaluate: DSolve[{2*y[x] + y'[x]^2 + a*x^b == 0}, y[x], x] During evaluation of In[26]:= Solve::ifun: Inverse functions are being used by Solve, so some ...
0
votes
1answer
127 views

Solving a nonlinear systems of coupled differential equations with boundary conditions

I am trying to solve the following systems of coupled differential equations with boundary conditions (BC) at $0$ and at $∞$ : ...
0
votes
1answer
110 views

Understanding the symbolic solutions of a system of differential equations

I'm trying to solve a simple system of differential equations but get outputs i just cant grasp, they seem way overkill. Someone said it might be because mathematica interprits the constants as ...
0
votes
2answers
93 views

Plotting the solution of a differential equation with its derivative [closed]

The question requires checking of numerical solutions to a differential equation by substituting in those solutions. The original equation is as follows; ...
0
votes
1answer
87 views

ODE system with “parameter” which is a (piecewise) function of one of the dependent variables

I have a 2-dimensional dynamical system where time (t) is the independent variable. My ultimate aim is to perform stability analysis, but I bumped into a problem related to the construction of the ...
1
vote
2answers
58 views

Export NDSolve data

I have calculated a solution of a system of differential equations using NDSolve and generated a plot. How could I extract the data of the solutions? The input is ...
1
vote
1answer
76 views

DSolve to make the path of a ball in a rotating turntable

I have two equations from a analisys of the following situation. A ball starts to move from the center of a Turntable (rotating disk). I tried to solve these equations on Mathematica but the ...
3
votes
0answers
177 views

How to solve the Tsunami model and animate the shallow water wave?

Recently when I was learning differential equations, I noticed there is a shallow water wave equation to model the tsunami propagation. How to establish and solve the initial and boundary conditions ...
3
votes
1answer
216 views

How to solve 3D heat equation with Neumann B.C

I'm trying to solve the 3D heat equation on a cuboid to know if all the perimetric surfaces of a Cuboid achieve the desired temperature of a 873K on deadline time of 2 hours. i have convection on all ...
1
vote
0answers
55 views

NIntegration and DSolve with same variable

I need to solve this problem: Then follows the code, with: ...
7
votes
2answers
333 views

Plotting the image of a curve under a flow

I have some explicit time-independent vector field on the plane, and I would like to study how points evolve under the flow generated by this vector field. The flow is rather complicated and cannot be ...
0
votes
0answers
62 views

Kovacic's algorithm

Kovacic’s algorithm is used to find out whether or not an ODE has a Liouvillian solution, and if it has, then what that solution is. My question is: Is there any code written in Mathematica for ...
0
votes
1answer
43 views

Getting error when trying to solve two coupled differential equations with NDSolve

Mathematica newbie here. I have been trying to use NDSolve to solve two coupled differential equations. Each equation on its own (minus the coupling terms) solves ...
0
votes
0answers
80 views

Coupled PDEs with discontinuous boundary conditions?

So I have a system of coupled PDEs that represent the movement of an interface of fluid. The interfaces are functions of time and space, ie. $h_1(x,t)$, $h_2(x,t)$, $h_3(x,t)$. The PDEs are horrendous ...
4
votes
1answer
81 views

Solving PDE involving Hilbert transform numerically

I'm trying to solve the following equation numerically: $$u_{tt}-\mathcal{H}(u_x)=A^2_{xx},$$ where $\mathcal{H}$ is the Hilbert transform and $A$ is a prescribed forcing function which we assume ...
1
vote
1answer
220 views

Using the Levenberg - Marquardt algorithm to minimize a user-defined function

I wrote a function in Matlab that optimizes another user defined function using lsqnonlin with 'levenberg-marquardt' option. Now, I'd like to use the first to ...
3
votes
1answer
87 views

Cylindrical GL Equations

I am trying to solve numerically the GL equations in cylindrical coordinates. I already know what the shape of the solutions to these equations should look like, but I am not able to get the correct ...
3
votes
0answers
83 views

Measuring Temp on a cubical using heat equation with Neumann and Dirichlet boundary conditions

I'm trying to solve 3D heat equation to understand how the temperature varies with time on cubical. i have mixed boundary conditions: Dirichlet on the bottom surface (constant Temp) Neumann on all the ...
5
votes
2answers
108 views

Why is NDSolve's StartingStepSize with ExplicitEuler not working? How do I set the step size?

Why do I get a smooth curve (that is obviously wrong) when I run this? ...
3
votes
1answer
54 views

How can I manage discrete variables that oscillate when using NDSolve?

I have a system of differential equations with two discrete variables (var1 and var2) that can be either 0 or 1. I'm trying to ...
0
votes
1answer
96 views

Problem with solving an optimal control system with Shooting Method

I try to solve a system of differential equation by shooting method. I looked at questions asked on the site but I could not find that could help me. My functions are ; ...
4
votes
1answer
302 views

Solving the Frenet Serret equations for non-constant curvature and torsion, obtaining parametric equations

I wish to solve for the curvature and torsion functions $k_1 = \dfrac{1}{1+s^2}, k_2 = \dfrac{s}{1+s^2}$ using the Frenet Serret system and obtain the parametric equations for the curve. I need the ...
3
votes
1answer
134 views

NDSolve for optimal control theory

For an optimal control problem, I need to solve a differential equation in which (a) the system is formulated according to a function and (b) initial/terminal conditions are discontinuous. Here is an ...
1
vote
1answer
134 views

NDSolve: Indeterminate expression 0. ComplexInfinity encountered? [closed]

The following code give errors that I don't understand. Could you please suggest a method to solve these differential equations? ...
0
votes
0answers
60 views

NDSolve: problem due to boundary conditions

I am trying to solve a 2nd order ordinary linear differential equation (DE==0). Here is the code: ...
10
votes
0answers
113 views

Differential equations with rational functions as solution

I have some families of nonlinear, first order differential equations. When I try to use DSolve, I usually get a mess (if anything at all) in terms of ...
1
vote
1answer
72 views

Partial differential equation with infinity limit

How do I get Mathematica to solve the following partial differential equation $$ \frac{\partial u(y,t)}{\partial t} = \nu \frac{\partial^2 u(y,t)}{\partial y^2}+\frac{\partial U_0(t)}{\partial t}$$ ...
3
votes
2answers
113 views

Reaction-diffusion PDE with NDSolve: either very slow or very inaccurate

I am struggling to have Mathematica 10.3 solve a system of PDE's (with periodic boundary conditions and random initial conditions), but either it produces a set of very noisy InterpolatingFunction ...
2
votes
1answer
109 views

Solve and plot differential equation

I want to solve this differential equation $\qquad y'(x)=\frac{x-y(x)}{1-y(x)-x}$ and plot it's solution. But DSolve doesn's work. ...
2
votes
1answer
129 views

How to handle a special Neumann-like boundary condition in NDSolve?

How to solve the following nonlinear ODE with two algebraic equations and one boundary condition? ...
1
vote
1answer
72 views

To find a limit cycle with period 1 in NDSolve

I am trying to find a limit cycle with period 1 of an ode system. The idea is to compare the values of the system during two periods, for example, ...
0
votes
1answer
55 views

Multiple stopping constraints in NDSolve

I need to numerically solve several differential equations, with several constraints, like this : ...
2
votes
0answers
160 views

Schrödinger equation for Hydrogen atom [closed]

I'm trying to solve Schrödinger 1D equation for hydrogen atom but I found several difficulties. To get in context I want to solve this equation For Z and l real and arbitraries. To start with I ...
0
votes
1answer
90 views

How to solve the following GR PDE?

Someone could help me for solving the following GR PDE? (this code do not work yet) ...
2
votes
1answer
76 views

How to solve two ODE differential equations, numerically, with four dependent boundary condition

I have two ODE differential equations(DE) and four dependent B.C. $$y_1''(x)=-\frac{1}{2} \left(x-y_1(x)\right)$$ $$y_2''(y)=-\frac{1}{4} \left(y-y_2(y)-5\right)$$ $$y_1(0)=y_2(1)$$ $$y_1(1)=y_2(0)$$ ...
0
votes
0answers
63 views

How to define an exclusion zone for NDSolve

I need to numerically integrate a differential equation and define two exclusion zones to stop the integration. The Mathematica code looks like this : ...
11
votes
3answers
517 views

Using NDSolve to find particle trajectory

I'm trying to simulate a particle in an electric and magnetic fields, but numerically instead of analytically. This is basically solving the equation $$q \cdot \left(p'\times B\right) + q\cdot E = m ...
2
votes
0answers
73 views

NDeigensystem returns error due to mesh discretization when calculating vibrations of a cantilever

There was an transcription error in the code I provided in a previous post: NDeigensystem returns complex non-hermitian error for the calculation of vibrations of a cantilever This resulted in it ...
4
votes
0answers
70 views

WhenEvent and Resetting of Variable in PDE when operation succeeds [closed]

I have had success in using WhenEvent to reset or change a variable within NDSolve with ordinary differential equations. My ...
0
votes
0answers
60 views

Fourier-style solutions to differential equations, not piecewise polynomials

NDSolve returns piecewise polynomials. Is there any way I can get a single (non piecewise) function consisting of sines and cosines instead? Sort of a "Fourier approximation" to the solution of my ...
2
votes
0answers
41 views

Changing from Cartesian to polar coordinates in partial differential equation [duplicate]

I have a partial differential equation: $$\left(x^2+y^2\right)\frac{{{\partial ^2}u(x,y)}}{{\partial {x^2}}} + x^2\frac{{{\partial ^2}u(x,y)}}{{\partial {y^2}}}=0$$ How to change from Cartesian to ...
4
votes
1answer
94 views

System of equations is solved by NDSolve over just a tiny domain

I'm solving numerically a system of differential equations with the use of NDSolve. The numerical integration works only a very small interval of the argument. I'm ...
1
vote
1answer
48 views

Issue plotting PDE solution with manipulate

I am trying to plot the solution to the following PDE with the help of mathematica, however, when trying to employ manipulate to animate the behavior, I find that if I try this: ...
4
votes
0answers
60 views

Numerical instabilities of a convection-(non-)diffusion equation when shrinking from a square to a triangular domain

I am trying to evaluate a parameter-dependent indefinite integral using a PDE-based scheme I described here, and I'm having some trouble when I try and cut down the domain from a square to a triangle. ...
4
votes
1answer
79 views

How can I numerically pre-compute an indefinite integral with a parameter?

Suppose I have a function $f(t)$, and I want to compute its indefinite integral $$F(t)=\int_0^tf(\tau)\mathrm d\tau.$$ Moreover, suppose that, for any of a number of reasons, I require this integral ...
3
votes
1answer
69 views

Split Boundary Value Problems win Algebraic Equations

Is it possible at all to solve with NDSolve (or other built in function) a split boundary value problem with algebraic equations? Please look at the following example: ...
0
votes
0answers
19 views

Solving a simple ODE [duplicate]

I am trying to solve $y(x)' = \sqrt{y(x)}$ with initial condition $y(0)=0$ for $x\in \mathbb{R}$. I have tried this: DSolve[{y'[x] == Sqrt[y[x]], y[0] == 0}, y[x], x] but this gives {{y[x] -> ...
4
votes
2answers
254 views

Rotating an Interpolating Function

I have used the following code to generate eigenfunctions of a PDE: ...
3
votes
0answers
58 views

Weird behaviour for a vector InterpolatingFunction inside an NDSolve

I have run into some weird behaviour on the part of NDSolve which I find pretty bizarre and which I would like to understand better. Suppose, for the sake of ...