A=M1+M2+M3+M4+M5
B=(M1*y1+M2*Y2+M3*Y3+M4*Y4+M5*Y5)/A

A = 0
do i=1,5
   A = A+M(i)
end do

B=0
do i=1,5
   B = B+ M(i)*y(i)
end do

B = B/A

A=M1+M2+M3+M4+M5
B=(M1*y1+M2*Y2+M3*Y3+M4*Y4+M5*Y5)/A
-----------------------------------------
能不用循环就不用:

A = sum(M(1:5))
B = sum(M(1:5)*y(1:5))/A

或者

B = dot_product(M(1:5),y(1:5))/A