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

learn more… | top users | synonyms

1
vote
2answers
48 views

Using manipulate with a defined funtion to solve a ODE

As many in this community I'm new in Mathematica and while exploring the numerical solver for differential equations, I tried: ...
6
votes
1answer
690 views

Schroedinger eigenvalue problem in two dimensions (Harmonic Oscillator)

I read here, the discussion about how to solve one dimensional eigenvalue problem. I am wondering, how can one generalize these methods to two dimensions. For example: ...
2
votes
1answer
93 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 ...
0
votes
3answers
88 views

Delay difference equation

How can I correctly specify and solve a system of delay difference equations in mathematica similar to NDSolve with regard to delay differential equations? For example if I want to "see" the values ...
2
votes
1answer
255 views

Finding the eigenfunctions of one and two dimensional Harmonic Oscillator

(Edited) For finding the ground state wave function of: $ H\psi(x) = (-1/2)d^2\psi(x)/dx^2 + (1/2)x^2\psi(x) = E \psi(x)$ I have written: ...
23
votes
3answers
574 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
2answers
56 views

Using parts of piecewise function

I define: f[x_] := Piecewise[{{5, 0 <= x < 10}, {g[x], 10 <= x < 15}, {h[x], x > 15}}] Then I try to solve ODE: ...
3
votes
1answer
128 views

How to solve Laplace equation in 3D?

Basically I want to solve Laplace equation for truncated octahedron in a cube matrix. The boundary condition is Concentration u=200 at surface of truncated octahedron and u=15 at boundary of cube. I ...
0
votes
1answer
99 views

Implicit differentiation and the Folium of Descartes [duplicate]

Given the function (Folium of Descartes) $x^3 + y^3 = 3xy$, how would I find the equations for tangent and normal lines at the point $\left( \frac{3}{2},\frac{3}{2} \right)$? I know that I must use ...
0
votes
1answer
73 views

Phase Portrait Trajectories

How do I make a phase portrait which looks like the ones demonnstrated here? Not the fancy version where I can use the controls to change it, just a screen-shotted version of say a spiral point? Any ...
-2
votes
0answers
39 views
0
votes
0answers
50 views

Solving a system of differential equations using NDSolve

I am trying to model the descent of a diver jumping from space, and have the following system of equations: ...
0
votes
1answer
51 views

Using InterpolatingFunction in equations

I tried to use InterpolatingFunction in another ODE, but it doesn't work, because it seems for Mathematica that 2 functions are unknown instead of one. ...
0
votes
0answers
30 views

DSolve of the system of two equations with known constants

In my previous post DSolve solution of two unknowns we conclude that Math 10 obtained results fast and that Math 9 did not give the solutions at all after one day. Now I have a problem if there are ...
2
votes
0answers
63 views

DSolve solution of two unknowns

I have one problem with the solution of the system of two differential equations. I waited hours, but Mathematica didn't give me a solution. Is there any way to obtain the results if we suppose that ...
-2
votes
1answer
41 views

Find constants of differential equations [closed]

I have the equation $y'' + 16y' + 64y = 0$ Repeated root of aux eqn $r = -8$ Given the initial conditions, how do I solve for the constants $A$ and $B$? ...
0
votes
0answers
16 views

Differential equation : y'(x)=(x+y(x)+1)/(x-y(x)+2) [migrated]

I have a question related to the differential equations.I have to solve the following equation: y'(x)=(x+y(x)+1)/(x-y(x)+2) However I don't see what should I ...
0
votes
1answer
48 views

Solving coupled ODEs with NDSolve

Given the system of differential equations: $i\frac{d}{dt}c_{1}(t)=H_{12} E(t)e^{(iw_{12}t)}c_2(t)+H_{13} E(t)e^{(iw_{13}t)}c_3(t)$, $i\frac{d}{dt}c_{2}(t)=H_{21} E(t)e^{(iw_{21}t)}c_2(t)+H_{13} ...
7
votes
3answers
300 views

Envelope for harmonic oscillator [duplicate]

I am able to take the expression $\cos 11t-\cos 12t$ and use the sum to product identity to write it in the form $$\left(2\sin\frac t2\right)\sin\frac{23}{2}t.$$ I can then plot the function and use ...
1
vote
1answer
104 views

error DSolve::dsfun: cannot be used as a function

I tried to solve system of ODEs, but something was wrong. ...
1
vote
1answer
41 views

WhenEvent error in NDSolve

I’m trying to solve a discontinuous differential equation where there are lost of energy when the trajectory cross the discontinuous manifold. The code I used was ...
0
votes
0answers
30 views

NDSolveValue:dsvar error in a system of non linear differential equations generated using a Table

I have a system of 4 non-linear differential equation, that I generate using Table command. On solving, I keep getting a dsvar error that is not going away. I debugged my code and figured it has to do ...
5
votes
2answers
76 views

Initial time as parameter in ParametricNDSolve

I need a help with the function ParametricNDSolve. My goal is solve the equation \begin{array} &&\dot{x}(t) = y(t) \\ &\dot{y}(t) = x(t)-1-\varepsilon Cos(\omega t) \\ &x(t_0) = x_0\\ ...
4
votes
1answer
361 views

