$SPAD/src/input quat1.input

Timothy Daly December 3, 2016

Abstract

Contents

2

1 License

— license —

--Copyright The Numerical Algorithms Group Limited 1991.

———-

— * —

)set break resume
)spool quat1.output
)set message test on
)set message auto off
)clear all
--S 1 of 11
q := quatern(2/11,-8,3/4,1)
--R
--R
--R         2        3
--R   (1)  -- - 8i + - j + k
--R        11        4
--R                                          Type: Quaternion(Fraction(Integer))
--E 1
--S 2 of 11
[real q, imagI q, imagJ q, imagK q]
--R
--R
--R          2     3
--R   (2)  [--,- 8,-,1]
--R         11     4
--R                                                Type: List(Fraction(Integer))
--E 2
--S 3 of 11
inv q
--R
--R
--R          352     15488      484       1936
--R   (3)  ------ + ------ i - ----- j - ------ k
--R        126993   126993     42331     126993
--R                                          Type: Quaternion(Fraction(Integer))
--E 3
--S 4 of 11
q**6
--R
--R
--R          2029490709319345   48251690851     144755072553     48251690851

--R   (4)  - ---------------- - ----------- i + ------------ j + ----------- k
--R             7256313856        1288408         41229056         10307264
--R                                          Type: Quaternion(Fraction(Integer))
--E 4
--S 5 of 11
r := quatern(-2,3,23/9,-89); q + r
--R
--R
--R          20        119
--R   (5)  - -- - 5i + --- j - 88k
--R          11         36
--R                                          Type: Quaternion(Fraction(Integer))
--E 5
--S 6 of 11
q * r - r * q
--R
--R
--R          2495             817
--R   (6)  - ---- i - 1418j - --- k
--R           18               18
--R                                          Type: Quaternion(Fraction(Integer))
--E 6
--S 7 of 11
i:=quatern(0,1,0,0); j:=quatern(0,0,1,0); k:=quatern(0,0,0,1)
--R
--R
--R   (7)  k
--R                                                    Type: Quaternion(Integer)
--E 7
--S 8 of 11
[i*i, j*j, k*k, i*j, j*k, k*i, q*i]
--R
--R
--R                                2         3
--R   (8)  [- 1,- 1,- 1,k,i,j,8 + -- i + j - - k]
--R                               11         4
--R                                    Type: List(Quaternion(Fraction(Integer)))
--E 8
--S 9 of 11
norm q
--R
--R
--R        126993
--R   (9)  ------
--R         1936

--R                                                      Type: Fraction(Integer)
--E 9
--S 10 of 11
conjugate q
--R
--R
--R          2        3
--R   (10)  -- + 8i - - j - k
--R         11        4
--R                                          Type: Quaternion(Fraction(Integer))
--E 10
--S 11 of 11
q * %
--R
--R
--R         126993
--R   (11)  ------
--R          1936
--R                                          Type: Quaternion(Fraction(Integer))
--E 11
)spool
)lisp (bye)

———-

References^{[1] nothing}

6