help links subscribe changes refresh edit

	SandBox
SandBoxEtc ←	↑	→ SandBoxExpr

SandBoxExpOfEnd

Edit detail for SandBoxExpOfEnd revision 7 of 18

	1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
Editor: Bill Page Time: 2014/09/19 23:21:48 GMT+0
Note:

added:
n:=2
eq2_1:= m[X]=(x-r[1])*(x-r[2])

removed:
-\end{axiom}
-\begin{axiom}
-n:=2

added:
  Example 2.1

added:
eval(eq2_1,[r[2]=-r[1]])

added:
eq3_1:= m[X]=(x-r[1])*(x-r[2])*(x-r[3])
eq3_2:= exp(X)=g[0]*id+g[1]*X+g[2]*X^2

removed:
-factor(numer rhs %)/factor(denom rhs %)

removed:
-factor(numer rhs %)/factor(denom rhs %)

removed:
-factor(numer rhs %)/factor(denom rhs %)

added:
  Example 3.1
\begin{axiom}
eval(eq3_1,[r[2]=-r[1],r[3]=0])
eq3_4:= eval(eval(eq3_2,[eq3_3a,eq3_3b,eq3_3c]),[r[2]=-r[1],r[3]=0])
htrigs rhs %
\end{axiom}

Choose n items from a list. Returns list of size binomial(#a,n) of lists.

fricas

choose(a,n) ==
  j:=[i for i in 1..n]
  r:=[[a(j(i)) for i in 1..n]]
  k:=n
  while k>0 and j(k)+n-k<#a repeat
    j(k):=j(k)+1
    for i in k..n-1 repeat j(i+1):=j(i)+1
    r:=concat(r,[a(j(i)) for i in 1..n])
    k:=n; while j(k)+n-k>=#a and k>1 repeat k:=k-1
  if binomial(#a,n)~=#r then error "error in choose"
  return r

Type: Void

The Main Result

fricas

sum(x)==reduce(+,x,0)

Type: Void

fricas

product(x)==reduce(*,x,1)

Type: Void

fricas

-- specify n
f(i,j) == sum [ product x for x in choose([r[q]::Expression Integer for q in 1..n|q~=j],n-i-1)]

Type: Void

fricas

groupPolyCoeff(i) == (-1)^(i+n+1)*reduce(+,[exp(r[j])/reduce(*,[r[j]-r[m] for m in 1..n | j~=m])*f(i,j) for j in 1..n])

Type: Void

Polynomial of degree 2

fricas

n:=2

$\label{eq1}2$

(1)

Type: PositiveInteger?

fricas

eq2_1:= m[X]=(x-r[1])*(x-r[2])

$\label{eq2}{m_{X}}={{{x}^{2}}+{{\left(-{r_{2}}-{r_{1}}\right)}\ x}+{{r_{1}}\ {r_{2}}}}$

(2)

Type: Equation(Polynomial(Integer))

fricas

eq2_2:= exp(X)=g[0]*id+g[1]*X

$\label{eq3}{{e}^{X}}={{{g_{0}}\ id}+{{g_{1}}\ X}}$

(3)

Type: Equation(Expression(Integer))

fricas

eq2_3a:= g[0]=groupPolyCoeff(0)

fricas

Compiling function choose with type (List(Expression(Integer)),
      Integer) -> List(List(Expression(Integer)))

fricas

Compiling function product with type List(Expression(Integer)) -> 
      Expression(Integer)

fricas

Compiling function sum with type List(Expression(Integer)) -> 
      Expression(Integer)

fricas

Compiling function f with type (NonNegativeInteger,PositiveInteger)
       -> Expression(Integer)

fricas

Compiling function groupPolyCoeff with type NonNegativeInteger -> 
      Expression(Integer)

$\label{eq4}{g_{0}}={{-{{r_{1}}\ {{e}^{r_{2}}}}+{{r_{2}}\ {{e}^{r_{1}}}}}\over{{r_{2}}-{r_{1}}}}$

(4)

Type: Equation(Expression(Integer))

fricas

eq2_3b:= g[1]=groupPolyCoeff(1)

fricas

Compiling function f with type (PositiveInteger,PositiveInteger) -> 
      Expression(Integer)

fricas

Compiling function groupPolyCoeff with type PositiveInteger -> 
      Expression(Integer)

$\label{eq5}{g_{1}}={{{{e}^{r_{2}}}-{{e}^{r_{1}}}}\over{{r_{2}}-{r_{1}}}}$

(5)

Type: Equation(Expression(Integer))

Example 2.1

fricas

eval(eq2_1,[r[2]=-r[1]])

$\label{eq6}{m_{X}}={{{x}^{2}}-{{r_{1}}^{2}}}$

(6)

Type: Equation(Polynomial(Integer))

fricas

eq2_4:= eval(eval(eq2_2,[eq2_3a,eq2_3b]),r[2]=-r[1])

$\label{eq7}{{e}^{X}}={{{{\left({{r_{1}}\ id}+ X \right)}\ {{e}^{r_{1}}}}+{{\left({{r_{1}}\ id}- X \right)}\ {{e}^{-{r_{1}}}}}}\over{2 \ {r_{1}}}}$

(7)

Type: Equation(Expression(Integer))

fricas

htrigs rhs %

$\label{eq8}{{X \ {\sinh \left({r_{1}}\right)}}+{{r_{1}}\ id \ {\cosh \left({r_{1}}\right)}}}\over{r_{1}}$

(8)

Type: Expression(Integer)

Polynomial of degree 3
fricas
```
n:=3
```
$\label{eq9}3$ (9)

Type: PositiveInteger?

fricas
```
eq3_1:= m[X]=(x-r[1])*(x-r[2])*(x-r[3])
```
$\label{eq10}\begin{array}{@{}l} \displaystyle {m_{X}}={ \begin{array}{@{}l} \displaystyle {{x}^{3}}+{{\left(-{r_{3}}-{r_{2}}-{r_{1}}\right)}\ {{x}^{2}}}+{{\left({{\left({r_{2}}+{r_{1}}\right)}\ {r_{3}}}+{{r_{1}}\ {r_{2}}}\right)}\ x}- \ \ \displaystyle {{r_{1}}\ {r_{2}}\ {r_{3}}}$ (10)

Type: Equation(Polynomial(Integer))

fricas
```
eq3_2:= exp(X)=g[0]*id+g[1]*X+g[2]*X^2
```
$\label{eq11}{{e}^{X}}={{{g_{0}}\ id}+{{g_{2}}\ {{X}^{2}}}+{{g_{1}}\ X}}$ (11)

Type: Equation(Expression(Integer))

fricas
```
eq3_3a:= g[0]=groupPolyCoeff(0)
```
$\label{eq12}\begin{array}{@{}l} \displaystyle {g_{0}}={{\left( \begin{array}{@{}l} \displaystyle {{\left({{r_{1}}\ {{r_{2}}^{2}}}-{{{r_{1}}^{2}}\ {r_{2}}}\right)}\ {{e}^{r_{3}}}}+ \ \ \displaystyle {{\left(-{{r_{1}}\ {{r_{3}}^{2}}}+{{{r_{1}}^{2}}\ {r_{3}}}\right)}\ {{e}^{r_{2}}}}+ \ \ \displaystyle {{\left({{r_{2}}\ {{r_{3}}^{2}}}-{{{r_{2}}^{2}}\ {r_{3}}}\right)}\ {{e}^{r_{1}}}}$ (12)

Type: Equation(Expression(Integer))

fricas
```
eq3_3b:= g[1]=groupPolyCoeff(1)
```
$\label{eq13}{g_{1}}={{{{\left(-{{r_{2}}^{2}}+{{r_{1}}^{2}}\right)}\ {{e}^{r_{3}}}}+{{\left({{r_{3}}^{2}}-{{r_{1}}^{2}}\right)}\ {{e}^{r_{2}}}}+{{\left(-{{r_{3}}^{2}}+{{r_{2}}^{2}}\right)}\ {{e}^{r_{1}}}}}\over{{{\left({r_{2}}-{r_{1}}\right)}\ {{r_{3}}^{2}}}+{{\left(-{{r_{2}}^{2}}+{{r_{1}}^{2}}\right)}\ {r_{3}}}+{{r_{1}}\ {{r_{2}}^{2}}}-{{{r_{1}}^{2}}\ {r_{2}}}}}$ (13)

Type: Equation(Expression(Integer))

fricas
```
eq3_3c:= g[2]=groupPolyCoeff(2)
```
$\label{eq14}{g_{2}}={{{{\left({r_{2}}-{r_{1}}\right)}\ {{e}^{r_{3}}}}+{{\left(-{r_{3}}+{r_{1}}\right)}\ {{e}^{r_{2}}}}+{{\left({r_{3}}-{r_{2}}\right)}\ {{e}^{r_{1}}}}}\over{{{\left({r_{2}}-{r_{1}}\right)}\ {{r_{3}}^{2}}}+{{\left(-{{r_{2}}^{2}}+{{r_{1}}^{2}}\right)}\ {r_{3}}}+{{r_{1}}\ {{r_{2}}^{2}}}-{{{r_{1}}^{2}}\ {r_{2}}}}}$ (14)

Type: Equation(Expression(Integer))
Example 3.1
fricas
```
eval(eq3_1,[r[2]=-r[1],r[3]=0])
```
$\label{eq15}{m_{X}}={{{x}^{3}}-{{{r_{1}}^{2}}\ x}}$ (15)

Type: Equation(Polynomial(Integer))

fricas
```
eq3_4:= eval(eval(eq3_2,[eq3_3a,eq3_3b,eq3_3c]),[r[2]=-r[1],r[3]=0])
```
$\label{eq16}{{e}^{X}}={{{{\left({{X}^{2}}+{{r_{1}}\ X}\right)}\ {{e}^{r_{1}}}}+{{\left({{X}^{2}}-{{r_{1}}\ X}\right)}\ {{e}^{-{r_{1}}}}}+{2 \ {{r_{1}}^{2}}\ id}-{2 \ {{X}^{2}}}}\over{2 \ {{r_{1}}^{2}}}}$ (16)

Type: Equation(Expression(Integer))

fricas
```
htrigs rhs %
```
$\label{eq17}{{{r_{1}}\ X \ {\sinh \left({r_{1}}\right)}}+{{{X}^{2}}\ {\cosh \left({r_{1}}\right)}}+{{{r_{1}}^{2}}\ id}-{{X}^{2}}}\over{{r_{1}}^{2}}$ (17)

Type: Expression(Integer)
Polynomial of degree 4
fricas
```
n:=4
```
$\label{eq18}4$ (18)

Type: PositiveInteger?

fricas
```
eq4_2:= g[0]=groupPolyCoeff(0)
```
$\label{eq19}\begin{array}{@{}l} \displaystyle {g_{0}}={{\left( \begin{array}{@{}l} \displaystyle {{\left({ \begin{array}{@{}l} \displaystyle {{\left(-{{r_{1}}\ {{r_{2}}^{2}}}+{{{r_{1}}^{2}}\ {r_{2}}}\right)}\ {{r_{3}}^{3}}}+ \ \ \displaystyle {{\left({{r_{1}}\ {{r_{2}}^{3}}}-{{{r_{1}}^{3}}\ {r_{2}}}\right)}\ {{r_{3}}^{2}}}+ \ \ \displaystyle {{\left(-{{{r_{1}}^{2}}\ {{r_{2}}^{3}}}+{{{r_{1}}^{3}}\ {{r_{2}}^{2}}}\right)}\ {r_{3}}}$ (19)

Type: Equation(Expression(Integer))

fricas
```
eq4_3:= g[1]=groupPolyCoeff(1)
```
$\label{eq20}\begin{array}{@{}l} \displaystyle {g_{1}}={{\left( \begin{array}{@{}l} \displaystyle {{\left({ \begin{array}{@{}l} \displaystyle {{\left({{r_{2}}^{2}}-{{r_{1}}^{2}}\right)}\ {{r_{3}}^{3}}}+ \ \ \displaystyle {{\left(-{{r_{2}}^{3}}+{{r_{1}}^{3}}\right)}\ {{r_{3}}^{2}}}+{{{r_{1}}^{2}}\ {{r_{2}}^{3}}}- \ \ \displaystyle {{{r_{1}}^{3}}\ {{r_{2}}^{2}}}$ (20)

Type: Equation(Expression(Integer))

fricas
```
eq4_4:= g[2]=groupPolyCoeff(2)
```
$\label{eq21}\begin{array}{@{}l} \displaystyle {g_{2}}={{\left( \begin{array}{@{}l} \displaystyle {{\left({ \begin{array}{@{}l} \displaystyle {{\left(-{r_{2}}+{r_{1}}\right)}\ {{r_{3}}^{3}}}+{{\left({{r_{2}}^{3}}-{{r_{1}}^{3}}\right)}\ {r_{3}}}- \ \ \displaystyle {{r_{1}}\ {{r_{2}}^{3}}}+{{{r_{1}}^{3}}\ {r_{2}}}$ (21)

Type: Equation(Expression(Integer))

fricas
```
eq4_5:= g[3]=groupPolyCoeff(3)
```
$\label{eq22}\begin{array}{@{}l} \displaystyle {g_{3}}={{\left( \begin{array}{@{}l} \displaystyle {{\left({ \begin{array}{@{}l} \displaystyle {{\left({r_{2}}-{r_{1}}\right)}\ {{r_{3}}^{2}}}+{{\left(-{{r_{2}}^{2}}+{{r_{1}}^{2}}\right)}\ {r_{3}}}+ \ \ \displaystyle {{r_{1}}\ {{r_{2}}^{2}}}-{{{r_{1}}^{2}}\ {r_{2}}}$ (22)

Type: Equation(Expression(Integer))

Old

fricas

eq42 := _
  -r2*r3*r4*exp(r1)/((r1-r2)*(r1-r3)*(r1-r4)) + _
  -r1*r3*r4*exp(r2)/((r2-r1)*(r2-r3)*(r2-r4)) + _
  -r1*r2*r4*exp(r3)/((r3-r1)*(r3-r2)*(r3-r4)) + _
  -r1*r2*r3*exp(r4)/((r4-r1)*(r4-r2)*(r4-r3))

$\label{eq23}{\left( \begin{array}{@{}l} \displaystyle {{\left({ \begin{array}{@{}l} \displaystyle {{\left(-{r 1 \ {{r 2}^{2}}}+{{{r 1}^{2}}\ r 2}\right)}\ {{r 3}^{3}}}+ \ \ \displaystyle {{\left({r 1 \ {{r 2}^{3}}}-{{{r 1}^{3}}\ r 2}\right)}\ {{r 3}^{2}}}+ \ \ \displaystyle {{\left(-{{{r 1}^{2}}\ {{r 2}^{3}}}+{{{r 1}^{3}}\ {{r 2}^{2}}}\right)}\ r 3}$

(23)

Type: Expression(Integer)

three

fricas

htrigs eval(eq42, [r3=-r1,r4=-r2])

$\label{eq24}{-{{{r 1}^{2}}\ {\cosh \left({r 2}\right)}}+{{{r 2}^{2}}\ {\cosh \left({r 1}\right)}}}\over{{{r 2}^{2}}-{{r 1}^{2}}}$

(24)

Type: Expression(Integer)

fricas

htrigs limit(%,r2=r1)

$\label{eq25}{-{r 1 \ {\sinh \left({r 1}\right)}}+{2 \ {\cosh \left({r 1}\right)}}}\over 2$

(25)

Type: Expression(Integer)

fricas

eq43 :=
  (r2*r3 + r3*r4 + r4*r2)*exp(r1)/((r1-r2)*(r1-r3)*(r1-r4)) + _
  (r1*r3 + r3*r4 + r4*r1)*exp(r2)/((r2-r1)*(r2-r3)*(r2-r4)) + _
  (r1*r2 + r2*r4 + r4*r1)*exp(r3)/((r3-r1)*(r3-r2)*(r3-r4)) + _
  (r1*r2 + r2*r3 + r3*r1)*exp(r4)/((r4-r1)*(r4-r2)*(r4-r3))

$\label{eq26}{\left( \begin{array}{@{}l} \displaystyle {{\left({ \begin{array}{@{}l} \displaystyle {{\left({{r 2}^{2}}-{{r 1}^{2}}\right)}\ {{r 3}^{3}}}+{{\left(-{{r 2}^{3}}+{{r 1}^{3}}\right)}\ {{r 3}^{2}}}+ \ \ \displaystyle {{{r 1}^{2}}\ {{r 2}^{3}}}-{{{r 1}^{3}}\ {{r 2}^{2}}}$

(26)

Type: Expression(Integer)

fricas

htrigs eval(eq43, [r3=-r1,r4=-r2])

$\label{eq27}{-{{{r 1}^{3}}\ {\sinh \left({r 2}\right)}}+{{{r 2}^{3}}\ {\sinh \left({r 1}\right)}}}\over{{r 1 \ {{r 2}^{3}}}-{{{r 1}^{3}}\ r 2}}$

(27)

Type: Expression(Integer)

fricas

htrigs limit(%,r2=r1)

$\label{eq28}{{3 \ {\sinh \left({r 1}\right)}}-{r 1 \ {\cosh \left({r 1}\right)}}}\over{2 \ r 1}$

(28)

Type: Expression(Integer)