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

learn more… | top users | synonyms (3)

2
votes
0answers
92 views

How to include the radiation from a surface in a heat flow PDE? [closed]

I have a PDE were the temperature distribution T[t,x,y] is to be determined. The region is a cross section of a (glass) plate. On one side I heat using a Neumann condition, on the opposite side I want ...
0
votes
1answer
637 views

What is the meaning of the error message “Set::write: Tag Plus in … is Protected”? [closed]

i'm trying to solve system of 9 coupled linear differential equations.I've this problem ... what is it's meaning?
0
votes
1answer
150 views

Iterating FindRoot to solve a differential equation

I am new to Mathematica style programming, better at traditional programming. I would like to know if I can create a table of results using FindRoot, where the ...
4
votes
1answer
145 views

Help nedeed for understanding NDSolve architecture -Vacuum Heating Simulation

I am a beginner in mathematica and I am facing a problem that is way bigger than my knowledge.. I am tryng to simulate a semplificated version of Brunel Effect, wich is a Laser-Plasma interaction ...
1
vote
0answers
56 views

Getting ALL solutions of the logistic ODE with DSolve?

In Mathematica 10.0, the command DSolve[y'[x] == y[x]*(1 - y[x]), y[x], x] gives {{y[x] -> E^x/(E^x + E^C[1])}} I take ...
3
votes
2answers
139 views

Coupled differential equations from Lagrangian — singularity being reached very quickly

I was wondering if there was a way to help the solutions to these differential equations not hit a singularity as quickly as they currently do. The equations are very complex (derived from a ...
7
votes
2answers
309 views

Solve a time-dependent matrix DE using NDSolve

I am trying to numerically solve a time-dependent matrix DE using Mathematica. However, it seems that Mathematica freezes even for two-dimensional matrix equations. For example, we have: ...
3
votes
0answers
96 views

Using NIntegrate and DiscretePlot to visualize pseudodifferential operators

In harmonic analysis, pseudodifferential operators are a way to generalize the notions of derivatives, through the use of Fourier transforms. The basic idea being, Let ...
1
vote
2answers
138 views

How to determine the parameter in a differential equation with some fixed solution

Assuming that I have a ODE system with undetermined parameter $$x''(t) == y(t) x(t)$$ $$y'(t) == 2 - a x(t)$$ and I have some fixed solution condition $$x(0)=0$$ $$x(10)=8$$ $$y(10)=3.5$$ Is there ...
2
votes
1answer
84 views

Using ouput of ParametericNDSolve as argument in NMinimize: NMinimize::nnum: function value not a number

I am trying to fit some experimental data and I am trying to use NMinimize function to calculate the parameters that fit the experimental dataset. I am also using ...
4
votes
1answer
171 views

EquationTrekker using Mathematica 10.2

I wonder whether the problem that arise with the following code is specific to my machine, as it seem to work with the documentation example. ...
3
votes
1answer
76 views

NDSolve: obtain time difference when different solutions get same distance to equilibrium value

I have a numerically solved a differential equation which has an equilibrium value and a starting point, both included as a parameter. I want to plot the difference between the times it takes to come ...
2
votes
1answer
131 views

Decouple 3 simultaneous 2nd order ODEs (of 3 dependent variables)

I have three simultaneous ODEs with coupled dependent variables. I was wondering whether there's a way for Mathematica to decouple these equations (into three ODEs of a single dependent variable ...
0
votes
0answers
190 views

How to solve the heat flow PDE for the constant source term

I have the two dimensional heat flow equation solved for a two dimensional reagion over time. That is an instance I found in the internet. What I am looking for is a method how to determine the heat ...
0
votes
2answers
169 views

How to solve this equation with NDSolve?

Here is the code: ...
0
votes
0answers
83 views

Piecewise NDSolveValue: Nested Variable Iteration

I am attempting to make the variable Rc a piecewise-function of [z] and have it update as NDsolve iterates through the range of 'z' such as {z,0,Le), but to no avail. The following code is an attempt ...
2
votes
1answer
86 views

Simplifying the result of DSolve

I have a system of differential equations that can be solved with DSolve. However, the coefficients are very long there are like terms which can be combined. I used //FullSimplify in my code but the ...
0
votes
1answer
141 views

Using NDSolve initial conditions as parameters for a function and then a replacement rule

I am trying to make a function based on two initial conditions of a set of differential equations: ...
12
votes
1answer
188 views

“ParametricSensitivity” in ParametricNDSolve

"ParametricSensitivity" is listed as a Method in the documentation for ParametricNDSolve, ...
1
vote
0answers
59 views

Solving an ODE using shooting method [duplicate]

I have been trying to solve the following nonlinear ordinary differential equation: $$-\Phi''-\frac{3}{r}\Phi'+\Phi-\frac{3}{2}\Phi^{2}+\frac{\alpha}{2}\Phi^{3}=0$$ with boundary conditions ...
2
votes
0answers
59 views

Solving a differential equation with a complex (independent) varible

I'm trying to plug in complex values in the numerical solutions of ODEs without success. For instance y[I] /. NDSolve[{y'[x] == 0, y[0] == 3}, y, {x, -10, 10}] ...
-1
votes
1answer
128 views

Problem with NDSolve and If

