This tag concerns computational problems central to mathematical and scientific computating. The scope includes algorithms, numerical analysis, optimization, and linear algebra, computational topology, computational geometry, symbolic methods, and inverse problems.

learn more… | top users | synonyms

1
vote
2answers
41 views

The eventual advantage of a primality test without known exceptions

The primality test of Fermat with base $2$ seems to be as secure as the computer hardware for testing numbers big enough. However, I think there are an infinite numbers of false primes using this ...
0
votes
0answers
20 views

Maple vs Mathematica [on hold]

Can you tell which one is a better software for evaluation of questions... in terms of finding their closed forms. And solving algebra equations. I don't know about these softwares but I only heard ...
1
vote
2answers
163 views

Sequences of a computable function

Is there any computable function $f(n)$, which given any integer $n$ has been proven to return either $0$ or $1$ in finite time, and for which the statement "$f(1), f(2), f(3),\ldots$ contains ...
0
votes
0answers
48 views

Is the cohomology ring of a CW complex computable?

There is a well-developed technology for computing the cohomology groups of a CW complex, cellular cohomology. It reduces the problem of computing cohomology to the two simpler problems of (1) ...
3
votes
4answers
97 views

Is there an alternative way to represent the $\operatorname{diag}$ function?

In optimization, it is common to see the so called $\operatorname{diag}$ function Given a vector $x \in \mathbb{R}^n$, $\operatorname{diag}(x)$ = $n \times n$ diagonal matrix with components of $x$ ...
2
votes
2answers
525 views

Space spanned by matrices

I have a set of $5$ by $5 $matrices, $M_1,M_2,...,M_{19} ,M_{20}$. I want to try to find a basis from this set and also to find relationships between these matrices. This is how I think I should ...
0
votes
0answers
16 views

How to define a variable which is an integral involving cauchy principal value inside?

How to define a variable which is an integral involving cauchy principal value inside in any computer programming language? I want to know how to break down the procedure step by step from a ...
0
votes
0answers
48 views

To solve multivariate polynomial equations

