くだすれFORTRAN（超初心者用）その4at TECHくだすれFORTRAN（超初心者用）その4 - 暇つぶし2ch■コピペモード□スレを通常表示□オプションモード□このｽﾚｯﾄﾞのURL■項目テキスト36:35 09/01/30 00:07:53 ! Searching the factor and its location which held the maximum number in a column j do j = 1,sizeN ValMax = 0.0d0 ImaxHeld = 0 do k = j, sizeN if( ValMax < abs( P(j,k) ) )then ImaxHeld = k ValMax = abs( P(j,k) ) end if end do write(*,*) 'ImaxHeld=', ImaxHeld, 'ValMax=', ValMax ! Two rows are exchanged so that the maximum factor is placed at a diagonal tmp(:) = P(:,j) P(:,j) = P(:,ImaxHeld) P(:,ImaxHeld) = tmp(:) ! Now the diagonal factor P(j,j) is a pibot P(:,j) = P(:,j) / P(j,j) ! <=(＊) do i = 1, sizeN if (i /= j) then P(:,i) = P(:,i) - P(j,i)*P(:,j) end if end do end do 37:デフォルトの名無しさん 09/01/30 01:55:43 >>35 一列全部０なら、シンギュラーな行列なんだから、そもそも解が無いだろ！ 線形代数の教科書を１００回読み直せ！ｗ 次ページ最新レス表示レスジャンプ類似スレ一覧スレッドの検索話題のニュースおまかせリストオプションしおりを挟むスレッドに書込スレッドの一覧暇つぶし2ch