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

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

http://telegram.me/joinchat/CivilStructure_Surveying

نقل قول :

کد فرترن روش ماتریس معکوس برای حل دستگاه معادلات جبری خطی...این کد اول تعداد معادلات یا n را میگیرد و ماتریس ضرایب n در n رو تشکیل میدهد و همچنین دو ماتریس n در 1 را تشکیل میدهد که یکی برای مجهولات و دیگری برای معلومات در طرف راست تساوی هست.بعد خط به خط ماتریس ضرایب و ماتریس معلوم را میگیرد و شروع به محاسبه ماتریس همسازه ها میکند.اول دترمینان ماتریس ضرایب را محاسبه میکنه که اگر 0 باشه یک پیام خطا نمایش میده که ماتریس معکوس پذیر نیست.اگر دترمینان 0 نباشد محاسبه کردن را ادامه میدهد.برای ماتریس همسازه ها باید با استفاده از کهاد سطر ها و ستون ها مرتب دترمینان محاسبه بشود که برای این مورد یک سابروتین نوشته شده است که این کار را میکند.سپس ترانهاده ماتریس همسازه ها را بدست می اورد و ضرب در دترمینان ماتریس ضرایب میکند.ماتریس بدست امده همان معکوس ماتریس ضرایب است.سپس این ماتریس معکوس شده را در ماتریس معلوم ضرب میکند که این کار را هم یک سابروتین انجام میدهد.جواب های بدست امده از ضرب این دو ماتریس همان جواب های مجهولات ما هستند و انها را نمایش میدهد.…

program matrix_inverse_method
implicit none
INTEGER::i2,i,j,k,p,i1,j1,n,j2,i3
REAL::s,detr,deta
REAL,ALLOCATABLE::a(:,:),b(:,:),c(:,:),f(:),x(:,:),t(:,:),g(:,:)
PRINT*
PRINT*,"   _                _     _  _     _  _     "
PRINT*,"  | a11  a12 ..a1n1  |   | x1 |   | b1 |    "
PRINT*,"  |                  |   |    |   |    |    "
PRINT*,"  | a21  a22 ..a2n1  | X | x2 | = | b2 |    "
PRINT*,"  | .                |   | .  |   | .  |    "
PRINT*,"  | .                |   | .  |   | .  |    "
PRINT*,"  |_am11 am12..am1n1_|   |_xn_|   |_bn_|    "
PRINT*,"                    nxn       nx1       nx1 "
PRINT*
PRINT*,"------------------------------------------------------------------------"
PRINT*

در ادامه مطلب