完全ビポッティング法 -完全ビポッティング法によるガウスの消去法のプログラ-

ガウスの消去

ガウスの消去法のプログラムを作ったのですがうまく動きませんどこが間違っているのでしょうちなみに連立方程式を解くプログラムです * the numerical solution of linear equation * by gauss method parameter(ll=10) dimension a(ll,ll+1), x(ll) read(5,100) n 100 format(i3) do 150 i = 1, n 150 read(5,200) (a(i,j), j=1, n+1) 200 format(11f4.1) write(6,300) 300 format(' ', 10x, 'coefficient') do 10 i = 1, n write(6,400) (a(i,j), j=1, n+1) 400 format(' ', 11 (5x, f4.1)) 10 continue esp=10.0e-19 call rgaule(a, x, ll, n, esp, ipivot) if (ipivot. eq. 1) go to 20 write(6,500) 500 format(' ', 10x, 'the pivot is little '/ 1 'so the solution is singular') go to 110 20 write(6,600) 600 format(/' ',2x, 8hsolution) write (6,700) (i, x(i), i = 1, n) 700 format (' ',5x, 'x (', i2, ' ) =', 2x, e14.7) 110 stop end subroutine rgaule(a, x, ll, n, esp, ipivot) dimension a(ll, ll+1), x(ll), ln(100) ipivot = 1 * the order of x do 5 i = 1, n ln(i) = i 5 continue do 100 m = 1, n-1 * the selection of pivot amax = 0.0 do 10 i = m, n do 20 j = m, n if (amax. ge. abs(a(i,j))) go to 20 amax = abs(a(i, j)) irow = i icolum = j 20 continue 10 continue if (amax. le. eps) go to 70 if (m. eq. irow) go to 22 * the exchange of row do 27 l = m, n+1 swap = a(irow, l) a(irow, l) = a(m, l) a(m, l) = swap 27 continue 22 if (m. eq. icolum) go to 30 * the exchange of colum do 25 l =1, n swap = a(l, icolum) a(l, icolum) = a(l,m) a(l, m) = swap 25 continue * the exchange of x iswap = ln(m) ln(m) = ln(icolum) ln(icolum) = iswap * gaussian elimination 30 do 35 i = m+1, n do 37 j = m+1, n+1 a(i,j) = a(i,j) - a(i,m) * a(m,j) / a(m,m) 37 continue 35 continue 100 continue if (abs(a(n,n)). le. eps) go to 70 * back subsititution x(n) = a(n, n+1) / a(n,n) kk = ln(n) a(n, kk) = x(n) do 50 i = n-1, 1, -1 k = n-i x(i) = 0.0 do 52 j = 1, k ll = i + j x(i) = a(i, ll) * x(ll) + x(i) 52 continue x(i) = (a(i, n+1) - x(i)) / a(i,i) kk = ln(i) a(n, kk) = x(i) 50 continue do 60 i = 1, n x(i) = a(n, i) 60 continue return 70 ipivot = 0 return end

連立方程式をSOR法を用いてコンピュータで解く

{ 2 1 1 }{ X1 }　{ 2 } { 2 3 1 }{ X2 }＝{ 4 } { 1 1 3 }{ X3 } {-1 } という問題をSOR法を用いてコンピュータ（fortran）で解け。ただし初期値は０とし許容誤差eps=10^-5とする。そこで私は次のようなプログラムで解きました。 ____dimension a(20,20),b(20),x(20),x0(20) ____open(5,file='senkei.d') ____open(6,file='senkei.r') ____open(7,file='senkei.k') ____read(5,*)n,max,eps,u ____read(5,*)((a(i,j),i=1,n),j=1,n) ____read(5,*)(b(i),i=1,n) ____read(5,*)(x0(i),i=1,n) ____do i=1,n ____ad=a(i,i) ____a(i,i)=0.0 ____b(i)=b(i)/ad ____do j=1,n ____a(i,j)=a(i,j)/ad ____end do ____end do ____do k=1,max ____do i=1,n ____w=0.0 ____do j=1,n ____if(j.lt.i)x0(j)=x(j) ____w=w+a(i,j)*x0(j) ____end do ____x(i)=u*(b(i)-w)+(1-u)*x(i) ____write(7,10)x(i) _10_format(3f10.5) ____end do ____w=0.0 ____do i=1,n ____w=w+(x(i)-x0(i))**2 ____end do ____w=sqrt(w/n) ____if(w.lt.eps)go to 20 ____do i=1,n ____x0(i)=x(i) ____end do ____end do ____write(6,*)'error' ____go to 30 _20_write(6,*)'sulution' ____write(6,40)(x(i),i=1,n) _40_format(3f15.5) _30_stop ____end max:最大繰り返し回数,u:緩和係数,a(i,j):係数行列, b(i):既知項,x(i):未知項としました。データではeps=0.00001,max=50,緩和係数uは0.8から0.1ずつ増やして1.5までそれぞれ答えを求めました。私のプログラムではu=1.1の時7回で収束しました。しかしu=1.5の時3回で収束したのですが答えはかなりの誤差ありました。許容誤差eps=10^-5としたのになぜ誤差が大きいまま、しかも三回で収束してしまったのかわかりません。誰かわかる方いらっしゃったら教えていただけませんか？長々と申し訳ありません。お願いします。

