Матрицы и векторы

[otto07]

Нужно написать программу которая считывает значения матриц(A,B) и векторов(x,y) с файла и вычисляет выражение A*x*((B*A)*y). Значение матриц и векторов записаны в файлах(.txt) и находятся в том же каталоге что и проект с программой. Заранее спасибо!

#include <stdio.h>
#include <conio.h>
#include <math.h>
#define n 5
void main(void)
{FILE *fp;
float x[n],y[n],a[n][n],b[n][n],ba[n][n],vr1[n],vr2[n],d[n],s;
int i,j,k;
fp=fopen("VectX.txt","r");
for(j=1;j<n;j++)
fscanf(fp,"%f",&x[j]);
fclose(fp);
printf("VectX=%f%f%f%f\n",x[j]);
fp=fopen("VectY.txt","r");
for(j=1;j<n;j++)
fscanf(fp,"%f",&y[j]);
fclose(fp);
printf("VectY=%f%f%f%f\n",y[j]);
fp=fopen("MatrA.txt","r");
for(i=1;i<n;i++)
for(j=1;j<n;j++)
fscanf(fp,"%f",&a[i][j]);
fclose(fp);
printf("MatrA=%f%f%f%f\n",a[i][j]);
fp=fopen("MatrB.txt","r");
for(i=1;i<n;i++)
for(j=1;j<n;j++)
fscanf(fp,"%f",&b[i][j]);
fclose(fp);
printf("MatrB=%f%f%f%f\n",b[i][j]);
{
s=0;
for(i=1;i<n;i++)
for(j=1;j<n;j++)
for(k=1;k<n;k++)
s=s+b[i][k]*a[k][j];
ba[i][j]=s;
}
{
s=0;
for(i=1;i<n;i++)
for(j=1;j<n;j++)
s=s+a[i][j]*x[j];
vr1[j]=s;
}
{
s=0;
for(i=1;i<n;i++)
for(j=1;j<n;j++)
s=s+ba[i][j]*y[j];
vr2[j]=s;
}
{
s=0;
for(j=1;j<n;j++)
s=s+vr1[j]*vr2[j];
d[j]=s;
printf("Result=%f%f%f%f",d[j]);
}
}