login  home  contents  what's new  discussion  bug reports     help  links  subscribe  changes  refresh  edit

Edit detail for SandBoxExpOfEnd revision 6 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 20:59:53 GMT+0
Note:

changed:
-1.3 The Main Result
1 The Main Result

changed:
-Polynomial of degree 2
2 Polynomial of degree 2
\begin{axiom}
eq2_2:= exp(X)=g[0]*id+g[1]*X
\end{axiom}

changed:
-Polynomial of degree 3
\begin{axiom}
eq2_4:= eval(eval(eq2_2,[eq2_3a,eq2_3b]),r[2]=-r[1])
htrigs rhs %
\end{axiom}

3 Polynomial of degree 3

changed:
-Polynomial of degree 4

4 Polynomial of degree 4

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

  1. 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
  2. Polynomial of degree 2
    fricas
    eq2_2:= exp(X)=g[0]*id+g[1]*X
    
\label{eq1}{{e}^{X}}={{{g_{0}}\  id}+{{g_{1}}\  X}}(1)
    Type: Equation(Expression(Integer))

fricas
n:=2

\label{eq2}2(2)
Type: PositiveInteger?
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{eq3}{g_{0}}={{-{{r_{1}}\ {{e}^{r_{2}}}}+{{r_{2}}\ {{e}^{r_{1}}}}}\over{{r_{2}}-{r_{1}}}}(3)
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{eq4}{g_{1}}={{{{e}^{r_{2}}}-{{e}^{r_{1}}}}\over{{r_{2}}-{r_{1}}}}(4)
Type: Equation(Expression(Integer))

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

\label{eq5}{{e}^{X}}={{{{\left({{r_{1}}\  id}+ X \right)}\ {{e}^{r_{1}}}}+{{\left({{r_{1}}\  id}- X \right)}\ {{e}^{-{r_{1}}}}}}\over{2 \ {r_{1}}}}(5)
Type: Equation(Expression(Integer))
fricas
htrigs rhs %

\label{eq6}{{X \ {\sinh \left({r_{1}}\right)}}+{{r_{1}}\  id \ {\cosh \left({r_{1}}\right)}}}\over{r_{1}}(6)
Type: Expression(Integer)

  1. Polynomial of degree 3
    fricas
    n:=3
    
\label{eq7}3(7)
    Type: PositiveInteger?
    fricas
    eq3_3a:= g[0]=groupPolyCoeff(0)
    
\label{eq8}\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}}}}
(8)
    Type: Equation(Expression(Integer))
    fricas
    factor(numer rhs %)/factor(denom rhs %)
    
\label{eq9}{\left(
\begin{array}{@{}l}
\displaystyle
{{\left({{r_{1}}\ {{r_{2}}^{2}}}-{{{r_{1}}^{2}}\ {r_{2}}}\right)}\ {{e}^{r_{3}}}}+{{\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}}}}
(9)
    Type: Fraction(Factored(SparseMultivariatePolynomial?(Integer,Kernel(Expression(Integer)))))
    fricas
    eq3_3b:= g[1]=groupPolyCoeff(1)
    
\label{eq10}{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}}}}}(10)
    Type: Equation(Expression(Integer))
    fricas
    factor(numer rhs %)/factor(denom rhs %)
    
\label{eq11}-{{{{\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)}\ {\left({r_{3}}-{r_{2}}\right)}\ {\left({r_{3}}-{r_{1}}\right)}}}(11)
    Type: Fraction(Factored(SparseMultivariatePolynomial?(Integer,Kernel(Expression(Integer)))))
    fricas
    eq3_3c:= g[2]=groupPolyCoeff(2)
    
\label{eq12}{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}}}}}(12)
    Type: Equation(Expression(Integer))
    fricas
    factor(numer rhs %)/factor(denom rhs %)
    
\label{eq13}{{{\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)}\ {\left({r_{3}}-{r_{2}}\right)}\ {\left({r_{3}}-{r_{1}}\right)}}(13)
    Type: Fraction(Factored(SparseMultivariatePolynomial?(Integer,Kernel(Expression(Integer)))))
  2. Polynomial of degree 4
    fricas
    n:=4
    
\label{eq14}4(14)
    Type: PositiveInteger?
    fricas
    eq4_2:= g[0]=groupPolyCoeff(0)
    
\label{eq15}\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}}}
(15)
    Type: Equation(Expression(Integer))
    fricas
    eq4_3:= g[1]=groupPolyCoeff(1)
    
\label{eq16}\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}}}
(16)
    Type: Equation(Expression(Integer))
    fricas
    eq4_4:= g[2]=groupPolyCoeff(2)
    
\label{eq17}\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}}}
(17)
    Type: Equation(Expression(Integer))
    fricas
    eq4_5:= g[3]=groupPolyCoeff(3)
    
\label{eq18}\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}}}
(18)
    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{eq19}{\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: Expression(Integer)

three

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

\label{eq20}{-{{{r 1}^{2}}\ {\cosh \left({r 2}\right)}}+{{{r 2}^{2}}\ {\cosh \left({r 1}\right)}}}\over{{{r 2}^{2}}-{{r 1}^{2}}}(20)
Type: Expression(Integer)
fricas
htrigs limit(%,r2=r1)

\label{eq21}{-{r 1 \ {\sinh \left({r 1}\right)}}+{2 \ {\cosh \left({r 1}\right)}}}\over 2(21)
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{eq22}{\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}}}
(22)
Type: Expression(Integer)

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

\label{eq23}{-{{{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}}(23)
Type: Expression(Integer)
fricas
htrigs limit(%,r2=r1)

\label{eq24}{{3 \ {\sinh \left({r 1}\right)}}-{r 1 \ {\cosh \left({r 1}\right)}}}\over{2 \  r 1}(24)
Type: Expression(Integer)