何処の会社のＢＡＳＩＣか教えて下さい。

何処の会社のＢＡＳＩＣか教えて下さい。下記のようなソースファイルです。宜しくお願いします。 DIMENSION EBK(20,20),ECK(20,20),SH(20),SL(20),CRM(20,20), 1CLM(20,20),Q(20),SOM(20),CDM(20,20),CUM(20,20),BML(20,20), 2BMR(20,20),BMU(20,20),BMD(20,20),DVL(20,20),DVR(20,20),DVU(20,20), 3DVD(20,20),RVU(20,20),RVD(20,20),TML(20,20),TMR(20,20), 4TMD(20,20),TMU(20,20),RMD(20,20),RMU(20,20),BBML(20,20), 5BBMR(20,20),BBMD(20,20),BBMU(20,20),TQL(20,20),TQR(20,20), 6TMQ(20,20),QLL(20,20),QLR(20,20),CMM(20,20),TMM(20,20) READ(5,10) M,L READ(5,10) N1,L1 N12=Nl+2 N13=N1+1 L12=Ll+2 L13=L1+1 10 FORMAT(2I5) DO 15 I=1,N12 DO 15 J=1,L12 EBK(I,J)=0.0 ECK(I,J)=0.0 15 CONTINUE READ(5,11)((EBK(I,J),J=2,L13),I=2,N12) READ(5,11)((ECK(I,J),J=2,L12),I=2,N13) 11 FORMAT(7F10.2) READ(5,11)(SH(I),I=2,N13) READ(5,11)(SL(I),I=2,L13) IF(L-1)20,30,30 20 DO 21 I=2,N12 DO 21 J=1,L12 READ(5,12)CLM(I,J),CRM(I,J),CMM(I,J),QLL(I,J),QLR(I,J) 21 CONTINUE 12 FORMAT(5F10.2) DO 40 I=2,N12 DO 40 J=1,L12 CLM(I,J)=-CRM(I,J) 40 CONTINUE

FORTRAN→Cに翻訳

　どなたか、次のFORTRANのプログラムを、Cに、翻訳して頂けないでしょうか。C++ではなく、Cです。ANSI準拠のCでお願いします。　プログラムの内容は、最小二乗法による計算プログラムです。MS-DOS Ver3.3～6.0の頃の、MS FORTRANコンパイラ仕様のものです。その頃持っていたFORTRANの本も処分してしまい、今からFORTRANを学びなおすのにも多大な労力と時間がかかりそうなので、Cに翻訳して頂ければ大変ありがたいです。よろしくお願いします。（“□”はタブ） ◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆ C□LEAST SQUARE APPROXIMATION □PROGRAM MAIN9 □DIMENSION X(100),Y(100),S(0:18),T(0:9),SM(10,10),TV(10),AV(10) □WRITE(*,*) 'N ?' □READ(*,*) N □WRITE(*,*) 'x1,x2,..,xn ?' □READ(*,*) ( X(I),I=1,N ) □WRITE(*,*) 'y1,y2,..,yn ?' □READ(*,*) ( Y(I),I=1,N ) □WRITE(*,*) 'M ?' □READ(*,*) M □DO 110 K=0,M*2 □□VS=0. □□DO 100 I=1,N □100□VS=VS+X(I)**K □□S(K)=VS □110□CONTINUE □□DO 130 K=0,M □□□VS=0. □□□DO 120 I=1,N □120□VS=VS+Y(I)*X(I)**K □□□T(K)=VS □130 CONTINUE □□DO 140 I=1,M+1 □□□DO 140 J=1,M+1 □□□□K=I+J-2 □□□□SM(I,J)=S(K) □140 CONTINUE □□DO 150 I=1,M+1 □150 TV(I)=T(I-1) □□CALL SIMULE( AV, SM, TV, M+1 ) □□DO 160 I=1,M+1 □160 WRITE(*,1000) I-1,AV(I) □1000 FORMAT(1H ,'A',I1,'=',F10.5) □□STOP □□END ◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆◆

