39 lines
1.7 KiB
Fortran
39 lines
1.7 KiB
Fortran
|
program test1a
|
||
|
|
implicit none
|
||
|
|
|
||
|
|
integer, parameter :: order = 10
|
||
|
|
integer i, pivot(order), ok
|
||
|
|
|
||
|
|
real, dimension(order,order) :: A = reshape( (/ &
|
||
|
|
196.356483, 255.783202, 254.495740, 248.699780, 382.610650, &
|
||
|
|
62.171758, 210.460418, 792.184862, 492.866147, 16.418566, &
|
||
|
|
465.590320, 419.182266, 276.934004, 92.750120, 633.946433, &
|
||
|
|
363.964274, 330.618290, 458.146049, 752.371999, 612.926015, &
|
||
|
|
232.746176, 47.023783, 46.899871, 780.529468, 40.868139, &
|
||
|
|
923.453047, 793.390477, 450.089661, 455.645932, 212.431431, &
|
||
|
|
56.827484, 66.928259, 977.022044, 688.254626, 618.276179, &
|
||
|
|
9.551896, 863.743164, 602.405740, 787.380819, 422.227709, &
|
||
|
|
374.758678, 634.315885, 432.129400, 787.574506, 573.625780, &
|
||
|
|
902.975285, 890.988258, 605.412751, 874.101587, 511.508295, &
|
||
|
|
4.822238, 522.909350, 256.688848, 122.516340, 900.301032, &
|
||
|
|
748.466705, 812.826296, 789.245172, 250.487399, 982.528843, &
|
||
|
|
378.805393, 602.230048, 364.264573, 654.350329, 622.863758, &
|
||
|
|
565.516688, 522.821921, 276.101853, 536.717666, 849.518816, &
|
||
|
|
476.899007, 895.548565, 398.298216, 45.213335, 535.759609, &
|
||
|
|
55.798948, 802.599529, 599.711219, 817.372219, 291.477778, &
|
||
|
|
486.338273, 596.013067, 442.259842, 149.362509, 592.750920, &
|
||
|
|
47.573334, 385.515622, 246.314796, 708.110189, 896.706744, &
|
||
|
|
676.303218, 896.507931, 523.151304, 582.015342, 609.873193, &
|
||
|
|
741.178150, 21.299662, 941.739146, 556.672861, 627.576228/), (/order, order/) )
|
||
|
|
real, dimension(order) :: b = (/955.191884, 900.038298, 166.121084, &
|
||
|
|
164.253073, 767.644681, 457.815496, 792.551618, 232.543770, 254.975493, 666.654631/)
|
||
|
|
|
||
|
|
|
||
|
|
call SGESV(order, 1, A, order, pivot, b, order, ok)
|
||
|
|
|
||
|
|
do i=1, order
|
||
|
|
write(*,*) b(i)
|
||
|
|
end do
|
||
|
|
|
||
|
|
end program test1a
|