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

learn more… | top users | synonyms (3)

2
votes
3answers
642 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. ...
4
votes
1answer
322 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: ...
2
votes
1answer
414 views

Different Poincaré maps for same Equation

I've seen the Wolfram Poincaré example. I'd like to modify this sample in order to have a list of Poincaré maps. In the example the event is captured only when ...
1
vote
1answer
343 views

Why my differential equations become True? [duplicate]

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 ...
1
vote
1answer
266 views

Get Cos[2 x] in “Solve”

My goal is to get an analytical solution of x ( $\lambda$ in the problem), when $\epsilon$ (f[x] in the problem) is maximum (use inverse form). I found that if I can use trigonometric identity ...
0
votes
0answers
65 views

Derive a PDE from other PDEs

I'm trying to do some symbolic calculation like this: Substitute (11.2) into the last two of equations (11.1) and, in turn, substitute these into the first of equations(11.1), I want to ...
0
votes
0answers
175 views

NDSolve: Couple ODE and PDE

I am trying to numerically solve an equation with NDSolve, where there is a ODE coupled to a PDE, like the following: ...
0
votes
1answer
256 views

How to work with Interpolating Function?

I'm starting using Mathematica. I have two coupled differential equations and after using NDSolve the output is in the form of a Interpolating Function. I want to Plot the square modulus of one of the ...
0
votes
1answer
78 views

Unable to configure a system of differential equations

My differential equations have been written as ...
1
vote
1answer
256 views

predator prey model with prey refuge

I want to plot the phase diagram of prey predator versus prey refuge to see how the prey refuge influences the population of prey and predator.And this is the system $x'=\alpha ...
5
votes
1answer
309 views

Solving a nonlinear, second order d.e. in Mathematica

I am trying to solve the following differential equation $$ s''(x) = \frac{8}{9}s(x) -2s^2(x) + s^3(x), $$ with boundary conditions $s(-30) = 4/3 $ and $s(30)=0$, in Mathematica by using the following ...
4
votes
3answers
739 views

Constraining function found by NDSolve to stay positive

I am trying to constrain my variables in NDSolve so that they stay positive for the integration interval. Is is possible to do that? Using Assumptions as in ...
0
votes
0answers
90 views

Solving systems of equations and derivatives

I am new to Mathematica, and I'm not sure how to solve a system of equations that include differentiation. Given that f and g are differentiable functions: ...
1
vote
1answer
193 views

Solution of Coupled Second Order ODEs with Boundary Conditions

D[Tu1[x], {x, 2}] qu - Tu1[x] (fu0 + I w p) == Tou[x] f1 D[T1[x], {x, 2}] q - T1[x] (f0 - b g0 + I w p)== To[x] (f1 - b g1) - g1 For T1 , ...
0
votes
0answers
294 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 ...
3
votes
1answer
637 views

Solving a differential equation with initial conditions

I want to get a solution to an ordinary differential equation with imposed initial conditions. I'm confused about why my following code doesn't work: ...
1
vote
1answer
916 views

Non-numerical value for a derivative on NDSolve function

I am trying to solve the following system of differential equations: ...
0
votes
1answer
236 views

DSolve does not give solution

No answer is given when I evaluate the following expression: ...
2
votes
1answer
240 views

Runge-Kutta-2 on System

After spending some time using the Mathematica documentation and this Mathematica.SE answer, I implemented the Runge-Kutta-2 routines. I am hoping someone can validate what I did and tell me that it ...
1
vote
1answer
219 views

Solution of an ODE in implicit form

