Skip to main content

Questions tagged [rootfinding]

For questions about the theory and process of finding the roots of a function (values where the function returns zero).

0 votes
0 answers
60 views

I am working on evaluating the integral: $ I(\vec{r}) = \int f(k) e^{i\vec{k} \cdot \vec{r}} d\vec{k}, $ for a system with a tight-binding dispersion relation given by: $ \epsilon_k = -2t \left[\cos(...
Lawerence's user avatar
0 votes
1 answer
116 views

Tried to solve this equation for a week with multiple initial guesses using Muller method ( In scilab ) . But no success as of now. Tested the algorithm by finding root for various different complex ...
Bailiwick3310's user avatar
2 votes
0 answers
31 views

I need to implement a root-searching algorithm for a function, that can possibly have some divergence points (DPs), or non-removable singular points, e.g. ...
GioK's user avatar
  • 21
5 votes
1 answer
238 views

I have a differentiable function $\mathbb{R}^n \to \mathbb{R} $ of several variables $f(x_1,\ldots,x_n)$, whose form I can write down and compute derivatives of. Typically $n = 8$. The function is ...
math_lover's user avatar
4 votes
1 answer
176 views

I have a smooth function $F: \mathbb{R}^n \to \mathbb{R}^{n-1}$ and points $x_0, y_0$ with $F(x_0) = y_0$. For theoretical reasons, I know that $y_0$ is a regular value of $F$, which means that the ...
Glenn Davis's user avatar
1 vote
0 answers
41 views

I want to solve the problem to determine the next intersection of a Keplerian orbit with the Sphere of Influence of a celestial body to find the next intersection within one future period of the ...
lamont's user avatar
  • 143
1 vote
2 answers
130 views

The definition of convergence for root finding algorithms is given in a few sources as: A sequence ${x^k}$ generated by a numerical method is said to converge to the root $\alpha$ with order $p\geq 1$ ...
Username_57's user avatar
1 vote
1 answer
761 views

I have been dealing with the standard libraries of scipy.optimize for rootfinding and optimization problems, but the problems i want to solve are very large, which makes the standard solvers run out ...
Klaus3's user avatar
  • 133
2 votes
0 answers
254 views

I've made a previous question here and also in SO wondering why only the fsolve solver converges for the simple one dimensional unsteady conduction problem $$ \frac{\partial T}{\partial t} = \alpha \...
Klaus3's user avatar
  • 133
0 votes
2 answers
849 views

I'm running a test problem to set up larger problems. Solving the simple unsteady heat equation via finite differences: $$ \frac{\partial T}{\partial t} = \alpha \frac{\partial^2T}{\partial x^2}$$ $\...
Klaus3's user avatar
  • 133
7 votes
1 answer
449 views

I'm using a numeric root-finder to find $k$ satisfying $\|A^k x\|=c$ where $A$ is a symmetric $d\times d$ diagonal + rank-1 matrix. How to compute $A^k x$ efficiently? For integer $k$, I can get the ...
Yaroslav Bulatov's user avatar
3 votes
1 answer
107 views

I have coded a simple program for a root finding problem using Halley's method. Here is the code: ...
Hosein Javanmardi's user avatar
4 votes
1 answer
155 views

Consider the following expression: $$f(t) = B+\sum_{k=1}^{N} A_k\cos(\omega_kt)$$ where $A$ and $B$ are known. the frequencies are also known but are not multiples of a fundamental frequency. However, ...
Hosein Javanmardi's user avatar
0 votes
3 answers
584 views

Let $f: \Bbb R \to \Bbb R$ be continuous. What are efficient algorithms to finding all the zeros in an interval $[a, b]$? I am actually only interested in the smallest zero in that interval, if there'...
Joe Shmo's user avatar
  • 121
3 votes
0 answers
198 views

I was working on a problem recently (calculating all flows in a network given input and output flows, basically what Hardy-Cross tries to achieve) which can be formulated as a well-determined system ...
nardi's user avatar
  • 131

15 30 50 per page