===================================================================
RCS file: /home/cvs/OpenXM/src/asir-port/cgi/r-fd.rr,v
retrieving revision 1.4
retrieving revision 1.7
diff -u -p -r1.4 -r1.7
--- OpenXM/src/asir-port/cgi/r-fd.rr	2014/08/06 07:39:47	1.4
+++ OpenXM/src/asir-port/cgi/r-fd.rr	2015/02/25 04:47:50	1.7
@@ -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.6 2014/12/12 08:29:53 takayama Exp $ */
 load("tk_fd.rr")$
 import("tk_r.rr")$
 import("oh_number.rr")$
@@ -18,7 +18,15 @@ 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);
+  Ans=a_ahvec(A,B,C,Y);
+  /*  Fans=map(rtostr,map(tk_fd.tk_number_rattofloat,Ans)); */
+  Fans=map(deval,Ans);
+  Fans = tk_r.asir2r_c(Fans);
+  return Fans;
+}
+
+def a_ahvec(A,B,C,Y) { 
+  R=tk_fd.ahvec_abc(A,B,C,Y|all=1);
   Gamma=R[1];
   Der=R[0];
   Z=R[2]*Gamma;
@@ -26,10 +34,7 @@ def r_ahvec(A,B,C,Y) { 
   for (I=0; I<length(Der); I++) Der2[I] = Der[I]*Gamma;
   Der2 = vtol(Der2);
   Ans=cons(Z,Der2);
-  /*  Fans=map(rtostr,map(tk_fd.tk_number_rattofloat,Ans)); */
-  Fans=map(deval,Ans);
-  Fans = tk_r.asir2r_c(Fans);
-  return Fans;
+  return(Ans);
 }
 
 /* temporary */
@@ -54,4 +59,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$