login  home  contents  what's new  discussion  bug reports     help  links  subscribe  changes  refresh  edit
FriCAS Problems
List To Matrix  ← ↑ →  Mutual Recursion

Manipulating Expressions
  
last edited 12 years ago by test1

Edit detail for Manipulating Expressions revision 1 of 6
1 2 3 4 5 6
Editor: 127.0.0.1
Time: 2007/11/15 20:15:15 GMT-8
Note: transferred from axiom-developer 
changed:
-
The domain InputForm can be quite useful for manipulating parts
of expressions. For example

\begin{axiom}
ex1:=integrate(log(x)+x, x)
%::InputForm
ex2:=interpret((%::InputForm).2.2)
\end{axiom}

If you would like to do this with a more common type of
expression and hide the details, you can define
\begin{axiom}
op(n,x) == interpret((x::InputForm).(n+1))
\end{axiom}

Then manipulating expressions looks like this:
\begin{axiom}
op(1,op(1,ex1))
(op(1,op(1,ex1))-op(2,op(1,ex1)))/op(2,ex1)

\end{axiom}

Other Operations

  The domain FunctionSpace includes the following operations::

    isExpt(p,f:Symbol) returns [x, n] if p = x**n and n <> 0 and x = f(a)
    isExpt(p,op:BasicOperator) returns [x, n] if p = x**n and n <> 0 and x = op(a)
    isExpt(p) returns [x, n] if p = x**n and n <> 0
    isMult(p) returns [n, x] if p = n * x and n <> 0
    isPlus(p) returns [m1,...,mn] if p = m1 +...+ mn and n > 1
    isPower(p) returns [x, n] if p = x**n and n <> 0
    isTimes(p) returns [a1,...,an] if p = a1*...*an and n > 1

  If these conditions are not met, then the above operations
  return "failed".

For example,
\begin{axiom}
isMult(3*x)
\end{axiom}
but
\begin{axiom}
isMult(x*y)
\end{axiom}

In the context of 'Expression Integer', or 'Polynomial Integer'
the parameter n must be an Integer. The Symbol 'y' is not an
'Integer'.

Not exactly analogously
\begin{axiom}
isPower(x**y)
\end{axiom}
whereas
\begin{axiom} 
isPower(x**10)
\end{axiom}
In the first case the 'Integer' is assume to be 1.

We have:
\begin{axiom} 
isTimes(x*y*z)
isPlus(x+y+z*y)
\end{axiom}

Whereas 
\begin{axiom}
isTimes((x+y)*z)
\end{axiom}
That is because the expression is internally treated as a
'MultivariatePolynomial' like this:
\begin{axiom}
((x+y)*z)::MPOLY([x,y,z],INT)
\end{axiom}

If you say:
\begin{axiom}
isPlus((x+y)*z)
\end{axiom}
perhaps the result makes sense?

For some of the details of these operations I consulted the
actual algebra code at:

http://wiki.axiom-developer.org/axiom--test--1/src/algebra/FspaceSpad

Click on 'pdf' or 'dvi' to see the documentation.

You can also enter expressions like 'isTimes' in the search box
on the upper right and see all the places in the algebra where
this operation is defined and used.



The domain InputForm can be quite useful for manipulating parts
of expressions. For example


fricas
(1) -> ex1:=integrate(log(x)+x, x)

\begin{equation}
\label{eq1}\frac{{2 \  x \ {\log \left({x}\right)}}+{{x}^{2}}-{2 \  x}}{2}\end{equation}

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


fricas
%::InputForm

\begin{equation}
\label{eq2}\left(/ \ {\left(+ \ {\left( \ {\left( \  2 \  x \right)}\ {\left(\log \  x \right)}\right)}\ {\left(+ \ {\left(^\  x \  2 \right)}\ {\left( \  - 2 \  x \right)}\right)}\right)}\  2 \right)\end{equation*}

Type: InputForm?


fricas
ex2:=interpret((%::InputForm).2.2)

\begin{equation}
\label{eq3}2 \  x \ {\log \left({x}\right)}\end{equation}

Type: Expression(Integer)


If you would like to do this with a more common type of
expression and hide the details, you can define

fricas
op(n,x) == interpret((x::InputForm).(n+1))


Type: Void


