The tag has no usage guidance.

learn more… | top users | synonyms

1
vote
2answers
86 views

Numerical Integration involving multiple singularities [closed]

I want to solve the following double integral numerically. This integral has singularities in two points. ...
3
votes
2answers
240 views

NIntegrate and Integrate giving different results for ill-behaved function

I'm trying to integrate the following function with Mathematica 8: $$ I(a,b)= \int_0^1 \mathrm{d}x\int_0^1\mathrm{d}y \,\theta(1-x-y) \frac{1}{x a^2-y(1-y)b^2},$$ where $\theta$ is the Heaviside ...
0
votes
1answer
49 views

a question about singularity in numerical calculation

Suppose I want to numerically calculate a PDE as following $$a\frac{\partial^2}{\partial x^2}f(x)+\frac{b}{x}\frac{\partial }{\partial x}f(x)+c=0, x\in (-1, 1)$$ where $a$, $b$, $c$ are constants. So ...
0
votes
0answers
89 views

Help numerically solve this non-linear PDE with singularity

I have trouble finding a numerical solution to the partial differential equation below. It seems to be a singularity in the solution somewhere, so I searched online and found that it is suggested to ...
4
votes
0answers
57 views

How can I make NIntegrate aware of a singularity along a curve, e.g. a circle in a 3D integral?

I am having some trouble trying to get Mathematica to do a numerical integral over three dimensions which contains a singularity of dimension 1, and I would like some pointers to solid resources on ...
1
vote
1answer
132 views

To deal with infinity in NDSolve

Edit: According to Michael E2's comments, I did a few trail and error, and modify the question so that it uses exactly the same expressions I am having difficulties....
1
vote
0answers
42 views
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 ...
4
votes
0answers
82 views

Why does Mathematica gives a solution which is not in the domain of the cot function? [closed]

I've tried to solve the equation $\cot \theta = 2 \cot 2\theta$ with the command Reduce and Mathematica gave $\theta = n\pi$ as solution, with $n$ an integer. But $\...
0
votes
1answer
119 views
7
votes
2answers
296 views

NIntegrate fails to converge under almost any PrecisionGoal, MinRecursion etc. How can I trust the result?

I have been getting some ideas by reading other related questions in the forum, but the integral I have to do is not converging in many cases. The integrand is of the form: ...
1
vote
1answer
175 views

Plotting singular integral

I need to plot the following singular double integration $u(x,t)$ for $t=0, t=0.2, t=0.4, t=1$ and $t=2$ in one figure as a line. That is, I need $u$-coordinate vertically and $x$-coordinate ...
1
vote
1answer
85 views

Is it possible to avoid the singularity problem of a 2nd order differential equations by restricting the solution range?

I would greatly appreciate your help in implementation of Mathematica for solving a 2nd order differential equation expressing some important engineering problem. Apparently, this equation has no ...
1
vote
1answer
107 views

Singularity problem of a 2nd order differential equation

I would greatly appreciate your help in implementation of Mathematica for solving a 2nd order differential equation expressing some important engineering problem. Apparently, this equation has no ...
12
votes
1answer
172 views

Why do singularities disappear when performing convolutions with the Dirac Delta function?

I am attempting to perform a convolution involving the Dirac Delta function: $\int_{-\infty}^{\infty} \frac{1}{t+1} \cdot \delta(t+1)\ dt$ I would expect that the result of this integral is ...
0
votes
0answers
195 views

Boundary value problem: complicated functionals

I am trying to solve a non-linear boundary value problem where the differential equation is a convoluted function of the variable I am trying to solve for. The equation below presents the final form ...
7
votes
2answers
246 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 integration,...
4
votes
1answer
584 views

Multi-dimensional integral in the complex plane with poles and essential singularity

I've passed the last week searching a way to numerically integrate this multi-dimensional integral in the complex plane at the poles and avoiding the singularity at z=0: $$ \oint_{C}\oint_{C\ auound\ ...
0
votes
1answer
122 views

How to compute DSolve or NDSolve for kind of piece wise function?

I need to Compute the expression for two functions. q[t] and q1[t] over certain range of time. The code is given below. ...
2
votes
1answer
247 views
0
votes
0answers
103 views

solving singular matrix [duplicate]

Hello I have a matrix that is singular, and i wanted to use linear solve to solve it but that doesn't work with mathematica. The matrix given as a list looks like this: ...
3
votes
0answers
233 views

How to do Fast Fourier transform (FFT) for singular functions?

I want to do a 3-dimensional FFT on this function $\frac{\cos (x) \cos (y) \cos (z)-\sin (x) \sin (y) \sin (z)}{\left((1.0001+\sin (y)+\cos (z))^2+(0.0001+\cos (x)+\sin (z))^2+(0.0001+\sin (x)+\cos (y)...
2
votes
1answer
978 views
0
votes
1answer
165 views

Multiple NIntegral with singularity

I need to compute the following integral FF[0]. When evaluating FF[0] Mathematica simply takes too much time. The problem may ...
10
votes
0answers
457 views

How does Mathematica find a series expansion of expressions containing logarithms when there is a singularity at the expansion point?

I am looking for a good approximation to a function containing logarithms, especially at values close to zero. When I used Mathematica's Series command I found ...
2
votes
2answers
870 views

How to remove a singularity with WhenEvent

Bellow I have a differential equation which hits a singularity at low values of t. What I want to do is somehow utilize the WhenEvent command in order to replace ...
4
votes
2answers
822 views

What can I do to eliminate the error FindFit::nrjnum:?

I am testing the "Power Law with finite-time singularity" hypothesis for world population growth for a project. The data I'm using (same behaviour should also be exhibited by the stock market, thats ...
1
vote
1answer
1k views

How to find local minima of a multivariate function with singularity?

I am working on a logistic regression problem which requires minimizing a cost function J[{theta0, theta1, theta2}, X, y] to find the optimal value for fitting ...