#include
#include
#define N 7 /* u array size */
#define M 2 /* response function dimension */
#define N2 N-M+1
int main() {
int i,j,n,m;
array double u[N] = {6,19,32,45,58,71,42},
v[M] = {6,7};
array double q[N2], r[N];
printf("u=%f\n",u);
printf("v=%f\n",v);
deconv(u,v,q); /* real v, not specify remaind result */
deconv(u,v,q,r); /* real v, specified remaind result */
printf("q=%f\n",q);
printf("r=%f\n",r);
}
u=6.000000 19.000000 32.000000 45.000000 58.000000 71.000000 42.000000
v=6.000000 7.000000
q=1.000000 2.000000 3.000000 4.000000 5.000000 6.000000
r=0.000000 5.000000 32.000000 45.000000 58.000000 71.000000 42.000000