Then manipulating expressions looks like this:

fricas
op(1,op(1,ex1))


fricas
Compiling function op with type (PositiveInteger, Expression(Integer
      )) -> Any


fricas
Compiling function op with type (PositiveInteger, Any) -> Any

\begin{equation}
\label{eq4}2 \  x \ {\log \left({x}\right)}\end{equation}

Type: Expression(Integer)


fricas
(op(1,op(1,ex1))-op(2,op(1,ex1)))/op(2,ex1)

\begin{equation}
\label{eq5}\frac{{2 \  x \ {\log \left({x}\right)}}-{{x}^{2}}+{2 \  x}}{2}\end{equation}

Type: Expression(Integer)


Other Operations


  The domain FunctionSpace? includes the following operations:

    isExpt(p,f:Symbol) returns [x, n] if p = x**n and n <> 0 and x = f(a)
    isExpt(p,op:BasicOperator) returns [x, n] if p = x**n and n <> 0 and x = op(a)
    isExpt(p) returns [x, n] if p = x**n and n <> 0
    isMult(p) returns [n, x] if p = n * x and n <> 0
    isPlus(p) returns [m1,...,mn] if p = m1 +...+ mn and n > 1
    isPower(p) returns [x, n] if p = x**n and n <> 0
    isTimes(p) returns [a1,...,an] if p = a1*...*an and n > 1





  If these conditions are not met, then the above operations
  return "failed".


For example,

fricas
isMult(3*x)

   There are 1 exposed and 0 unexposed library operations named isMult 
      having 1 argument(s) but none was determined to be applicable. 
      Use HyperDoc Browse, or issue
                             )display op isMult
      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 
      isMult with argument type(s) 
                             Polynomial(Integer)

      Perhaps you should use "@" to indicate the required return type, 
      or "$" to specify which version of the function you need.


but

fricas
isMult(x*y)

   There are 1 exposed and 0 unexposed library operations named isMult 
      having 1 argument(s) but none was determined to be applicable. 
      Use HyperDoc Browse, or issue
                             )display op isMult
      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 
      isMult with argument type(s) 
                             Polynomial(Integer)

      Perhaps you should use "@" to indicate the required return type, 
      or "$" to specify which version of the function you need.


In the context of Expression Integer, or Polynomial Integer
the parameter n must be an Integer. The Symbol y is not an
Integer.


Not exactly analogously

fricas
isPower(x**y)

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

   Cannot find a definition or applicable library operation named ** 
      with argument type(s) 
                                 Variable(x)
                                 Variable(y)

      Perhaps you should use "@" to indicate the required return type, 
      or "$" to specify which version of the function you need.


whereas

fricas
isPower(x**10)

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

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

      Perhaps you should use "@" to indicate the required return type, 
      or "$" to specify which version of the function you need.


In the first case the Integer is assume to be 1.


We have:

fricas
isTimes(x*y*z)

\begin{equation*}
\label{eq6}\left[ z , \: y , \: x \right]?\end{equation*}

Type: Union(List(Polynomial(Integer)),...)


fricas
isPlus(x+y+z*y)

\begin{equation*}
\label{eq7}\left[{y \  z}, \: y , \: x \right]?\end{equation*}

Type: Union(List(Polynomial(Integer)),...)


Whereas 

fricas
isTimes((x+y)*z)

\begin{equation}
\label{eq8}\verb#"failed"#\end{equation}

Type: Union("failed",...)


That is because the expression is internally treated as a
MultivariatePolynomial like this:

fricas
((x+y)*z)::MPOLY([x,y,z],INT)

\begin{equation}
\label{eq9}{z \  x}+{z \  y}\end{equation}

Type: MultivariatePolynomial?([x,y,z],Integer)


If you say:

fricas
isPlus((x+y)*z)

\begin{equation*}
\label{eq10}\left[{y \  z}, \:{x \  z}\right]?\end{equation*}

Type: Union(List(Polynomial(Integer)),...)


perhaps the result makes sense?


For some of the details of these operations I consulted the
actual algebra code at:


http://wiki.axiom-developer.org/axiom--test--1/src/algebra/FspaceSpad


Click on pdf or dvi to see the documentation.