If T=T1, program can calculate. But, when T = T1*If[Abs[q2[t] - q2d] < e1 && Abs[q2'[t]] < e1, 0, 1], program can ...
0
votes
1answer
102 views

DSolve Error 3D Navier Stokes (z,r,t)

I am essentially wanting to solve Navier Stokes in the [z,r,t] dimensions (2D+Unsteady Flow) for pressure driven flow. However, I keep getting an error I am unable to resolve in my Dsolve such as: ...
2
votes
1answer
124 views

ODE of an eccentric circle

How can we simplify the solution in: DSolve[ {Sqrt[ 1 - (r[t] - T^2/r[t])^2/(2 a)^2]/(r[t] - T^2/r[t])/(2 a) == r'[t]/ r[t], r[0] == T}, r, t] (r,t) are polar ...
1
vote
2answers
72 views

Condition on Inside of Region

How do I handle boundary conditions present in the region. Mathematica complains about the condition v [1 , t ] == u [1,t]. The condition should be given to this ...
4
votes
3answers
600 views

Nonlinear differential equation: numerical solution

I have to find and plot a numerical solution for this second order differential equation: u''[x] + (u'[x]/x) - (u[x]/(x^2)) + u[x] - u[x]^3 = 0 where $0\leq x ...
1
vote
1answer
78 views

How do I restrict the range of an interpolating function in NDSolve?

I have a numerical solution for a differential equation. I know that the gray dotted solution is unphysical. Is there a way I can tell mathematica to ignore the gray dotted solution (e.g. restrict the ...
0
votes
1answer
129 views

How to solve time dependent coupled differential equations? [closed]

I am trying to solve time dependent coupled differential equations with know initial conditions. I am just a beginner in terms of coding so would really appreciate if someone can give me an example of ...
1
vote
1answer
111 views

Getting the message NDSolve::underdet:

I am solving the wave equation with initial position and velocity function: ...
0
votes
1answer
69 views

NDSolve does not satisfy initial conditions

For the problem I'm facing, I have three coupled diff eqns that I'm trying to solve ...
0
votes
2answers
238 views

How to solve a pair of coupled 2nd order equations using shooting method

I have a pair of equations for x[t] and y[t]. Both are 2nd order and they are coupled. Now say that I know the following ...
0
votes
0answers
39 views

Boundary condition [duplicate]

I'm having the following issue that can not find a solution , by using this boundary condition v[0, t] == u[1, t] that Mathematica does not understand, someone ...
0
votes
1answer
183 views

Improving NDSolve computation speed for 3D stationary PDE

Continuing my last question, I am trying to solve 3D Schrodinger equation (with mass equals to 1) and with no potential (I will add a potential later). The analytical solution for this equation is ...
2
votes
1answer
142 views

Equal flow boundary conditions

I'm having the following issue that can not find a solution , by using this boundary condition v[0, t] == u[1, t] that Mathematica does not understand, someone ...
7
votes
2answers
230 views

Extending NDSolve beyond a singularity

The $\tan$ function satisfies the following IVP: $$y'=1+y^2 ,\quad y(0)=0 $$ and has simple poles at the points $x=\pi/2+ \pi n$ for integer $n$. When trying to get $\tan$ via numerical ...
0
votes
2answers
138 views

Convert DE to first order?

I have the following set of second-order differential equations: ...
5
votes
4answers
238 views

How should I deal with messages being produced by my Manipulate?

I have made a knockoff of John Polking's dfield for a web-embedded CDF (via Enterprise Mathematica). My code extends that from Belisarius who last year showed me how to use ...
0
votes
1answer
84 views

Problems with solving two coupled lorenz equation with NDsolve

I have a system of two coupled lorenz equations and I'm to integrate it: ...
1
vote
0answers
50 views

NDSolve error with a given system of equation [closed]

I know that my question will be very similar to NDSolve giving NDSolve::deqn: message when given my system of equations . Sorry for that, I am so new on Mathematica. I have tried to use the ...
0
votes
1answer
146 views

How to speed up NonlinearModelFit applied on a NDSolve solution?

I'm trying to perform fittings of a model defined through a system of ODEs to data consisting of time courses of a measured quantity. As described in a previous question, the system of ODEs is as ...
2
votes
0answers
252 views

NDSolve fixed step problem

Working example here: Drive folder (have both files in the same directory! Notice: the line <<variables' in the file seems to throw an error for me, but ...
3
votes
1answer
87 views

DSolve: Why does Mathematica does a separate integration from 0 instead of 1 and then from 1 to upper limit?

I am trying DSolve in Mathematica for an ODE. Here is the command I am trying and the result I get DSolve[{x'[t] + y[t] x[t] == 0, x[0] == x0}, x[t], t] $$ ...
5
votes
1answer
137 views

List of unknowns in a differential equation

I often find myself wanting to see what variables/constants are undefined in a large set of differential equations, is there a function that will do this for me? As an example, for the following ODE ...
2
votes
1answer
106 views

Fails to apply NonlinearModelFit on a numerically evaluated model defined as a system of ODEs

I call for your help with because I'm getting troubles to fit a model defined through a system of ODEs. The system of ODEs is as follows: ...
2
votes
1answer
147 views

solving coupled differential equations with parameters subject to optimization problem

I have an equation, but I don’t know how to program it in Mathematica: x[0]==1; y[0]==1; x'[t]==(a+b) x[t] y[t]; y'[t]==a b y[t]/x[t]; where ...
4
votes
2answers
278 views

NMinimize Constraints on Interpolating Functions

I have the following piece of code and am trying to use NMinimize (the last line in the following chunk of code) to minimize ...
1
vote
1answer
39 views

Reverse order in evaluated list

I solve the following set of ODE and plot solution of S1 vs k2 (Figure 1). ...
7
votes
1answer
198 views

Why should the spatial derivative order of the ODE *not* exceed two?

Following this question I came across this strange behaviour. Let me define a 1 D interval implicitely ...
7
votes
3answers
382 views

FEM Solution desired for “Plate with orifice” deflection: Application of Boundary Conditions and use of Regions

I found the deflection of an orifice plate (circular plate with a hole) subject to uniform pressure using Mathematica's NDSolve functionality. The plate is fixed at its outer circumference and is ...