Page 12 - Demo

Planar Truss Example for Comrel Add-on RCP Consult, 2023-2026 Page 12end doreturnendsubroutine MATMLT(a,b,c,m,n,l)!.....Multiplies the matrixes A(M,N) and B(N,L)!.....Result will be stored in matrix C(M,L).!.....Calling: dimension a(m,n),b(n,l),c(m,l)!..... call MATMLT(a,b,c,m,n,l)implicit real*8(a-h,o-z)dimension a(m,n),b(n,l),c(m,l)do k=1,l do i=1,m d=0. do j=1,n d=d+a(i,j)*b(j,k) end do c(i,k)=d end doend doreturnendsubroutine LINSOL(a,b,c,m)!.....Solves a system of linear equations A*B=C!.....Result will be stored in vector B(M).!.....Calling: dimension a(m,m),b(m),c(m)!..... call LINSOL(a,b,c,m)implicit real*8(a-h,o-z)dimension a(m,m),b(m),c(m)!.....https://stackoverflow.com/questions/37701747/program-to-solve-a-system-of-linear-equations-in-c!.....Triangularizationdo i=1,m-1 do k=i+1,m t=a(k,i)/a(i,i) do j=1,m a(k,j)=a(k,j)-t*a(i,j) end do c(k)=c(k)-t*c(i) end doend do!.....Resolutiondo i=m,1,-1 b(i)=c(i) do j=m,i+1,-1 b(i)=b(i)-a(i,j)*b(j) end do b(i)=b(i)/a(i,i)end doreturnend