A computer algebra system (CAS) is a program which is able to carry out various symbolic manipulations with mathematical expressions. Some well-known computer-algebra systems: Mathematica, Maple, Wolfram Alpha, GAP, SAGE. For questions about Mathematica please see the ...

learn more… | top users | synonyms

1
vote
4answers
76 views

How would you represent an abstract set on a computer?

I'm trying to come up with a simple language for specifying mathematics. We can't simply define a set as a collection of objects because that's meant for human understanding and that would also lead ...
1
vote
2answers
67 views

How can I plot $y^x$?

How can I plot $y^x$? To keep things simple and to not have another $z$ variable on the other end of the equation, let's assume $y^x=10$. As long as that value is not $0$, the curve we get should look ...
1
vote
1answer
78 views

Chinese remainder theorem for polynomial evaluation

Let $R$ be a euclidean domain, $m_0,\ldots ,m_{k-1}\in R$ be pairwise coprime and $m:=m_0\cdots m_{k-1}$. The Chinese remainder theorem states: $$\varphi:R\to R/(m_0)\times\cdots \times ...
1
vote
3answers
59 views

If there is an $a\in\mathbb{Z}$ with $a^{n-1}\equiv 1\mod n$ but $a^{\frac{n-1}p}\not\equiv 1$ for all primes $p\mid n-1$, then $n$ is a prime

Let $n\in\mathbb{N}$ with $n\ge 3$ and $a\in\mathbb{Z}$ such that $$a^{n-1}\equiv1\text{ mod } n\;\;\;\wedge\;\;\;a^{\frac{n-1}{p}}\not\equiv1\text{ mod }n\;\;\;\forall p\in\mathbb{P}:p\mid n-1$$ ...
1
vote
1answer
80 views

Comparing expressions exactly

Suppose I have two expressions; call them $A$ and $B$. The following values of $A$ and $B$ are good examples for my question... $$ A = \pi e^2\\ B = \pi^2 e $$ Is there a method to determine the ...
1
vote
1answer
25 views

How to create a matrix in pari/gp [closed]

How does one create a matrix in $\texttt{pari}$? I currently understand how to create a vector: ? vector(6,n,n^2) %1 = [1, 4, 9, 16, 25, 36] I also know how ...
1
vote
1answer
318 views

MuPAD - rearrange an equation

I have a huge expression, say $f(x_1,x_2,x_3) = 0,$ and I want MuPAD to try to express $x_1 $ in terms of $x_2$ and $x_3$ Is this possible?
1
vote
1answer
106 views

how to use hilbert function in gap system with loadpackage singular

how to calculate hilbert function as it do in singular code singular code ...
1
vote
1answer
232 views

Programming PARI/GP to do a sum

I'm trying to compute the following sum in PARI/GP $C=\sum_{n=1}^{\infty} \frac{g(n)}{n^2}$ where $g(n)$ defined as $$g(n)=(-1)^r, \qquad r=\text{number of even indexed prime factors of $n$}$$ By ...
1
vote
1answer
257 views

Implementing a function in PARI/GP

I want to define a function: $$g(n)= \begin{cases} +1 & \text{if $n=1$},\\ +1 & \text{if $n$ is an odd indexed prime}, \\ -1 & \text{if $n$ is an even indexed prime},\\ (-1)^r & ...
1
vote
1answer
47 views

How many $\overline{a}\in\left(\mathbb{Z}/91\mathbb{Z}\right)^\times$ pass the Fermat and Miller-Rabin primability tests?

Let $$\text{F}_{91}:=\left\{\overline{a}\in\left(\mathbb{Z}/n\mathbb{Z}\right)^\times:91\text { passes the Fermat primality test to base }a\right\}$$ and ...
1
vote
1answer
252 views

Knuth-Bendix completion algorithm: word problem

Can someone explain me how to set up an algorithm to find the 12 normal forms of the group $A_4$ by making use of the Knuth-Bendix completion algorithm? So we have that $RRR=1, SSS=1$ and $RSRS=1$. ...
1
vote
1answer
43 views

Specifying if a function has an elementary integral

In Algorithms for Computer Algebra in the last chapter about Risch algorithm, the Rothstein-Trager method is applied to see if an elementary function has an elementary integral. For this, the ...
1
vote
2answers
200 views

wolfram mathematica, numerical integration, precision of a function/expression

I want to obtain the best numerical approximation (up to 10 decimal place would be ok for me) to an integral: $$ \int^{\infty}_{0} f(r)r^2dr $$ I am using the function $f(r)$, which is related to ...
1
vote
1answer
1k views

Existence of non-trivial solution of Sylvester equation.

I'm trying to solve a special case of Sylvester equation in my case it looks like $$A*X=X*B$$ so it can be written in form $$A*X+X*(-B)=C$$ where C consist of all 0 items. I tried to solve it in ...
1
vote
1answer
947 views

