عضویت در کانال و سوپر گروه مهندسین راه و ساختمان

http://telegram.me/joinchat/BULTFT3mvgx-z6Vq3rLs3Q

http://telegram.me/joinchat/CivilStructure_Surveying

نقل قول :

این کد اول m1 و n1 یا همان تعداد سطر و ستون ماتریس اول را میگیرد و سپس m2 و n2 مربوط به ماتریس دوم را میگیرد.بعد بررسی میکند که n1 و m2 برابر هستند یانه.برای ضرب ماتریس ها این دو عبارت باید برابر باشند.سپس یک ماتریس m1 در n2 تشکیل میدهد که همان ماتریس جواب ماست و شروع به محاسبه میکند و در آخر هم ماتریس بدست امده را نشان میدهد.…

program matrix
implicit none
INTEGER::m1,n1,m2,n2,t,k1,k2,i,k,j
REAL,allocatable::a(:,:),b(:,:),c(:,:),c1(:),c2(:)
REAL::s
PRINT*
PRINT*,"   _                _     _                _     _                _     "
PRINT*,"  | a11  a12 ..a1n1  |   | b11  b12 ..b1n2  |   | c11  c12 ..c1n2  |    "
PRINT*,"  |                  |   |                  |   |                  |    "
PRINT*,"  | a21  a22 ..a2n1  | X | b21  b22 ..b2n2  | = | c21  c22 ..c2n2  |    "
PRINT*,"  | .                |   | .     .     .    |   | .     .     .    |    "
PRINT*,"  | .                |   | .     .     .    |   | .     .     .    |    "
PRINT*,"  |_am11 am12..am1n1_|   |_bm21 bm22..bm2n2_|   |_cm11 cm12..cm1n2_|    "
PRINT*,"                    m1xn1                  m2xn2                  m1xn2 "
PRINT*
PRINT*,"------------------------------------------------------------------------"
PRINT*
PRINT*,"input m1 and n1 : "
READ*,m1,n1
PRINT*,"input m1 and n2 : "
READ*,m2,n2
if (n1/=m2) then
PRINT*," n1 must be equal m2"
PRINT*
stop
end if
ALLOCATE(a(m1,n1),b(m2,n2),c(m1,n2),c1(n1),c2(m2))
PRINT*
do k1=1,m1
PRINT*,"input line",k1,"of matrix a :"
READ*,a(k1,:)
end do
do k2=1,m2
PRINT*,"input line",k2,"of matrix b :"
READ*,b(k2,:)
end do
do i=1,m1
do j=1,n2
c1=a(i,:)
c2=b(:,j)
s=0.
do k=1,m2
s=s+(c1(k)*c2(k))
end do
PRINT*
c(i,j)=s
end do
end do
PRINT*," C(",m1,",",n2,") :"
PRINT*
do t=1,m1
PRINT*,"  ",c(t,:)
end do
PRINT*
end program