start Risa/Asir with the asir-contrib library. load("sf.rr"); bynd(); myprune(); A=mypf()$ /* Pfaffian system is saved in the file tmp-pf-sf32.ab It is a list of 6 matrices A[0], ..., A[5], which are coefficients of the Pfaffian system. V=dsf(3,2)[1] a list of variables. Let F be a solution vector. Then, dF/dV[i] = A[i]F holds. */
The function dsf(3,2) returns the system of differential equations for c(\Theta).
Start R. source("etrcalc.r"); SO.init(Ts.Xm) # Initial value for HGD for commets data. The output is embedded in sf2.rr SO.init(TsAst,XmAst) # Initial value for HGD for asteroids data. SO.int2(Ts) # value of c(\Theta) for Ts.
load("sdo.rr"); bynd()$
Start Risa/Asir Pfs=bload("pf.ab")$ /* Pfs[0] is a list of 9 matrices, which are the coefficients of the Pfaffian system. A hashed data structure is used to save polynomials. */The data is generated by mydso9.rr
Start Risa/Asir load("hgd-test.rr"); load("hgd-test16.rr");