PDE with Stefan Conditions, a.k.a variable boundary

I want to solve the one-dimensional one-phase Stefan problem, but I don't know how to make Mathematica understand the conditions. If you are not familiar with what I'm asking please refer to this ...
7
votes
0answers
2k views

Integro-differential equation

I have to numerically solve a nonlinear partial integro-differential equation using Mathematica. This is my equation, $$\frac{\partial y(x,t)}{\partial t}=\int_{-\infty}^\infty K_0(|x-u|) ...
17
votes
4answers
1k views

How to splice together several instances of InterpolatingFunction?

I have a set of InterpolatingFunction returned by NDSolve which are valid over different (but overall continuous) domains. How ...
0
votes
0answers
34 views

Fredholm integral equations [migrated]

I'm having problem obtaining the solution of the homogeneous Fredholm integral equation of the 2nd kind with a separable kernel. I always get a zero if I use the normal method I was taught for the non ...
0
votes
1answer
156 views

Defining a curved region for NDSolve

I have an equation to solve using NDSolve which involves two parameters: $en$ and $k$. The equation produces a set of curves when $en$ is plotted against $k$. In ...
0
votes
1answer
93 views

How to solve and plot this 2nd order non linear ODE? [closed]

Is this the correct way to plot the solution of this ODE from $(-8, 8)$? $$y'' - (1 - (y')^2)y' + y = 0$$ ...
1
vote
1answer
703 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? ...
0
votes
1answer
42 views

NDSolve emitting error message NDSolve::ndnum [closed]

I am trying to solve this system equn1 = ψ'[x] == - Elec[x] equn2 = Elec'[x] == - ((di/ a)^2)[(Exp[ψ[x]] - (x/(xoi*Sqrt[-ψ[x]])))] For that I am using ...
2
votes
1answer
44 views

How to write ODE to standard form?

Package like odepack needs the ODE written in standard form, which means write the high order ODE to first order ODE equations. The steps of converting ODE to ...
0
votes
0answers
34 views

Graphics Plotfield

A number of older books use the command: <<Graphics`PlotField` But Mathematica 10.0.2 says it cannot find it. Now, I know VectorPlot, but I am wondering if ...
1
vote
1answer
71 views

Failed to create a special random-disturbed initial condition for partial differential equation

I try to use BSplineFunction to creat a initial condition for NDSolve which has periodic boundary condition, that is realized by ...
0
votes
1answer
103 views

Generate a small random disturbance on a flat surface for NDSlove

I want to generate a small random perturbation on a flat surface defined on a square domain, e.g. 1 + 0.05 rand(x,y), which can be used as initial condition with periodic boundary conditions in ...
9
votes
1answer
576 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: ...
1
vote
2answers
103 views

Linearization of differential equations

I was wondering if one could define an operator such that, when I give a certain number of (differential) equations as an output, and an "equilibrium" value for each of the variables, it returns the ...
-1
votes
0answers
81 views

Linearization of differential equations at first order

Currently, I migrate my Maple Code to Mathematica code. On a part linked to linearization, I meet some difficulties. In fact, i try to linearize some differential equations around an equilibrium ...
0
votes
1answer
101 views

Solving Integro-Differential equations

I need help solving this equation. Is there a built in function that solves this type of equations? DSolve wouldn't work. Updated equation: ...
0
votes
0answers
855 views

PDE problem with initial and boundary conditions

I'm new to mathematica and I don't know what I do wrong, I have written the following PDE problem, with both initial and boundary (four in total) ...
0
votes
1answer
48 views

Unable to break out of two Do loops

I wrote a mathematica program to solve one dimensional schrodinger equation for particle in box model with zero potential. I used two do loops to achieve the task. The inner do loop will solve ...
-1
votes
1answer
71 views

amplitude-phase modulation equations solution

I am trying to solve the following equation ...
0
votes
1answer
66 views

Use ParametricNDSolve with Plotting

I am attempting to solve a differential equation of the form $$ \frac{dx}{dt}=ax+by \\ \frac{dy}{dt}=cx $$ with parameters $0<a<1$ and $a,b,c \in \mathbb{R}$. I am able to use NDParametric ...
2
votes
2answers
117 views

How to change the form $A\cos \omega t + B\sin \omega t$ into the form $\sqrt{A^2+B^2}\cos(\omega t-\phi)$ [duplicate]

I'm preparing notes for teaching next week's class in differential equations. My hand calculation solution of $$y''+192y=0, y(0)=1/6, y'(0)=-1$$ is $$y=\frac16\cos 8\sqrt3 t-\frac{\sqrt3}{24}\sin ...
-2
votes
1answer
66 views

Transform differential equation to get equilibrium position

I would like to determine the equation in the equilibrium position. eq = 1/2 g l m Cos[θ[t]] == J (θ^′′)[t] ...
1
vote
0answers
48 views

ParametricNDSolveValue over a region

I'm getting the following behavior I don't understand when using ParametricNDSolveValue over a region. You can reproduce the behavior considering Poisson's equation ...
0
votes
0answers
40 views

NDSolve — using a PDE solution in another PDE

My problem is to use the numerical solution of a PDE as a function appearing in another PDE. Here is the code solving the first PDE (heat equation in a domain) : ...