Items tagged with algebra algebra Tagged Items Feed

See the attached worksheet: myDet.mw

Why does the determinant of A  contain  x, x^3, x^5, x^7 and x^9 ?

Maple15/Linux/Fedora16_x86_64

Hello

I am using Maple 15 on a Macbook Pro, and I was using this simple code

   restart;with(LinearAlgebra):
     P := Matrix(4, 4,    [[1, s, t, 0], [0, 1, s, 0], [0, 0, 1, s], [0, 0, 0, s]]);
DiagonalMatrix(Column(P, 2)).P;

which should be easy but then Maple gives me a pop-up message:

"Kernel connection has been lost.  The kernel returned the following error message:

I'm trying to write a program that solves sudoku's using a Groebner basis. I introduced 81 variables x₁ to x₈₁, this is a linearisation of the sudoku board.

The space of valid sudokus is defined by:

for i=1,…,81 : F_i=(x_i−1)(x_i−2)⋯(x_i−9) This represents the fact that all squares have integer values between 1 and 9.

for all x_i and x_j which...

I am trying to compute commutator like this [d_t +J_x,d_x+J_t], where d_t and d_x are differential operators, J_x and J_t, which contains generators of group, are some currents on group manifold. terms like this [J_x,J_t] could sovled by physics package, but how to manipulate [d_t,J_x]? How to define operator d_t in order to make the commutator works? and if it works we might have another problem as [d_t,J_x]=[d_t, J^a_x*T_a], where T_a are generators, the result is T_a *d_t(J^a_x...

Hello, I'm Tam and I had a problem with Conditional of Jacobi Iteration.

So, the proplem is:

>restart;

>with(LinearAlgebra):

A:=Matrix([[5,1,-1],[0,4,2],[-1,2,-7]]);

And I want to check the condition

with the maple program.

Fisrt, I typed

b:=sum(abs('A'[1,j]),j=1..ColumnDimension(A)); and I got b:=7

Hello,

I would like to write some functions and algebra that work with dual numbers.

http://en.wikipedia.org/wiki/Dual_number 

http://en.wikipedia.org/wiki/Dual_quaternion 

I have not found a library that supports this.

The basis of dual numbers it is epsilon^2=0, similar to i,j,k in complex algebra where i^2=-1,i*j=-1 etc.

Does anyone happen to have a worksheet that implements LinearAlgebra[Generic][SmithForm] 

for Euclidean domains like Z[w], w^3 = 1, or Z[sqrt(d)] where d = −11, −7, −3, −2, −1, 2, 3, 5, 6, 7, 11, 13, 17, 19, 21, 29, 33, 37, 41, 57, 73?

If so I would greatly appreciate receiving a copy of it.

--Edwin Clark

i find that modulo is (Im(M)+Im(N))/Im(N) but where is Im command in maple?

after tried the following package, also not correct, i guess i should use quotient and Im command

with(MTM):
with(LinearAlgebra[Modular]):

Hom := proc(M, N)
Mdim := Size(M);
M_n := Mdim[1];
M_m := Mdim[2];
Ndim := Size(N);
N_n := Ndim[1];
N_m := Ndim[2];
F := kontraHom(M,N_n);
B := kohom(N,M_m);
C := kohom(N,M_n);
D := Mod(F, B);
E := Mod(D, C);

u = -1/2*ln(-2*t*x+x^2+1)/x;

what is command to find x in terms of u of above formula? i.e. x = u^2+u+1

Hello mathematicians :)

I have a problem finding the inverse of the following 12x12 matrix:

C:=Matrix([

[ 1, x1, y1,    x1^2, x1*y1, y1^2,      x1^3,    x1^2*y1, x1*y1^2,   y1^3,      x1^3*y1,   x1*y1^3],
[ 0,  0,  1,       0,    x1, 2*y1,         0,      ...

Hi all,

 I am trying to solve a system of equations- ode's and nonliear algebraic equations numerically. The maple dsolve actually solves this type of system but the problem is it is taking lot of time to solve. I think it is not solving the algebraic equation numerically.

i am using dsolve as: dsolve({eqs,ic},type=numeric,output=listprocedure). 

Is there any other way to specify maple dsolve to solve all equations including algebraic numerically?  

Answering to that question, I posted several procedures finding minimal polynomials for the elements of finite fields. The best one was the following,

alias(a=RootOf(T^100+T^97+T^96+T^93+T^91+T^89+T^87+T^86+T^82+T^81+T^71+T^70+T^67+T^61+
T^60+T^57+T^54+T^53+T^52+T^49+T^48+T^45+T^44+T^42+T^39+T^36+T^33+T^32+T^31+T^29+T^28+T^27+
T^26+T^24+T^23+T^22+T^18+T^17+T^16+T^14+T^13+T^12+T^10+T^8+T^7+T^6+T^3+T+1)):

F:=GF(2,100,op(a)):
z:=F:-input(2):

MinPolyGF:=proc(x,y:=_X)
local A, i;
A:=Matrix(100,...