For a system of multivariate polynomial equations like this: $$ \left( {\begin{array}{*{20}c} {\frac{{124}} {3}} & { - 24} & {\frac{{ - 68}} {3}} & {\frac{{68}} {3}} \\ {32} & {...
0
votes
1answer
27 views

Analyzing Accuracy of a method to compute summation

How does one analyze the accuracy of the following method to evaluate the summation? Also, what is the difference between analyzing the accuracy and proving that it is backward stable? I think I ...
1
vote
1answer
30 views

How to find all positive integer solutions of a Diophantine equation?

Here is the equation $$ 6a+9b+20c=16 $$ To solve this, i follow the below steps : $\gcd(6,9)(2a+3b)+20c = 16$ let, $w = 2a+3b$ So, $3w+20c =16$ then, specific solution of $w = 112+20n$, $c = -16-...
0
votes
0answers
18 views

Composite integration rule derivation

How does one derive a composite integration rule $\tilde{I}$ for $I=\int\limits_a^b f(x)dx$ using the integration rule $\hat{I} = w_0 y_0+w_1y_1+w_2 y_1'$ and $N+1$ equidistant points $x_i = a+i(b-a)$,...
3
votes
0answers
52 views

Why using primes as base in the Rabin-Miller test?

I have done some computer tests with the Rabin-Miller primality test: To test an odd number $n$, write $n=2^r\cdot s + 1$, where $s$ is odd. Given a number $a$ such that $1<a<n-1$, if $...
2
votes
0answers
34 views

terms of taylor expansions of multiple variables at the origin

By the fundamental theorem of symmetric polynomials, $X_1,X_2,\cdots,X_n$ are polynomials of $ e_1,\cdots,e_n$ and $$ \mathbb{Z}[ e_1,\cdots,e_n]=\mathbb{Z}[X_1,X_2,\cdots,X_n]. $$ We define a ...
0
votes
1answer
27 views

Finite difference formula to approximate second derivative

I have one question which asks to derive a finite difference formula to approximate $f''(x)$ in the form of $$f''(x)\approx Af(x+2h)+Bf(x+h)+Cf(x)$$ with the method of undetermined coefficients. ...
0
votes
1answer
1k views

Fast methods to check linearity of differentials? Generalizing linearity?

The L1 Mat-1.1010 -course here has taught me the linearity conditions $f(a x)=a f(x)$ and $f(a+b)=f(a)+f(b)$. I want to generalize it, some quite irrelevant slow investigation here. It requires time ...
4
votes
0answers
116 views

How to keep up when converting between bases?

Here is a schematized binary channel that neatly conveys a decimal number. $ \require{begingroup}\begingroup \def\T {{ \cal T }} \def \Ti {{ \T \raise5mu{ \text- \scriptsize 1 } }} \def\Bx #1{{ ~ ...
0
votes
0answers
21 views

First Intersection Of Periodically Repeating Intervals

I have a set of coupled tasks, let's say $M$ of them. The $ith$ coupled task is represented as the following 3-tuple $\{A_i,D_i,B_i\}$ where $A_i$ represents the time it takes to perform the first ...
0
votes
1answer
19 views

Combination with Repetitions Including Duplication of N Value

Is a Combination with Repetition the correct term for the following problem N - Letters a, b, c R - 2 Example Result should equal aa ab ac bb ba bc cc ca cb Total ...
3
votes
1answer
217 views

Conjecture about Rabin-Miller pseudo prime test

I tested the Rabin-Miller pseudo prime algorithm using a single test value and found that the number of false calls depends on the size of the number to test, reducing to a (conjectured) negligible ...
7
votes
2answers
282 views

A conjecture about the prime function $p_n$

While testing my system Zet for computational mathematics I find possible relations now and then. The latest is: Conjecture: For all $(m,n)\in\mathbb Z_+^2$ except $(3,4),(4,3) \text{ and } (4,4)$...
0
votes
1answer
20 views

Product of all Square Roots, taken only Decimal Digits

How and where could I compute the decimal reminder of a product of square roots times ten: $$Dr\left( \prod_{x=1}^{k}x^\frac{1}{2} \right) \times 10$$ Where $k$ is a power of $10$. I would like to ...
0
votes
0answers
68 views

What is the fastest known algorithm for finding eigenvalues?

What is the fastest known algorithm for finding eigenvalues? Second and third fast are also of interest if they are simpler, basically anything better than the standard solve characteristic polynomial ...
5
votes
1answer
622 views

Cylinder-ray intersections equation

I found an article involving infinite cylinder-ray intersections, and I don't know how they develop this equation: $$(q - p_a - (v_a, q - p_a)v_a)^2 - r^2 = 0$$ In the end of the first page I quote: ...
2
votes
1answer
66 views

Is this polynomial time for greatest prime factor of odd numbers?

For natural numbers $n$ and $x,$ the number of $n^{th}$ roots that have $x$ in the whole numbers place can be represented as $(x+1)^{n}-x^{n}.$ For $p$ prime, $(x+1)^{n}-x^{n}-1\equiv0\bmod p$ iff $n=...
0
votes
2answers
24 views

Newton Raphson Example help

For h:= $\mathbb{R} \rightarrow \mathbb{R}, x \rightarrow e^{x}-x^2+1$ I know the formula as $$X_{n+1}=X_{n}-\frac{f(X_{n})}{f'(X_{n})}$$ so this would give me: $$ X_{n+1}=X_{n}-\frac{e^{x}-x^{2}+1}{e^...
1
vote
0answers
32 views

R $\subseteq \omega$ recursive iff $\exists m \in \omega$ such that $R=\{n \ | \ \bar{\omega} \models \phi[m,n] \}$.

The queston I'm trying to solve is use Kleene's enumeration theorem to show R $\subseteq \omega$ recursive iff $\exists m$ such that $R=\{n \ | \ \bar{\omega} \models \phi[m,n] \}$ for some $m \in \...
1
vote
0answers
38 views

Normalizing an elliptic curve to find integer solutions

I have an elliptic curve $$ c_1y^2 + a_1xy + a_3 = c_2x^3 + a_2x^2 + a_4x + a_6 $$ with integers $a_1,a_2,a_3,a_4,a_6,c_1,c_2$ and I would like to find all integer solutions of this elliptic curve. I ...
4
votes
2answers
689 views

Software for numerical solution of a non-linear ODE system?

I have been given a nonlinear system of ODEs which has arisen out of a colleague's engineering research: $$\begin{array}{rcl} \dot{x}_0&=&x_1\\ \dot{x}_1&=&-\frac{\lambda}{(x_2)^n-k^2\...
0
votes
1answer
26 views

Write column form elementary matrix in terms of element form elementary matrices

Recall that any unit lower triangular matrix $L\in\mathbb{R}^{n\times n}$ can be written in factored form as \begin{equation} L = M_1 M_2\ldots M_{n-1} \end{equation} where $M_i = I + l_i e_i^{T}$ ...
2
votes
5answers
113 views

Minimize $a^5+b^5+c^5+d^5+e^5 = p^4+q^4+r^4+s^4 = x^3+y^3+z^3 = m^2 + n^2$ with distinct positive integers

Find the minimum value of the following: $$a^5+b^5+c^5+d^5+e^5 = p^4+q^4+r^4+s^4 = x^3+y^3+z^3 = m^2 + n^2$$ where all numbers are different/distinct positive integers. I know the answer (see ...
10
votes
6answers
574 views

Calculating logs and fractional exponents by hand

In view of what we can compute by hand, on a piece of paper, without having to use a computer or a calculator, how far can we go with the evaluation of $\log$-functions and fractional powers? More ...
2
votes
1answer
42 views

For what value of k onwards is it pointless for a computer to compute the probability mass function of the Poisson distribution

I am asking my computer to compute the probability mass function of $X \sim \text{Pois}(\lambda)$, a Poisson random variable. The function is: $$\Pr(X = k) = \frac{{e^{-\lambda} \lambda^k }}{k!}$$ ...
1
vote
2answers
35 views

Could we find an element on finite field? [closed]

Let $F$ be a finite field. Given an element $a^x$ in $F\setminus\{0\}$, could we find $a$?? I know that finding an integer $x$ is very hard problem (Discrete Logarithm Problem). However, I don't know ...
8
votes
1answer
98 views

Given two algebraic conjugates $\alpha,\beta$ and their minimal polynomial, find a polynomial that vanishes at $\alpha\beta$ in a efficient way

Inspired by this question, I was wondering about the following problem. $\alpha,\beta,\gamma,\ldots$ are the roots of an irreducible polynomial over $\mathbb{Q}$. How to compute the coefficients ...
1
vote
1answer
53 views

Generating all prime powers $\leq N$

Some very good algorithms exist to generate all primes $p$ up to some bound $N$, like the sieve of Erastothenes and the sieve of Atkin. However, suppose I want to generate a (sorted) list of all prime ...
2
votes
1answer
45 views

Mysterious functions

I originally asked the following question in stackoverflow, but the question is closed because some members meant that the question is about math(see the following thread) So I will give a try here: ...
0
votes
1answer
19 views

Number system question

Is it true that if we are given a number system with base 3, mantissa 2, $-1 \le p \le 1$, determined by $$\pm 0.d_1 d_2 \times 3^p$$ where each number is normalized, unless it is zero, then the ...
3
votes
1answer
51 views

Conjecture about odd primes

For each odd prime $p$ there exist $n\in\mathbb{N}$ such that $p\equiv n^2 \text{ (mod }\varphi(n^2))$, where $\varphi$ is Euler's totient function. I'm developing my Forth based computational ...
7
votes
5answers
384 views

What is a nice way to compute $f(x) = x / (\exp(x) - 1)$?

I want it to be stable near $f(0) = 1$. Is there a nice function that does this already, like maybe a hyperbolic trig function or something like expm1, or should I just check if $x$ is near zero and ...
1
vote
0answers
25 views

Ordering of elements in the base of a group

In section 4.6.7 of HANDBOOK OF COMPUTATIONAL GROUP THEORY, the authors use an ordering $\prec$ for the elements in a coset. That ordering, $\prec$, was defined in section 4.6 as follows. Throughout ...
0
votes
2answers
41 views

A matrix of a single 1 in each row and 0 elsewhere

Is there a particular name given to a matrix of m rows and n columns such that it must have one and only one 1 in each row and 0 elsewhere? For instance: ...
0
votes
0answers
14 views

Nature of state - recursion?

I always wondered how mathematicians define state (or rather: where it comes from?). This is tricky, because I always thought that in math there is only one "thing" - a pure, stateless function. Few ...
0
votes
0answers
14 views

Getting the most weekend for a specific date

Is there a formula that can find all the date(s) that fall on weekend the most times given a range of years. Example: Let say that it is given that the min year as 2000 and max year as 2005. and ...
0
votes
0answers
35 views

Positivity of the last component of non negative least squares based on active set method

I have followed the instructions given in Lawson and Hanson book for non-negative least squares using active set method. I am having a trouble in justifying one of the statements they have made about ...
0
votes
1answer
185 views

Multi layer perceptron activation function

How can you show that the Fourier series approximation of a function (so $f(x)=\sum\limits_{n=0}^{\infty} (a_n cos(nx) + b_n sin(nx))$ can be approximated to arbitrary precision by a feedforward ...
4
votes
0answers
537 views

Two-layer Perceptron for XOR

I'm reading Neural Networks for Pattern Recognition by Christopher M. Bishop. It's for a physics class, but I think the problem is closer to mathematics so I'm asking here instead of PSE. Chapter 4 of ...
0
votes
0answers
12 views

Finding a quadratic function the bounds any function below.

Given a function $f(x)$ I need to find another function $g(x)$, where $g(x) = Ax^2 + C$, where $A$ and $C$ are constants, such that $g(x) \le f(x)$, for all $x$. I can use a computer to do this for me ...
3
votes
1answer
1k views

Radial Basis Function and Neural Networks

I need a simple explanation about what is the radial basis function? And what is the relationship between the radial basis function and neural networks? And are there any simple examples to explain ...
0
votes
0answers
14 views

Help with Hidden Markov model and SMC methods

So its quite a long background i don't really know where to start but here goes. The background is as follows: Background Observation model As the target is moving, it measures the signal (RSSI) ...
2
votes
1answer
26 views

What is the algorithm used by Matlab for computing the Bessel function?

I am quite curious about the algorithm behind. It is definitely not the power series expansion, right? So, what is the trick? I cannot find it in the help file of Matlab.