help links subscribe changes refresh edit

	SandBox
SandBoxProp ←	↑	→ SandBoxRealSpace

SandBoxREN

Edit detail for SandBoxREN revision 1 of 2

	1 2
Editor: Time: 2007/11/18 18:33:04 GMT-8
Note:

changed:
-
\begin{axiom}
As := matrix([ [-3,1,1,1], [1,1,1,1], [1,1,1,1], [1,1,1,1]])
A := subMatrix(As, 2,4,2,4)
ob := orthonormalBasis(A)
P : Matrix(Expression Integer) := new(3,3,0)
setsubMatrix!(P,1,1,ob.3) 
setsubMatrix!(P,1,2,ob.1) 
setsubMatrix!(P,1,3,ob.2)
Pt := transpose(P)
Ps : Matrix(Expression Integer) := new(4,4,0)
Ps(1,1) := 1
setsubMatrix!(Ps,2,2,P)
PsT := transpose(Ps)
PsTAsPs := PsT * As * Ps
b1 := PsTAsPs(2,1)
l1 := PsTAsPs(2,2)
Us : Matrix(Expression Integer) := new(4,4,0)
Us(1,1) := 1
Us(2,2) := 1
Us(3,3) := 1
Us(4,4) := 1
Us(2,1) := -b1 / l1
Us
PsUs := Ps * Us
PsUsT := transpose(PsUs)
PsUsTAsPsUs := PsUsT * As * PsUs
C := inverse(PsUs) 
c := PsUsTAsPsUs(1,1)
gQ := PsUsTAsPsUs / c 
x1 := transpose(matrix([[1,2,3,4]]))
v1 := transpose(x1) * As * x1
x2 := C * x1
v2 := transpose(x2) * PsUsTAsPsUs * x2
\end{axiom}


From ren Fri Jul 7 14:33:03 -0500 2006
From: ren
Date: Fri, 07 Jul 2006 14:33:03 -0500
Subject: 
Message-ID: <20060707143303-0500@wiki.axiom-developer.org>

\begin{axiom}
As := matrix([ [-c,-1,-2,-1], [-1,3,-1,0], [-2,-1,3,0], [-1,0,0,-6]])
A := subMatrix(As, 2,4,2,4)
ob := orthonormalBasis(A)
P : Matrix(Expression Integer) := new(3,3,0)
setsubMatrix!(P,1,1,ob.3) 
setsubMatrix!(P,1,2,ob.1) 
setsubMatrix!(P,1,3,ob.2)
Pt := transpose(P)
Ps : Matrix(Expression Integer) := new(4,4,0)
Ps(1,1) := 1
setsubMatrix!(Ps,2,2,P)
PsT := transpose(Ps)
PsTAsPs := PsT * As * Ps
Us : Matrix(Expression Integer) := new(4,4,0)
Us(1,1) := 1
Us(2,2) := 1
Us(3,3) := 1
Us(4,4) := 1
Us(2,1) := -PsTAsPs(2,1) / PsTAsPs(2,2)
Us(3,1) := -PsTAsPs(3,1) / PsTAsPs(3,3)
Us(4,1) := -PsTAsPs(4,1) / PsTAsPs(4,4)
Us
PsUs := Ps * Us
PsUsT := transpose(PsUs)
PsUsTAsPsUs := PsUsT * As * PsUs 
cc := PsUsTAsPsUs(1,1)
so := solve(cc = 0, c)
c0 := rhs so.1
gQ := PsUsTAsPsUs / cc
eval(PsUsTAsPsUs, c = c0)
\end{axiom}


From unknown Fri Jul 7 14:38:35 -0500 2006
From: unknown
Date: Fri, 07 Jul 2006 14:38:35 -0500
Subject: 
Message-ID: <20060707143835-0500@wiki.axiom-developer.org>

7.5

axiom

As := matrix([ [-3,1,1,1], [1,1,1,1], [1,1,1,1], [1,1,1,1]])

