==========================================================
Fortran90 Source Code:
==========================================================

implicit real*8 (a-h,o-z)
dimension a(5,7),b(7,5),c(5,7),d(5),result(7,7)

data ((a(i,j),j=1,7),i=1,5)/ &
0.7577,0.706, 0.8235,0.4387,0.4898,0.276, 0.4984,&
0.7431,0.0318,0.6948,0.3816,0.4456,0.6797,0.9597,&
0.3922,0.2769,0.3171,0.7655,0.6463,0.6551,0.3404,&
0.6555,0.0462,0.9502,0.7952,0.7094,0.1626,0.5853,&
0.1712,0.0971,0.0344,0.1869,0.7547,0.119,0.2238/,&
((b(i,j),j=1,5),i=1,7)/ &
0.4314,0.5797,0.4018,0.0497,0.3692,&
0.9106,0.5499,0.076, 0.9027,0.1112,&
0.1818,0.145, 0.2399,0.9448,0.7803,&
0.2638,0.853, 0.1233,0.4909,0.3897,&
0.1455,0.6221,0.1839,0.4893,0.2417,&
0.1361,0.351, 0.24,  0.3377,0.4039,&
0.8693,0.5132,0.4173,0.9001,0.0965/,&
((c(i,j),j=1,7),i=1,5)/ &
0.1299,0.1622,0.602,0.4505,0.8258,0.1067,0.8687,&
0.5688,0.7943,0.263,0.0838,0.5383,0.9619,0.0844,&
0.4694,0.3112,0.6541,0.229,0.9961,0.0046,0.3998,&
0.0119,0.5285,0.6892,0.9133,0.0782,0.7749,0.2599,&
0.3371,0.1656,0.7482,0.1524,0.4427,0.8173,0.8001/,&
(d(i),i=1,5)/0.1320,0.9421,0.9561,0.5752,0.0598/

print *,''
print *,'Matrix A='
do i=1,5
   write(*,'(1x,7f8.4)') (a(i,j),j=1,7)
enddo

print *,''
print *,'Matrix B='
do i=1,7
   write(*,'(1x,7f8.4)') (b(i,j),j=1,5)
enddo

print *,''
print *,'Matrix C='
do i=1,5
   write(*,'(1x,7f8.4)') (c(i,j),j=1,7)
enddo

print *,''
print *,'Vector d='
do i=1,5
   write(*,'(1x,f8.4)') d(i)
enddo

print *,''
print *,'A - C='
do i=1,5
   do j=1,7
      result(i,j)=a(i,j)-c(i,j)
   enddo
   write(*,'(1x,7f8.4)') (result(i,j),j=1,7)
enddo

print *,''
print *,'Bd='
do i=1,7
   result(i,1)=0
   do j=1,5
      result(i,1)=result(i,1)+b(i,j)*d(j)
   enddo
   write(*,'(1x,f8.4)') result(i,1)
enddo

print *,''
print *,'AB='
do i=1,5
   do j=1,5
      result(i,j)=0
      do k=1,7
	 result(i,j)=result(i,j)+a(i,k)*b(k,j)
      enddo
   enddo
   write(*,'(1x,7f8.4)') (result(i,j),j=1,5)
enddo

print *,''
print *,'BC='
do i=1,7
   do j=1,7
      result(i,j)=0
      do k=1,5
         result(i,j)=result(i,j)+b(i,k)*c(k,j)
      enddo
   enddo
   write(*,'(1x,7f8.4)') (result(i,j),j=1,7)
enddo
	
end

==========================================================
C++ Source Code:
==========================================================

	int i;
	int j;
	int k;

	double a[5][7]={0.7577, 0.706,  0.8235, 0.4387, 0.4898, 0.276,  0.4984,
					0.7431, 0.0318, 0.6948, 0.3816, 0.4456, 0.6797, 0.9597,
					0.3922, 0.2769, 0.3171, 0.7655, 0.6463, 0.6551, 0.3404,
					0.6555, 0.0462, 0.9502, 0.7952, 0.7094, 0.1626, 0.5853,
					0.1712, 0.0971, 0.0344, 0.1869, 0.7547, 0.119,  0.2238};

	double b[7][5]={0.4314, 0.5797, 0.4018, 0.0497, 0.3692,
					0.9106, 0.5499, 0.076,  0.9027, 0.1112,
					0.1818, 0.145,  0.2399, 0.9448, 0.7803,
					0.2638, 0.853,  0.1233, 0.4909, 0.3897,
					0.1455, 0.6221, 0.1839, 0.4893, 0.2417,
					0.1361, 0.351,  0.24,   0.3377, 0.4039,
					0.8693, 0.5132, 0.4173, 0.9001, 0.0965};

	double c[5][7]={0.1299, 0.1622, 0.602,  0.4505, 0.8258, 0.1067, 0.8687,
					0.5688, 0.7943, 0.263,  0.0838, 0.5383, 0.9619, 0.0844,
					0.4694, 0.3112, 0.6541, 0.229,  0.9961, 0.0046, 0.3998,
					0.0119, 0.5285, 0.6892, 0.9133, 0.0782, 0.7749, 0.2599,
					0.3371, 0.1656, 0.7482, 0.1524, 0.4427, 0.8173, 0.8001};

	double    d[5]={0.1320,
					0.9421,
					0.9561,
					0.5752,
					0.0598};

	double x[5][7];
	double y[7];
	double z[5][5];
	double w[7][7];

