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

learn more… | top users | synonyms

1
vote
1answer
105 views

When plotting vector fields what is the difference between the following?

I am a new Mathematica user, and I have run into a question. What is the difference between the output of the commands: ...
1
vote
1answer
90 views

Using spherical derivatives with NDSolve

I am trying to plot the path of an object around a centre of mass. Obviously this needs spherical coordinates. The problem I encounter is with the derivatives. The Derivatives in spherical coordinates ...
5
votes
2answers
78 views

Automatically detect largest interval over which NDSolve can find a solution

Question: Consider the following numerical resolution: NDSolve[eqn, {x1[t], x2[t], y[t]}, {t, tmin, tmax} where eqn, ...
2
votes
1answer
163 views

Solve ODE $d^2u/dx^2 + u/A = 0$

How can I solve following ODE with Mathematical: $$d^2u/dx^2 + u/A = 0 \quad (\text{or } \ C),$$ with the conditions: $$ \left.\frac{d^2u}{dx^2}\right|_{x=0} = 0, $$ $$u(x=0) = B$$ and ...
4
votes
1answer
117 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] ...
25
votes
3answers
815 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 ...
1
vote
1answer
543 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? ...
1
vote
1answer
36 views

WhenEvent “StopIntegration” problem

The first example on the WhenEvent help page is ...
0
votes
0answers
36 views

NDSolve with Table and Boundary Conditions

I am trying to implement a 1D model of coupled field equations where I split up space into pieces to get a system of coupled ODEs which I call Eeqns and Ieqns. The problem I can't figure out is how to ...
0
votes
1answer
161 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 ...
1
vote
0answers
103 views

Use Runga-Kutta method to solve a system of many ODEs [closed]

For my simulation purpose, I need to solve the following system of DEs using 4th-order Runga-Kutta (RK) method in Mathematica: $$ I_{k,0}'(t) = -\beta_1~I_{k,0}(t) + k~\lambda~ \Theta(t) + \gamma~ ...
0
votes
1answer
60 views

Solution of a set of ODEs

How can I write code to solve a set of ODEs in two-coordimate (i,j) system? Example: ...
1
vote
2answers
228 views

How do I plot the differential of this function?

Suppose I have this function: $$z(x,y) = \left| \frac{ \frac{1}{3x +iy} -2x}{iy + \frac{1}{x}} \right| $$ I want the contour plot of $\frac{\partial z}{\partial x}$ with axes $(x,y)$. Tried this ...
2
votes
0answers
85 views

Edge finite elements in Mathematica [closed]

I am wondering if it's possible to extend the Finite Element capabilities of Mathematica so as to include other types of elements such as edge vector elements (Nedelec) which can be used in ...
2
votes
1answer
142 views

Understanding of method for NDSolve

I used automatic method for NDSolve. Then I asked myself - which method Mathemathica prefered? I got an answer for this question on this forum, that I need to use: ...
1
vote
2answers
72 views
3
votes
1answer
182 views

Indexing of Large Autonomous System of Equations for Use in NDSolve

I currently have a list of ~150 autonomous DAE's (differential-algebraic equations) that I need to be able to work with using NDSolve. Fortunately I already have ...
1
vote
0answers
28 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
31 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
50 views

NDSolve solution methods issue

I´m modelling a biological reaction thorugh a set of differential equations. ...
5
votes
3answers
151 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
130 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: ...
2
votes
1answer
372 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
109 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
53 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 ...
0
votes
0answers
67 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
55 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
43 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
88 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
0answers
58 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
53 views

solution of differential equation with complex roots [closed]

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 ...
8
votes
3answers
198 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 ...
4
votes
1answer
75 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
75 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
60 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
49 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; ...
2
votes
1answer
65 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
68 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
147 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
53 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
106 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
115 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
54 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
287 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 ...
12
votes
4answers
829 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
29 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
59 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: ...