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
0answers
7 views

stability of FTCS scheme for parabolic equation

Can you suggest any method for stability analysis of FTCS scheme for the the following parabolic equation ? D.E: $u_{t}=a(x,t)u_{xx}+f(x,t,u)$, $0<x<1$, $0<t<T$, $T>0$ BCs: ...
0
votes
0answers
16 views

Proof of theorem about iterative methods

How do I prove that if $A$ is a tridiagonal (or block tridiagonal) matrix then the corresponding $P_J$ and $P_G$ iteration matrices for the Jacobi and Gauss-Seidel methods satisfy that if $\lambda$ is ...
1
vote
1answer
13 views

Bounding volume of catmull-rom splines

I need to compute a 2D "spherical" bounding volume for the part of a catmull-rom spline $S(t)$ with four control points $P1$, $P2$, $P3$ and $P4$ in the domain $0 \le t \le1$. The purpose is to reduce ...
0
votes
0answers
25 views

The condition number of a function

Compute the condition number of $f(x_1, x_2)=x_1 + x_2$ I really don't see where is my mistake. Here is my solution: The condition number of a function $f$ is given by $\frac{\left\Vert x ...
2
votes
1answer
32 views

Solve system of first order ODEs when values are a matrix

I have a set of first order ODEs such that $$ \frac{d \Psi}{ds} = H(s)\Psi$$ H(s) is an 8x8 matrix. I would like to find $\Psi(b)$ where $b$ is 10, for instance. The initial value for $\Psi(a)$ is ...
1
vote
1answer
49 views

Find the coefficients in quadrature formula on $[0,1]$ with the nodes at $1/4$, $1/2$, $3/4$

In my worksheet I was given a question about numerical integration that says: Find the formula for $\int_{0}^{1}f(x)dx=A_{0}f(\frac{1}{4})+A_1f(\frac{1}{2})+A_2f(\frac{3}{4})$ I suppose the goal ...
0
votes
0answers
16 views

Gap between numbers in fixed-point vs. floating-point arithmetic.

If $r$ is a machine-representable number and $f(r)$ is the next larger machine representable number, are the following true or false? (1) In fixed-point arithmetic, the distance between $r$ and ...
1
vote
0answers
29 views

Why do we interpolate - no guarantee of success

this is somewhat of a general question about interpolation, I don't fully understand how can we be confident that our approximation is good, even if we know a lot of points. An example would be: ...
1
vote
1answer
35 views

Newtons Method for approximation of absolute maximum [on hold]

Use Newton's method (with $x_0 = 2.0$) to approximate the absolute maximum (correct to $5$ decimal points) of $f(x)=x\sin(x)$ on the interval $[0, \pi]$. Using Newton's method, the error does not ...
-1
votes
0answers
29 views

Runge's Phenomenon [on hold]

" Runge's phenomenon is a problem of oscillation at the edges of an interval that occurs when using polynomial interpolation with polynomials of high degree over a set of equispaced interpolation ...
1
vote
0answers
20 views

Will Gauss quadrature numerical integration work with a variable dx

The question kind of says it all, but I'm reading about Gauss quadrature from here: http://www.damtp.cam.ac.uk/lab/people/sd/lectures/nummeth98/integration.htm which gives an equation of this form: ...
1
vote
3answers
28 views

Book for ODEs and numerical solution

I would like to ask you information for a book. I want to (self) study ordinary differential equation and their numerical solution (with MATLAB). I am not a math student (life science) so I want a ...
1
vote
3answers
41 views

How to calculate power of a number with decimal exponent programatically?

I am trying to code and algorithm that can allow me to calculate power of a function with decimal exponent. The language that I am using to code in doesn't has any predefined power functions. I ...
0
votes
0answers
12 views

Solving implicit theta-method function numerically faster than using fixed point iteration