/* (ii)A-C */

	for(i=0; i<5; i++){
		for(j=0; j<7; j++){
			x[i][j]=a[i][j]-c[i][j];
		}
	}

/* (iv)Bb */

	for(i=0; i<7; i++){
		y[i]=0;
		for (k = 0; k < 5; k++) {
			y[i]+=b[i][k]*d[k];
		}
	}

/* (v)AB */

	for(i=0; i<5; i++){
		for(j=0; j<5; j++){
			z[i][j]=0;
			for (k = 0; k < 7; k++) {
				z[i][j]+=a[i][k]*b[k][j];
            		}
		}
	}

/* (vii)BC */

	for(i=0; i<7; i++){
		for(j=0; j<7; j++){
			w[i][j]=0;
			for (k = 0; k < 5; k++) {
				w[i][j]+=b[i][k]*c[k][j];
            		}
		}
	}

	printf("Matrix A-C: \n");
	for(i=0; i<5; i++){
		for(j=0; j<7; j++){
			printf(" %7.4f", x[i][j]);
		}
		printf("\n");
	}

	printf("Matrix Bb: \n");
	for(i=0; i<7; i++){
		for(j=0; j<1; j++){
			printf(" %7.4f", y[i]);
		}
		printf("\n");
	}

	printf("Matrix AB: \n");
	for(i=0; i<5; i++){
		for(j=0; j<5; j++){
			printf(" %7.4f", z[i][j]);
		}
		printf("\n");
	}

	printf("Matrix BC: \n");
	for(i=0; i<7; i++){
		for(j=0; j<7; j++){
			printf(" %7.4f", w[i][j]);
		}
		printf("\n");
	}

	return 0;		

==========================================================
Solutions
==========================================================

Matrix A=
   .7577   .7060   .8235   .4387   .4898   .2760   .4984
   .7431   .0318   .6948   .3816   .4456   .6797   .9597
   .3922   .2769   .3171   .7655   .6463   .6551   .3404
   .6555   .0462   .9502   .7952   .7094   .1626   .5853
   .1712   .0971   .0344   .1869   .7547   .1190   .2238

Matrix B=
   .4314   .5797   .4018   .0497   .3692
   .9106   .5499   .0760   .9027   .1112
   .1818   .1450   .2399   .9448   .7803
   .2638   .8530   .1233   .4909   .3897
   .1455   .6221   .1839   .4893   .2417
   .1361   .3510   .2400   .3377   .4039
   .8693   .5132   .4173   .9001   .0965

Matrix C=
   .1299   .1622   .6020   .4505   .8258   .1067   .8687
   .5688   .7943   .2630   .0838   .5383   .9619   .0844
   .4694   .3112   .6541   .2290   .9961   .0046   .3998
   .0119   .5285   .6892   .9133   .0782   .7749   .2599
   .3371   .1656   .7482   .1524   .4427   .8173   .8001

Vector d=
   .1320
   .9421
   .9561
   .5752
   .0598

A - C=
   .6278   .5438   .2215  -.0118  -.3360   .1693  -.3703
   .1743  -.7625   .4318   .2978  -.0927  -.2822   .8753
  -.0772  -.0343  -.3370   .5365  -.3498   .6505  -.0594
   .6436  -.4823   .2610  -.1181   .6312  -.6123   .3254
  -.1659  -.0685  -.7138   .0345   .3120  -.6983  -.5763

Bd=
  1.0379
  1.2368
   .9801
  1.2620
  1.0770
   .7965
  1.5207

AB=
  1.7773  1.9784   .9740  2.4498  1.4497
  1.5681  1.8828  1.1603  2.2208  1.4436
  1.1600  1.8853   .7672  1.7887  1.1750
  1.3415  2.0202  1.0066  2.2912  1.5921
   .5384   .9432   .3682   .8313   .4258

BC=
   .6994   .7429   .9855   .4366  1.2359   .9458   .8926
   .5150  1.1036  1.4479  1.3151  1.2435  1.4168  1.1914
   .4930   .8479  1.5395  1.1308   .8865  1.5298  1.1360
   .7145  1.0827  1.0937   .7263  1.0107  1.5481   .7898
   .5464   .8736   .8896   .6435   .7835  1.1915   .5730
   .4702   .6209   .8662   .5157   .7456   .9450   .6547
   .6440  1.1702  1.6238  1.3670  1.5229  1.3647  1.2765