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

learn more… | top users | synonyms

4
votes
5answers
129 views

Plot all curves in a table with a single color

I want all of my solutions in a table to be the same color. My attempt: ...
0
votes
3answers
49 views

Trouble using Manipulate and NDsolve for coupled ODE's

I'm aware that there are a few threads on the topic of solving coupled ODEs using NSolve and Manipulate. Based on those, I wrote the following code, solving my system of equations: ...
1
vote
0answers
48 views

System of differential equations [on hold]

I have to solve the following system of two differential equations with certain parameters inside which I just put them together for simplicity and called them t: ...
0
votes
0answers
38 views

Selecting values over lists for discrete inputs in NDSolve

I'm currently trying to simulate a step response in a simple tank system with NDSolve. I would like it to start at say time 2 (the precise time is not important). ...
4
votes
1answer
232 views

Differentiating ParametricNDSolve solutions

Is there any way to differentiate a solution obtained by ParametricNDSolve? For instance, I have the position $\phi_\gamma(t)$ as a function of time, parametrized ...
3
votes
1answer
130 views

Boundary sphere partial differential equation

I am trying to solve partial differential equation in spherical coordinates $(\theta,\phi)$,but I don't know how to properly include boundary conditions of $\theta$. For $\phi$ it is periodic, but ...
2
votes
0answers
47 views

NDSolve is running an extremely long time: how can I save the existing data?

I am trying to solve a PDE by the following code. It takes 1 hour or so to reach t=25.72 but about 20 hours to reach t=25.72404031638060174049337306126853310997. Actually, the time step is extremely ...
1
vote
1answer
35 views

Variables within functions within DSolve

Hi I am just beginning to learn Mathematica and this is my first time I have been exposed to any type of coding. I am encountering a problem in a basic physics problem. For example we are always ...
2
votes
2answers
125 views

Too high differential order in boundary conditions?

The following boundary value problem has a unique solution: $$ \begin{cases} x''-x=0\\ x''(0)=0\\ x(1)=1 \end{cases}\\ \text{General solution: } x=c_1 \sinh{t}+c_2\cosh{t}\\ \text{Solution for BVP: } ...
2
votes
2answers
54 views
1
vote
1answer
74 views

Preventing complex term in a differential equation solution

I am solving the initial value problem: $$y'=\frac{t+1}{t(t+4)},\quad y(-1)=0$$ I believe the correct answer is $$y=\frac14\left[\ln(-t)+3\ln(4+t)-3\ln 3\right]$$ with an interval of existence ...
3
votes
1answer
114 views

Fick's Law over Implicit Regions

I'm trying to solve 2D Fick's law of diffusion from the boundaries of a triangle (future shapes will include more complex implicit regions). I'm able to model diffusion without time in complex ...
3
votes
3answers
104 views

Differential equation with a Heavisidetheta function running for ever

So I'm trying to get the solution of the following differential equation: ...
20
votes
2answers
373 views

Can Mathematica solve Plateau's problem (finding a minimal surface with specified boundary)?

Given a closed curve $\mathcal C$ in three dimensions, is it possible to use Mathematica's built-in functionality to compute a minimal surface whose boundary is $\mathcal C$? For simplicity, let us ...
0
votes
1answer
181 views

Help solving the following non-linear differential equation

After spending several hours exploring different options, I still couldn't not figure out how to solve the following differential equations. Any suggestions would be highly appreciated. ...
1
vote
1answer
81 views

RSolve works slow and does not give an analytical solution for solving high order linear difference equations

RSolve@@{{a[n]-a[n-10]==101}~Join~Table[a[i+1]==101+10i,{i,0,9}],a[n],n} I want to get the exact form of general formulas, but ...
4
votes
0answers
41 views

Extraneous reaped data with NDSolve and WhenEvent

In exploring What is the range of values $x$ for which $f(x)$ is higher than $k$ over a given domain?, I came across the following strange behavior: ...
1
vote
1answer
69 views

Initial Value Problem with initial conditions as closed region

Does Mathematica have functionality to solve an Initial Value Problem when the initial conditions consist of a closed region (rectangle) instead of a point? That is, instead of initial conditions as ...
1
vote
2answers
64 views

Plotting solutions to ODE for continuously varying parameters

I am trying to plot the solution of an ODE with continuously varying one of its three parameters. (I just fixed two of them for simplicity.) Any suggestions? ...
0
votes
0answers
27 views

What means the definition F[t_?NumericQ] :=…? [duplicate]

I've come along a function defined in this way: F[t_?NumericQ] := NDSolve[{-u''[x] - u[x]/x == -t^2 u[x], u[x0] == x0, u'[x0] == 1}, u, {x, x0, x1}] What does ...
3
votes
2answers
66 views

Saving a plot as a pdf to use in Latex

