Questions on numerical analysis/numerical methods; methods for approximately solving various problems that often do not admit exact solutions.

learn more… | top users | synonyms (2)

0
votes
1answer
33 views

Jacobi Method for the linear systems (for first two iterations)

I have the following 3 linear systems: $$4x_1 + x_2 - x_3 = 5$$ $$-x_1 + 3x_2 + x_3 = -4$$ $$2x_1 + 2x_2 + 5x_3 = 1$$ for x^(0) = 0 Then, I write it as following: $$x_1 = -1/4x_2 + 1/4x_3 + ...
3
votes
1answer
31 views

Newton Rhapson Algorithm Accuracy

I read somewhere that the NR algorithm in general (given an appropriate initial value) increases in accuracy by roughly two decimal places per iteration. Is this something that can be proven, or is ...
3
votes
1answer
76 views

I want to study Numerical linear algebra [closed]

Would you like to recommend a book to me? the proof is explicit and easy to understand is preferred.
0
votes
1answer
41 views

Parity of number of factors up to a bound?

Consider $b,n\in\mathbb{N}$ where $b\leq n$. We want to find the parity (ie. odd or even) of the number of divisors of $n$ that are $\leq b$. The question is to find a fast algorithm to find that ...
3
votes
1answer
111 views

Numerical integration of $\sin(p_{m})$ and $\cos(p_{m})$ for a polynomial $p_{m}$

I was wondering if anyone knew about any numerical methods specifically designed for integrating functions of the form $\sin(p_{m})$ and $\cos(p_{m})$ where $p_{m}$ is a polynomial of degree $m$. I ...
1
vote
1answer
55 views

Computing a solution of the Laplace-Eigenvalueproblem with Neumann-b.c.

Good day! I was considering the Laplace-Eigenvalueproblem with Neumann b.c., i.e. find $u \in H^1(\Omega) \setminus \{0\}$ and $\lambda \in \mathbb{R}$, such that: \begin{eqnarray} -\Delta u \ ...
0
votes
0answers
19 views

Runge-Kutta for newton's law with dependency

I'm trying to determine the changes (position and velocity) on a mechanical system during a step of time. I have a mobile mass whose position (everything is 1D-only) is denoted $x(t)$, velocity $v(t)$ ...
2
votes
0answers
10 views

Remainder of the minimax approximation polynomial - number of extrema

Recall some definitions. Let $f \in C [a,b]$. The minimax polynomial $p_n$ is the polynomial $p_n (x) $ of degree $\leq n$ that minimizes $||f-p_n||_\infty $. It can be proved that this polynomial ...
1
vote
0answers
13 views

Order of Romberg's method

We call a method(numerical integration) of $n-$th order, if it can integrate any polynomial of degree $n-1$ without any error. In this sense: The simpson rule is of $4$-th order and the trapezium ...
0
votes
1answer
9 views

Change of variables from intinite to bounded support.

I may be missing something simple, but I am stuck. My question: I am solving a system of partial differential equations numerically, but one of the variables can take on any value, ie $x \in ...
0
votes
1answer
37 views

Show that Newton method converges for every choice of $x_0$

Let $f\in\mathcal C^{2}(\mathbb R)$ with $f'(x)>0$ and $f''(x)<0$. If $f(x^*)=0$, show that the Newton method converges to $x^*$ for every choice of $x_0$ Wlog, choose a start point $x_0$, ...
2
votes
3answers
94 views

Intuitive Numerical Analysis Texts

Steven Strogatz has a great informal textbook on Nonlinear Dynamics and Chaos. I have found it to be incredibly helpful to get an intuitive sense of what is going on and has been a great supplement ...
1
vote
1answer
64 views

Function Approximation

