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

learn more… | top users | synonyms (3)

0
votes
1answer
63 views

Sum the several solution functions returned by DSolve [duplicate]

Could you please explain why this code does not work. ...
0
votes
0answers
44 views

fit multiple dataset with a single ODE to extract parameters

I found my problem very similar to the one discussed in here. But I'm having difficulties to fit multiple dataset (with varied temperatures) to a single ODE. The hurdle for me now is since the rate ...
9
votes
1answer
147 views

Beam deformation due to localized force distribution

I would like to simulate the deformation of a beam which is fixed on one end and has some localized forces applied along the top. I have been inspired by this Mathematica example: http://www....
1
vote
0answers
67 views

Solving two PDEs (for f[x,y] and g[x,y]) coupled through a non-linear function |Grad[f[x,y]]|^2

I am new to Mathematica, trying to solve two PDEs (for f[x,y] and g[x,y]) coupled through Grad[f[x,y]]. In the following code when n is set to 1, NDSolve seems to be able to solve for f[x,y] and g[x,y]...
2
votes
1answer
76 views

Problem with 3rd order ODE solution from DSolve

I Dsolved the following 3rd order ODE: x'[t]==(3x''[t]^2-Sqrt[9x''[t]^4+4x'''[t]^2])/(2x'''[t]) Mathematica (version 8) gave me three solutions including the ...
1
vote
1answer
59 views

How to rescale a NDSolve curve to fit in a box?

I'm numerically solving some differential equations and the 3D curve I get has a very complicated shape and size that cannot be predicted from the initial conditions. I need to rescale the whole ...
0
votes
1answer
75 views

Solve the n ODE in the loop, add the solutions to create one function

I need to solve the n ODE and create a function which is a sum of all the solutions. I wanted to use the code below, but only part (solve) and (create solution) works well. How can I store the ...
2
votes
1answer
167 views

The numerical solution of a nonlinear ODE: boundary value problem

I want to solve the following nonlinear ordinary differential equation with boundary values $f(0)=1$ and $f(1)=1$: $${\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}f \left( x \right) ={\dfrac { \left( 51-51\,x ...
2
votes
2answers
74 views
1
vote
2answers
85 views

A way to show all the 3D curve, with view still centered on {0, 0, 0}

I frequently have this graphics output problem : I draw a 3D curve using NDSolve, and depending on the initial data, there's no way I could predict the shape and ...
4
votes
3answers
660 views

Does it mean that there no analytical solution if Mathematica cannot find analytical solution?

I'm solving simple but coupled ODEs recently. I use both MATLAB symbolic computation and Mathematica. For example, my coupled ODE is the following \begin{align*} &\dot{x}(t)=y(t)-\rho b\frac{x(t)...
1
vote
1answer
73 views

Variable coefficient Differential Equation

Can anyone help with the following differential equation. The errors i am getting are also given. ...
1
vote
1answer
101 views

Solve a system of nonlinear ODE

I have a set of coupled non-linear ODEs, they look symmetric. But when I use matlab symbolic calculation to solve, matlab told me that no explicit solution found. I was recommended by my friend to try ...
0
votes
1answer
66 views

NDSolve and a vector function [duplicate]

I'm trying to solve a set of differential equations which include a numerical function that takes a vector and returns a vector. ...
2
votes
1answer
121 views

Improving the speed of evaluating numerical solution to a sytem of PDE

I am running a simulation for a kinetic reaction inside a pellet modeled by a system fo reaction-diffusion equations. ...
1
vote
1answer
57 views

How to define “DiffusionCoefficients” in InitializePDECoefficients that depend on ElementMarker?

Modelling a Poisson Equation on piecewise constant material, ElementMarkers on the PointElements are used for the boundary ...
3
votes
2answers
255 views

Why can't FindFit get the proper result in this problem?

Cross-posted in http://tieba.baidu.com/p/4425634583 The mathematica helps me solve one differential equation in numerial method and got the data of the answer,but the taylor series fit is not good in ...
0
votes
1answer
54 views

Reduce in NDSolve [duplicate]

Second order ode; neither the trig RHS simplifies nor is the solution obtained readily: ...
1
vote
0answers
61 views

How to use interpolating functions in NSolve with a piecewise function?

I have a system of two non-linear differential equations that are piecewise, the solutions of which I would like to use to obtain the values of a third variable that is also piecewise. I have trouble ...
1
vote
1answer
115 views

Problem with solving a system of differential equations [closed]

I am trying to reproduce a result that is part of a derivation of the flow due to a rotating disk. I have as given this system of differential equations. (In these equations primes indicate ...
0
votes
0answers
59 views

ndsz : step size is effectively zero; singularity or stiff system suspected

This is the first time I ask a question. I have seen many solutions to this error and tried but they are not working. Here is the code. Needs["DifferentialEquations...
0
votes
0answers
50 views

Output of NDSolve is useless when the system gets bigger

I have a bridge, which I want to discretize on $n$ stiff rods. To find a solution I wrote the following code in mathematica: ...
3
votes
0answers
59 views

EvaluationMonitor slows down NDSolve considerably

Hi guys I'm solving a nonlinear system of ODEs with NDSolve. Since the dimension of the system ranges from 2700 to 20000 equations, depending on the parameters, the time required to solve the entire ...
-3
votes
1answer
86 views

How to implement Bifurcation of delayed Lorenz system? [closed]

The delayed Lorenz system is as follows: (Bifurcation parameter as $\tau$) $$\begin{align*} D^\alpha x_1(t) = & a_1(x_2(t-\tau)-x_1(t)) \\ D^\alpha x_2(t) = & a_2x_1(t)-x_2(t)-x_1(t)x_3(t) \...
5
votes
1answer
92 views

Passing starting values to FindRoot when called by the Shooting Method of NDSolve

According to the documentation of the NDSolve Method -> "Shooting" option, it is possible to pass to ...
0
votes
0answers
93 views

NDSolve with default settings solves much faster than classic Runge-Kutta 4

NDSolve with default parameters takes 0.25 s for the first run. With Runge Kutta 4 specified, it takes about 4.6 s, ...
0
votes
0answers
67 views

How to systematically approach when solving a large system of differential equations?

I defined several functions in order to get to the final system of a partial differential equations I have to solve. Here is how I did those definitions: ...
0
votes
1answer
93 views

Proper way to produce StreamPlot

I have faced with problem of not completely correct phase trajectories of ODE, which was produced by StreamPlot. Next following code ends up with this plot: ...
3
votes
0answers
56 views

Singularity treatment in a simple problem [duplicate]

With a simple search, several posts can be found on treating singularities encountered while using NDSolvefamily of numerical methods. However, most of them involve ...
1
vote
2answers
120 views

Discontinuous boundary condition in NDSolve

Assume we have any PDE to be solved on the rectangular domain $0<X<4$ and $0<y<2$ How do we tell Mathematica to impose the following boundary conditions? $\cases{ U[0,y,t]=0,& if $0&...
1
vote
2answers
72 views

If statement, differential equation

I've got problem with If statement. It doesn't work. I've never used Mathematica before, probably it's a problem with syntax. ...
1
vote
0answers
37 views

wrong sign with NeumannValue

I want to use NDSolveValue with Neumann boundary condition. The documentation says this https://reference.wolfram.com/language/ref/NeumannValue.html . In the ...
1
vote
2answers
241 views

Differential equation and For loop in Mathematica

I'm trying to optimize some kind of process in Mathematica. I did it without any problems in Matlab and now I'm trying to do the same in Mathematica which I don't know at all. I'm trying to find the ...
1
vote
1answer
83 views

How to plot an Ordinary differential equations

Hey guys I'm using mathematica 5 now, and so new to this software thing. May I know how to plot the direction field using PlotVectorField[] and also finding ...
1
vote
0answers
46 views

NDSolve with (approximately continuous) Piecewise input

I am finding some troubles when using NDSolve with this piece of code: ...
0
votes
1answer
35 views

WhenEvent StopIntegration not working [closed]

I am trying to use a WhenEvent to StopIntegration when Tmelt=1356 is reached. Although the NDSolve appears to be working, "StopIntegration" is not stopping the integration and no value is given for ...
3
votes
1answer
81 views

NDSolve with an interpolating function in Piecewise

Again I have question I hope someone can answer. I made a very simple example illustrating my more difficult problem. I want to solve a differential equation numerically with a conditional in the ...
3
votes
4answers
225 views

Stiff coupled system of ODE's

I am having trouble with a possible stiff nonlinear system of ODE's. ...
0
votes
1answer
108 views

NDSolve coupled Non Linear differential equation, Obtain solution

Trying to solve a coupled differential equations of motion. When there is no Working Precision mention it runs and give following plot for θ[t] but this is wrong, θ must be symmetric about 0. ...
0
votes
0answers
25 views

Discretization of a region for NDSolve using [duplicate]

I would like to mesh the region below in order to use it for a calculation by means of NDSolve. Could anyone help me to discretize this region for NDSolve using. Thank you very much. Needs["NDSolve<...
2
votes
1answer
87 views

System of differential algebraic equations with one of the variables as an argument of the other varaibles

This is my first question here and I don't even know what to search for to find a solution. So I'd appreciate it even if someone can point me in the right direction. EDIT: My post had many problems, ...
6
votes
3answers
148 views

Force WhenEvent to respect default stepsize

I have to numerically integrate an equation system and monitor the accumulating datapoints. For example, I fit a line to a subsample of the points and terminate via ...
6
votes
0answers
236 views

Problems with NDSolve for solving coupled Schrödinger-like PDEs with perfectly matched layer boundary condition

I'm trying to solve a full-vectorial wave equation for an arbitrarily shaped wave guide, by using NDSolve and perfectly matched layer (PML) conditions. The PML ...
0
votes
1answer
82 views

Difficulty using NDSolve

I am tying to solve two differential equations with two delays numerically. The result didn't come out and it was written that “NDSolve::ndsz:...step size is effectively zero; singularity or stiff ...
0
votes
0answers
22 views

The result didn't come out and it was written that “step size is effectively zero; singularity or stiff system suspected” [duplicate]

I am tying to solve two differential equations with two delays numerically. The input is given by the picture below. I have tried the methods that have been posted in the web but it didn't work. ...
0
votes
1answer
99 views

How to prevent extrapolation in Manipulate? [duplicate]

How can I limit the time slider to the NDSolve solution range in a Manipulate box, so no extrapolation is used ? More specifically, here's a stripped down version of my Manipulate code which shows ...
2
votes
0answers
31 views

exe file creation [duplicate]

HOW CAN I MAKE IT exe FILE? PLEASE HELP... ...
4
votes
2answers
116 views

Integrating an interpolating function

I have been using the following NDEigensystem command to generate the eigenfunctions of the laplacian with a square domain: ...