FORTRAN 初心者です

以下の連立一次方程式をＳＯＲ法で解く問題です。初心者なりにガウスザイデル法を応用してプログラムしたつもりですが、やはり難しいです（答えは違います）。どこをどうすれば良いのか分かりませんので、よろしければヒントや助言をいただきたいです。 PROGRAM SOR REAL A(10,10),B(10),X(10),X0(10) INTEGER N,I,J,K,Kmax,w N=3 A(1,1)=4 ;A(1,2)=1 ;A(1,3)=2 A(2,1)=1 ;A(2,2)=3 ;A(2,3)=1 A(3,1)=1 ;A(3,2)=2 ;A(3,3)=5 B(1)=16 B(2)=10 B(3)=12 X0(1)=1 ;X0(2)=1 ;X0(3)=2 w=1.2 Kmax=50 EPS=1.D-5 DO I=1,N D=A(I,I) S=B(I) B(I)=B(I)/D END DO DO K=1,Kmax DO I=1,N DO J=1,N if(J<I) X0(J)=X(J) S=S-A(I,J)*X0(J) END DO X(I)=(1-w)*X(I)+w*S END DO DO I=1,N S=S-(X(I)-X0(I))**2 END DO IF(S<EPS) GOTO 10 DO I=1,N X0(I)=X(I) END DO END DO 10 WRITE(*,*) K DO I=1,N WRITE(*,*) 'SOR法で求めた解は' WRITE(*,*) 'X(',I,')=',X(I) END DO END PROGRAM SOR !------------------------------------ ※wは緩和係数です

計算速度

a(m,n) m=20000 n=10000 のとき do j=1,n do i=1,m x=x+a(i,j) end do end do の計算速度と do j=1,n do i=1,m x=x+a(j,i) end do end do の計算速度が違うのはなぜですか？

fortran　errorについて

fortranを勉強していたのですがエラーがでてしまい、何時間かけても理解できなかったので質問させてください。以下プログラム program test !ここからメインルーチン !前準備　配列の用意 implicit none integer N integer,dimension(0:N,0:N) :: A integer :: i,j,k read * ,N !初期状態の代入 do i=0,N do j=0,N A(i,j)=0 end do end do do i=N/2,N-1 A(N/2,i)=1 end do do i=N/2,N-1 A(N/2+1,i)=-1 end do !ループ　50回ループさせる do k=0,50 !状態の表示 call visualize !サブルーチン visualize subroutine visualize do i=0,N do j=0,N if(A(i,j)== 1) write(*,'(A1)',advance='NO') "*" if(A(i,j)== 0) write(*,'(A1)',advance='NO') " " if(A(i,j)==-1) write(*,'(A1)',advance='NO') "+" end do write(*,*) end do !end subroutine visualize call insert !サブルーチン insert subroutine insert do i=0,N do j=0,N if(A(i,j)== 1) A(i,j)=-1 if(A(i,j)== 0) A(i,j)=max(0,A(i-1,j),A(i,j-1),A(i,j+1),A(i+1,j)) if(A(i,j)==-1) A(i,j)=0 end do end do !end subroutine insert end do end program test これでコンパイラすると In file test.f90:48 subroutine visualize 1 Error: Unclassifiable statement at (1) In file test.f90:69 subroutine insert 1 Error: Unclassifiable statement at (1) とでますいろいろ調べたのですが全くわかりませんでしたできればよろしくお願いします

Delphi6で、実行後、一時停止をしたい。

