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

learn more… | top users | synonyms

1
vote
1answer
45 views
0
votes
1answer
167 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. ...
3
votes
1answer
132 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 ...
0
votes
0answers
31 views

Customizing display of partial differential equations

I am manipulating partial differential equations symbolically, and would like to get the easily readable form $\rho \frac{\partial v}{\partial t}$, leaving variables implicit. Based on suggestions ...
0
votes
0answers
14 views

How can I calculate matrix differentiation? [on hold]

I am studying about the Matrix Differentiation. I don't know if this red box differential metric, which is how it is calculated.
1
vote
1answer
86 views

Need help solving a system of two 1st order nonlinear differential equations

The original system of equations reads: $\begin{cases} f'(r) + f(r) \left(a(r) - \frac{1}{r}\right) = 0,\\ f^2(r) + a'(r) + \frac{a(r)}{r} - 1 = 0\,, \end{cases}$ with boundary conditions $f(0) = ...
3
votes
1answer
253 views

Transform recursion for coefficients into differential equation for generating function

Assume, one is given a linear recursion with polynomial coefficients for a sequence $(a_i)_i$, such as a[i] == i a[i-1] I would like to convert this recursion ...
16
votes
4answers
728 views

How to splice together several instances of InterpolatingFunction?

I have a set of InterpolatingFunction returned by NDSolve which are valid over different (but overall continuous) domains. How ...
0
votes
0answers
37 views

DSolve solves PDE only without boundary condition, but fails otherwise

I'm trying to solve for Te[w, Pprobe, t] in a partial differential equation. What's surprising is that it manages to solve it when I don't put in any initial ...
0
votes
1answer
48 views

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

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 ...
0
votes
1answer
137 views

Solving System of Nonlinear with Three Differential Equations

I've been trying to solve a system of nonlinear differential equation, but the conditions are a bit weird. Two of the differentials equate to the same equation, but have different boundary ...
0
votes
1answer
133 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 ...
2
votes
0answers
44 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
86 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 ...
0
votes
1answer
131 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 ...
1
vote
1answer
167 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: ...
34
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 ...
1
vote
1answer
41 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: ...
-1
votes
0answers
68 views

Can Mathematica solve my system of differential equations? [closed]

I have a system of differential equations about two unknown functions $f(x)$ and $g(x)$: $$\phi_1(x)=\alpha_1 g(x) \frac{d}{dx}f\left(x+\frac{g(x)}{g'(x)}\right)$$ $$\phi_2(x)=\alpha_2 [1-f(x)] ...
0
votes
2answers
64 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]$ ? ...
2
votes
1answer
139 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
152 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
97 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
140 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: ...
3
votes
1answer
96 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} ...
9
votes
1answer
1k views

How to solve a system of partial differential equations?

Edit: since the upgrade to Mathematica 10, this problem seems solved I just want to solve a system of partial differential equations, for example: $$ \left\{ \begin{array}{l} ...
1
vote
1answer
103 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
87 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
73 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
160 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
109 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] ...
24
votes
3answers
800 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
537 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
34 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
159 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
102 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
59 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
224 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
84 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
139 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
179 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
24 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
30 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
45 views

NDSolve solution methods issue

I´m modelling a biological reaction thorugh a set of differential equations. ...
5
votes
3answers
149 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
127 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
370 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
103 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 ...