You can also enter expressions like isTimes in the search box
on the upper right and see all the places in the algebra where
this operation is defined and used.
Some or all expressions may not have rendered properly, because Latex returned the following error:
This is pdfTeX, Version 3.141592653-2.6-1.40.24 (TeX Live 2022/Debian) (preloaded format=latex)
 restricted \write18 enabled.
entering extended mode
(./7270444483375878928-16.0px.tex
LaTeX2e <2022-11-01> patch level 1
L3 programming layer <2023-01-16>
(/usr/share/texlive/texmf-dist/tex/latex/base/article.cls
Document Class: article 2022/07/02 v1.4n Standard LaTeX document class
(/usr/share/texlive/texmf-dist/tex/latex/base/size12.clo))
(/usr/share/texlive/texmf-dist/tex/latex/ucs/ucs.sty
(/usr/share/texlive/texmf-dist/tex/latex/ucs/data/uni-global.def))
(/usr/share/texlive/texmf-dist/tex/latex/base/inputenc.sty
(/usr/share/texlive/texmf-dist/tex/latex/ucs/utf8x.def))
(/usr/share/texlive/texmf-dist/tex/latex/bbm-macros/bbm.sty)
(/usr/share/texlive/texmf-dist/tex/latex/jknapltx/mathrsfs.sty)
(/usr/share/texlive/texmf-dist/tex/latex/base/fontenc.sty)
(/usr/share/texlive/texmf-dist/tex/latex/pstricks/pstricks.sty
(/usr/share/texlive/texmf-dist/tex/generic/iftex/iftex.sty)
(/usr/share/texlive/texmf-dist/tex/latex/xcolor/xcolor.sty
(/usr/share/texlive/texmf-dist/tex/latex/graphics-cfg/color.cfg)
(/usr/share/texlive/texmf-dist/tex/latex/graphics-def/dvips.def)
(/usr/share/texlive/texmf-dist/tex/latex/graphics/mathcolor.ltx)
(/usr/share/texlive/texmf-dist/tex/latex/graphics/dvipsnam.def))
(/usr/share/texlive/texmf-dist/tex/generic/xkeyval/pst-xkey.tex
(/usr/share/texlive/texmf-dist/tex/latex/xkeyval/xkeyval.sty
(/usr/share/texlive/texmf-dist/tex/generic/xkeyval/xkeyval.tex
(/usr/share/texlive/texmf-dist/tex/generic/xkeyval/xkvutils.tex
(/usr/share/texlive/texmf-dist/tex/generic/xkeyval/keyval.tex)))))
(/usr/share/texlive/texmf-dist/tex/generic/pstricks/pstricks.tex
(/usr/share/texlive/texmf-dist/tex/generic/pstricks/pst-fp.tex
`pst-fp' v0.06, 2020/11/20 (hv))
(/usr/share/texlive/texmf-dist/tex/latex/pgf/utilities/pgffor.sty
(/usr/share/texlive/texmf-dist/tex/latex/pgf/utilities/pgfrcs.sty
(/usr/share/texlive/texmf-dist/tex/generic/pgf/utilities/pgfutil-common.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/utilities/pgfutil-latex.def)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/utilities/pgfrcs.code.tex
(/usr/share/texlive/texmf-dist/tex/generic/pgf/pgf.revision.tex)))
(/usr/share/texlive/texmf-dist/tex/latex/pgf/utilities/pgfkeys.sty
(/usr/share/texlive/texmf-dist/tex/generic/pgf/utilities/pgfkeys.code.tex
(/usr/share/texlive/texmf-dist/tex/generic/pgf/utilities/pgfkeyslibraryfiltered
.code.tex))) (/usr/share/texlive/texmf-dist/tex/latex/pgf/math/pgfmath.sty
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmath.code.tex
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathutil.code.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathparser.code.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.code.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.basic.code
.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.trigonomet
ric.code.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.random.cod
e.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.comparison
.code.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.base.code.
tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.round.code
.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.misc.code.
tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.integerari
thmetics.code.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathcalc.code.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pgf/math/pgfmathfloat.code.tex)))
(/usr/share/texlive/texmf-dist/tex/generic/pgf/utilities/pgffor.code.tex))
`PSTricks' v3.18  <2022/11/28> (tvz,hv)
--- We are running latex or xelatex ---
(/usr/share/texlive/texmf-dist/tex/generic/pstricks/pstricks.con)
(/usr/share/texlive/texmf-dist/tex/generic/pstricks/pstricks-color.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pstricks/pstricks-arrows.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pstricks/pstricks-dots.tex)
(/usr/share/texlive/texmf-dist/tex/generic/pstricks/pstricks.con))
(/usr/share/texlive/texmf-dist/tex/generic/pstricks/pst-fp.tex
`pst-fp' v0.06, 2020/11/20 (hv)))
(/usr/share/texlive/texmf-dist/tex/latex/graphics/epsfig.sty
(/usr/share/texlive/texmf-dist/tex/latex/graphics/graphicx.sty
(/usr/share/texlive/texmf-dist/tex/latex/graphics/graphics.sty
(/usr/share/texlive/texmf-dist/tex/latex/graphics/trig.sty)
(/usr/share/texlive/texmf-dist/tex/latex/graphics-cfg/graphics.cfg))))
(/usr/share/texlive/texmf-dist/tex/latex/pst-grad/pst-grad.sty
(/usr/share/texlive/texmf-dist/tex/generic/pst-grad/pst-grad.tex
`pst-grad' v1.06, 2006/11/27 (tvz,dg,hv)))
(/usr/share/texlive/texmf-dist/tex/latex/pst-plot/pst-plot.sty
(/usr/share/texlive/texmf-dist/tex/latex/xkeyval/pst-xkey.sty)
(/usr/share/texlive/texmf-dist/tex/latex/multido/multido.sty
(/usr/share/texlive/texmf-dist/tex/generic/multido/multido.tex
 v1.42, 2010/05/14 <tvz>))
(/usr/share/texlive/texmf-dist/tex/generic/pst-plot/pst-plot.tex
(/usr/share/texlive/texmf-dist/tex/generic/pst-tools/pst-tools.tex
`PST-tools' v0.12, 2021/09/23 (hv))
(/usr/share/texlive/texmf-dist/tex/generic/pstricks-add/pstricks-add.tex
(/usr/share/texlive/texmf-dist/tex/generic/pst-node/pst-node.tex
 v1.43, 2022/01/31)
(/usr/share/texlive/texmf-dist/tex/generic/pst-arrow/pst-arrow.tex
`pst-arrow' v0.05, 2021/11/16 (dr,hv))
(/usr/share/texlive/texmf-dist/tex/generic/pst-3d/pst-3d.tex
`PST-3d' v1.11, 2010/02/14 (tvz))
(/usr/share/texlive/texmf-dist/tex/generic/pst-math/pst-math.tex
`pst-math' v0.66 , (CJ,hv)) `pstricks-add' v3.93, 2022/11/21 (dr,hv))
 v1.94, 2022/11/21 (tvz,hv)))
(/usr/share/texlive/texmf-dist/tex/latex/geometry/geometry.sty
(/usr/share/texlive/texmf-dist/tex/generic/iftex/ifvtex.sty)


Package geometry Warning: `lmargin' and `rmargin' result in NEGATIVE (-108.405p
t).
    `width' should be shortened in length.


) (/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsmath.sty
For additional information on amsmath, use the `?' option.
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amstext.sty
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsgen.sty))
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsbsy.sty)
(/usr/share/texlive/texmf-dist/tex/latex/amsmath/amsopn.sty))
(/usr/share/texlive/texmf-dist/tex/latex/amsfonts/amsfonts.sty)
(/usr/share/texlive/texmf-dist/tex/latex/amsfonts/amssymb.sty)
(/usr/share/texlive/texmf-dist/tex/latex/amscls/amsthm.sty)
(/usr/share/texlive/texmf-dist/tex/latex/setspace/setspace.sty)
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xy.sty
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xy.tex Bootstrap'ing:
catcodes, docmode, (/usr/share/texlive/texmf-dist/tex/generic/xypic/xyrecat.tex
) (/usr/share/texlive/texmf-dist/tex/generic/xypic/xyidioms.tex)
 Xy-pic version 3.8.9 <2013/10/06>
 Copyright (c) 1991-2013 by Kristoffer H. Rose <krisrose@tug.org> and others
 Xy-pic is free software: see the User's Guide for details.




Loading kernel: messages; fonts; allocations: state, direction,
utility macros; pictures: \xy, positions, objects, decorations;
kernel objects: directionals, circles, text; options; algorithms: directions,
edges, connections;  Xy-pic loaded)
(/usr/share/texlive/texmf-dist/tex/generic/iftex/ifpdf.sty))
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xyall.tex
 Xy-pic option: All features v.3.8
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xycurve.tex
 Xy-pic option: Curve and Spline extension v.3.12 curve, circles, loaded)
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xyframe.tex
 Xy-pic option: Frame and Bracket extension v.3.14 loaded)
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xycmtip.tex
 Xy-pic option: Computer Modern tip extension v.3.7
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xytips.tex
 Xy-pic option: More Tips extension v.3.11 loaded) loaded)
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xyline.tex
 Xy-pic option: Line styles extension v.3.10 loaded)
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xyrotate.tex
 Xy-pic option: Rotate and Scale extension v.3.8 loaded)
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xycolor.tex
 Xy-pic option: Colour extension v.3.11 loaded)
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xymatrix.tex
 Xy-pic option: Matrix feature v.3.14 loaded)
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xyarrow.tex
 Xy-pic option: Arrow and Path feature v.3.9 path, \ar, loaded)
