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

learn more… | top users | synonyms (3)

4
votes
2answers
96 views

Problem with ParametricNDSolve

i have a strange problem with ParametricNDSolve. Given the sistem of differential euqations ...
0
votes
1answer
35 views

Numerical error in DSolve [duplicate]

I'm solving the 1D Poisson equation and am getting some funky numerical artifacts in the solution. ...
1
vote
1answer
37 views

DSolve with Absolute Value

When playing around with DSolve, I enter the following: DSolve[{Derivative[1][y][x] == x + Abs[x]*y[x]}, y[x], x] The output is: $$y(x)\to e^{\frac{1}{2} (x ...
0
votes
1answer
82 views

Method that combine LinearSolve and Dsolve for a set of differential equations?

I have a matrix for example: $\begin{bmatrix} 0 & 5 & 0 & 3 & 3\\ 5 & 5 & 5 & 3 & 1\\ 0 & 3 & 0 & 3 & 10\\ 1 & 3 & 2 & 3 & 10 \\ 1 ...
1
vote
1answer
78 views

System of differential equation with delay

I try to solve this differential system with delay: ...
0
votes
0answers
42 views

System of ODEs with Parameters

FYI: I am new to Mathematica so my questions are probably pretty amateur. I am attempting to write a code to numerically solve the system of equations: {dp/dt= ap-bpx {dx/dt= spx-yx where t is the ...
0
votes
2answers
138 views

Lotka Volterra Equation Plotting

So I have a basic Mathematica assignment where we need to numerically solve a Lotka-Volterra predator-prey population equation. So far I have managed to come up with this ...
3
votes
5answers
272 views

Taylor series without expanding factorial in denominator

A Taylor series is produced with the following code: ...
0
votes
0answers
46 views

Problem of WhenEvent in using variable updated in Event

equation1: value[inside RotationTransform] is a constant ...
2
votes
1answer
129 views

How can I invoke the solution of NDSolve to determine a parameter in my equation just inside NDSlove?

I am trying to solve a differential equation by NDSlove for $h(x,t)$. It reads $$h_t=h_{xx}-V_h-\lambda(t)$$ where $V_h$ is a given function of $h(x,t)$ denoted by ...
3
votes
1answer
112 views

What boundary is added when NDSolve::bcart pops up?

When insufficient boundary conditions are given to NDSolve for solving PDE, the warning NDSolve::bcart pops up: ...
2
votes
1answer
66 views

Solving PDE with Condition on Inside of Region

I'm trying to figure out how to use NDSolve to numerically solve PDEs where some internal values may be specified on the region of the solution. An example, which I picked from a textbook, is solving ...
0
votes
1answer
45 views

What function is the best for solving differential equations? [closed]

dy/dt = -y^2, what is the best function to solve this equation?
3
votes
1answer
121 views

NDSolve grid refinement for PDEs

I am experiencing trouble when trying to solve a PDE for various parameter values within mathematica: The PDE in principle can easily be solved numerically, however the stepsize of the spatial ...
1
vote
1answer
122 views

Integral Curve Question using Contour Plot

The ODE $$(xy-2)+(x^2-xy)y'=0$$ is NOT an exact differential equation. However, if I multiply both sides by an integrating factor $\mu=1/x$, the resulting equation ...
0
votes
0answers
51 views

DSolve gives back input (coupled ODEs)

I'm trying to solve a system of coupled ODEs, however I keep getting back the input whenever I try to solve with DSolve. The variables I'm solving for are complex. Here's what I have : ...
1
vote
2answers
43 views

Filling a region defined by a reduce command

I'd like to shade the region defined in the xy-plane produced by this code: ...
1
vote
0answers
118 views

Simulate a polygon bouncing with collision detection inside a circle/polygon

Simulate a polygon bouncing with collision detection inside a circle/polygon. the example can be see in help page by searching WhenEvent; ...
0
votes
1answer
62 views

PDE raises NDSolve::ntdvdae, then kernel quits

In a related question, the SolveDelayed->True option seemed to solve the problem. SolveDelayed is not a valid option in M10 ...
9
votes
2answers
215 views

Solve a PDE over a region defined by a Bezier patch

I am using NDSolve to find the solution to a PDE over an arbitrary domain. The domain is specified by a Bezier patch. ...
-1
votes
1answer
139 views

Small mismatch in theoretical vs MMA solutions of a damped harmonic oscillator differential equation

I solved the differential equation $$ \ddot y + 0.43 \dot y + y =0.1 \cos(2x), \quad y(0)=1,\ \dot y(0)=1 \quad(0) $$ using both the theoretical results from wikipedia and using MMA (I also tested ...
1
vote
1answer
71 views

How can I simplify the result of a differential equation by using my own defined expressions? [duplicate]

I want to solve the following differential equation with mathematica: α β * w''''''[ξ] + ( 1 + α - p ) * w''''[ξ] + p/β * w''[ξ] = 0 The answer seems at ...
1
vote
1answer
151 views

Solving a nasty partial differential equation

I have a differential equation that I would like to solve numerically in the region $z \in [0,L]$ and $t \in [0,t_{max}]$: $$ \partial_t S(z,t) = f(z)S(z,t) + \int_0^L \text{d} z'g(z,z') S(z',t), $$ ...
1
vote
1answer
47 views

How to get the raw data of solutions of partial differential equation in mathematica

I am using NDSolve command of mathematica solving a very complicated partial differential equation (two independents "t" and "x", three master variable "u", "v", "w"). Each time I solve it, it takes ...
0
votes
1answer
64 views

Third-order ODE with NDSolve

I would like to use NDSolve to draw the integral curve of a third-order ODE ...
0
votes
1answer
116 views

NDSolve with initial condition: initial conditions did not evaluate to an array of numbers of depth 1(for 1D) or 2 (for 2D) on spatial grid

I have the following naive code to solve a PDE in two spatial dimensions (x,y) with periodic boundary conditions: ...
1
vote
1answer
77 views

DSolve for Second Order Differential

I have an equation y''[t] + w^2*Sin[y[t]] == 0 So I'm using DSolve like this: ...
0
votes
1answer
56 views

Simplifying an answer returned by DSolveValue

Consider the IVP: $$x'=\frac92-\frac{x}{300+2t},\quad x(0)=0$$ Working the solution by hand, the answer is: $$x(t)=450+3t-\frac{4500\sqrt3}{\sqrt{300+2t}}$$ Entering code and then attempting to ...
3
votes
0answers
73 views
4
votes
5answers
256 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: ...
1
vote
3answers
101 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: ...
5
votes
2answers
187 views
0
votes
0answers
70 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). ...
2
votes
0answers
70 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
42 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
149 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: } ...
0
votes
0answers
40 views
2
votes
2answers
62 views

Strange behavior with Show command

I enter this: ...
1
vote
1answer
84 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
229 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 ...
3
votes
1answer
166 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 ...
6
votes
1answer
132 views

Long running ToElementMesh with very “large” domains

I'm trying to solve a system of PDE over a large domain. This doesn't means I need to have a huge amount of mesh points and mesh elements to discretize the domain. Just that the domains has a big ...
4
votes
0answers
53 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: ...
26
votes
3answers
638 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 ...
1
vote
1answer
101 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
98 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
1answer
50 views

EllipticF reduction

Is it possible to reduce the following Z to Legendre form Elliptic integral of the First kind? $$ \int \dfrac {\sec u\; du } {\sqrt{ (1- {(\nu \tan u)}^2 }} ..(1*)$$ With ...
0
votes
0answers
28 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
83 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
1answer
165 views

Runge-Kutta Butcher tables

I like to have the Butcher's table for Explicit (or implicit as well) Runge-Kutta method of a fixed order. I do not understand reading ...