Algorithm for finding limits of compositions of simple functions?

There are two questions: Define the set $S$. Compute the limit of functions $f/g$ for functions $f,g\in S$, where $S$ is defined in the following way. All constant function are in $S$, $f(n) = ...
1
vote
1answer
21 views

Ideal intersection in Boolean polynomial ring

I'm trying to calculate $ \langle x,y,z\rangle \cap \langle x+1,y+1,z+1\rangle $ in the ring of Boolean polynomials. In CoCoA initially I set the ring as ...
1
vote
1answer
21 views

Computing Extensions of an Ideal in Singular or Macaulay2

Does Macaulay2 or Singular compute extensions of ideals under ring homomorphisms? Specifically, if $\phi:R\rightarrow S$ is a ring homomorphism (say polynomial rings over $\mathbb{Q}$ which can ...
1
vote
1answer
51 views

Higher Level Group Actions beyond OnSets

Is there a way to produce an orbit of a chain in a poset (whose elements are sets of sets) (that is, the group action is a group acting on the order complex of a poset)? I really want something that ...
1
vote
1answer
70 views

Math software to calculate in different rings.

I was doing computational experiment by using rings of polynomial like this let $f(x)=x^6+x\in\mathbb{Z_n[x]}$ for any given $n$. Is there any software which help me to calculate $f(a)$ where ...
1
vote
1answer
45 views

Variance of discrete probability distribution

I was wondering how I should calculate the variance of the following discrete probability distribution: $$P(y = 0|X) = w + (1-w)e^{-\mu}$$ $$P(y = j|X) = (1-w)e^{-\mu}\mu^{y}/y! \qquad j=1,2...$$ ...
1
vote
1answer
43 views

Number conversion in decimal fraction

980.85D convert to hexadecimal number = 3D4 . ?? how to solve the answer after the decimal point? Thank you in advanced.
1
vote
1answer
52 views

If $A\ne 0$ is a square matrix over a commutative ring with $\det A=0$, then its null space contains an element whose components are minors of $A$

Let $R$ denote a commutative ring and $A\ne 0$ a $n\times n$ matrix over $R$ with $\det A=0$. Then there exists a $x\in\ker A\setminus\left\{0\right\}$ such that all components of $x$ are minors of ...
1
vote
1answer
90 views

Number of monic irreducible polynomials over a finite field

Let $\mathbb{K}=\mathbb{F}_q$ and $\nu_n$ denote the number of monic irreducible polynomials over $\mathbb{K}$. It holds $$\nu_n=\frac{1}{n}\sum_{d\mid n}\mu\left(\frac{n}{d}\right)q^d$$ What I need ...
1
vote
2answers
156 views

Some questions about similar matrices

Two matrices A and B are similar, if and only if there exists an invertible matrix C with $A=C^{-1}BC$. A necessary condition for the similarity is, that the characteristic polynomials coincide. I ...
1
vote
1answer
137 views

Prevent Maple to evaluate before simplify a function

I have been trying to find a domain of $f(x)=\frac{x}{\frac{(x+2)}{(x-3)}}$ using different kind of software ( its clear the domain of this function is $\mathbb{R}\backslash \{-2,3\}$ ). When I tried ...
1
vote
1answer
54 views

Computational complexity of expanding a MacLaurin/Taylor Series

What methods exist to computationally determine the first $k$ coefficients of a function (possibly polynomial or rational polynomial function)? How do Mathematica/MatLab/Maple/etc. solve this ...
1
vote
1answer
57 views

Checking if a polynomial expression is constant in SAGE

I have a huge fractional-polynomial expression in SAGE that I have good reasons to believe is the constant function. Is there a command in SAGE like "== constant function" that I could use to check ...
1
vote
1answer
207 views

Computing univariate resultant via modified Euclidean algorithm

In an answer to the question Resultant of Two Univariate Polynomials, a PDF of course slides was linked which describes a modification of Euclid's algorithm for computing univariate polynomial ...
1
vote
1answer
177 views

How to calculate resultant of two polynomials without knowing the roots.

So in Rothstein - Trager's Method of evaluating logarithmic part they need resultant of two polynomial as shown in the image. My question is that how do they calculate the resultant without knowing ...
1
vote
1answer
87 views

Help in understanding risch algorithm

This is an excerpt from the paper that explains risch algorithm for symbolic integration. I do not understand how do we get to the residue directly and also how does it help us in the integration. ...
1
vote
1answer
156 views

Resultant of Two Univariate Polynomials

I am trying to implement an algorithm for computing Res(f(x),g(x),x) where f(x) and g(x) uni variate polynomials with integer coefficients. Could any one list the various algorithms for computing ...
1
vote
2answers
656 views

Plotting parametric equations in Maple

I have the following problem. Let $\alpha:=\alpha(v)$ and $\beta:=\beta(v)$, where $v\in\mathbb{R}$. Using MAPLE, I would like to plot $\alpha$ and $\beta$ in the $(\alpha,\beta)$ space as $v$ varies ...
1
vote
1answer
138 views

Computing Resultant

The resultant of two polynomials is defined as the determinant of the Sylvester matrix. If the polynomials are of degree $n$ and $m$, than the Sylvester matrix will be of dimension $(m+n)\times ...
1
vote
1answer
60 views

What are the explicit expression for this interpolation problem

We want to fit $f(x) = a_0 + a_1 *x + a_2 * x^2 + ... + a_n * x^n$ to the data $(x_i,f(x_i))$ for $i = 0 ... n.$ It will give rise to the following system $ A a = b $ Here $ a = [ a_1 a_2 a_3 ...
1
vote
1answer
97 views

How to test a formula online?

I have an equation I have devised, it is: $$\frac{H+O}{25}+\frac{D}{12970}+\frac{E}{363}+110 = X$$ The equation doesn't matter, I just posted it as an example. I need to be able to plug in values ...
1
vote
1answer
308 views

How to effectively compute a periodic function?

I'm writing a program to compute a value of periodic function for any arbitrary large argument: $f(k) = (\sum_{n=1}^{2^k} n)\mod\ (10^9 + 7)$, where $n,k \in \mathbb{N} $ I know that $ f(k + P) = ...
1
vote
1answer
126 views

Is there CAS design wisdom regarding categorizing expressions?

Is there CAS design wisdom regarding categorizing expressions? For instance. A user types in the expression $X^2 + X + 1$. Would it be best to classify it as a member of $R[X]$ (for some R, yet to ...
1
vote
0answers
9 views

Algorithm for computing the inverse limit of a finite inverse system

Let $k$ be a field (finite if you'd like), let $(I,\le)$ be a finite directed poset with $|I|=n$, and let $(A_i,f_{ij})_{i\le j\in I}$ be an inverse system of finitely generated, graded, commutative ...
1
vote
0answers
14 views

Algorithm for computing an inverse image

Let $k$ be a field (finite if you'd like), and let $f:A\to B$ be a map of graded, commutative $k$-algebras. Suppose further that $A$ is finitely generated and choose a presentation ...
1
vote
1answer
25 views

Simplify fog type function chain using sympy

I need to simplify some lengthy $fog(x)$ where both $f(x)$ and $g(x)$ are known rational polynomial in one variable. For e.g. $f(x) = (7x-9)/12$ and $g(x) = (2x-5)/7$. Then I would like to get a ...
1
vote
0answers
26 views

Determining equations for Lie symmetries

Am I right that using methods of the 1-forms (that usually implemented in system of computer algebra for Lie Symmetries), we can always generate determining equations for ODE, that solved for highest ...
1
vote
0answers
57 views

Maxima CAS and programmatically defining a function with a variable number of arguments. -?

This is a very simplified question of what I had asked. In Maxima, how can I include a for-loop counter in the left-hand side of an assignment, e.g. ...
1
vote
1answer
49 views

Generating a 2D coordinate from a 1D Index

Imagine you have a 5x5 matrix (zero indexed), where only the upper right coordinates are valid. I am unsure if this is the terminology to use, so the following coordinates should be the valid ones: ...
1
vote
0answers
84 views

Understanding BlowUp Computation in Singular

Many of us might know that "Singular" is a computer algebra system for Algebraic Geometry, Commutative Algebra and Non-commutative algebra. This is a procedure in "Singular" for computing blowups. ...
1
vote
0answers
51 views

Lower bounds on possible integer relations from the PSLQ algorithm

For the equation: $$ \sum_{i=1}^na_ix_i=0 $$ where all $x_i$ are real numbers and all $a_i$ are integers, the PSLQ algorithm can either find an integer relation or give lower bounds on the norm of ...
1
vote
0answers
75 views

Is there any efficient progam or software to calculate the fractional chromatic number?

The fractional chromatic number $\chi_f(G)$ is a generation of the chromatic number of a graph $G$. It can be formulated as a linear programming question: Let $\mathcal{I}(G)$ be the set of all ...
1
vote
1answer
102 views

How to detect carry and overflow?

Let's say that 0011 + 0111 = 1010 How to detect whether the operation generate carry and overflow?
1
vote
0answers
48 views

Recent benchmarks of CAS programs? (HNF, Groebner)

I've been trying to find recent benchmarks of CAS systems for computing the Hermite Normal Form as well as Groebner basis over both $\mathbb{Q}$ and $\mathbb{C}$. I've been unable to find anything ...
1
vote
2answers
109 views

Boolean Expression simplification help

Hi I am new to the board. Taking a computer architecture course and I am having trouble understanding further simplification on a question I got on a previous quiz. When I type in the expression ...