(/usr/share/texlive/texmf-dist/tex/generic/xypic/xygraph.tex
 Xy-pic option: Graph feature v.3.11 loaded) loaded)
(/usr/share/texlive/texmf-dist/tex/latex/tools/verbatim.sty)
(/usr/share/texlive/texmf-dist/tex/latex/graphviz/graphviz.sty
(/usr/share/texlive/texmf-dist/tex/latex/psfrag/psfrag.sty))
(/usr/share/texmf/tex/latex/sagetex.sty
Writing sage input file 7270444483375878928-16.0px.sage
) (/usr/share/texlive/texmf-dist/tex/latex/gnuplottex/gnuplottex.sty
(/usr/share/texlive/texmf-dist/tex/latex/moreverb/moreverb.sty)
(/usr/share/texlive/texmf-dist/tex/latex/base/ifthen.sty)
(/usr/share/texlive/texmf-dist/tex/generic/catchfile/catchfile.sty
(/usr/share/texlive/texmf-dist/tex/generic/infwarerr/infwarerr.sty)
(/usr/share/texlive/texmf-dist/tex/generic/ltxcmds/ltxcmds.sty)
(/usr/share/texlive/texmf-dist/tex/generic/etexcmds/etexcmds.sty))


Package gnuplottex Warning: Shell escape not enabled.
(gnuplottex)                You'll need to convert the graphs yourself.


) (/usr/share/texlive/texmf-dist/tex/latex/l3backend/l3backend-dvips.def)
No file 7270444483375878928-16.0px.aux.
(/usr/share/texlive/texmf-dist/tex/latex/ucs/ucsencs.def)
geometry driver: auto-detecting
geometry detected driver: dvips
(/usr/share/texlive/texmf-dist/tex/latex/jknapltx/ursfs.fd)
(/usr/share/texlive/texmf-dist/tex/latex/amsfonts/umsa.fd)
(/usr/share/texlive/texmf-dist/tex/latex/amsfonts/umsb.fd)


Package amsmath Warning: Foreign command \over;
(amsmath)                \frac or \genfrac should be used instead
(amsmath)                 on input line 123.


[1]
 Missing { inserted.
<to be read again> 
                   \ 
l.125 ...x \right)}\right)}\ {\left(+ \ {\left(^\ 
                                                   x \  2 \right)}\ {\left(*...
 Missing } inserted.
<inserted text> 
                }
l.125 ...ht)}\ {\left(+ \ {\left(^\  x \  2 \right
                                                  )}\ {\left(* \  - 2 \  x \...




[2] [3] [4] [5] [6] [7] (/usr/share/texlive/texmf-dist/tex/latex/base/t1cmtt.fd
) [8] [9] [10] (./7270444483375878928-16.0px.aux)


LaTeX Warning: Label(s) may have changed. Rerun to get cross-references right.


 )
(see the transcript file for additional information)
Output written on 7270444483375878928-16.0px.dvi (10 pages, 3544 bytes).
Transcript written on 7270444483375878928-16.0px.log.