#include <stdio.h>
#include <math.h>
void main()
{void u0(float *p0);
void u1(float *p1,float x1[]);
void u2(float *p2,float x2[]);
void u3(float *p3,float x3[]);
void juzh(float *p4,float b[][5],float f1[],int n);
float x[5]={0.1,0.4,0.5,0.7,0.9};
float f[5]={0.61,0.92,0.99,1.52,2.03};
float a[4][5], a1[4][5];
float (*p)[5];
int i,j;
p=a;
u0(*p);
u1(*(p+1),x);
u2(*(p+2),x);
u3(*(p+3),x);
for(i=0;i<=3;i++)
{for(j=0;i<=4;j++)
    a1[i][j]=a[i][j];
}
p=a;
for(j=0;j<=3;j++,p++)
{juzh(*p,a1,f,j);
}
for(i=0;i<=3;i++)
{for(j=0;i<=4;j++)
if(j%4==0)
   printf("\n");
printf("%.6f",a[i][j]);
}
}
void u0(float *p0)
{int i;
for(i=0;i<=4;i++,p0++)
   *p0=1;
}
void u1(float *p1,float x1[])
{int i;
for(i=0;i<=4;i++,p1++)
   *p1=x1[i];
}
void u2(float *p2,float x2[])
{int i;
float x;
for(i=0;i<=4;i++,p2++)
{x=x2[i];
 *p2=sin(x);
}
   
}
void u3(float *p3,float x3[])
{int i;
float x;
for(i=0;i<=4;i++,p3++)
{x=x3[i];
 *p3=exp(x);
}
}
void juzh(float *p4,float b[][5],float f1[],int n)
{int i,j;
float sum1=0,sum2=0;
for(i=0;i<=3;i++,p4++)
{for(j=0;j<=4;j++)
{sum1=sum1+b[n][j]*b[i][j];

}
*p4=sum1;

}
for(i=0;i<=4;i++)
{sum2=sum2+b[n][i]*f1[i];
}
p4++;
*p4=sum2;
}