Tagged Questions

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

learn more… | top users | synonyms

1
vote
1answer
366 views

Non-numerical value for a derivative on NDSolve function

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

DSolve does not give solution

No answer is given when I evaluate the following expression: ...
1
vote
0answers
141 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
146 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
133 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
66 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
70 views

Locate Blow-up in NDSolve with Whenevent

At some point I get this error: ...
0
votes
0answers
67 views

fitting system of ODEs - “part2 of {} does not exist” errod

I'm trying to fit my data to a system of ODEs: x'[t]==a*x[t]+b*y[t] y'[t]==c*y[t]+d*x[t]+e*z[t] z'[t]==f*z[t]+g*y[t] having t,x[t],y[t],z[t]. I followed the code ...
0
votes
1answer
77 views

NDSolve accuracy related problem

I want to solve this equation: ...
1
vote
0answers
126 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
101 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
114 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: ...
3
votes
1answer
210 views
0
votes
1answer
283 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
36 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
1answer
248 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
96 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
305 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
0answers
69 views

Simplify DSolve solutions (all roots are known)

I have problem how to present the solution of the system of two differential equations in the form as in the case of one differential equation if I know the nature of roots. For example ...
0
votes
1answer
111 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
117 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
65 views

Shooting method problem

I am trying to solve the BVP problem: $ u''-u'=-6u^5+(8+4a)u^3-(2+4a)u, u(0.001)=1, u(10)=-1, u'(0.001)=-0.001$ ...
0
votes
0answers
131 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
155 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
292 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 ...
0
votes
1answer
116 views

NDSolve Initial Condition Returned as True [closed]

I am attempting to use an NDSolve solution as an initial condition for another NDSolve as suggested in this post: Using NDSolve solution as initial condition for another NDSolve ; however, Mathematica ...
1
vote
2answers
134 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
60 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
78 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
0answers
106 views

Question about NDSolve

Here is my code: ...
0
votes
1answer
269 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 ...
2
votes
1answer
111 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
468 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
152 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
0answers
205 views

NDSolve::ndsz:At t==0.000039576, step size is effectively zero

I adjusted my code accordinly to the advice I got. However I still have a problem with this code which now looks like this: ...
0
votes
1answer
98 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
198 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
40 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
111 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: ...
0
votes
0answers
42 views

How do you explicitly define Box Operator (D'Alembert operator)? [duplicate]

I wanted to see how one can define a box operator in mathematica? $$ \Box\phi=(-\frac{\partial^{2}}{\partial t^{2}}+\frac{\partial^{2}}{\partial x^{2}}+\frac{\partial^{2}}{\partial ...
0
votes
0answers
186 views

Shooting method for solving 3rd Oder ODE with RK method

This is my 3rd Order IVP. D[y[x], {x, 3}] == (m/2)*((D[y[x], x])^2 - 1) - ((m + 1)/2)* y*(D[y[x], {x, 2}])^2. Initial conditions given are as follows ...
0
votes
3answers
311 views

Find an equations of motion using Euler Eqution [closed]

I am a completely new to Mathematica, so my question and code may look silly. I need to find a equation of motion of a point particle of mass m moving in a potential V(x) = A abs[x]^n. (n=1,2,3,4,5,6. ...
3
votes
1answer
123 views

NDSolve with varying PrecisionGoal and WorkingPrecision

Sometimes we need higher numerical precision to deal with large number cancellation in an equation. But if this cancellation happens only in a small (and known) parameter space, would it be possible ...
6
votes
2answers
581 views

StreamPlot for Bifurcation Diagram

When we do a StreamPlot, I want to show the bifurcation when $a = 0$ transitions to $a >0$, but do not see a better way to do this than the following. ...
2
votes
0answers
106 views

Internal Shooting Method of NDSolve in combination with NDSolve`Reinitialize?

To explain my problem, I am trying to extend the BVP problem example from the help that illustrates how to use the shooting method of NDSolve: ...
10
votes
3answers
850 views

Plot Matlab icon

I started to explore this on a whim and hasn't succeeded yet… Some introduction for the icon is found here but I can't understand it very well. (I admit that, though playing with ...
0
votes
1answer
113 views

Setting plot size within vector; dynamic slider (two part question)

I have a question about the size of a plot in mathematica. I'm working on some differential equations (a simplified model of heat transfer in a building). I've programmed it so that you can set ...
0
votes
0answers
32 views

How to Invert a Multivariate Interpolating Function? [duplicate]

I used NSolve on a partial differential equation with boundary conditions to get a numerical solution that Mathematica says is of the form InterpolatingFunction. I then defined a function f(x,y) which ...
0
votes
1answer
162 views

Infinite expression 1/0.^(1/4) encountered in NDSolve

I am trying to solve very simple system of coupled differential equations using the NDSolve function, but I am getting ...
0
votes
0answers
204 views

How to apply Robin boundary condition in Mathematica

Robin condition is the boundary condition which combines the Dirichlet and Neumann condition at the same boundary. For example: At the boundary x = 0 to 2: I want to apply ...