Delphi6を使ってPascalを勉強しています。 Delphi6でプログラムを記述し実行すると、MS－DOSが出てきて実行します。しかし、実行した後すぐに終了してしまうので困っています。合っているのか確認できません。そのために一時停止用に型宣言のときに、　x:string;　 end.の前に　　readln(x)　　を付けて止めるようにしているのですが、プログラムによって止まるのと止まらないのがあります。どうしてなのでしょうか。 ↓は止まりません。どうしたら止まってくれますか？配列にデータを追加するプログラムです。 program ng; {$AppType Console} uses SysUtils; var 　i, j, n, dat : integer; 　a : array[1..100] of integer; 　{一時停止のため変数名} 　x:string; begin {配列にデータを入力する} 　n := 6; 　a[1] := 20; a[2] := 30; a[3] := 60; 　a[4] := 70; a[5] := 80; a[6] := 90; 　writeln('（１）現在の配列の内容'); 　for i := 1 to n do write(a[i]:3); 　writeln; 　writeln('（２）挿入するデータ？'); 　read(dat); 　repeat 　　　writeln('（３）挿入する位置？'); 　　　read(i); 　until(i < n) and (i > 0); {配列への挿入} 　for j := n downto i do a[j+1] := a[j]; 　a[i] := dat; 　n := n+1; {挿入後の配列を表示} 　writeln('（４）挿入後の配列の内容'); 　for i := 1 to n do write(a[i]:3); 　{一時停止のためのread文} 　readln(x) end.

並べ替え

次のようなプログラムを作れ。 1.数値を次々に読み込んで配列に1番目から順に格納する。 2.配列に格納されている数値を逆の順番に並べ換える (最後に読み込んだ数値が配列の1番目に入る)。 3.その配列の内容を1番目から順番に画面に出力する。読み込み、逆転、出力の部分はそれぞれ関数にすること。という課題をやっています。下のように作ってみたんですが上手くいかなくて・・・何処が悪いか教えて下さい。 #include <stdio.h> #define MAX 100 void reve(int *, int *); void read(int *, int *); void write(int *,int *); int main () { int x[MAX]; int n=0; read(&n, x); reve(&n, x); write(&n, x); return; } void read(int *n, int x[]) { while(scanf("%d", &x[n]) == 1 ) { n += 1; } return; } void reve(int *n, int x[]) { int i, j, y; for(i = 0, j = n-1; i <j; i++, j--) { y = x[j]; x[j] = x[i]; x[i] = y; } return; } void write(int *n, int x) { int i; for(i = n; i = 0; i--) { printf("x[%d] = %d\n", i, x[i]); } return; }

fscanfの書式について

fscanfの書式について a[n][n]が、 a11,a12,a13---a1n, a21,----------a2n, . . an1,----------ann, という風に並んでいるファイルを読み込む場合、どのように記述すれば良いのでしょうか？例えばFortranでは、 DO I=1,N READ(11,'(N(F8.4,A1))')(A(I,J),KAMA, I=1,N) END DO と書けば良いのですが、Cでの書き方が分りません。元々Fortranを使いでCは初心者です。ちなみに、 for(i=0;i<=n+1;i++){ for(j=0;j<=n+1;j++){ fscanf(fp,"%f7.6 %1s ",&a[i][j],kama); } } と書いて見たものの、駄目でした。改行してしまうのだと思うのですが、避ける方法も、そもそも、それが出来るのかどうかも分りません。。よろしくお願いいたします。

fortran 行列ベクトル積

行列ベクトル積を計算するプログラムを下のように書いたのですが、実行した結果の答えが実際計算した答えと異なります。初期の要素の設定がおかしいのでしょうか？教えて下さい。よろしくお願いします。 program list2_14 implicit none integer , parameter :: n = 2 real(8) a(n,n), x(n), y(n) integer i, j, k, l a(1,1:2) = (/1.2d0,3.4d0/) a(2,1:2) = (/5.6d0,7.8d0/) x(:) = (/9.0d0,10.0d0/) do i = 1, n y(i) = 0.0d0 do j = 1, n y(i) = y(i) + a(i,j) * x(j) enddo enddo do k = 1, n write(*,*) (a(k,l), l = 1, n) enddo write(*,*) x(:) write(*,*) y(n) end program list2_14 実行結果 1.2 3.4 5.6 7.8 9. 10. 128.4

Fortranについて教えてください！

