<![CDATA[
#include <stdio.h>
#include <numeric.h>
#define N 3 /* data array size */
#define M 6 /* response array size */
int main() {
array double u[N] = {1,2,3};
array double v[M] = {1, 5, 9, 7, 2, 0};
array double r[M-1],p[M-1];
array double k[1];
residue(u,v,r,p,k);
printf("coef. of u=%f\n",u);
printf("coef. of v=%f\n",v);
printf("coef. of r=%f\n",r);
printf("coef. of p=%f\n",p);
printf("coef. of k=%f\n",k);
}
]]>
<![CDATA[
coef. of u=1.000000 2.000000 3.000000
coef. of v=1.000000 5.000000 9.000000 7.000000 2.000000 0.000000
coef. of r=1.500000 0.000000 0.000000 -2.000000 1.500000
coef. of p=-2.000000 -1.000000 -1.000000 -1.000000 0.000000
coef. of k=0.000000
]]>