$\label{eq1}\left[ \begin{array}{cccc} - 3 & 1 & 1 & 1 \ 1 & 1 & 1 & 1 \ 1 & 1 & 1 & 1 \ 1 & 1 & 1 & 1$

(1)

Type: Matrix(Integer)

axiom

A := subMatrix(As, 2,4,2,4)

$\label{eq2}\left[ \begin{array}{ccc} 1 & 1 & 1 \ 1 & 1 & 1 \ 1 & 1 & 1$

(2)

Type: Matrix(Integer)

axiom

ob := orthonormalBasis(A)

$\label{eq3}\begin{array}{@{}l} \displaystyle \left[{\left[ \begin{array}{c} -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}} \ {{\sqrt{2}}\over{\sqrt{3}}} \ -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}$

(3)

Type: List(Matrix(Expression(Integer)))

axiom

P : Matrix(Expression Integer) := new(3,3,0)

$\label{eq4}\left[ \begin{array}{ccc} 0 & 0 & 0 \ 0 & 0 & 0 \ 0 & 0 & 0$

(4)

Type: Matrix(Expression(Integer))

axiom

setsubMatrix!(P,1,1,ob.3)

$\label{eq5}\left[ \begin{array}{ccc} {1 \over{\sqrt{3}}}& 0 & 0 \ {1 \over{\sqrt{3}}}& 0 & 0 \ {1 \over{\sqrt{3}}}& 0 & 0$

(5)

Type: Matrix(Expression(Integer))

axiom

setsubMatrix!(P,1,2,ob.1)

$\label{eq6}\left[ \begin{array}{ccc} {1 \over{\sqrt{3}}}& -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}& 0 \ {1 \over{\sqrt{3}}}&{{\sqrt{2}}\over{\sqrt{3}}}& 0 \ {1 \over{\sqrt{3}}}& -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}& 0$

(6)

Type: Matrix(Expression(Integer))

axiom

setsubMatrix!(P,1,3,ob.2)

$\label{eq7}\left[ \begin{array}{ccc} {1 \over{\sqrt{3}}}& -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}& -{1 \over{\sqrt{2}}} \ {1 \over{\sqrt{3}}}&{{\sqrt{2}}\over{\sqrt{3}}}& 0 \ {1 \over{\sqrt{3}}}& -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}&{1 \over{\sqrt{2}}}$

(7)

Type: Matrix(Expression(Integer))

axiom

Pt := transpose(P)

$\label{eq8}\left[ \begin{array}{ccc} {1 \over{\sqrt{3}}}&{1 \over{\sqrt{3}}}&{1 \over{\sqrt{3}}} \ -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}&{{\sqrt{2}}\over{\sqrt{3}}}& -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}} \ -{1 \over{\sqrt{2}}}& 0 &{1 \over{\sqrt{2}}}$

(8)

Type: Matrix(Expression(Integer))

axiom

Ps : Matrix(Expression Integer) := new(4,4,0)

$\label{eq9}\left[ \begin{array}{cccc} 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0$

(9)

Type: Matrix(Expression(Integer))

axiom

Ps(1,1) := 1

$\label{eq10}1$

(10)

Type: Expression(Integer)

axiom

setsubMatrix!(Ps,2,2,P)

$\label{eq11}\left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \ 0 &{1 \over{\sqrt{3}}}& -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}& -{1 \over{\sqrt{2}}} \ 0 &{1 \over{\sqrt{3}}}&{{\sqrt{2}}\over{\sqrt{3}}}& 0 \ 0 &{1 \over{\sqrt{3}}}& -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}&{1 \over{\sqrt{2}}}$

(11)

Type: Matrix(Expression(Integer))

axiom

PsT := transpose(Ps)

$\label{eq12}\left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \ 0 &{1 \over{\sqrt{3}}}&{1 \over{\sqrt{3}}}&{1 \over{\sqrt{3}}} \ 0 & -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}&{{\sqrt{2}}\over{\sqrt{3}}}& -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}} \ 0 & -{1 \over{\sqrt{2}}}& 0 &{1 \over{\sqrt{2}}}$

(12)

Type: Matrix(Expression(Integer))

axiom

PsTAsPs := PsT * As * Ps