こんにちわ。ヤマといいます。以下について教えて頂けないでしょうか？ ------------------- 途中省略 ------------------- DO 31 I = 1,NCARDS LC = LC + 1 READ(8,FMAT) (XR(J), J=1,NPL) IF(I.NE.NCARDS) GO TO 6 IF(JL.EQ.0) GO TO 6 JL=NPL+1-JL DO 5 J=JL,NPL 5 XR(J)=0. C--------------------------------------------------------------- C ONLY PRINT OUT A FEW OF THE FIRST AND THE LAST LINES of Input Motion C--------------------------------------------------------------- 6 ICHECK = NCARDS - I IF (I .LE. 5 .OR. ICHECK .LT. 5) WRITE(6,2008) I,(XR(J), J=1,NPL) IF (I .EQ. 10) WRITE (6,2009) 2009 FORMAT(3X,'........ INPUT MOTION READ NOT ECHOED...........') C ENDIF C C FIND MAX. INPUT ACC. (XMAX) C 311 DO 31 J = 1,NPL,2 N = N + 1 X(N) = CMPLX(XR(J),XR(J+1)) 31 CONTINUE ........................................... とあります。いま、XRは、 ----------------データ------------------------ -0.001694 -0.001668 -0.000086 -0.001356 -0.000678 0.000700 -0.001209 -0.000604 ・・・・・・・・・・・・・・・・・・・・・ ----------------データ------------------------ と1行に8列分のデータが存在します。ソース中に「CMPLX(XR(J),XR(J+1))」とあります。この意味がよく分かりません。通常だと実部にXR(J),虚部にXR(J+1)を入れるというように解釈できるのですが、本当でしょうか？ご存知の方ご意見よろしくお願いいたします。

Ｃ言語への翻訳お願いします

以下の文はある本から抜粋してきたガウス・ジョルダン法のサブルーチンです。Ｃに翻訳してくださいって言ったら怒りますか？（笑 c Gauss-Jordan method subroutine gj(a,n1,nn,m,epsl,isw) double precision a,epsl,p,w dimension a(n1,*) nplsm=nn+m do 100 n=1,nn p=0.0 do 10 i=n,nn if(p.lt.abs(a(i,n))) then p=abs(a(i,n)) ip=i end if 10 continue if(p.le.epsl) then isw=1 write(6,*) 'error' return end if do 20 j=n,nplsm w=a(n,j) a(n,j)=a(ip,j) 20 a(ip,j)=w do 30 j=n+1,nplsm 30 a(n,j)=a(n,j)/a(n,n) do 40 i=1,nn if(i.ne.n) then do 45 j=n+1,nplsm 45 a(i,j)=a(i,j)-a(i,n)*a(n,j) end if 40 continue 100 continue isw=0 return end だいたい自分でもＣに書き換えることはできたんですが、次のところがどのように書き換えたらよいのか分かりません。 subroutine gj(a,n1,nn,m,epsl,isw), dimension a(n1,*), do 20 j=n,nplsm w=a(n,j) a(n,j)=a(ip,j) 20 a(ip,j)=w どうかよろしくお願いします。

プログラムに内容と計算の質問です。

こんにちは。補間多項式についての、コンピュータのプログラムの解読に困っています。内容は、「For the numerical experiments suggested in the computer problems, the following two procedures should be satisfactory. The first is called Coef. It requires as input the number n and tabular values in the array {Xi} and {Yi}. Remember that the number of points in the table is n＋1. The procedure then computes the　coefficients required in the Newton interpolating polynomial, storing them in the array{Ai}. -------------------------------------------------------- procedure; Coef(n,{Xi},{Yi},{Ai}) real array; {Xi}0:n, {Yi}0:n, {Ai}0:n integer; i,j,n for i=0 to n do {Ai}←{Yi} end for for j=1 to n do for i=n to j step -1 do Ai←({Ai}－{Ai-1})/({Xi}－{Xi-j}) end for end for end procedure Coef --------------------------------------------------------- このプログラムのn=3の時を考えるとき、　(1)j=1のとき、i=3,2,1 <j=1,i=1> {A1}=({A1}-{A0})/({X1}-{X0}) =({Y1}-{Y0})/({X1}-{X0}) <i=1,i=2> {A2}=({A2}-{A1})/({X2}-{X1}) =({Y2}-{Y1})/({X2}-{X1}) <i=1,i=3> {A3}=({A3}-{A2})/({X3}-{X2}) =({Y3}-{Y2})/({X3}-{X2}) (2)j=2のとき、i=3,2 <j=2,i=1> A1=({A2}-{A1})/({X2}-{X0}) 　　　　　　　　 ={[({Y2}-{Y1})/({X2}-{X1})]-[({Y1}-{Y0})/({X1}-{X0})]}/({X2}-{X0}) =この式変形をしたいのですが、どのように　　　　　　　　すれば良いのかわかりません。ラグランジェ　　　　　　　　型になりそうでなりません（泣）　　　　　　　　(1)で求めた{A1},{A2},{A3}を使って求めな　　　　　　　　いといけないみたいです。見にくい表し方で申し訳ありません。アドバイスお願いします！！

