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

learn more… | top users | synonyms

0
votes
1answer
50 views

How to find the value of a solution of an ODE at a particular point [closed]

I'm trying to find the value of $z(y)$ at a particular time, say $y=0.5$. I've solved the differential equation. I tried using z3 = z[3] /. sol, but it's not ...
2
votes
0answers
48 views

EquationTrekker-like behavior for state space?

EquationTrekker is great for phase space plots, however I want to plot the results of $$\phi '(t)=-b \sin (\phi (t))+g \sin (\Phi (t)-\phi (t))+1\\\Phi '(t)=g y ...
3
votes
1answer
152 views

Issue in ParallelTable after evaluating another function using NDSolve and FindRoot. What is wrong with this inverse?

I am trying to find the inverse of a function which is defined through NDSolve and NIntegrate. The question is pretty similar ...
3
votes
1answer
92 views

Already solved DE, now I need to rearrange and plot

I have solved for $z(w,x,y)$ in a differential equation: $$ 3\frac{\partial z}{\partial y} = 2(z-1) + (1-wy^2 )x $$ And I obtained the general solution: $z = f(w,x,y)$ Now putting in $x=0$, we have ...
1
vote
1answer
47 views

Solved the differential equation, how to rearrange to find solution? [closed]

I have solved for values of $y(x)$ from $1 < x < 2$ and plotted them: ...
0
votes
2answers
67 views

Solve simple differential equation - Error?

I get the error "The function Te appears with no arguments." when running this code. I'm not sure why. Is it because of the $Abs[\Gamma]$ ? ...
5
votes
2answers
221 views

how to solve ODE with boundary at infinity

y''[x]-x y[x]==0 y[0]==AiryAi[0], y[infinity]==0 the analytic solution to this ODE is the Airy function y[x]=AiryAi[x] if I ...
0
votes
2answers
171 views

Plotting a solution of the Allee Effect

I am doing a project on the Allee Effect. I am able to successfully create a stability analysis. That is, I can find the relevant equilibrium points ($y_e = 0,\ y_e = \alpha$ and $y_e=k$) and draw ...
0
votes
1answer
109 views

DSolve problem with system of linear ODEs

I encounter a rather strange problem in Mathematica when trying to solve the following system of linear differential equations: ...
1
vote
0answers
148 views

Problem with NDSolve in Mathematica 9 / 10

I'm having trouble by solving the following differential equation in Mathematica 9 and 10, where the code works fine in version 7: ...
1
vote
1answer
99 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
80 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, ...
3
votes
1answer
97 views

Is it possible to obtain explicit symbolic solutions to such linear ordinary differential equations?

The ordinary differential equations to solve have symbolic parameters $k_1,k_2,k_3,k_4,k_5,k_6 \in \mathbb{R}$. $$ \left\{ \begin{array}{l} {y_1}'(t)=-{k_1} {y_1}(t)-{k_2} {y_1}(t),\\ {y_2}'(t)={k_2} ...
1
vote
1answer
211 views

Solving a nonlinear PDE with Mathematica10 FEM Solver

I am trying to solve a system of coupled nonlinear PDEs in a rectangular region with the new FEM solver in Mathematica 10. However, I come across an error stating NDSolveValue::femnonlinear: ...
1
vote
1answer
107 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: ...
2
votes
1answer
171 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 ...
1
vote
1answer
37 views

WhenEvent “StopIntegration” problem

The first example on the WhenEvent help page is ...
0
votes
0answers
38 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 ...
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~ ...
1
vote
2answers
229 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
86 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 ...
1
vote
2answers
75 views
1
vote
0answers
33 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 ...
35
votes
2answers
1k 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 ...
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
53 views

NDSolve solution methods issue

I´m modelling a biological reaction thorugh a set of differential equations. ...
4
votes
2answers
140 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: ...
5
votes
3answers
163 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: ...
6
votes
0answers
115 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
39 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 ...
4
votes
1answer
119 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
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
57 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 ...
25
votes
3answers
826 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
0answers
44 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
94 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
60 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
204 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
77 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
87 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
64 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
67 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
1answer
136 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 ...
0
votes
2answers
120 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
69 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 ...
0
votes
0answers
56 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 ...