For some non-linear ODEs there is only implicit form of solution using DSolve. For example DSolve[(y[x] + x - 1)*y'[x] - y[x] + 2 x + 3 == 0, y[x], x] gives ...
3
votes
1answer
270 views

Polar Differential Equation

I have a polar differential equation (this is the reduced variant) as: $$r' = 0, \theta' = 1$$ I figured out (from previous answers) that I can nicely convert this to Cartesian and use StreamPlot ...
1
vote
1answer
69 views

Variable location for intial value in ParametericNDSolve

I am trying to figure out if Mathematica allows a variable location for a known initial value for ParametricNDSolve. I haven't been able to find an example where ...
2
votes
0answers
89 views

Locate Blow-up in NDSolve with Whenevent

At some point I get this error: ...
0
votes
1answer
119 views

NDSolve accuracy related problem

I want to solve this equation: ...
1
vote
0answers
209 views

BVP system of nonlinear coupled ODEs

Here I am, trying to solve this system of coupled ODEs (up to a minus sign): $ u''=6u^5-(8+4a)u^3+(2+4a)u+\frac{2u((w^2-s)^2+bw)}{(u^2+c)^2}$ $ w''=\frac{4w^3-4bw+b}{u^2+c}$ with the boundary ...
1
vote
1answer
204 views

Discrete sampling of interpolating function returned by NDSolve

When solving an ODE with NDSolve, Mathematica returns an interpolation function. I need a discrete sampling of this function however. Naively, I can write this as (example): ...
0
votes
1answer
134 views

Replacement rules in combination with pure functions to make a change of variables

I have defined my function in this way w = f[r,Θ] ; After some calculations i obtained my results with respect to the previous function. For instance: ...
5
votes
1answer
485 views
0
votes
1answer
457 views

Heat Transfer equation by numerical methods

I want to solve the following heat conduction equation using numerical methods: D[u[x, t], t] -alpha*D[u[x, t], {x, 2}] == 0 u[x, 0] == 1/(1 + x^2)^0.25, u[-10, t] == u[10, t] == 0, ...
1
vote
0answers
47 views

How to find discretezation error of NDSolve

Is there a way to find out how large the truncation, round-off, and other errors that occur from discretizing a differential equation are while using the default settings in NDSolve? Or would I have ...
2
votes
2answers
576 views

Trying to solve a differential equation with a piecewise initial condition

I am trying to solve $$u_t=\frac{1}{4}u_{xx}$$ $-\infty<x<\infty,\: t>0$ With the initial condition $u(x,0)=\phi(x)$ where $$ \phi(x)= \left\{ \begin{array}{lr} 1 & ...
1
vote
0answers
141 views

ExpToTrig transforms solution to 4th order ODE into unwanted form

Mathematica gives the solution of the second order differential equation DSolve[a y''[x] + b*y[x] == 0, y[x], x] in trigonometric form ...
-2
votes
2answers
540 views

How do we solve third order nonlinear differential equation f’’’+ff’-f’2-Re2 f’=0

How do we solve third order nonlinear differential equation f’’’+ff’-f’2-Re2 f’=0 f(0)=0,f'(0)=1,f(Infinity)=0 From the OP's comment: ...
0
votes
1answer
127 views

NDSolve in exact points - nonlinear equation (Degrees and Radians confusion)

I found solutions of two nonlinear equations, but I have three small questions. First I am in confusion because of initial value. The system is mechanical and unknown ϕ is an angle, which I set to ...
1
vote
1answer
143 views

Numerically solving an ODE for a parameter

My equation is as follows: sol = ParametricNDSolve[{(f'[r]^2 - 1) f'[r] r == 6.2 (f'[1])^2/1000, f[1] == a}, {f}, {r, 1, 3}, {a}] where the function ...
0
votes
0answers
278 views

Boundary Value Problem- using NDSolve or another method

I am trying to solve a set of coupled partial differential equations, with defined boundary conditions using mathematica. Here are the equations and the boundary conditions. ...
2
votes
2answers
269 views

WhenEvent for several variables

I want to solve a system of nonlinear two differential equations (say, in $\theta(x,t)$ and $\phi(x,t)$ with NDSolve, but I want to stop the evaluation when one of ...
4
votes
2answers
478 views

How to transform transfer functions into differential equations?

is there a way with Mathematica to transform transferfunctions (Laplace) into differential equations? Let's say I have the transfer function $\frac{Y(s)}{U(s)}=\text{Kp} \left(\frac{1}{s ...
1
vote
2answers
176 views

NDsolve with variable end-point

I want to plot the solution a coupled ODEs as a function of the end point only. Mathematica code: ...
1
vote
1answer
77 views

Invariant error plot for an arbitray system of ODE

i try to use InvariantErrorPlot command for a generic system of ODE but i find only example with pre-building equations from Mathematica Packages. This is one of ...
2
votes
0answers
93 views

Are there some other ways to solve a second PDE except DSolve?

I have a partial differential equation as follows: $$\frac{\partial p(x,t)}{\partial t}=\text{Dp} \frac{\partial ^2p(x,t)}{\partial x^2}-\frac{p(x,t)-\text{p0}}{\tau }$$ What I try to do was to get ...
0
votes
1answer
468 views

NDSolve example - analytical vs numerical solution. How to specify initial conditions?

So, I have the following equation: y[t]^-3 == w^2 y[t]/t^2 + y''[t] where w^2 = (t0 k0)^2. Analytic solution is easy to ...
1
vote
1answer
924 views

Euler's method for system of differential equation

I need to program Euler's method to solve a system of two diffferential equations of first order. Fist, I have programmed the Euler's method for just one differential equation: ...
2
votes
1answer
200 views

The output is the same as Input when I add the boundary conditions to a PDE

I am trying to solve a PDE in the first order with specific boundary conditions. When I solve use DSolve without the boundary conditions, Mathematica gives me an answer in an arbitrary function. When ...
2
votes
2answers
564 views

How to visualize slope fields of differential equations without vectors?

I'm looking to visualize slope fields of differential equations for my differential equations course. Every example I see draws them as vectors, adding unnecessary "arrows" that, to me, are visually ...
0
votes
1answer
264 views

Initial and boundary value errors

NDSolve::ivone: Boundary values may only be specified for one independent variable. Initial values may only be specified at one value of the other independent variable. >> I keep getting this, ...
0
votes
1answer
133 views

Plotting the discrete solution to a PDE

So I went through all the grunt work of solving a PDE in discrete form for a research project. Now I have a Temperature solution as a function of time and space: $$ T_{i}^{n+1}=\left( \dfrac ...
0
votes
0answers
294 views

How to obtain differential equations of motion using Lagrangian dynamics?

I'm confused to find a set of differential equations of motion of a pair of masses, m1 and m2 joined by a spring of constant ...
0
votes
1answer
52 views

How do I pull from a data list for parameter values in a system of ODEs and then solve and plot?

I have a parameter set list and I would like to solve a system of ODEs multiple times using each parameter set in the list. Then I would like each set of solutions to be graphed on the same plot. ...
0
votes
0answers
146 views

Boundary Condition Problem with Green's Function

My temperature distribution looks like: Mean wire temperature is : To[x_] := 2 c1 Cosh[x a1] + a2 And it is continuous at two points. In positive part: ...