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

learn more… | top users | synonyms

0
votes
0answers
25 views
20
votes
2answers
401 views

Numerically solving Helmholtz equation in 2D for arbitrary shapes

I would like to solve the Helmholtz equation with dirichlet boundary conditions in 2 dimensions for an arbitrary shape. (for a qualitative comparison of the eigenstates to periodic orbits in the ...
1
vote
0answers
21 views

How to fit data with numerical solution of system of parametric ODE? [duplicate]

I need to find the parameters (k1,k2,k3,k4,k1r,k2r,k3r,k4r) that fit my data (list of [Intensity, time]) using the function ...
0
votes
0answers
29 views

How can Mathematica solve 7 equations and want to find 7 variables [duplicate]

I have 7 equations as below * e1 := (1 - a) p Subscript[A, 0] Subscript[s, K ] ! *SubsuperscriptBox[(s), (L), (-a)] = (1 - b) Subscript[B, 0] (1 - Subscript[s, K]) (1 - Subscript[s, ...
-1
votes
1answer
40 views

NDSolve solution methods issue

I´m modelling a biological reaction thorugh a set of differential equations. ...
5
votes
3answers
142 views

Can NDSolve handle discountinuos data?

It is possible to numerically solve a differential equation if not-smooth data are involved? For example the following instruction return the error NDSolve::bvdisc: ...
4
votes
2answers
116 views

NDSolve fails to solve vector-valued function x'[t] == a . x[t] + u . b

I get no solution when I run the following command: ...
21
votes
3answers
719 views

Symbolic solution(s) to generalized Heat equation

Symbolic solution(s) to Heat equation? or more generally,(eventually) Green functions to known PDEs I am interested in variations of the heat equation: or more generally or even more generally ...
2
votes
1answer
362 views

How can I solve my system of differential equations?

I have a question pertaining to my code below: I require to find the values of the roots E2 that will set the function ...
6
votes
0answers
98 views

Partial Differential Equation in Parallel

is there any native way to implement multi-core parallel solving of PDE in Wolfram Mathematica? WM 10 now supports Finite Elements Method, but it is actually useless without parallelization. Usually ...
2
votes
1answer
45 views

NDSolve: EventAction list specification using Table fails

I have an NDSolve problem that needs a list of a large number of events that are generated programmatically. Here is a simple example that demonstrates the problem on Mathematica 7 and 8 (the versions ...
0
votes
0answers
38 views

Complex infinities in an Euler-Lagrange equation

I now wish to solve for $L(t)$ and $H(t)$ some ugly Lagrangian (in the code below) in configuration space, respectively the length and the height of a system. I tried running the code below, with the ...
-1
votes
0answers
94 views

Poincare section for double pendulum [closed]

I would like to reproduce the problem here, for a specified initial condition. I almost have copied the code. But I do not know why it does not work: ...
4
votes
0answers
88 views

1/0 encountered when solving an ODE

What can cause this error to show up? Clear[lambda, a, b, x, y]; ode = lambda y[x] + a b y'[x] + a (-1 + b x) y''[x] + y''''[x] == 0; DSolve[ode, y[x], x] ...
0
votes
0answers
57 views

Plotting Geodesics in Kerr

I'm interested in plotting the trajectories of null geodesics near a rotating black hole (given by the Kerr solution) which involves a system of first order differential equations. Some Context (not ...
3
votes
0answers
52 views

Error solving third order ODE in version 10. No error in V9. DSolve`DSolveKovacicDump`

This is an ode which is solved with no error in V9, but gives a very strange error in V10. Is this a regression bug? On version 9, the ode is solved with no error ...
2
votes
0answers
39 views

DSolve breaks when the ordering of independent variables aren't proper?

I encountered this when trying to solve this problem with DSolve: ...
1
vote
1answer
77 views

Use NDSolve to solve a PDE

I have the following PDE (a master equation, and $P$ is probability density, $0\le x\le1$ and $0\le y\le1$): $$ \partial_t P(x,y,t)=x\partial_xP(x,y,t)+(1-y)\partial_yP(x,y,t)+2P(x,y,t) $$ The ...
0
votes
1answer
143 views

DSolve will not apply assumption `m ∈ Integers`

I am trying to solve a linear second order ODE using DSolve which involves an arbitrary integer m. ...
0
votes
0answers
54 views

Trouble outputting Plot

I am trying to solve a system of ODEs and plot them with Manipulate but I keep getting blank outputs. Below is my code: ...
1
vote
0answers
44 views

solution of differential equation with complex roots

How to find the solution of differential equation with complex roots For example of [{y^2 n[t]+ 2 k y n'[t] + n''[t]= M p sin(o t), n[0]==0, n'[0]==0}, n,t] given k lies between 0 to 1 which makes the ...
6
votes
3answers
177 views

Should DSolve always return solution with constant of integration?

Version 10.0 Clear[y,x]; DSolve[D[y[x], x] - y[x]^2 + y[x]*Sin[x] - Cos[x] == 0, y[x], x, GeneratedParameters -> C] or ...
0
votes
1answer
84 views

how to solve second order nonlinear coupled differential equations using NDSolve with hyperbolic function

i have to solve some solitons scattering through this coupled equations. i need to get two different graph, but still the graph did not come out. and also the equations quite complicated containing ...
3
votes
1answer
64 views

Factoring a differential equation into an operator product

Following the fundamental theorem of algebra, I can (as stated here in Sec. 1.1.4) factor an $n$th-order linear ordinary differential equation $$ a_n \frac{d^n x}{dt^n} + a_{n-1} \frac{d^{n-1} ...
1
vote
0answers
53 views

Speeding up NDSolve for system of differential equations

I am wondering if there is a way to speed up this function that solves a system of ordinary differential equations with NDSolve? Thus far I've tried specifying a few different methods such as LSODA, ...
1
vote
1answer
50 views

Numerical solution NDSolve for SIR model [closed]

I'm trying to numerically solve the SIR model. This should be simple but I'm can't understand why this is not working: ...
2
votes
0answers
39 views

Apply IC and BCs on Second - Order Linear PDE

I am now trying to solve second - order linear partial differential equation in that interested eq. has been separated into variables to simplify the procedure for Mathematica. Here is my equation; ...
0
votes
1answer
102 views

How do I subtract two contours?

Suppose I have this contour described by the equation the root $z$ of this equation $ \frac{1}{x^2 + y^2} + \frac{1}{xz} = 2y $ Now suppose the equation is tweaked slightly, with an addition of $w$, ...
2
votes
1answer
62 views

Why DSolve giving Inconsistent or redundant transcendental equation on this problem

Kamke's differential equation #574 has a solution, but Mathematica generates an error message as it tries to solve it. It still gives the solution. My question is: What could have caused Mathematica ...
0
votes
0answers
25 views

How DSolve with a multisymbol function argument

I used DSolve on a differential equation e got the following result: You may notice that C[ 1 ] is a function of the variables r and z in a defined form. If I substitute this result in my ...
0
votes
0answers
45 views

Poincaré Section

I have encountered somewhat the same problem as here. But with the equations, $x'(t) = p(t), p'(t) = - x(t) - y(t), y'(t) = q(t), q'(t) = - y(t) - x(t)$ My code is, ...
0
votes
0answers
63 views

Optimizing a functional using variational calculus

Here I am, trying to get the trajectory $(x(t),y(t))$ that will minimize the following Lagrangian (i.e. the integrand of the functional) between $(-1,-1)$ and $(-1,y_{eff})$ where $x_{eff}$ is defined ...
3
votes
2answers
130 views

Problem with Poincare section

these are equations for some double oscillators: $x'(t)=p, p'(t)=-x-3y, y'(t)=q, q'(t)=-y-3x$ I would like to plot the Poincare section for collection of $(y,q)$ when $x=0$ and $x'(t)>0$ (or for ...
0
votes
0answers
47 views

Help to plot Poincare section for a system of delay differential equations

I'am trying to plot Poincare section for a 4-dimentional system of Delay differential equations. I don't know how to specify initial conditions. Can I use each of the code in the following ? Also how ...
0
votes
2answers
92 views

Solve differential equation with 3 variables and plotting contour

Suppose we have this equation: $$ 2 - x g(z) f \left(x, \frac{\partial y}{\partial x}\right) = 3y $$ using initial condition $y = 2$ where $g = \left| \frac{3}{2-iz} \right|$ and $f = ...
2
votes
1answer
104 views

Simple differential equation to solve

Suppose we have this differential equation: $$ x^2 + y^2 + z^2 = \frac{\frac{\partial y}{\partial x}}{x+y+z}$$ I want to find $\frac{\partial y}{\partial x}$ at x=1,y=1,z=1. I tried this code but ...
0
votes
0answers
45 views

Energy Equation: Boundary and initial condition

I'm trying to solve the energy equation and NDSolve is able to solve the problem. Even though the result is acceptable, I have always the warning "boundary and initial conditions are inconsistent". ...
7
votes
2answers
252 views

Nonrectangular region for NDSolve

I have a PDE with mixed boundaries (Neumann and Dirichlet on some sides) in the region $(t,x,y) \in \left( 0, T\right) \times\left\{ -L \leq x \leq L, 0 \leq y \leq h(x) \right\}$ where $h(x)$ is ...
11
votes
4answers
779 views

How do I solve a PDE with a strange boundary condition?

How do I solve the PDE with boundary value like this $$u(t,x,y)=0, \textrm{when } F(x,y)=0$$ using DSolve? As a specific example, I want to solve heat equation $$\frac{\partial u}{\partial ...
12
votes
1answer
1k views

Numerically solving an inhomogeneous partial differential equation

I'm trying to solve a cylindrical partial differential equation with boundary conditions. But I got an error message saying ...
0
votes
0answers
25 views

Ndsolve, transmission line transients

I am trying to use NDSolve to analyze transients in transmission lines. I started with this simple one: ...
0
votes
2answers
58 views

Plot integration curve of a system of ODEs

I'm really n00b in Mathematica, so please bear with me, as this seems to be my only option to learn how to do what I wany to do. I have a system of two differential equations: ...
2
votes
3answers
190 views

Implicitly differentiate an equation, then solve the resulting equation

Suppose I have an extremely tedious equation to differentiate and want mathematica to help do the differentiation and solve. Consider a less tedious equation: $$y (x,z) = sin \left(\frac{1}{x} ...
1
vote
1answer
520 views

Integro-differential eqn with double integral

I am looking at the following variation of a integro-differential, with y[0]=1. The output is not great, any solutions to this? ...
7
votes
0answers
2k views

Integro-differential equation

I have to numerically solve a nonlinear partial integro-differential equation using Mathematica. This is my equation, $$\frac{\partial y(x,t)}{\partial t}=\int_{-\infty}^\infty K_0(|x-u|) ...
0
votes
1answer
133 views

Numerical solution of the hyperbolic equation

I am trying to solve the following hyperbolic equation with given boundary conditions: I choose as initial condition $u=1$, and evolve the above hyperbolic equation until reaching a stationary ...
0
votes
0answers
22 views

Pass model parameters to NDSolve'Reinitialize

I am wondering if there is a way to pass parameter values during the NDSolve'Reinitialize step of the NDSolve process? I've read that one way to possibly speed up repeated calls to NDSolve is to ...
4
votes
1answer
146 views

Using NDSolve to solve Equation of Motion in cylindrical coordinates

I have a set of coupled differential equations which represents the equation of motion of a particle in cylindrical coordinates with the following Hamiltonian: $$ H=\frac{1}{2m} \left( p_r^2 + ...
0
votes
1answer
143 views

Problem with Plot and ParametricNDSolve

This is my first question on this site, I hope it is written in an understandable way. Let's start. My aim is to plot a certain parametric region, with $x$ and $y$ coordinates given by the following ...
0
votes
0answers
97 views

Im trying to model balls in a box with NDsolve, but the balls keep escaping the box

This Code produces a box filled with balls at different positions with random charges on them and then calculates their motion according to Coulombs law. rendering the electric field as well as the ...