/* Copyright (C) 1994 by Thomas Glen Smith. All Rights Reserved. */ /* residuex APL2 V1.0.0 ************************************************ * Residue of complex numbers. * ***********************************************************************/ #define INCLUDES 0 #include "includes.h" void residuex(left,rite,ret) double *left,*rite,*ret; { Dividex; Floorx; Minusx; Timesx; double wrka[2],wrkb[2]; if (*left == 0e0 && *(left+1) == 0e0) { *ret = *rite; *(ret+1) = *(rite+1); } else { dividex(rite,left,wrka); floorx(wrka,wrkb); timesx(left,wrkb,wrka); minusx(rite,wrka,ret); } }