#define INCLUDES APLCB+APLMAIN+STDIO
#include "includes.h"
main()
{
     Dyadic; Endoper; Idyadic; Imod; Iresidu;
     Litvect; Quadout; Residue; Residuep; Residuex;
     double left[2],rite[2],ret[2];
     #include "aplinit.h"

     endoper(quadout(dyadic(residue,
          litvect("0  1     5 5 -5"),
          litvect("-3 2.385 2 -2 2"))));
     endoper(quadout(idyadic(imod,
       litvect(" 5  5 -5  0  2  2  2 -2 -2"),
       litvect(" 2 -2  2 -3  8  3 -3  3 -3"))));
     endoper(quadout(idyadic(iresidu,
       litvect(" 2 -2  2 -3  8  3 -3  3 -3"),
       litvect(" 5  5 -5  0  2  2  2 -2 -2"))));
     *left = 2.1;
     *rite = 5.6;
     residuep(left,rite,ret);
     printf("residuep returned %f\n",*ret);
     *left = 0;
     *(left + 1) = 12.0;
     *rite = 4.0;
     *(rite + 1) = 0;
     residuex(left,rite,ret);
     printf("residuex returned %fj%f\n",*ret,*(ret+1));
}