[BACK]Return to appell.sm1 CVS log [TXT][DIR] Up to [local] / OpenXM / src / kan96xx / Doc

Diff for /OpenXM/src/kan96xx/Doc/appell.sm1 between version 1.1 and 1.4

version 1.1, 1999/10/08 02:12:02 version 1.4, 2003/08/18 06:36:49
Line 1 
Line 1 
 %% appell.sm1, 1998,  11/8  %% appell.sm1, 1998,  11/8
   % $OpenXM: OpenXM/src/kan96xx/Doc/appell.sm1,v 1.3 2003/07/29 08:37:16 takayama Exp $
 /appell.version (2.981108) def  /appell.version (2.981108) def
 appell.version [(Version)] system_variable gt  appell.version [(Version)] system_variable gt
 { (This package requires the latest version of kan/sm1) message  { (This package requires the latest version of kan/sm1) message
Line 161  $appell.sm1 generates Appell hypergeometric differenti
Line 162  $appell.sm1 generates Appell hypergeometric differenti
  ]   ]
 ] putUsages  ] putUsages
   
   /appell4r {
     /arg1 set
     [/in-appell4r /typev /setarg /b /n /vv /i /a /c /bb /ans
     ] pushVariables
     [(CurrentRingp) (KanGBmessage)] pushEnv  %% push current global environment.
     [
       /aa arg1 def
       aa isArray { } { (array appell) message (appell4) usage error } ifelse
       /setarg 0 def
       aa { tag } map /typev set
       typev [ ArrayP ] eq
       {  /b aa 0 get def
          /setarg 1 def
       } { } ifelse
       typev [ ] eq
       {
          /b appell.b def
          /setarg 1 def
       } { } ifelse
       setarg { } { (Argument mismatch) message (appell4) usage error } ifelse
   
       [(KanGBmessage) appell.verbose] system_variable
   
       /n b length 2 sub def   %% Lauricella F_C^n
   
       %% vv = [(x1) (x2)]
       [
         1 1 n {
           /i set
           (x) i gensym
         } for
       ] /vv set
   
       [vv from_records ring_of_differential_operators 0] define_ring
   
       %% b = [a  b  c_1 ... c_n ]
       /a b 0 get def
       /c b 1 get def
       /bb b rest rest def
   
       [ 1 1 n {
           /i set
              [@@@.Dsymbol (x)] cat i gensym  .
              i i appellr.euler . bb i 1 sub get (1)..  sub add
              mul
   
              1 n appellr.euler . (0).. a add add
              1 n appellr.euler . (0).. c add add
              mul
   
              sub
              (numerator) dc cancelCoeff dehomogenize
              toString
         } for
       ] /ans set
       /arg1 [ans vv] def
     ] pop
     popEnv
     popVariables
     arg1
   } def
   %% [ [(1).. (2).. div -4 -2 5 6] ] appell4r
   
   
 /appell.euler {  /appell.euler {
   /arg2 set    /arg2 set
   /arg1 set    /arg1 set
Line 177  $appell.sm1 generates Appell hypergeometric differenti
Line 240  $appell.sm1 generates Appell hypergeometric differenti
   arg1    arg1
 } def  } def
   
   /appellr.euler {
     /arg2 set
     /arg1 set
     [/n /i /n0] pushVariables
     [
       /n0 arg1 def
       /n arg2 def
       [ n0 1 n { /i set (x) i gensym  ( ) [@@@.Dsymbol (x)] cat i gensym
                  i n eq not { ( + ) } {  } ifelse } for ]  cat
       /arg1 set
     ] pop
     popVariables
     arg1
   } def
   
   /appell2 {
     /arg1 set
     [/in-appell2 /typev /setarg /b /n /vv /i /a /c /bb /ans
     ] pushVariables
     [(CurrentRingp) (KanGBmessage)] pushEnv  %% push current global environment.
     [
       /aa arg1 def
       aa isArray { } { (array appell) message (appell2) usage error } ifelse
       /setarg 0 def
       aa { tag } map /typev set
       typev [ ArrayP ] eq
       {  /b aa 0 get def
          /setarg 1 def
       } { } ifelse
       typev [ ] eq
       {
          /b [1 [2 3] [4 5]] def
          /setarg 1 def
       } { } ifelse
       setarg { } { (Argument mismatch) message (appell2) usage error } ifelse
   
       [(KanGBmessage) appell.verbose] system_variable
   
       /n b 1 get length def   %% Lauricella F_A^n
   
       %% vv = [(x1) (x2)]
       [
         1 1 n {
           /i set
           (x) i gensym
         } for
       ] /vv set
   
       %% b = [a  [b_1 ... b_n]  [c_1 ... c_n] ]
       /a b 0 get def
       /c b 2 get def
       /bb b 1 get def
   
       [ 1 1 n {
           /i set
           [  [@@@.Dsymbol (x)] cat i gensym
              $ ($  i i appell.euler c i 1 sub get 1 sub $) - ( $
              1 n appell.euler a $) ($ i i appell.euler bb i 1 sub get $ ) $] cat
         } for
       ] /ans set
       /arg1 [ans vv] def
     ] pop
     popEnv
     popVariables
     arg1
   } def
   (appell2 ) messagen-quiet
   [(appell2)
    [(param appell2 c)
     (array param; array c;)
     (appell2 returns the Lauricella function F_A(a,b_1, ..., b_n, c_1, ..., c_n; x_1,...,x_n))
     (for the parameter << param >> = [a, [b_1, ..., b_n],[c_1, ..., c_n]].)
     (In case of n=2, the function is called the Appell function F_2.)
     (c = [ generators,  variables ])
     (Example: [ [1 [-4 -2] [5 6]] ] appell2 rank ::)
    ]
   ] putUsages
   
   /appell2r {
     /arg1 set
     [/in-appell2r /typev /setarg /b /n /vv /i /a /c /bb /ans
     ] pushVariables
     [(CurrentRingp) (KanGBmessage)] pushEnv  %% push current global environment.
     [
       /aa arg1 def
       aa isArray { } { (array appell) message (appell2) usage error } ifelse
       /setarg 0 def
       aa { tag } map /typev set
       typev [ ArrayP ] eq
       {  /b aa 0 get def
          /setarg 1 def
       } { } ifelse
       typev [ ] eq
       {
          /b [1 [2 3] [4 5]] def
          /setarg 1 def
       } { } ifelse
       setarg { } { (Argument mismatch) message (appell2) usage error } ifelse
   
       [(KanGBmessage) appell.verbose] system_variable
   
       /n b 1 get length def   %% Lauricella F_A^n
   
       %% vv = [(x1) (x2)]
       [
         1 1 n {
           /i set
           (x) i gensym
         } for
       ] /vv set
   
       [vv from_records ring_of_differential_operators 0] define_ring
   
       %% b = [a  [b_1 ... b_n]  [c_1 ... c_n] ]
       /a b 0 get def
       /c b 2 get def
       /bb b 1 get def
   
       [ 1 1 n {
           /i set
              [@@@.Dsymbol (x)] cat i gensym .
              i i appellr.euler . c i 1 sub get (1).. sub add
              mul
   
              1 n appellr.euler . (0).. a add add
              i i appellr.euler . (0).. bb i 1 sub get add add
              mul
   
              sub
              (numerator) dc cancelCoeff dehomogenize
              toString
         } for
       ] /ans set
       /arg1 [ans vv] def
     ] pop
     popEnv
     popVariables
     arg1
   } def
   %%[[(1).. (2).. div [(1).. (2).. div (1).. (2).. div] [1 1]] ] appell2r rank ::
   
 ( ) message-quiet ;  
   
   ( ) message-quiet ;
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>