I use the following code: sol = DSolveValue[{y'[t] == 1/(2 y[t] + 3), y[0] == 1}, y, t]; Plot[sol[t], {t, -10, 10}, ImageSize -> Small] Then I Ctrl+click the ...
1
vote
2answers
114 views

Solving a second order non-linear differential equation

I am trying to solve the following equation DSolve[{u''[t] + 4 u[t] + 0.1 u[t]^3 u'[t] ^2 == 0, u[0] == 1, u'[0] == 1}, u, t] Unfortunately Mathematica is ...
14
votes
2answers
182 views

Least effort to handle a point source inside the domain of PDE(s)

By point source I mean a constrained condition at one point inside the domain of PDE(s). For example: $$\frac{\partial ^2u(t,x,y)}{\partial t^2}=\frac{\partial ^2u(t,x,y)}{\partial ...
2
votes
1answer
96 views

PDE: Solving Burgers' equation with initial value given by a self consistency equation

I would like to solve in Mathematica the well known inviscid Burgers' equation \begin{align} \begin{cases} u_t(x,t) + u(x,t)u_x(x,t)= 0 \\ u(x,0) = m(x) \end{cases} \end{align} where $$m(x) ...
2
votes
1answer
100 views

DSolve not returning a rule

I am trying to solve a PDE, but the DSolve is not returning a rule. How exactly can I get Mathematica to yield a solution? ...
8
votes
1answer
315 views

Modelling the effect of a structure on a “tsunami” (hyperbolic wave equation)

So, the hyperbolic wave equation can be quite easily solved in Mathematica like this: ...
0
votes
2answers
71 views

NDSolve ensemble of initial points

I want to illustrate how the differential equation depends on the initial conditions. First if all differential equations: ...
0
votes
0answers
56 views

Flux continuity with NDSolve

When NDSolve is applied to model mass diffusion (or heat flow) through different materials, is the continuity of flux condition automatically satisfied at the ...
1
vote
2answers
216 views

Trying to model Heat flow trough different materials with NDsolve

What I'm trying to achieve is model of the heat flow, in this case for the simplest 1D case,its relatively easy to do for the steady state case, but when I try to do it with NDsolve so I get the ...
7
votes
1answer
167 views

mathematica 10 not showing numerical solution of differential equations?

I just got the new mathematica version 10 and tried to solve the following system of differential equation. $$r^2\frac{d^2f}{dr^2} = 2f(1-f)(1-2f)+\frac{r^2}{4}h^2(f-1)$$ ...
-1
votes
1answer
88 views
2
votes
2answers
367 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 & ...
-2
votes
1answer
73 views

Ordinary Differential Equation Question [closed]

I was wondering whether someone might be able to point me in the right direction for this ODE that I am trying to solve? It is: $$\frac{du}{dx} = \frac{x}{s}-\frac{1}{2}-\frac{u}{x}$$ with $s$ ...
0
votes
0answers
51 views

Solving for a coordinate transformation for a rank 2 tensor

I am working on a project in general relativity, and I have a metric or rank 2 tensor $g_{ab}$ of the form, $$g =\left( \begin{array}{ccccc} f(x^a) & 0 & 0 & 0 & 0\\ 0 & -f(x^a) ...
6
votes
2answers
129 views

How to use WhenEvent with a vector ODE in NDSolve

I have an ODE system I'd like to specify as a vector equation in NDSolve. I'm not clear on how to use WhenEvent for a system ...
3
votes
1answer
191 views

Using WhenEvent for derivative of discontinuous function

I have a discontinuous function ($u(t)$, a square wave) and I would like WhenEvent to trigger when the signal goes high/low, i.e. when the value of $u(t)$ changes. ...
4
votes
2answers
2k views

How can I solve this PDE with initial conditions?

I wish to solve the following second order partial differential equation $$u_{tt}=c^2 u_{xx}$$ with the initial conditions $$u(x,0)=e^x,\quad u_t(x,0)=\sin x$$ In Mathematica, I evaluated the ...
2
votes
0answers
129 views

How to control pedestrians to enter the simulation region with random functions?

I am working on simulating crowds of a cross-typed region, given that pedestrians entering from 4 gates of north, east, west and south, ...
0
votes
0answers
81 views

Solving a system of ODEs with DSolve

I am trying to solve a system of ODEs using DSolve ...
1
vote
1answer
660 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? ...
8
votes
1answer
495 views

NDSolve::ndcf: Repeated convergence test failure. How to solve?

I am trying to simulate a system of $n$ pendulums with some friction in Mathematica 9. This is the code I am using: ...
0
votes
2answers
49 views

trouble with adding a legend to a plot

I am having trouble with adding a legend to a plot. this is my code,I tryed to use serval things,like ...
1
vote
2answers
171 views

NDSolve with Explicit, Implicit Euler and Trapezoidal method

I am using Mathematica 9.0 both on Linux and Windows and I would like to integrate the Van der Pol equation numerically using various techniques such as Explicit and Implicit Euler and Trapezoidal ...
7
votes
3answers
241 views
1
vote
1answer
122 views

Plot multiple solutions of a system of differential equations with different values of constant c by using NDsolve

I am traying to solve system of differential equations with different values of constant c then plot the solution with different values of the constant c. where c= 2.2758, 2.04822, 1.82064, 1.59306, ...
0
votes
0answers
64 views

Issue in NDSolve and unknown error

I try to solve a simple BVP: $A {v}''(t)+B v(t)=f(t)$ where $A, B$ — $5\times 5$ matrices, $f(t)$ — 5-dimentional vector-function with the use of NDSolve. Each component of $f(t)$ is a combination ...
2
votes
0answers
77 views

DSolve versus DSolveValue [closed]

I've just discovered the new command DSolveValue in Mathematica 10. Is this new command now the preferred instead of DSolve? Is it the better command to teach to students? Same question for ...