=================================================================== RCS file: /home/cvs/OpenXM/src/asir-port/cgi/r-fd.rr,v retrieving revision 1.4 retrieving revision 1.6 diff -u -p -r1.4 -r1.6 --- OpenXM/src/asir-port/cgi/r-fd.rr 2014/08/06 07:39:47 1.4 +++ OpenXM/src/asir-port/cgi/r-fd.rr 2014/12/12 08:29:53 1.6 @@ -1,4 +1,4 @@ -/* $OpenXM: OpenXM/src/asir-port/cgi/r-fd.rr,v 1.3 2014/07/31 07:00:41 takayama Exp $ */ +/* $OpenXM: OpenXM/src/asir-port/cgi/r-fd.rr,v 1.5 2014/08/27 02:14:40 takayama Exp $ */ load("tk_fd.rr")$ import("tk_r.rr")$ import("oh_number.rr")$ @@ -18,7 +18,7 @@ def r_d2rat(Y) { } def r_ahvec(A,B,C,Y) { Y = r_d2rat(Y); - R=tk_fd.ygahvec(A,B,C,Y|all=1); + R=tk_fd.ahvec_abc(A,B,C,Y|all=1); Gamma=R[1]; Der=R[0]; Z=R[2]*Gamma; @@ -54,4 +54,17 @@ def checkrats2() { } } } + +def a_expect(A,B,C,Y) { + E=tk_fd.expectation_abc(A,B,C,Y); + return(E); +} +def r_expect(A,B,C,Y) { + Y = r_d2rat(Y); + E=a_expect(A,B,C,Y); + Fans=map(deval,E); + Fans = tk_r.asir2r_c(Fans); + return Fans; +} + end$