$\label{eq13}\left[ \begin{array}{cccc} - 3 &{3 \over{\sqrt{3}}}& 0 & 0 \ {3 \over{\sqrt{3}}}& 3 & 0 & 0 \ 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0$

(13)

Type: Matrix(Expression(Integer))

axiom

b1 := PsTAsPs(2,1)

$\label{eq14}3 \over{\sqrt{3}}$

(14)

Type: Expression(Integer)

axiom

l1 := PsTAsPs(2,2)

$\label{eq15}3$

(15)

Type: Expression(Integer)

axiom

Us : Matrix(Expression Integer) := new(4,4,0)

$\label{eq16}\left[ \begin{array}{cccc} 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0$

(16)

Type: Matrix(Expression(Integer))

axiom

Us(1,1) := 1

$\label{eq17}1$

(17)

Type: Expression(Integer)

axiom

Us(2,2) := 1

$\label{eq18}1$

(18)

Type: Expression(Integer)

axiom

Us(3,3) := 1

$\label{eq19}1$

(19)

Type: Expression(Integer)

axiom

Us(4,4) := 1

$\label{eq20}1$

(20)

Type: Expression(Integer)

axiom

Us(2,1) := -b1 / l1

$\label{eq21}-{1 \over{\sqrt{3}}}$

(21)

Type: Expression(Integer)

axiom

Us

$\label{eq22}\left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \ -{1 \over{\sqrt{3}}}& 1 & 0 & 0 \ 0 & 0 & 1 & 0 \ 0 & 0 & 0 & 1$

(22)

Type: Matrix(Expression(Integer))

axiom

PsUs := Ps * Us

$\label{eq23}\left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \ -{1 \over 3}&{1 \over{\sqrt{3}}}& -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}& -{1 \over{\sqrt{2}}} \ -{1 \over 3}&{1 \over{\sqrt{3}}}&{{\sqrt{2}}\over{\sqrt{3}}}& 0 \ -{1 \over 3}&{1 \over{\sqrt{3}}}& -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}&{1 \over{\sqrt{2}}}$

(23)

Type: Matrix(Expression(Integer))

axiom

PsUsT := transpose(PsUs)

$\label{eq24}\left[ \begin{array}{cccc} 1 & -{1 \over 3}& -{1 \over 3}& -{1 \over 3} \ 0 &{1 \over{\sqrt{3}}}&{1 \over{\sqrt{3}}}&{1 \over{\sqrt{3}}} \ 0 & -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}}&{{\sqrt{2}}\over{\sqrt{3}}}& -{{\sqrt{2}}\over{2 \ {\sqrt{3}}}} \ 0 & -{1 \over{\sqrt{2}}}& 0 &{1 \over{\sqrt{2}}}$

(24)

Type: Matrix(Expression(Integer))

axiom

PsUsTAsPsUs := PsUsT * As * PsUs

