=================================================================== RCS file: /home/cvs/OpenXM/src/kan96xx/Doc/appell.sm1,v retrieving revision 1.1.1.1 retrieving revision 1.3 diff -u -p -r1.1.1.1 -r1.3 --- OpenXM/src/kan96xx/Doc/appell.sm1 1999/10/08 02:12:02 1.1.1.1 +++ OpenXM/src/kan96xx/Doc/appell.sm1 2003/07/29 08:37:16 1.3 @@ -1,4 +1,5 @@ %% appell.sm1, 1998, 11/8 +% $OpenXM$ /appell.version (2.981108) def appell.version [(Version)] system_variable gt { (This package requires the latest version of kan/sm1) message @@ -176,6 +177,69 @@ $appell.sm1 generates Appell hypergeometric differenti 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 -( ) message-quiet ; \ No newline at end of file +( ) message-quiet ;