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

learn more… | top users | synonyms (3)

1
vote
2answers
217 views

ParametricNDSolve and plotting problem

I have a problem with ParametricNDSolve: ...
4
votes
1answer
569 views

Combining Gravity Turn and Orbit Models

I have a mathematical model for the motion of an orbiting spacecraft about Earth: ...
0
votes
2answers
81 views

Derivative in function form

I have a function in mathematica defined as b[n_, x_] := x^n/n!/\!\(\*UnderoverscriptBox[\(\[Sum]\), \(k = n\), \(\[Infinity]\)]\*FractionBox[\(x^k\), \(k!\)]\) ...
1
vote
0answers
77 views

Updating procedural fitting algorithm to more efficient style?

For some time now, I have used this procedural programming approach to fit data in my research that is too cumbersome for the built in Mathematica functions. (Although this simple example works fine ...
0
votes
0answers
27 views

Identification of a parameter in a first order differential equation [duplicate]

I wan't to identify a parameter in a first oder differential equation. My goal is proof of the calculation principle i use. First I calculate a response where the paramater is set to 0.3. Later this ...
1
vote
1answer
71 views

NDSolve with WhenEvent, resetting system for a prolonged period

I am currently working on a complex system where I would like to (for the lack of a better description) reset (part) of the system for a certain period. Say I five ODE's: ...
1
vote
0answers
71 views

A Stupid Question Maybe - EigenNDSolve

I would like to find out how accurate and how useful the chebyshev (collocation?) method is in finding many eigenvalues to a second-order ODE in one go. Specifically, I used the Mathematica package ...
5
votes
2answers
129 views

How to add (energy) constraint when using NDSolve to Equation of Motion

To simplify my problem, I will try and solve the Equation of Motion for a particle in a 1D Harmonic Potential. energy[x_, p_, m_, ω_] := p^2/(2 m) + (m ω^2)/2 x^2 ...
1
vote
1answer
91 views

Fit to inverse function

I am trying to fit my data to an inverse function from solving a DGL: ...
0
votes
1answer
64 views

Why does DSolve fail here?

I have the following differential equation: a*L'[t] == 40000*L[t]^(-3/2) - 1.4 - extf I want to fit the solution to data in order to acquire a. I tried DSolve ...
0
votes
1answer
115 views

Tough Calculation, novice mathematica user

I have an equation, that I've been calling $b_N(x)$ that satisfies the following identity: $$-Nb_N(x)^2=(x-N)b_N(x)+xb_N'(x)$$ where $b_N'(x)$ is the first derivative. I take the derivative then of ...
0
votes
0answers
122 views

NDSolve error message: Cannot find starting value for the variable y'

I'm attempting to solve a differential equation using NDSolve. This is the code: ...
2
votes
1answer
161 views

Analytic approximation of NDSOLVE

I have to solve the following differential equation L'[t] == L[t]^(-3/2) - 0.1 I tried DSolve: ...
0
votes
1answer
202 views

StreamPlot in a system x´ = f[x] + g[t] [closed]

StreamPlot can be used in the O.D.E system: x' = 2 x - 3 y + 3 t y' = 5 x + y - t How? ...
6
votes
1answer
414 views

Poisson PDE over a irregular region with FDM

The problem modeled by the Poisson PDE is related to the torsion of prismatic beams and I use the Finite Differences Method (FDM). I've managed to solve the equation over a rectangle region with ...
1
vote
1answer
97 views

Series expansion of InterpolatingFunction obtained from NDSolve

I am trying to obtain a series expansion of the numerical solution of a differential equation. I encounter difficulties going beyond first-order expansions which I believe might be due to my inability ...
0
votes
1answer
119 views

Problem with DSolve [closed]

I am trying to use DSolve: DSolve[{y'[x] == x/y, y'[8] == 2}, y[x], x] But I get an error message: DSolve::dvnoarg: The ...
0
votes
0answers
56 views

How to tell NDSolve to ignore small values in choosing step size

I have a very large system of first order differential equations, which can be written (on paper) as dX/dt=F(X), where X is a vector and F is a vector function. All elements of X are strictly positive ...
0
votes
2answers
783 views

recursive depth of 1024 exceeded [closed]

Can some tell me what is wrong in the following code. ...
1
vote
1answer
236 views

Solving system of differential equations using loops

I have $F$ system of differential equations. Out of those $F$ equations except for first and last I have general form for the remaining equations (say $ dP_{i}/dt)$. Let $dP_0/dt,dP_F/dt$ denotes the ...
1
vote
1answer
76 views

NDSolve with a constant

I have a simple differential equation like this s = NDSolve[{y'[x] == A y[x] Cos[x + y[x]], y[0] == 1}, y, {x, 0, 30}] where $A$ is a constant. Mathematica gives ...
6
votes
1answer
236 views

Second order differential equation with boundary conditions solving repeatedly

I want to solve a second order differential equation in the interval[-1:1], which does not have a analytic solution, \begin{eqnarray} y''(x) &=& k \phi^2(x)y(x) \\ \phi(x) &=& ...
2
votes
0answers
59 views

Can't find the limit of this complicated expression

Here is the limit I am trying to calculate ...
2
votes
3answers
127 views

Solving a differential equation with NDSolve and plotting it WITH the use of Manipulate

I want to solve, let's say, this differential equation: y'[x]=y[x]^n Given any initial condition (which I might want to manipulate too). I want to plot the ...
0
votes
0answers
159 views

How to solve system of non linear Differential Equations with 1 independent variable and 3 dependent variables?

I have been trying to solve this system of differential equations using mathematica. $$F(t) B''(t)+B'(t) F'(t)+B(t) F''(t)=0,\\ F(t) A''(t)+A'(t) F'(t)+A(t) F''(t)=0, \\ B(t) A'(t)+A(t) B'(t)=0$$ Here ...
0
votes
1answer
82 views

How to evaluate a “pure function”? [closed]

I need to evaluate the result of DSolve, e.g. {{y -> Function[{x}, x^2 - x^3]}} at some specific point. I tried y[1] and y /. x -> 1. Hope someone knows how to do it
1
vote
2answers
172 views

Solving Differential Algebraic Equations as BVP

I am trying to solve a set of DAEs. \begin{equation} -4 \nu (\lambda(s))^{(-1 - 4 \nu)} \theta'(s) \lambda'(s) + (\lambda(s))^{(-4 \nu)} \theta''(s) = -\alpha_y \cos\theta(s) + \alpha_x ...
0
votes
1answer
179 views

NDSolveValue with a multiple variable equation

I'm trying to plot this equation of motion for a pendulum with a periodically applied torque for 20 cycles, where b, m, L, g are given constant values. ...
7
votes
1answer
255 views

Problem with Neumann condition in quarter disc

Bug introduced in 10.0 and fixed in 10.2 So I'm following the available examples in version 10 for FEM, The plane stress operator is shown as this ...
0
votes
1answer
168 views

How can I stop the integration of NDSolve with a condition?

I try to solve a partial differential equation by NDSolve. At some point, I want to stop the integrating by a condition that compares the min value of the function ...
3
votes
3answers
273 views

Higher-order, nonlinear differential equation with Initial Values

I tried to solve for an non-Hookean spring's motion, but the output from Mathematica is weird. It seems that there is inverse functions involved. ...
2
votes
2answers
121 views

DSolve returning { } when given my ODE

Why does this code return {}? ...
2
votes
1answer
200 views

Problems with NDSolve and partial differential equations of several variables

Suppose we have the following partial differential equation: $$ 0 = \frac{ \partial w }{ \partial \tau } + \left( w + \sqrt{ h + \beta } \right) \frac{ \partial h }{ \partial \chi } $$ where $w$ = ...
1
vote
0answers
106 views

Numerical solution to two non-linear coupled differential equations [closed]

I am trying to solve two differential equations representing the position of an object in space. I have specified arbitrary initial conditions. ...
4
votes
1answer
240 views

Using NDSolve to solve a system of coupled PDEs

I am trying to solve the Gross-Neveu model in one dimension for a specific soliton initial condition. I am trying ...
5
votes
2answers
346 views

Simplifying general solutions of differential equations (driven harmonic oscillator)

Solving general differential equations in Mathematica usually leads to somewhat unsightly results. As an example, consider the solution of the driven, damped harmonic oscillator: ...
1
vote
0answers
114 views

PDE with Integral constraint

I am trying to solve the Non-linear Schrodinger equation $-\Delta \psi(r) + \psi(r) - |\psi(r)|^2\psi(r) = 0$ where $r \in \Omega$ In a square domain ($(x,y) \in \Omega$ where $\Omega=[0,1]\times ...
2
votes
1answer
95 views

StateResponse is non-deterministic

I observed non-deterministic behaviour in StateResponse. Let's look at an example. ...
5
votes
1answer
128 views

Drawing disk with coordinates from NDSolve

I have simple Manipulation expression, but continuously getting error Coordinate {x[0.], y[0.]} should be a pair of numbers, or a Scaled or Offset form. I assume that something is wrong with ...
0
votes
0answers
77 views

Why NDSolve With Orthogonal-Projection Method On Orr-Sommerfeld Equation Does Not Work(?)

I am attempting to solve the Orr-Sommerfeld equation for plane Poiseuille flow with the Orthogonal Projection method within NDSolve. The Orr-Sommerfeld equation is (a "stiff" problem); $\psi''''(x) ...
1
vote
1answer
33 views

Error following a variable in an ODE [closed]

I am trying to solve a 4th order differential equation and this is my code thus far ...
0
votes
0answers
76 views

Boundary condition is not specified on a single edge of the boundary of the computational domain

When I want to solve this PDE I get this error "NDSolve::bcedge: "Boundary condition u[0,0.05]==0 is not specified on a single edge of the boundary of the computational domain"" ...
0
votes
0answers
51 views

Resources about Geodynamics

I'm about to start a small project on numerical geodynamic. If possible, I whish to handle the computational work in Mathematica. The project will be small but probably not so small I can do all by ...
0
votes
0answers
384 views

Poincare Section of an Hamiltonian

I'm in desire to plot the Poincare Section of a differential equation defined by a hamiltonian system. The hamiltonian is as follow: ...
17
votes
1answer
701 views

Has this implementation of FDM touched the speed limit of Mathematica?

Still, I'll use the implementation of the 1D FDTD method (you can simply understand it as a kind of explicit finite difference scheme for the Maxwell's equation) as the example. Just for completeness, ...
0
votes
1answer
121 views

Attempt to solve for Differential Equation - Acc, Vel, Pos in (x,y,z) [closed]

In the problem I'm working on, there is a tennis ball that is subjected to a headwind, tailwind, and crosswind. I am trying to use NDSolve in order solve for the position functions and to eventually ...
0
votes
1answer
84 views

Formulating a second boundary condition to get an alternative solution to a ODE [closed]

I have the ODE $c'(t) = t^2c^3$ with the initial condition $c(1) = 20$. The differential equation $c'(t) = t^2c^3$ (without boundary conditions) has two branch solutions. I want to formulate a new ...
1
vote
0answers
115 views

Taking a derivative, then squaring a two dimensional numerical function

Context I am interested in first finding an interpolating function of the solution to the linearly damped wave equation. Here is the solution to the LDWE with smooth square inital function : ...
11
votes
2answers
492 views

How to create subregions for the NDSolve FEM Solver

I am trying to create a 2d region consisting of two subregions. The inner region has several holes, where boundary conditions are applied. The figure shows the idea. I have tried to create this ...
0
votes
0answers
218 views

Using NDSolve to solve a PDE with a Dirac Delta function

I'm trying to solve the following equation (in Mathematica 10): ...