$\label{eq25}\left[ \begin{array}{cccc} - 4 & 0 & 0 & 0 \ 0 & 3 & 0 & 0 \ 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0$

(25)

Type: Matrix(Expression(Integer))

axiom

C := inverse(PsUs)

$\label{eq26}\left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \ {{\sqrt{3}}\over 3}&{{\sqrt{3}}\over 3}&{{\sqrt{3}}\over 3}&{{\sqrt{3}}\over 3} \ 0 & -{{\sqrt{3}}\over{3 \ {\sqrt{2}}}}&{{2 \ {\sqrt{3}}}\over{3 \ {\sqrt{2}}}}& -{{{\sqrt{2}}\ {\sqrt{3}}}\over 6} \ 0 & -{{\sqrt{2}}\over 2}& 0 &{{\sqrt{2}}\over 2}$

(26)

Type: Union(Matrix(Expression(Integer)),...)

axiom

c := PsUsTAsPsUs(1,1)

$\label{eq27}- 4$

(27)

Type: Expression(Integer)

axiom

gQ := PsUsTAsPsUs / c

$\label{eq28}\left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \ 0 & -{3 \over 4}& 0 & 0 \ 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0$

(28)

Type: Matrix(Expression(Integer))

axiom

x1 := transpose(matrix([[1,2,3,4]]))

$\label{eq29}\left[ \begin{array}{c} 1 \ 2 \ 3 \ 4$

(29)

Type: Matrix(Integer)

axiom

v1 := transpose(x1) * As * x1

$\label{eq30}\left[ \begin{array}{c} {96}$

(30)

Type: Matrix(Integer)

axiom

x2 := C * x1

$\label{eq31}\left[ \begin{array}{c} 1 \ {{{10}\ {\sqrt{3}}}\over 3} \ 0 \ {\sqrt{2}}$

(31)

Type: Matrix(Expression(Integer))

axiom

v2 := transpose(x2) * PsUsTAsPsUs * x2

$\label{eq32}\left[ \begin{array}{c} {96}$

(32)

Type: Matrix(Expression(Integer))

... --ren, Fri, 07 Jul 2006 14:33:03 -0500 reply

axiom

As := matrix([ [-c,-1,-2,-1], [-1,3,-1,0], [-2,-1,3,0], [-1,0,0,-6]])

$\label{eq33}\left[ \begin{array}{cccc} 4 & - 1 & - 2 & - 1 \ - 1 & 3 & - 1 & 0 \ - 2 & - 1 & 3 & 0 \ - 1 & 0 & 0 & - 6$

(33)

Type: Matrix(Expression(Integer))

axiom

A := subMatrix(As, 2,4,2,4)

$\label{eq34}\left[ \begin{array}{ccc} 3 & - 1 & 0 \ - 1 & 3 & 0 \ 0 & 0 & - 6$

(34)

Type: Matrix(Expression(Integer))

axiom

ob := orthonormalBasis(A)

$\label{eq35}\begin{array}{@{}l} \displaystyle \left[{\left[ \begin{array}{c} 0 \ 0 \ 1$

(35)

Type: List(Matrix(Expression(Integer)))

axiom

P : Matrix(Expression Integer) := new(3,3,0)

$\label{eq36}\left[ \begin{array}{ccc} 0 & 0 & 0 \ 0 & 0 & 0 \ 0 & 0 & 0$

(36)

Type: Matrix(Expression(Integer))

axiom

setsubMatrix!(P,1,1,ob.3)

$\label{eq37}\left[ \begin{array}{ccc} -{1 \over{\sqrt{2}}}& 0 & 0 \ {1 \over{\sqrt{2}}}& 0 & 0 \ 0 & 0 & 0$

(37)

Type: Matrix(Expression(Integer))

axiom

setsubMatrix!(P,1,2,ob.1)

$\label{eq38}\left[ \begin{array}{ccc} -{1 \over{\sqrt{2}}}& 0 & 0 \ {1 \over{\sqrt{2}}}& 0 & 0 \ 0 & 1 & 0$

(38)

Type: Matrix(Expression(Integer))

axiom

setsubMatrix!(P,1,3,ob.2)

$\label{eq39}\left[ \begin{array}{ccc} -{1 \over{\sqrt{2}}}& 0 &{1 \over{\sqrt{2}}} \ {1 \over{\sqrt{2}}}& 0 &{1 \over{\sqrt{2}}} \ 0 & 1 & 0$

(39)

Type: Matrix(Expression(Integer))

axiom

Pt := transpose(P)

$\label{eq40}\left[ \begin{array}{ccc} -{1 \over{\sqrt{2}}}&{1 \over{\sqrt{2}}}& 0 \ 0 & 0 & 1 \ {1 \over{\sqrt{2}}}&{1 \over{\sqrt{2}}}& 0$

(40)

Type: Matrix(Expression(Integer))

axiom

Ps : Matrix(Expression Integer) := new(4,4,0)

$\label{eq41}\left[ \begin{array}{cccc} 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0$

(41)

Type: Matrix(Expression(Integer))

axiom

Ps(1,1) := 1

$\label{eq42}1$

(42)

Type: Expression(Integer)

axiom

setsubMatrix!(Ps,2,2,P)

$\label{eq43}\left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \ 0 & -{1 \over{\sqrt{2}}}& 0 &{1 \over{\sqrt{2}}} \ 0 &{1 \over{\sqrt{2}}}& 0 &{1 \over{\sqrt{2}}} \ 0 & 0 & 1 & 0$

(43)

Type: Matrix(Expression(Integer))

axiom

PsT := transpose(Ps)

$\label{eq44}\left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \ 0 & -{1 \over{\sqrt{2}}}&{1 \over{\sqrt{2}}}& 0 \ 0 & 0 & 0 & 1 \ 0 &{1 \over{\sqrt{2}}}&{1 \over{\sqrt{2}}}& 0$

(44)

Type: Matrix(Expression(Integer))

axiom

PsTAsPs := PsT * As * Ps

$\label{eq45}\left[ \begin{array}{cccc} 4 & -{1 \over{\sqrt{2}}}& - 1 & -{3 \over{\sqrt{2}}} \ -{1 \over{\sqrt{2}}}& 4 & 0 & 0 \ - 1 & 0 & - 6 & 0 \ -{3 \over{\sqrt{2}}}& 0 & 0 & 2$

(45)

Type: Matrix(Expression(Integer))

axiom

Us : Matrix(Expression Integer) := new(4,4,0)

$\label{eq46}\left[ \begin{array}{cccc} 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0 \ 0 & 0 & 0 & 0$

(46)

Type: Matrix(Expression(Integer))

axiom

Us(1,1) := 1

$\label{eq47}1$

(47)

Type: Expression(Integer)

axiom

Us(2,2) := 1

$\label{eq48}1$

(48)

Type: Expression(Integer)

axiom

Us(3,3) := 1

$\label{eq49}1$

(49)

Type: Expression(Integer)

axiom

Us(4,4) := 1

$\label{eq50}1$

(50)

Type: Expression(Integer)

axiom

Us(2,1) := -PsTAsPs(2,1) / PsTAsPs(2,2)

$\label{eq51}1 \over{4 \ {\sqrt{2}}}$

(51)

Type: Expression(Integer)

axiom

Us(3,1) := -PsTAsPs(3,1) / PsTAsPs(3,3)

$\label{eq52}-{1 \over 6}$

(52)

Type: Expression(Integer)

axiom

Us(4,1) := -PsTAsPs(4,1) / PsTAsPs(4,4)

$\label{eq53}3 \over{2 \ {\sqrt{2}}}$

(53)

Type: Expression(Integer)

axiom

Us

$\label{eq54}\left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \ {1 \over{4 \ {\sqrt{2}}}}& 1 & 0 & 0 \ -{1 \over 6}& 0 & 1 & 0 \ {3 \over{2 \ {\sqrt{2}}}}& 0 & 0 & 1$

(54)

Type: Matrix(Expression(Integer))

axiom

PsUs := Ps * Us

$\label{eq55}\left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \ {5 \over 8}& -{1 \over{\sqrt{2}}}& 0 &{1 \over{\sqrt{2}}} \ {7 \over 8}&{1 \over{\sqrt{2}}}& 0 &{1 \over{\sqrt{2}}} \ -{1 \over 6}& 0 & 1 & 0$

(55)

Type: Matrix(Expression(Integer))

axiom

PsUsT := transpose(PsUs)

$\label{eq56}\left[ \begin{array}{cccc} 1 &{5 \over 8}&{7 \over 8}& -{1 \over 6} \ 0 & -{1 \over{\sqrt{2}}}&{1 \over{\sqrt{2}}}& 0 \ 0 & 0 & 0 & 1 \ 0 &{1 \over{\sqrt{2}}}&{1 \over{\sqrt{2}}}& 0$

(56)

Type: Matrix(Expression(Integer))

axiom

PsUsTAsPsUs := PsUsT * As * PsUs

$\label{eq57}\left[ \begin{array}{cccc} {{43}\over{24}}& 0 & 0 & 0 \ 0 & 4 & 0 & 0 \ 0 & 0 & - 6 & 0 \ 0 & 0 & 0 & 2$

(57)

Type: Matrix(Expression(Integer))

axiom

cc := PsUsTAsPsUs(1,1)

$\label{eq58}{43}\over{24}$

(58)

Type: Expression(Integer)

axiom

so := solve(cc = 0, c)

   There are 18 exposed and 3 unexposed library operations named solve 
      having 2 argument(s) but none was determined to be applicable. 
      Use HyperDoc Browse, or issue
                              )display op solve
      to learn more about the available operations. Perhaps 
      package-calling the operation or using coercions on the arguments
      will allow you to apply the operation.

   Cannot find a definition or applicable library operation named solve
      with argument type(s) 
                        Equation(Expression(Integer))
                             Expression(Integer)

      Perhaps you should use "@" to indicate the required return type, 
      or "$" to specify which version of the function you need.
c0 := rhs so.1

   There are no library operations named so 
      Use HyperDoc Browse or issue
                                 )what op so
      to learn if there is any operation containing " so " in its name.

   Cannot find a definition or applicable library operation named so 
      with argument type(s) 
                               PositiveInteger

      Perhaps you should use "@" to indicate the required return type, 
      or "$" to specify which version of the function you need.
gQ := PsUsTAsPsUs / cc

$\label{eq59}\left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \ 0 &{{96}\over{43}}& 0 & 0 \ 0 & 0 & -{{144}\over{43}}& 0 \ 0 & 0 & 0 &{{48}\over{43}}$

(59)

Type: Matrix(Expression(Integer))

axiom

eval(PsUsTAsPsUs, c = c0)

   >> Error detected within library code:
   left hand side must be a single kernel

... --unknown, Fri, 07 Jul 2006 14:38:35 -0500 reply

7.5