![]() ![]() | ![]() |
version 1.4, 2014/08/06 07:39:47 | version 1.6, 2014/12/12 08:29:53 | ||
---|---|---|---|
|
|
||
/* $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")$ | load("tk_fd.rr")$ | ||
import("tk_r.rr")$ | import("tk_r.rr")$ | ||
import("oh_number.rr")$ | import("oh_number.rr")$ | ||
|
|
||
} | } | ||
def r_ahvec(A,B,C,Y) { | def r_ahvec(A,B,C,Y) { | ||
Y = r_d2rat(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]; | Gamma=R[1]; | ||
Der=R[0]; | Der=R[0]; | ||
Z=R[2]*Gamma; | Z=R[2]*Gamma; | ||
|
|
||
} | } | ||
} | } | ||
} | } | ||
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$ | end$ |