(Σa_n・x^n)^m

mを自然数として(Σ[n=0↑∞]a_n・x^n)^mが収束する場合にこれをべき級数で表した時のx^kの係数の計算の仕方がよくわかりません。a_nやxは実数とします。 Σ[n=0↑∞]Σ[n=n_1+n_2+…+n_m]a_n_1・a_n_2・…・a_n_m・x^nとして a_n_1・a_n_2・…・a_n_m=a_0^i_0・a_1^i_1・…・a_j^i_j・… と表すと有限個のjについてi_j>0でΣ[j=0↑∞]i_j=mであってnを固定するとこの係数をもつ項がm!/(i_1!・i_2!・…・i_n!)個あると考えればいいのかと思ったのですがこの推論は間違っているようです。別のやり方としてx=0でのk次微分係数を計算してk!で割ればいいと思ったのですが具体的な計算ができませんでした。

fortran エラーについて

fortranで、副プログラムを使ってデータを昇順または降順に並べ替えるプログラムを入力して実行しようとしたところ、・Unexpected junk in formal argument list at (1) ・Two main PROGRAMs at (1) and (2) という2つのエラーが出ました。これらの改善方法を教えて頂きたいです。初心者ですので簡単なところで間違えている可能性もありますが、ご指摘いただければ幸いです。以下、実際に入力したプログラムです。 ------------------------------ implicit none integer::i,n real::x(1000),a(1000),b(1000) n=1000 open(10,file='input-data-1.txt') do i=1,n read(10,*) x(i) end do close(10) open(10,file='output-data-1.txt') do i=1,n call koukan(i,x(i),a(i),b(i)) write(10,'(i4,2f10.3)') i,a(i),b(i) end do close(10) stop end subroutine koukan(i,x(i),shoujun,koujun) implicit none integer::i,n,made real::x(1000),w,shoujun,koujun do made=n-1,1,-1 do i=1,made if(x(i)>x(i+1)) then w=x(i) x(i)=x(i+1) x(i+1)=w end if end do end do shoujun=x(i) do made=n-1,1,-1 do i=1,made if(x(i)<x(i+1)) then w=x(i) x(i)=x(i+1) x(i+1)=w end if end do end do koujun=x(i) return end ------------------------------

FORTRAN…これってどんなプログラムになりますか？？

DO 10 I=1,47 CALL SUB1 10 CONTINUE STOP END SUBROUTINE SUB1 DIMENSION B1(3),B2(3),C(3),L(3),P(3) CHARACTER*12 A READ(5,50) A,B1,B2 50 FORMAT(A12,3F8.1,3F7.1) X=1.0 DO 11 K=1,300 Y1=(-1.0) Y2=0.0 DO 12 J=1,3 L(J)=(-NINT(B1(J)*10.0/B2(J))) Y1=Y1+X**L(J) Y2=Y2+L(J)*X**(L(J)-1) 12 CONTINUE W=X-Y1/Y2 IF(ABS(W-X).LT.1E-10) GO TO 13 X=W 11 CONTINUE 13 WO=W DO 14 J=1,3 C(J)=WO**L(J) 14 CONTINUE R1=0.0 DO 15 J=1,3 R2=R1+B2(J) R1=R2 15 CONTINUE D=0.0 DO 16 J=1,3 P(J)=B2(J)/R1 D=D+P(J)*ALOG(P(J)/C(J)) 16 CONTINUE E=0.0 DO 17 J=1,3 E=E+(B1(J)/B2(J)*P(J)) 17 CONTINUE WRITE(*,200) A,B2,P,E,D 200 FORMAT(1H,2X,A12,3X,3(F7.1,2X),4X,3(F9.6,X),4X,F9.6,2X,F9.6) RETURN END

プログラムです。