When we are using the theta-method to solve an IVP. We have the equation: $$y(x_{n+1}) \approx y(x_n) + h[(1-\theta) f(x_n, y(x_n)) + \theta f(x_{n+1}, y(x_{n+1})]$$ where $f(x,y)=y'(x)$ The only ...
0
votes
1answer
23 views

Initial value problems with known solutions?

I'm trying to find a list of IVPs with known solutions to test my implementation of some numerical techniques. The only one I know of is: $$f(x,y)= y' =-\lambda y\;,\;\;\; y(0)=1$$ with the ...
0
votes
0answers
8 views

Zero-stability for numerical methods, only applies to LMM's?

I'm trying to get a better grasp of the notion of zero-stability. Mainly I'm using a book by Leveque (Finite Difference Methods for Ordinary and Partial Differential Equations). Anywho, Leveque ...
0
votes
1answer
37 views

Where is the symmetry of Fourier transform in its implementation in Maxima and Wolfram Alpha?

From Wikipedia I saw that there is a symmetry of the Fourier transformation $F(F(f))(x) = f(-x)$ This matches the graphical explanation of the (German) Youtube video (9:15 to 9:45). I tried to see ...
0
votes
0answers
11 views

Does p-order consistency for a linear multistep method imply lower-order consistency?

In some worksheet for college it says: A multistep method $$\sum_{m=0}^s a_s y_{n+m} = h \sum_{m=0}^s b_mf(t_{n+m}, y_{n+m})$$ with $a_s = 1$ is consistent of order 1 iff the following algebraic ...
1
vote
1answer
11 views

Condition number of a matrix bounded from below and above?

Is condition number of an invertible matrix bounded from below? And is condition number always bounded from above for an invertible matrix?
-1
votes
0answers
39 views
+50

Different method for QR decomposition - is it possible

This method could also possibly be applicable to matrices of higher dimension, but for the simplicity of my question i will only ask it for $2$x$2$matrices. Suppose $A=\begin{pmatrix} a_{11} & ...
5
votes
1answer
98 views

Can I turn $Ax=b$ into $Ax=0$?

For a system of equations $$ \begin{bmatrix}d_1 & d_2 & \dots & d_n \end{bmatrix} \begin{bmatrix}u_1\\u_2\\ \vdots \\ u_n \end{bmatrix} = d_{n+1} $$ where each $d$ is a column of ...
1
vote
1answer
33 views

How to generate a point cloud with known symmetry?

So I would like to know if there are any published algorithms to generate point clouds with known symmetry groups, such as $D_{3h}$ or $O_h$ and stuff like that. I know lots and lots of point clouds ...
0
votes
0answers
11 views

Why does QUADPACK only enforce the least strict error boundary?

According to this reference (which is in agreement with my own numerical experiments), QUADPACK tries to fulfill the following accuracy requirement on the approximation error: |RESULT - I| $\le$ ...
1
vote
1answer
30 views

Algorithm for retrieving all the permutations (randomized) for a vector sequence 1…N with only unique values

Here is the problem: I have a vector of $N$ elements long (containing only unique values from $1...N$). I am searching for an algorithm to obtain all the (randomized) combinations possible, where ...
-2
votes
0answers
42 views

interpolating polynomial [on hold]

Let the function f(x) = ln (1/|x-1|) be approximated by an interpolating polynomial of degree 15 with the nodes equally spaced in the interval [-1, 0]. What bound can be placed on the absolute ...
1
vote
1answer
71 views

loss of significance errors in mathematical expressions

Find a way of computing $$f(x) = \dfrac{e^{3x} - \cos (2x)}{3x}$$ without serious loss of significance near $x=0$.
0
votes
0answers
29 views

Finite Difference - Forward Difference with 2nd order Accuracy: What to do at the boundary?

I implemented a BVP using a first-order finite difference scheme after the shooting method did not work reliably. Its the first time I have worked with this. The code works but I would like to move to ...
0
votes
0answers
15 views

traslate information-technology [closed]

i want to translate this text in frensh, any one can help me please On the one hand, DuMu x modules can be freely combined, on the other hand, dynamic polymorphism is avoided for reasons of ...
0
votes
1answer
27 views

finite element-question

can any any one help me to understand this paragraph First, the model domain G is discretized with a FE mesh consisting of nodes i and corresponding elements $E_ k$ . Then, a secondary FV mesh is ...
0
votes
0answers
10 views

fehlberg mthod and control step size

want to use Fehlberg method to solve a system of ode. But I don't how can I control step size? for example I use h=.5 and then approximate $y_{n+1}$ by RK4 and RK5 then what am I doing?
0
votes
0answers
13 views

Mapping from $\left(-\infty,\infty\right)$ to $[a,b]$ to reduce numerical error

Suppose $A = \left[\begin{array}{cc} \exp\left(x_1\right)&\exp\left(x_2\right)\\ \exp\left(x_3\right)&\exp\left(x_4\right) \end{array} \right]$, where each of $x_i\in\left(-1000,1000\right),$ ...
1
vote
1answer
34 views

Simpson's rule and Trapezoid Rule?

Let $S(n)$ and $T(n)$ be the approximations of a function using $n$ intervals by using Simpson's rule and the Trapezoid rule respectfully. My book then states: $$S(2n) = \frac{4T(2n) - T(n)}{3}$$ ...
0
votes
0answers
29 views

Numerical integration tolerance pitfalls

Consider that we want to estimate $$\int_{\pi/2}^{\pi/2+8\pi}sin(x)dx$$ (the value of this integrate is obviously zero) with the Midpoint rule. We start with the endpoints $a=\pi/2$ and $b=\pi/2+8\pi$ ...
0
votes
0answers
14 views

Is there a way to characterize the range of a chebyshev series through its coefficients?

Let $f$ be a Chebyshev series of order $n$ $$ f(x) = \sum_{i=0}^n a_i \cos\left( i \arccos\left(x\right)\right), x \in \left[-1, 1\right]. $$ Is it possible to characterize all the $\lbrace a_i ...
1
vote
1answer
59 views

Is Frobenius norm induced up to a scalar factor?

I know that the Frobenius norm is not induced since $||I||_F=\sqrt n\neq 1$. But what if we consider the norm $\frac 1 {\sqrt n} ||\cdot ||_F$? Thank you!
2
votes
1answer
94 views

How to solve $\tan x=x^2$ in radians?

How to solve $\tan x=x^2$ with $x \in [0, 2\pi]$? I try with trigonometry and many ways but the numerical solutions seems to be difficult.
0
votes
1answer
19 views

interpolate a function in a little interval

I want to interpolate a function that I have amount of it in this points: $0:.001:1$ and I have interpolate it in this points: $0:.0005:.001$. When I use lagrange method, I have a very bad result. I ...
2
votes
2answers
74 views

Solution of an equation involving even integers

If $x$ is any positive even integer $> 1$. I compute: $$ c = x + x! $$ Now assume instead $c$ (even integer) is given, and I want to get back the value $x$. Is it possible to find a simple ...
1
vote
1answer
24 views

Finding polynomial optimal in terms of least squares approximation

Find polynomial $w$ of degree at most $2$ optimal in terms least squares approximation for a function $f(x)=x^3$ in the norm $\|g\|=\sqrt{(g,g)}$, given that: $$ (f,g) = \int\limits^1_0 f(x)g(x)dx. ...
0
votes
1answer
44 views

Compute $\lim\limits_{x \to 1}[f(x)]$ and $\lim\limits_{x \to -1}[f(x)]$ for $f(x)=\frac{1}{1-x}-\frac{1}{1+x}$

Is it possible to rewrite expression $\frac{1}{1-x}-\frac{1}{1+x}$ in order to be able to find its values near $x=1$ and $x=-1$ more precisely? This is a question in a numerical methods course. Is the ...
0
votes
3answers
42 views

Definition of Global Convergence

I am confused by the notion of "global convergence" as used in numerical optimization literature, and did not find an exact definition for that yet. Now I try to double check my understanding here. ...
1
vote
1answer
42 views

Constrained Newton-Raphson method

Peace be upon you, I want to solve a system of two equations in which the existence of $ln\left(\frac{\alpha}{\alpha+\beta}\right)$ function makes some limitations in iterations of the Newton-Raphson ...
2
votes
0answers
32 views

What is the (currently) optimal root finding algorithm for multivariate functions? [closed]

Let's say we wish to find the roots of the function: $f(x,y,\cdots) = 0 \;,$ so, for a minimal example: $xy - 1 = 0 \; .$ I know there are different methods to solve this problem for the ...
1
vote
1answer
30 views

Solving second order differential equation numerically with values given at intermediate points.

I need to numerically solve the equation, \begin{equation} y''(x) + p(x)y(x) = 1 \end{equation} in the range [a,b] with conditions \begin{eqnarray} y'(\alpha) &=& 1\\ y(\beta) &=& 0 ...
0
votes
0answers
21 views

numerical (script) fit of function with 2 arguments

I would like to find the least-square fit for a 1D-function that takes two arguments. m(x,y) = d * (x-x0)^2 / (y-y0)^2 I would like to write a c++ routine to ...
2
votes
1answer
33 views

Comparison of Adams-Bashforth and Runge-Kutta methods of order 4

I have a system of ODE, that must to solve with numerical methods. I solve it by Adams_Bashforth with order4 and Runge-Kutta with order4 methods. Do you know with same length step which methods ...
0
votes
1answer
17 views

discretization of mixed boundary conditions in the advection diffusion reaction equation, the Crank Nicolson method

Sorry, I have some doubts regarding the discretization of mixed boundary conditions in a PDE. I have discretized my equation, but I doubt about the boundary conditions. I dont know if you have to ...
0
votes
2answers
48 views

Solving a differential equation numerically to plot particle path

I'm trying to plot the evolution of a particle in an accretion disk by solving the equation $$2X\frac{\partial X}{\partial\tau}=V_R(X,\tau)$$ where I have found $V_R$ numerically to be ...
2
votes
1answer
28 views

Online resources to learn numerical methods for PDEs?

I would like to get into a career that uses alot of applied math. I took a numerical analysis course in undergrad and liked it, so I plan to self-learn numerical methods for PDEs. Other than the MIT ...
1
vote
1answer
34 views

Approximation of values in [0,1] by sums of unit fractions

Let $U = \{(-1)^k\cdot\frac{1}{n}: n\in\mathbb{N}, k\in\{-1,1\}\}$ be the set of positive and negative unit fractions. For a positive integer $m\in\mathbb{N}$ and a real $x\in [0,1]$ we set $d(m,x) = ...