I need to solve the following equation $$-\frac{\partial S(x,y,t)}{\partial t}=ax^2+bx\frac{\partial S(x,y,t)}{\partial x}+c\Big[\frac{\partial S(x,y,t)}{\partial ...
1
vote
0answers
37 views

Error bound by the Simpson's rule

My lecture notes have a little exercise. Two functions are given: $$ f(x) = \cos(x) \ \text{and} \ g(x)=\sqrt{x+1} $$ And we're asked about the error bound of the Simpson's rule to estimate the ...
0
votes
1answer
15 views

Converting x number of petaFLOPS into a base 2 number

I would like a few different formulas or methods for doing a couple of conversions and calculations: 1) How can I convert petaFLOPS into a base $2$ number representing how many operations per second ...
0
votes
0answers
19 views
1
vote
2answers
135 views

What is the sum that the square root button on calculator does so I can do it without the calculator button [duplicate]

I am not very good when it comes to Maths but the current work I am doing means I need to get better and quick. I have been teaching myself about areas, diagonals and square roots. However I am ...
0
votes
0answers
32 views

Lagrange interpolation polynomial and error estimation

Given is a function $f(x)$ with $f(0)=1$, $f(\frac{1}{2})=2$ and $f(1)=-1$. Additionally is given that $max_{x\in \left [ 0,1 \right ]}f''(x)=1$. Find its Lagrange interpolation polynomial $P$ and ...
0
votes
0answers
24 views

How to design Boundary condition for Euler equations (CFD)?

I'm developing on the calculation of the euler equations using the finte volume method. As you may know each cell is calculated by the incoming and outgoing flux. That means I need in a 1D System the ...
0
votes
1answer
42 views

An integration formula using Gauss-Laguerre method

Using Gauss-Laguerre method show that: $ \int_{0}^{\infty}\frac{e^{-x}}{x+a}dx=\frac{a+3}{a^2+4a+2}+\frac{4}{\theta^5} $ where: $ 0<\theta<1, a>0 $
0
votes
1answer
37 views

Degree of Precision Effect on Quadrature Accuracy

For an $n$ point Gaussian quadrature, one can show that it has degree of precision $2n - 1$ meaning it will exactly integrate polynomials of that degree or lower. Is it always true that a quadrature ...
2
votes
2answers
77 views

An analytical proof that the sequence from the chord's iterational method is monotonic

Assume $f:[a,b]\to\mathbb R$ is twice differentiable and $f'f''\not=0$ on the interval. Assume $f(a)<0<f(b)$. Let $x_0$ be the point of intersection of the $x$-axis with the line through the ...
2
votes
1answer
40 views

How to compute the pade approximation?

Like $\log(1+x)$? Is there any algorithms? I have read many materials but doesn't have an idea
0
votes
0answers
16 views

Solving a nonlinear volterra integral equation with two integrals each with a non separable kernel

I am trying to solve the nonlinear volterra integral equation ...
1
vote
0answers
17 views

Why do we get oscillations in Euler's method of integration and what is the period?

When using Euler's method of integration, applied on a stochastic differential eq. : For example - given $$\dot v = -\gamma v \Delta t + \sqrt{\epsilon \cdot \Delta t }\Gamma (t) $$ we loop over $$ ...
0
votes
0answers
30 views

Calculating log and trigonometric functions using only +,-,/,*

How to calculate logarithm and trigonometric functions (sin, cos etc.) on base n with using +,-,/,* ? Is there any way to do it?
1
vote
0answers
39 views

implicit non-linear equations with complex variables

I am trying to understand a methodology for solving implicit non-linear equations with complex variables. I would like to solve for z1 below where z2 is known. Also both z1 and z2 are complex ...
2
votes
2answers
59 views

Determinant of an ill conditioned matrix

I have the following ill conditioned matrix. I want to find its determinant. How is it possible to calculate it without much error \begin{equation} \left[\begin{array}{cccccc} ...
6
votes
1answer
72 views

Solution of $\exp(z)=z$ in $\Bbb{C}$.

I have posted a related question here. I thinkg this one is more interesting: What about the solution of $\exp(z)=z$ in $\Bbb{C}$? My try : $z \mapsto e^z - z$ is entire non-constant. Perhaps ...
1
vote
0answers
28 views

Numerical Computation for K smallest eigenvalues of a large Real Symmetric Matrix with restricted methods

I'm writing some code on a distributed platform, using some programming language like Hadoop, and now I need to calculate the K smallest eigenvalues for a Large Matrix. K is a small constant at most ...
0
votes
0answers
9 views

How to determine the upper bound of the global error when I calculate a ODE using Euler method

I have an ODE say, $u^\prime(t) = tan^{-1}u(t)$ $(0 \le t \le T)$, $u(0) = a$. Now I would like to decide the upper bound of the global error when using Euler method to solve this ODE. I know that the ...
0
votes
1answer
21 views

Help inverting a non-linear system of equations

I have a set of two equations like this $$ \gamma_3=\left(\frac{1}{\sqrt{1+2c_3^2+6c_4^2}}\right)^3 \left( \alpha_1\,c_3^3 + \alpha_2\,c_3c_4^2 + \alpha_3\,c_3c_4 + \alpha_4\,c_4\right)\\ ...
0
votes
1answer
20 views

Quadratic Optimization Problem with Box Constraints

I want to solve a problem of form $$\min_x x'Ax + b'x \;\;\mbox{ s.t. } l\leq x \leq u$$ where $A$ is a positive semidefinite matrix, thus the function I'm optimizing should be convex. However the ...
0
votes
0answers
33 views

System of linrar equations and condition number

The relative error of the solution of a system of linear equation $Ax=b$, for any natural norm $\|\cdot\|$ is bounded by $$ \frac{1}{\| A\| \|A^{-1} \|} \frac{\|r\|}{\|b\|} \le \frac{\|e\|}{\|x\|} \le ...
1
vote
0answers
27 views

Binary search (bisection method) - is it worth checking continuity

I am implementing a rather simple matlab code, that gets a function $f$ and 3 real numbers $a,b, \epsilon$ where $\epsilon >0$ is a very small positive number (for instance, no larger than ...
0
votes
0answers
36 views

On an interesting boundary condition

So I am tackling an interesting boundary condition, where $B(Du)=0$, for $x\in\Omega$, where $B$ is the signed distance function to $\Omega^*$ (where $\Omega,\Omega^*$ are convex domains in $\Bbb ...
1
vote
1answer
26 views

How to show that the one step method can't have consistency $p=3$?

I was looking at some exercises from last years' of my Intro to numerical math class, and found this: Consider the following explicit one step method: $$\psi^h x=x+h \gamma_1 f(x)+ h \gamma_2 ...
1
vote
0answers
34 views

Fastest way of finding eigenvectors from eigenvalues

Given the eigenvalue of a matrix of large dimensions, I want to know if there is a fast way of finding the corresponded eigenvectors?
0
votes
2answers
32 views

Convergence of the newton method

If f(x) has m'th root at r. consider the newton method $x_{i+1}=x_i-\frac{mf(x_i)}{f'(x_i)}$. why it has quadratic convergence?
0
votes
2answers
63 views

Newton Raphson Method for double roots

I am currently working on Newton Raphson Method. I am kind of facing a problem how Newton Raphson Method work on more than second order quadratic functions with double roots. I have googled and ...
2
votes
2answers
46 views

Initial value of Newton Raphson Method

I am currently studying Newton-Raphson Method. I feel that I understand the concept of it. Somehow, I am facing some question in my head about how to actually apply it. The questions that I have are ...
2
votes
1answer
43 views

Determine the coefficients of a polynomial knowing its roots

My prof. gave this problem as a bonus in an exam, and I couldn't figure out a solution. Some hints or a general method of solving it would be very nice. Given the following polynomial: ...
3
votes
1answer
54 views

What is the Most Efficient Way to Calculate the Internal Rate of Return?

I have built a program that prices financial assets and it does this in part by calculating the IRR. The problem is that it does not run as quickly as I would like it to. I currently use the ...
1
vote
1answer
68 views

How do I construct such a numerical method for solving ODE?

I am asked to expand $x(t+h)$ and $x(t+2h)$ around $t$ up to the rest term of the third order, find $A, B, C \in \mathbb R$ such that $$x'(t)=\frac{Ax(t)+Bx(t+h)+Cx(t+2h)}{h} + O(h^2)$$ and based on ...
0
votes
1answer
36 views

Anyone recognize this pattern? Plotting relationship between two parameters and their response.

First time asking a question here so hopefully I can provide enough information to you guys without explaining more than necessary. I'm doing some amateurish numerical analysis in MATLAB on some ...
2
votes
2answers
58 views

Integral over implicit function

Suppose the equation is given: $$y^5+2y^4-7y^3+y-x=0$$ (or any other equation that cannot be expressed explicitly). Let the solutions be implicitly given as $y=g(x)$. Is there an approach to solve ...
0
votes
0answers
29 views

Derivation of composite Gaussian quadrature error formula

I am working on studying for the Numerical Analysis qualifying exams. One of the questions I am stuck on is the following: Derive the error term for the composite Gaussian quadrature rule with ...
0
votes
0answers
35 views

condition number of orthogonal matrix

Let $A\in M_n(\mathbb R)$ be an orthogonal matrix. Then: $cond (A) =1$. I tryed to use facts about the eigenvalues but is did not help. In 2-norm it is easy to prove it since $||A||_2 = \sqrt{\rho ...
2
votes
2answers
83 views

What do mathematicians mean by “analytical solution of an equation”?

Given a PDE equations of the form: $\dfrac{\partial}{\partial t} u(t,x) = \left(\hat{L}+\hat{N_u}\right)u(t,x) \;\;\;\;\;\;\hspace{10mm}(**)$ where $\hat{L}$ is a linear operator and ...
1
vote
1answer
50 views

Which numerical method to use for ODE?

In practice what is the most common way to numerically estimate $y(t)$ (possibly using a series expansion) in the ODE with initial conditions, $$ y'(t) = f(t,y(t)), \qquad y(t_0)=y_0 $$ Wikipedia has ...