　次のプログラムは何をするものでしょうか？ただし、自然数が入力されるそうです。　　program test(input,output) var i,j,n:integer; begin read(n); s:=0; i:=1; j:=1; while i <=n do begin s:=s+j*i i:=i+1; j:=j*(-1) end writeln('n=' ,n:3,' のときs=',s:5) end 宜しくお願いします。

Jacobi法

Jacobi法のアルゴリズム x[i]_{m+1} = ( b[i] - Σ{j≠i} a[i][j] * x[j]_{m} ) i=1,2,…,n m=0,1,2,… を基にプログラム（JAVA）を組んでみました。（質問に直接関係ないところは省いています） for( m=0; m<100; m++ ) { for( i=0; i<n; i++ ) { for( j=0; j<i; j++) { if( i != j ) { tmp = tmp + a[i][j] * x2[j]; } } x1[i] = (b[i] - tmp) / a[i][i]; } } もう１つ。 for( k=0; k<100; k++ ) { for( i=0; i<n; i++ ) { xnew[i] = b[i]; for( j=0; j<n; j++ ) { if( j != i ) { xnew[i] = xnew[i] - a[i][j] * xold[j]; } } xnew[i] = xnew[i] / a[i][i]; } 2つ目のプログラムの方が正確です。 1つ目のプログラムだと間違いでしょうか？私には2つ目のプログラムの xnew[i]=b[i]; の部分がどうしても理解できません。これを行ってしまうと、アルゴリズムが x[i]_{m+1} = ( Σ{j≠i}b[i] - a[i][j] * x[j]_{m} ) になってしまわないでしょうか・・・？どちらでも良いのでしょうか？納得いきません。よろしくお願いします。

fortran 上三角行列

一様乱数を要素とする上三角行列を設定するプログラミングを行ったのですが以下のプログラムで call random_number(a(1 : j , j)) a(j+1 : n , j ) = 0.0d0 の部分を理解した上で実行したのですが、実行結果がそれぞれの列に同じ一様乱数が表示されました。プログラム内容と実行結果が納得のいく物ではないのですが、果たして以下のプログラムは上三角行列を設定する正しく実行される物となっているのでしょうか？教えて下さい。よろしくお願いします。 program list2_8 implicit none real(8), allocatable :: a(:,:) integer n, i, j write(*,'(a)', advance ='no') ' input n (1<=n<=100) :' read(*,*) n if (n < 1 .or. 100 < n) stop 'stop, n is invalid' allocate (a(n,n)) call random_seed do j = 1, n call random_number(a(1:j,j)) a(j+1:n,j) =0.0d0 enddo do i = 1, n write(*,'(100e12.4)') a(1, 1:n) enddo end program list2_8

完全ビポッティング法

みんなの回答

関連するQ&A

ガウスの消去

連立方程式をSOR法を用いてコンピュータで解く

何処の会社のＢＡＳＩＣか教えて下さい。

FORTRAN→Cに翻訳

FORTRAN 初心者です

計算速度

fortran　errorについて

Delphi6で、実行後、一時停止をしたい。

並べ替え

fscanfの書式について

fortran 行列ベクトル積

Fortranについて教えてください！

Ｃ言語への翻訳お願いします

プログラムに内容と計算の質問です。

(Σa_n・x^n)^m

fortran エラーについて

FORTRAN…これってどんなプログラムになりますか？？

プログラムです。

Jacobi法

fortran 上三角行列

注目のQ&A

カテゴリ

あなたにピッタリな商品が見つかる！ OKWAVE セレクト

専門家に質問してみよう

完全ビポッティング法

みんなの回答

関連するQ&A

ガウスの消去

連立方程式をSOR法を用いてコンピュータで解く

何処の会社のＢＡＳＩＣか教えて下さい。

FORTRAN→Cに翻訳

FORTRAN 初心者です

計算速度

fortran errorについて

Delphi6で、実行後、一時停止をしたい。

並べ替え

fscanfの書式について

fortran 行列ベクトル積

Fortranについて教えてください！

Ｃ言語への翻訳お願いします

プログラムに内容と計算の質問です。

(Σa_n・x^n)^m

fortran エラーについて

FORTRAN…これってどんなプログラムになりますか？？

プログラムです。

Jacobi法

fortran 上三角行列

注目のQ&A

カテゴリ

あなたにピッタリな商品が見つかる！ OKWAVE セレクト

専門家に質問してみよう

fortran　errorについて