int main() {
array double x[6]={1,2,3,4,5,6}, y[6];
array double complex zx[2][3]={complex(1,1),2,3,complex(2,2),5,6}, zy[2][3];
array double x1[2][3][4]={{1,2,3,4,
5,6,7,8,
5,6,7,8},
{10,11,12,13,
10,11,12,13,
14,15,16,17}};
array double y1[2][3][4];
cumprod(y,x);
printf("x=%f",x);
printf("y=%f",y);
printf("\n");
cumprod(zy,zx);
printf("zx=%5.2f",zx);
printf("zy=%5.2f",zy);
printf("\n");
cumprod(y1,x1);
printf("x1=%f",x1);
printf("y1=%f",y1);
}
x=1.000000 2.000000 3.000000 4.000000 5.000000 6.000000
y=1.000000 2.000000 6.000000 24.000000 120.000000 720.000000
zx=complex( 1.00, 1.00) complex( 2.00, 0.00) complex( 3.00, 0.00)
complex( 2.00, 2.00) complex( 5.00, 0.00) complex( 6.00, 0.00)
zy=complex( 1.00, 1.00) complex( 2.00, 2.00) complex( 6.00, 6.00)
complex( 2.00, 2.00) complex(10.00,10.00) complex(60.00,60.00)
x1=1.000000 2.000000 3.000000 4.000000
5.000000 6.000000 7.000000 8.000000
5.000000 6.000000 7.000000 8.000000
10.000000 11.000000 12.000000 13.000000
10.000000 11.000000 12.000000 13.000000
14.000000 15.000000 16.000000 17.000000
y1=1.000000 2.000000 6.000000 24.000000
5.000000 30.000000 210.000000 1680.000000
5.000000 30.000000 210.000000 1680.000000
10.000000 110.000000 1320.000000 17160.000000
10.000000 110.000000 1320.000000 17160.000000
14.000000 210.000000 3360.000000 57120.000000