くだすれFORTRAN（超初心者用）

このスレッドは、他のスレッドでは書き込めない超低レベル、 
もしくは質問者自身何が何だが分からない質問を勇気を持って書き込むスレッドです。 
FORTRAN使いが優しくコメントを返しますが、 
お礼はFORTRANの布教と初心者の救済をお願いします。 

　　　　　　　ﾛｰﾘﾝｸﾞ!!　　　 ∧∧　　
　　　　　　　　　　　　　　　(ﾟДﾟ,,)　　
　　　　　　　　　　　　　　　⊂⊂,,ヽ　　
　　　　　　　　　　　　　　　　(_ ( 　）ノ
　　　　　ｸﾙﾝ　　　　　　　
　　　　　　　　　　　 　　/⌒⌒ヽノ 　））)
　　　　　　　　　　　 　 （ 　　）て ）　　
　　　　　 　　　 (（（　　 ∨∨⊂ノ
ｽﾞｻｷﾞｺ!!　
　　　　　　　　　　　　　　　　　　　　(´´
　　　∧∧　　）　　≡≡≡≡≡(´⌒(´≡≡　
　⊂(ﾟДﾟ⊂⌒つ　≡≡≡≡(´⌒;;≡　
　　　　　　　　　　　(´⌒(´⌒;;
2げっとーーーーーーーーー！！！

なにがなんだが　×
なにがなんだか　○

すみません、私が悪かったです。コピペするとき
添削してくれたらよかったのですが・・・。

任意の生年月日と今日の日付を入力すれば、
産まれてから今日まで生きてきた日数を出力するプログラムを作れ。
さらに、任意の整数ｎを入力して、産まれてからｎ日目が、
いつかを出力させるプログラムを作れ。
ちなみに、閏年は以下のようにして決められる。

１）基本的に西暦が4で割り切れたら閏年。
２）しかし、西暦が100で割り切れたら閏年ではない。
３）しかし、西暦が400で割り切れたら閏年とする。

レポート作成
予備的考察、プログラム、プログラムの説明、計算結果、考察

フォートラン使用･･・


一応書き込んでみました。

少しあきらめられない自分がいました・・・。

>>4
Rosenmaidenを見ていたおれ様が来ましたよｗ
まあ少し待たれよｗ

MODULE mod_subs
IMPLICIT NONE
INTEGER, PARAMETER :: ndays_of_month(12, 2) = &
RESHAPE( (/31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31, &
31, 29, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31 /), (/12, 2/) )
CONTAINS
LOGICAL FUNCTION qleapyear(iyear)
IMPLICIT NONE
INTEGER, INTENT(IN) :: iyear
qleapyear = .FALSE.
IF ( MOD(iyear, 4) == 0 ) qleapyear = .TRUE.
IF ( MOD(iyear, 100) == 0 ) qleapyear = .FALSE.
IF ( MOD(iyear, 1000) == 0 ) qleapyear = .TRUE.
RETURN
END FUNCTION qleapyear
!
INTEGER FUNCTION idays_from_jan1(iyear, imonth, iday)
IMPLICIT NONE
INTEGER, INTENT(IN) :: iyear, imonth, iday
INTEGER :: i, ileap, n
ileap = 1
IF ( qleapyear(iyear) ) ileap = 2
n = 0
DO i = 1, imonth - 1
n = n + ndays_of_month(i, ileap)
END DO
idays_from_jan1 = n + iday
RETURN
END FUNCTION idays_from_jan1
END MODULE mod_subs

PROGRAM AliceGame
USE mod_subs
IMPLICIT NONE
INTEGER :: iyear0, imonth0, iday0, iyear1, imonth1, iday1
INTEGER :: iyear, ndays, ndays0, ndays1
!
WRITE(*, *) "input birth date : Year, Month, Day"
READ(*, *) iyear0, imonth0, iday0
!
WRITE(*, *) "input today's date : Year, Month, Day"
READ(*, *) iyear1, imonth1, iday1
!
ndays0 = idays_from_jan1(iyear0, imonth0, iday0)
ndays1 = idays_from_jan1(iyear1, imonth1, iday1)
ndays = -ndays0
DO iyear = iyear0, iyear1 - 1
IF ( qleapyear(iyear) ) THEN
ndays = ndays + 366
ELSE
ndays = ndays + 365
END IF
END DO
ndays = ndays + ndays1
WRITE(*, *) "The number of days from birth =", ndays
STOP
END PROGRAM AliceGame

専ブラ　インデント付き用
>>6-7

>>4
結果を余りチェックしていないので、間違ってたらごめんｗｗｗｗ謝罪も賠償もしないｗｗ

とりあえず、課題の前半だ。後半は眠いから作らないｗ

紀元前の日付には対応していないが、入力値はチェックしていないので危ない。

出力例：
input birth date : Year, Month, Day
1900 1 1
input today's date : Year, Month, Day
2000 1 1
The number of days from birth = 36524
Press any key to continue

１世紀で１年の平均日数が３６５．２４日になっているのでグレゴリオ暦としては正しいと
思うが、間違ってたらごめんｗ

>>4
日数計算のアルゴリズムはこんな感じ↓
5/10の10日先は、10+10=20なので、5/20。
5/10の20日先は、10+20=30なので、5/30。
5/10の30日先は、10+30=40なので、5/40、
しかし、5月は31日までなので、日は40-31=9、月は5+1=6。よって、6/9。
これを繰り返せば、任意の日付のx日後を求められる。
3/23 + 100 = 3/123 = 4/92 = 5/62 = 6/31 = 7/1 といった要領で。

閏年の判定は
LOGICAL FUNCTION CheckLeap(year)
とかいう関数を作って、
入力：西暦(整数型year)、出力：フラグ(trueなら閏年)
みたいにしておけばいいと思う。

やる気があれば、これくらいのヒントで十分作れるだろう。
あとは頑張れ。

プログラムの解説だ。
関数　qleapyear　はうるう年なら　.TRUE. を、普通の年なら .FALSE. を返す。
関数　idays_from_jan1　は、1月１日から入力日までの日数を返す。
　　（　1月１日を入れた場合は０日、1月２日を入れたら１日と数える。　）
主プログラムのアルゴリズム
最初の年は、３６５（または３６６）−誕生日までの日数。
最後の年は、今日までの日数。
途中の年は、うるう年を考慮しつつ３６５（３６６）を足してゆく。

料金の１００万円はマカオ銀行に金正日あてで振り込んでおいてくれ。

ありゃあ作っちゃってるよ。
投稿するの遅かったな。
でも、module文とか使っていいの？
Fortran77じゃ動かないけど…。

>>11
ごめｗ　
練習のため９０で書いてしまったｗ　
ともかく後半を任せた　もう寝るｗ

次の宿題もってこい！ｗ

誤差関数E(x)を次のように定義する。
E(x)=2/√π∫（範囲：０からｘ）e(-yの２乗)乗dy
この時、E(x)=2x(2)乗-x(3)乗
の実数解を「すべて」求めよ。（解はひとつとは限らないので解の個数についても議論せよ。）

このプログラムがどうしてもできません…fortranで作ってほしいですm(--)m
if文とループがうまく使えなくてエラーばっかりでてしまいます↓時間がないので助けてほしいですM(--)M×∞

>>4　後半。だるいので適当ｗ　　マイナス入れると暴走、入力チェックしてないｗ

MODULE mod_subs
IMPLICIT NONE
INTEGER, PARAMETER :: ndays_of_month(12, 2) = &
RESHAPE( (/31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31, &
31, 29, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31 /), (/12, 2/) )
CONTAINS
LOGICAL FUNCTION qleapyear(iyear)
IMPLICIT NONE
INTEGER, INTENT(IN) :: iyear
qleapyear = .FALSE.
IF ( MOD(iyear, 4) == 0 ) qleapyear = .TRUE.
IF ( MOD(iyear, 100) == 0 ) qleapyear = .FALSE.
IF ( MOD(iyear, 1000) == 0 ) qleapyear = .TRUE.
RETURN
END FUNCTION qleapyear
END MODULE mod_subs

PROGRAM AliceGame2
USE mod_subs
IMPLICIT NONE
INTEGER :: iyear0, imonth0, iday0, iyear, imonth, iday
INTEGER :: ileap, ndays, idays
!
WRITE(*, *) "input birth date : Year, Month, Day"
READ(*, *) iyear0, imonth0, iday0
!
WRITE(*, *) "input days "
READ(*, *) ndays
!
ndays = ndays + iday0
iyear = iyear0
DO
IF ( qleapyear(iyear) ) THEN
ileap = 2
idays = 366
ELSE
ileap = 1
idays = 365
END IF
IF ( ndays - idays > 0 ) THEN
iyear = iyear + 1
ndays = ndays - idays
ELSE
EXIT
END IF
END DO

!
imonth = imonth0
DO
IF ( ndays - ndays_of_month(imonth, ileap) > 0 ) THEN
imonth = imonth + 1
ndays = ndays - ndays_of_month(imonth, ileap)
ELSE
EXIT
END IF
END DO
iday = ndays
WRITE(*, *) "The date is ", iyear, imonth, iday
STOP
END PROGRAM AliceGame2

センブラインデント用
>>15-17

結果チェックしてないｗ　自分でチェックして書き直すですよ！

ごめｗ　バグってたｗｗｗｗ

IF ( ndays - ndays_of_month(imonth, ileap) > 0 ) THEN
　imonth = imonth + 1
　ndays = ndays - ndays_of_month(imonth, ileap)
ELSE
↓
IF ( ndays - ndays_of_month(imonth, ileap) > 0 ) THEN
　ndays = ndays - ndays_of_month(imonth, ileap)
　imonth = imonth + 1
ELSE

順番入れ替えてくれｗ

>>14
ｆ（ｘ）=ERF(x) - 2x**2 + x**3　
とおく。求めたいｘはｆ（ｘ）＝０を満たすｘ。

１階及び２階の微係数を求めると、
f'(x)=exp(-x**2) - 4x + 3x**2
f"(x)=-2x exp(-2x**2) - 4 + 6x

また、それぞれのｘ＝０とｘ→∞の値は、
f(0) = 0,　lim f(x) >0　
f'(0) ＞ 0,　lim f'(x) >0　
f"(0) < 0, lim f"(x) >0

ここでf"(x)の概形を調べてみると
f"(x) = 2x( 3 - exp(-2x**2) ) - 4 > 2x(3 - 1) - 4 = 4x - 4　〜　１次直線

これらのグラフを書けば、f'(x)はｘ〜１付近で最小値を持つ。
f'(1) < 0 なので、f'(0) ＞ 0,　lim f'(x) >0より　f'(x) は２回０点を切ることがわかる。
つまり、ｆ（ｘ）は２個極値を持つ。

さらにこのことと、f(0) = 0,　lim f(x) >0　また　f(1) < 1 - 2 + 1 = 0 より　fはx＞０で
２回０点を切ることが分かる。

ゆえに自明な解ｘ＝０も含めて３個解を持つ。
（その解はx=1付近に前後に１こづづあるはずｗｗ）

FORTRANの宿題っていうか授業してるところってまだあるんだね。
機械とか建築関係ならまだ多いのかな？

理学部はFORTRANしかやらねーな。もしくはFORTRANすらやらねーな。

>>14
とりあえず、誤差関数ERF(x)は持っていると仮定した。
表式が解析的に表されているので、微係数の式を求めた上で、ニュートン・ラフソン法でといた。

MODULE mod_sub
IMPLICIT NONE
CONTAINS
REAL FUNCTION f(x)
IMPLICIT NONE
REAL, INTENT(IN) :: x
REAL :: ERF
EXTERNAL :: ERF
f = ERF(x) - 2.0 * x**2 + x**3
END FUNCTION f
!
REAL FUNCTION fderiv(x)
IMPLICIT NONE
REAL, INTENT(IN) :: x
REAL :: pi
pi = 4.0 * ATAN(1.0)
fderiv = 2.0 / SQRT(pi) * EXP(-x**2) - 4.0 * x + 3.0 * x**2
END FUNCTION fderiv
END MODULE mod_sub
!

PROGRAM erf0
USE mod_sub
IMPLICIT NONE
INTEGER :: i
REAL :: x
WRITE(*, *) ' input start value x0 :'
READ(*, *) x
DO i = 1, 100
x = x - f(x) / fderiv(x) ! Newton-Raphson method
PRINT *, x, f(x) ! debug info
IF ( ABS(f(x)) < EPSILON(x) ) THEN
WRITE(*, *) ' Root found x= ', x, ' : f(x) =', f(x)
STOP
END IF
END DO
WRITE(*, *) ' Not converged in 100 times iteration. Try another start value.'
STOP
END PROGRAM erf0

実行例
input start value x0 :
0.1
-2.5092300E-02 -2.9582733E-02
-1.0398688E-03 -1.1755297E-03
-1.9108115E-06 -2.1561270E-06
-6.5146001E-12 -7.3509384E-12
Root found x= -6.5146001E-12 : f(x) = -7.3509384E-12
Press any key to continue

input start value x0 :
0.6
0.7874328 -1.7304393E-02
0.7620823 -9.1377224E-05
0.7619469 1.8604648E-08
Root found x= 0.7619469 : f(x) = 1.8604648E-08
Press any key to continue

input start value x0 :
1.3
4.483723 50.93219
3.281823 14.80570
2.510043 4.213046
2.034683 1.139575
1.769606 0.2661909
1.657068 3.9245199E-02
1.633732 1.5385470E-03
1.632740 2.6429243E-06
1.632738 7.6136075E-08
Root found x= 1.632738 : f(x) = 7.6136075E-08
Press any key to continue

解いてくれてありがとうございます！！これはfortranでプログラム組む事できるんですか？？

うわ！！ありがとうございます！！！

つまり、求める解は、
x = 0.0, x = 0.7619469, x = 1.632738
の３つ。

ニュートン法の収束条件判定はいい加減につくった。
自分で工夫してくれ。
単精度実数では、1.0e-7~8　位が数値誤差だ。




この位の宿題は自分でやるです！
間違ってても、謝罪も賠償もしないです！

>>26
FORTRAN90で書いたから、理解したうえでFORTRAN77に直してくれｗ

なお、ここではERF（ｘ）は、IMSLのライブラリ外部関数を使っている。
その辺の特殊性がちょっとある。

fderiv は　ｆの微分だ。

ｆ（ｘ）のグラフを書いて理解を深めてくれ。

もう寝るです。

ありがとうございましたｍｍ

ごめんｗ　ちょっと間違ってたです。

>>20
ERF(x)の微分から　２／√π　が抜けてるｗｗｗ
微妙に式が変わってくるが結論はあんまり変わらないから許せｗ

＜デフォルトの名無しさん＞＜５＞さんへ
閏年の問題ありがとうございます！！
参考になりました！！

>>32
おまえ、いろんなとこにマルチしまくっているな。

Yahooの方にもいやがった。
トップ > コンピュータとインターネット > プログラミング言語 > 全般 > ◆◆◆プログラミング質問箱◆◆◆

YahooID検索したら、ここにもいたｗ

トップ > コンピュータとインターネット > ソフトウェア > オペレーティングシステム > Windows > Windows XP > fortranわかる人いますか？

何故XPトピ？

すいません、ちょっとお力をお貸しください。

曲線y=xlog(a+x)と直線y=axとで囲まれた部分の面積の近似値を
連立方程式の a+3b=1 2a-b=2 を行列計算を用いて求めたa,bを使い答えを求めよ。

というプログラミングなのですがうまくできません。
行列計算においてdo文の使い方がいまいち…
もしよろしければどなたかご助言を。

>>33
向こうでの回答はFORTRAN６６風味だったです。
ついでにだれかFORTRAN77風に書いてみやがれです。

>>32もマルチポストする気力があるなら、自分で解きやがれです。
サブルーチン２個とメインルーチン１個、それぞれ核心部分は１０行程度なんだから、
たいしたこと無いだろです。よく学んでおくです。

>>35
そもそも、ｂが式に現れていない。問題を写し間違えていないか？
それとも単に、連立方程式を解いて出たａだけを使うという問題のなのか？

問題の背景とここまでやった自分の考えを述べなさい。

>>37
ええと、問題はこれで合っています。
bは使わないで、問題を解いて出てきたaのみで計算するようです。

行列計算ということなのでガウス消去法を用いて解こうとし
それぞれの座標に
a(1,1)=1
というように値を当てはめていったのですが、その後計算するにあたって
fortranのプログラムにおいてガウス消去法がどのようにやるのか分からず、今に至ります。

とりあえずフォートランに多相型っていつはいるんですか？

>>39
polymorphismならFORTRAN2003でそれなりに入っている。

>>35
まず、両方の曲線の交点が２つあって、その間の面積は解析的に求めることは簡単ですね。
で、その解析的に出た面積はaが分かれば数値的に出てくると。
そのaは二元連立方程式の解の一つであると。

二元連立方程式はプログラミングしないでも、中学レベルの頭ならすぐ解ける。

どこにプログラミングの必要性があるのか疑問。

>>38
いまいち問題設定が見えてこないなｗｗｗ　これが問題のすべてか？背景は無いのか？

ガウスの消去法は、例外処理を考えなければ、手でやるとおりに組めｗ
例外処理はめんどくさい。

２＊２限定なら、高校で習った公式に突っ込んだほうが早い。例外がdet=0で排除できるし。
というか手でといたほうが早いわけだが。

はっきり言って２＊２をガウスの消去法で解く意味が分からんｗ
というわけだから張り切って一般のｎ次行列について作れｗ

>>41
これをプログラミングで解いてレポートして提出しろという課題なんですよ。
一応実際に計算して解くことは出来ました。
しかしプログラミングのほうが分からず…
>>42
ガウスは使わない方がいいのでしょうか。
となると逆行列を用いての計算のほうがやりやすいのでしょうか？

>>43
じゃあさ、面積を求める部分をSimpson法やいろんな数値積分法で組んでしまうのはどう？
面積を求める部分をサブルーチンか関数にしてaを渡せばそれなりにプログラミングの形になるし。

>>43
ガウスの消去法はアルゴリズムが単純なので、小さめの連立方程式を解くための
練習問題としてはよく使われる。行列がsingularな時や誤差が溜まりやすい状況を、
あまり考慮しなければ、比較的楽に作れるだろう。

素朴には、対角要素が全部０で無いなら、１列目からｎ列目まで、順に各列の対角要素以外を
０にするように行ごとひいていけばいい。手でとくのと同じだ。
対角要素が０や０付近だと例外処理をすることになる。このへんがまんどくさい。

この辺は、問題製作者の意図がどの辺にあるのかで、こだわり度が決まってくるのだが、
この問題設定では、いまいち理解しにくい。

２＊２での一般形なら単純に公式を使うほうが早いが、n*nの一般サブルーチンを求めているなら、
ガウス法なり何らかの方法が必要になる。

後半では数値積分をすることになるわけだが、問題として全く意図が読めないので、
これだけの情報では、どの方法をとるべきか答えがたい。

>>44
うーん、すいません、またｔくの初心者ですのでよくわかりません…
助言は非常に嬉しいのですが…すいません。
>>45
今回はどうやら2*2行列のようです。
それでは今回は公式の方を用いて計算することにします。
まぁ、普通の行列のプログラミングもよくわからないのですが…
後半部分はあれだけです。
タダ単に前半で求めた値を代入して面積を求めるだけでして、背景は無いと思います。

>>35
まあ、適当にやれｗ

a = 1.0, b = 0.0

交点
x = 0.0, x = e - 1

積分　ｆ（ｘ） = ax - x log(a + x) (a = 1.0　積分範囲　[0.0, e-1])

解析的な値
-e^2 / 4 + e - 5 / 4

台形公式　短精度　分割数 n=1000 で　こんな感じだ。　
誤差は( (e-1)/1000 )^2〜1e-6位なので、まぁもっともな気がする。

numerical integral 0.378981858
analytic value　　　0.378982186
Program Completed
Press Enter to Continue.

PROGRAM vip
IMPLICIT NONE
INTEGER, PARAMETER :: n = 1000
REAL :: x0, x1, s
x0 = 0.0
x1 = EXP(1.0) - 1.0
s = trapez_integ(n, x0, x1)
PRINT *, 'numerical integral', s
PRINT *, 'analytical value ', 5.0 / 4.0 + EXP(1.0)**2 / 4.0 - EXP(1.0)
STOP
CONTAINS
REAL FUNCTION trapez_integ(n, x0, x1)
INTEGER, INTENT(IN) :: n
REAL , INTENT(IN) :: x0, x1
INTEGER :: i
REAL :: s, h
h = (x1 - x0) / REAL(n)
s = 0.0
DO i = 1, n
s = s + 0.5 * ( func(x0 + (i - 1) * h) + func(x0 + i * h) )
END DO
trapez_integ = s * h
RETURN
END FUNCTION trapez_integ
!
REAL FUNCTION func(x)
REAL, INTENT(IN) :: x
REAL, PARAMETER :: a = 1.0
func = a * x - x * LOG(a + x)
RETURN
END FUNCTION func
END PROGRAM vip

訂正

解析的な値
-e^2 / 4 + e - 5 / 4
↓
e^2 / 4 - e + 5 / 4

正負逆だったｗ

交点が　0.0 と　e-1　なのは　func(0.0) と func(e-1) を書かせて確かめるのもよし。

前半は、２＊２の逆行列の公式でいいだろう。
Ax＝ｂ　
ｘ＝（A^-1）ｂ
で一応行列の積も入るし。

ひょっとして、数値計算プログラム作るなら
Cよりフォートランの方が高速コードを生成しますか？

>>50
昔はそうだった。最近はそれほどで差が無い。
コンパイラ自体がフロントエンド部とバックエンド部に分けて作られるようになり、
中間言語→機械語生成のバックエンド部は共通になったから。

しかしF90以降ならどう考えてもFORTRANの方がバグが少なく、
簡潔明快なソースでやりたいことを表現できる。

速度は？

スパコンのレベルになると、事実上FORTRANしか選択肢がなかったりする

でもそこらのパソコンでやるなら、Cで書こうがFORTRANで書こうが大して変わらんね。
g77とかはほとんど内部でf2cを動かして、C言語に変換してからgccでコンパイルしてるし

>>52
コンパイラメーカーが同じなら、速度もたいして変わらない。
要するに、ソースを中間言語に直すところでの最適化の違いだけ。
中間言語から機械語へは共通化している。
ただ、FORTRANの方が制約がきつい分コンパイラーの判断できる最適化は
Cより高いとされている。

>>53
スパコンに関して言えば、９０年代はじめごろまではベクトル化コンパイラーはFORTRAN
だけで、Cコンパイラーがあっても非ベクトル型だった。しかし、その後Cコンパイラーも
ベクトル化コードを吐くようになってきた。

FORTRAN９０／９５はコードを書いていて気持ちいい。
柔軟性と制約のバランスが取れていてぴったりくる手袋や靴下のよう。

イラつくのは宣言部が重くなりがちな点くらいか。

みなにお勧めする。

じゃあ、特に速度を要求する部分を、書けるならアセンブリで、書けないならフォートランで、
んで普段はC++で書くとかって感じですかね

MPIのルーチンはCとFortranしか提供されていないが、微妙にCより。

C++は動的バインディングのせいでCにくらべて１０〜２０％遅くなるという話だ。
次期FORTRANにOOP機能が付け加わったらどうなることか。

IMPLICIT REAL*8 (A-H,O-Z)
WRITE(*,*) 'あ'
N=3
COMPLEX*16 A
END


でコンパイルしようとして

$ g77 test3.f
test3.f: In program `MAIN__':
test3.f:2:
WRITE(*,*) 'あ'
1
test3.f:4: (continued):
COMPLEX*16 A
2
Statement at (2) invalid in context established by statement at (1)

って言われるんですけど、何がいけないんでしょう？
教えてください

>>58
COMLEX文は宣言文だから、最初の実行文よりも前に置かれねばならない。
このばあい、WRITE文の前。

それで怒られている。

2006年の1月から3月までのカレンダーを作れという課題がでました。

・・・・・・・・・・・
do 20 t=1,3
if(t.eq.2)then
b=28
else
b=31
end if
・・・・・・・・・
・・・・・・・・・
write(6,*)" ",t,"月　　　　　　　　"
write(6,*)"　日　月　火　水　木　金　土　　"
・・・・・・・・・・
・・・・・・・・・・
write(6,100)(a(i),i=6,b)
100 format(7i5)
・・・・・・・・・・
・・・・・・・・・・

・・・・の部分がわかりません。よろしくお願いします。

いまどきフォートラン使ってる人って、どういう立場だなんだ。

趣味？　仕事？　学校？　環境は？　作ってるアプリの種類は？

ちと興味ある。

>>61
失敬な！
FORTRANは現役ですよ！いまやモジュール型言語として完成。次はOOPｗ
INTELが直接開発・販売しているコンパイラはC＋＋とFORTRANだけですよん。

いやだから何作ってるのか教えてよ。現役なのは存じてるから

FORTRANを科学技術計算以外で使う人なんているの？

ということはここに居る人みんな科学者！？

>>61, >>63
俺は固体の量子計算。学校の研究です。
アプリっていうとGUIを想像してしまうけど、そうじゃなくて、
コマンドラインで実行するようなプログラム。

>>65
たぶんそうだと思う。
あとガッコウの宿題ができない人々…。

Fortran使って構造計算してる土建屋も科学者か？

FORTRAN　ｍｐ３　ｐｌａｙｅｒ
http://www.pbase.com/gadget_square/fortran

ﾜﾛｽ

fortran mp3 encoder
ttp://members.at.infoseek.co.jp/kitaurawa/cgi-bin/wiki.cgi?page=UZURA

ところで、fortranに高階関数はいつ入るんですか？

>>70
高階かんすうってなんですか？おしえてかださい

　　　　　　　　　　　 ,,.　 -‐―-　､
　　　　　　　　　 ／ _..-‐―- 、　　.ヽ
　　　　　　　　 / ／-‐―-..__..ヽ　_..、ヽ
　　　　　　　　.!三三三´　　　 ｀´　　:　i
　　　　　　　　|.!三三　　 __　　￣ _,. .;. i　暗黙の型宣言を守らない
　　　　　　　　|ii　i!"　　´ ｪｪ｀　l´rｪ､ i..!　　　　　　>>60は地獄へ逝く
　　　　　　　　 !ｉ､ i!:::..... -￣r　　ヽ　　.!
　　　　　　　　　!ｉｉ!.|:::::::　　/ ｀−´ヽ　/
　 _,,..　===＝"/%;;;::::.　　　'∠ﾆﾆ」　 /、　　　　　　　　　　　　_　-‐‐-､
/´　　　　　　 | .%　:::;;;;:　,,,.,,,:. ⌒: ./´　￣＝=-- 　_,,.. -‐‐'"´　_,..-‐ ´
　　　　　　　 |　 %　　::::::;;;;;;;;;:.""%.／　　　　　　ｒ"　(<二￣￣>
　　　　　　　.|　　%｡ﾟ+。　　。+% ／　　　　　　　｀ヽiﾉ<二.￣￣>
　　　　　　　|　　　%+*ﾟ+。*ﾟ%／　　　　　　　　　（　<二￣￣>

>>60　特別にFORTRAN77で書いてやったです。　感謝しやがれです。
PROGRAM calndr
C calender for Jan to Mar 2006
INTEGER ia(31)
CHARACTER*30 fmt
mon1 = 1
C Jan 1 2006 is Sunday (= 1)
DO 10 i = 1, 31
ia(i) = i
10 CONTINUE
DO 20 it = 1, 3
IF (it .eq. 2) THEN
ib=28
ELSE
ib=31
END IF
WRITE(fmt, '(a, i2, a, i1, a)')
1 '(', 4 * (mon1 - 1), 'x, ', (8 - mon1), 'i4, / (7i4))'
mon1 = MOD(ib + mon1 - 1, 7) + 1
WRITE(6, *) ' ', it, '月　　　　　　　　'
WRITE(6, *) '　日　月　火　水　木　金　土　　'
WRITE(6, fmt) (ia(i), i = 1, ib)
WRITE(6, *)
20 CONTINUE
STOP
END

>>73

実行結果です。

1 月　　　　　　
　日　月　火　水　木　金　土
1 2 3 4 5 6 7
8 9 10 11 12 13 14
15 16 17 18 19 20 21
22 23 24 25 26 27 28
29 30 31

2 月　　　　　　
　日　月　火　水　木　金　土
1 2 3 4
5 6 7 8 9 10 11
12 13 14 15 16 17 18
19 20 21 22 23 24 25
26 27 28

3 月　　　　　　
　日　月　火　水　木　金　土
1 2 3 4
5 6 7 8 9 10 11
12 13 14 15 16 17 18
19 20 21 22 23 24 25
26 27 28 29 30 31

Press any key to continue

画面が崩れたので、専ブラ用のアンカーを書いてやるです。　
>>73,75

FORMATの技巧が高度すぎて、>>60はもちろん>>60に出題した
どんくさいチビ人間には理解できないだろです！

専ブラ用のアンカー間違ったです。　
>>73,74

>>71 関数を値として引数に取れる関数だよ。

関数のポインタ？

>>77
単にexternalなら今でも出来るが、関数ポインターならFORTRAN2003で可能に。

今日の宿題はまだか？ｗ

ところで型推論っていつFortranにはいるんですか？

>>81
型推論ってなんですか？

すみません、質問なのですが。
フォートランで作った実行ファイルとソースファイルから
そのソースでcallから呼び出されてる関数
のライブラリ(モジュール？）の場所を探す方法はありますか？

というのも、プログラムを使っていたＵＮＩＸマシンから別のＵＮＩＸマシンにプログラムを移したところ
マルチバイト文字が正しくありません。とエラーが出てしまいます。
なのでソースがあるので別のマシンでコンパイル、リンクしようと思ったのです。
そうするとソースの中にcallで呼び出されたサブルーチン記述されていない名前の関数が
あったのでリンクのときにこのモジュールを指定しなければいけないんですよね？
その場所がわからずコンパイルでエラーが出てしまいます。

長文で申し訳ないのですがよろしくお願いします。

>>83
エラーメッセージを見ていないので状況がいまいちつかめないが、
MAP出力とかのコンパイラ（リンカ、ローダー）オプションで、
必要な情報が得られるかもしれない。

>>84
レスありがとうです。

いろいろ悩みましたが
ｺﾏﾝﾄﾞ　nm を使ってライブラリを全部grepに渡して関数名を検索したら出てきました。
新しいマシンに移してコンパイルリンクでき、実行も出来ました。

今日は宿題無いのか？

日付を2つ入力して，その間に何日あるかを書き出すプログラム(Nannichi)を作れ。実行例はつぎの通りである． 日付を２つ入力して下さい．
何年？ 2004
何月？ 7
何日？ 1
何年？ 2005
何月？ 7
何日？ 1
日数は 365日です．
これの答えを教えてください　お願いします

>>87
>>5-6をみろ。すでにより一般的なうるう年を考慮したサブルーチンが存在している。
それをちょっといじればよい。

質問があれば、日本語でおｋｗｗｗｗｗ

>>82 型を宣言しなくても勝手に正しく推論してくれるしくみだよ。

いつはいるの？

>>89
率直に言ってそれは下らないアイデアだと思う。
書きなぐるスクリプト言語ではよくあるが、なんんというかだらしなさのあらわれだ。

どっちにしろ暗黙の型宣言があるし。

ROZEN　MAIDEN ・・・　orz

ところで、Fortranにproof-carring-codeはいつ入るの？

>>92
proof-carring-codeってなんでちゅか？ばぶー

だれかこのプログラミングを教えてください。お願いします。
ｘ=0.1 , 0.2 , 0.3 , 0.4 , 0.5
ｙ=1.228 , 1.005 , 0.823 , 0.674 , 0.552
近似式　ｙ=a*exp(bx)
上の関数とデータで最小自乗法で近似した式を求めよ。
この問題のプログラムがわかりません・・・。お願いします。

>>94
PROGRAM vip
IMPLICIT NONE
INTEGER, PARAMETER :: n = 5
REAL, PARAMETER :: x0(n) = (/0.1 , 0.2 , 0.3 , 0.4 , 0.5 /), &
y0(n) = (/1.228, 1.005, 0.823, 0.674, 0.552/)
REAL :: y(n), t(n, 2), z(2, 2), zinv(2, 2), v(2)
INTEGER :: i
y = LOG(y0)
t(:, 1) = 1.0
t(:, 2) = x0
z = MATMUL(TRANSPOSE(t), t)
v = MATMUL(TRANSPOSE(t), y)
zinv = ainverse(z)
v = MATMUL(zinv, v)
WRITE(*, *) ' least square fit a=', EXP(v(1)), ' b=', v(2)
STOP
CONTAINS
FUNCTION ainverse(x)
IMPLICIT NONE
REAL, INTENT(IN) :: x(2, 2)
REAL :: ainverse(2, 2)
REAL :: determinant
determinant = x(1, 1) * x(2, 2) - x(1, 2) * x(2, 1)
IF (ABS(determinant) < 1.0e-7) STOP ' singular matrix '
ainverse(1, 1) = x(2, 2) / determinant
ainverse(2, 2) = x(1, 1) / determinant
ainverse(1, 2) = -x(1, 2) / determinant
ainverse(2, 1) = -x(2, 1) / determinant
RETURN
END FUNCTION ainverse
END PROGRAM vip

>>95
least square fit a= 1.499271 b= -1.998701
1.227659 1.005252 0.8231379 0.6740159 0.5519093
1.228000 1.005000 0.8230000 0.6740000 0.5520000

>>95さん　>>96さん
ほんとありがとうございます。
月曜にでも学校でやってみます。
また、わからないことがあれば質問させてください。
ありがとうございました。

>>95について。
y=a*exp(bx)を次の式に変形し
log(y)=log(a)+bx
yとxの各値を上の式に代入すると
次の5本の式となる。
log(y1)=log(a)*1+b*x1
log(y2)=log(a)*1+b*x2
log(y3)=log(a)*1+b*x3
log(y4)=log(a)*1+b*x4
log(y5)=log(a)*1+b*x5
ここで、
log(yi)がi番目の要素である列ベクトルをy、
i行1列を1、i行2列をxiとする行列をt、
1番目の要素をlog(a)、2番目の要素をbとする
列ベクトルをvとすると、上の5式は、
y=tv
と表せる。
tが1次独立な列から成るので、
vの最小二乗推定は、次の式
(t't)**(-1)t'y
となる定理を>>95は使用していると思います。
(注)t'はtの転置を表す。(t't)**(-1)はt'tの逆行列を表す。）

ところでＦｏｒｔｒａｎに配列のバウンダリチェックはいつ入るのでしょうか？

>>99
subscript check は、debug 用にoptionで入るのが普通。1960年代からある。
しかし、遅くなるから実際の実行時に使う奴はいない。
この辺はスピード優先の数値計算の価値観を反映している。

オッス！おれ極右！
>>98氏解説�d

最近寒いので湯たんぽ代わりにノートパソコンを腹の上において寝ているお。
昨夜そういう寝そべった状態で布団の中で作ったのが>>95だお。
オラはVipperだから間違っていても謝罪も賠償もしないお！

最小二乗法なんて普段使わないから、きょうの昼おっきしてからNumericalRecipesで
式を確かめといたおｗ　方式はちょっと違うが式おｋｗｗｗ　（15.2.4~6)あたり。
ttp://www.library.cornell.edu/nr/bookfpdf/f15-2.pdf

FORTRAN77で書くと３倍くらいの長さになるとおもうお。
MATMULやTRANSPOSEはF90にしかない行列積と転置だお。
LOGを取ったのは線形化するためだお。
普通、最小二乗法は線形の式にしか使わないと思ったから、そうしてみたおｗ

NumericalRecipesは、旧版のC言語版なら訳本が出ているお。
ttp://www.amazon.co.jp/exec/obidos/ASIN/4874085601/qid%3D1139049291/249-3009705-2437946
本のコードにはバグがちりばめられているので、本を打ち込むのはお勧めしないお。
金出して買うコードのほうはこまごまと修正されているお。

（＾ω＾）このくらいソースコードが短いと勉強する気にもなるだろ？

>>87
フォートランはわかんないからベーシック
あってるかどうかもわかんない
こんにちは

10 gosub 60
20 nannichi0=nannichi
30 gosub 60
40 print "日数は";nannichi-nannichi0
50 end
60 input "年,月,日",year,month,day
70 if month<3 then year=year-1:month=month+12
80 nannichi=year+int(year/4)-int(year/100)+int(year/400)+int((153*month-457)/5)+day
90 return

>>102
解説頼む

>>103
まちがいた。ちょと訂正。

80 nannichi=365*year+int(year/4)-int(year/100)+int(year/400)+int(306*(month+1)/10)+day-428

「これは Fairfield公式 といって、西暦1年1月1日からの経過日数を算出する式らしい。
　ただしグレゴリオ暦が採用されたのは1582年10月15日からだが。
　要するにこれを使って2つの日付の差を求めればよい。」
っておばあちゃんがいってた。

すいませんがこないだのプログラムを
ＤＯ文で作っていただけませんか。
何回もお願いして、すいません。

コンピューターおばあちゃん♪
コンピューターおばあちゃん♪
イエーイ　イエーイ　僕は大好きさ〜♪

>>104
把握したｗ
２月の２８日の例外を２月をお尻にすることで回避。
大の月と小の月の入れ替えを１ヶ月の平均日数を30.6日にすることで実現。
かなり絶妙な公式だ。ｵﾓｼﾛｽｗ


だが最近流行の明示的なプログラミングから離れているがｗｗｗｗｗ

>>105　世話が焼けるッス！　FORT７７で、素朴に書いといたからｗ
PROGRAM vipper
　　　PARAMETER (n = 5)
REAL x(n), y(n)
　　　DATA x /0.1 , 0.2 , 0.3 , 0.4 , 0.5 /
DATA y /1.228, 1.005, 0.823, 0.674, 0.552/
　　　sumx = 0.0
　　　sumy = 0.0
　　　sumxx = 0.0
　　　sumxy = 0.0
sum = REAL(n)
DO 10 i = 1, n
　　　 sumx = sumx + x(i)
sumy = sumy + LOG(y(i))
sumxx = sumxx + x(i) * x(i)
sumxy = sumxy + x(i) * LOG(y(i))
10 CONTINUE
det = sum * sumxx - sumx * sumx
alog = ( sumxx * sumy - sumx * sumxy ) / det
b = ( sum * sumxy - sumx * sumy ) / det
WRITE(6, *) ' a = ', EXP(alog), ' b = ', b
　　　END

>>108　インデント

PROGRAM vipper
PARAMETER (n = 5)
REAL x(n), y(n)
DATA x /0.1 , 0.2 , 0.3 , 0.4 , 0.5 /
DATA y /1.228, 1.005, 0.823, 0.674, 0.552/
sumx = 0.0
sumy = 0.0
sumxx = 0.0
sumxy = 0.0
sum = REAL(n)
DO 10 i = 1, n
sumx = sumx + x(i)
sumy = sumy + LOG(y(i))
sumxx = sumxx + x(i) * x(i)
sumxy = sumxy + x(i) * LOG(y(i))
10 CONTINUE
det = sum * sumxx - sumx * sumx
alog = ( sumxx * sumy - sumx * sumxy ) / det
b = ( sum * sumxy - sumx * sumy ) / det
WRITE(6, *) ' a = ', EXP(alog), ' b = ', b
C check
WRITE(6, *) ' i x y',
1 ' a * EXP(bx) (dy)^2'
DO 20 i = 1, n
yy = EXP(alog + b * x(i))
WRITE(6, *) i, x(i), y(i), yy, (y(i) - yy)**2
20 CONTINUE
STOP
END

>>109　インデント

a = 1.499271 b = -1.998701
i x y a * EXP(bx) (dy)^2
1 0.1000000 1.228000 1.227659 1.1664589E-07
2 0.2000000 1.005000 1.005252 6.3688631E-08
3 0.3000000 0.8230000 0.8231379 1.9023346E-08
4 0.4000000 0.6740000 0.6740159 2.5326941E-10
5 0.5000000 0.5520000 0.5519093 8.2298044E-09
Press any key to continue
>>110　インデント

>>94
死ぬほど易しく作ったから理解しろｗ
>>101の本の公式通りに作り、変数名も大体あわせておいた。
ここで誤差のσは各要素で共通一定と仮定した。この場合キャンセルしあうから式から無視してよし。

本を見れば分かるが、二乗誤差が極値をとる条件を、a,bに関してそれぞれ微分を取ることで求めている。
二乗誤差は極値を１個だけもちかつそれは最小だからこれでおｋとなる。
（誤差は外せばいくらでも大きく出来るが、合わせるほうはに限界がある。最高に合いまくって０。）

プログラム解説：
１０番LOOPで総和を取る。
その後、２＊２の逆行列をかけるのと同じ理屈で、２元連立方程式を解いている。
２０番LOOPでどのくらい近似できたか見ている。

求めていないが、一番右側の列の二乗差の総和が最小になっているはずである。

a,bの値を手動で変えて、求めた二乗差の総和が最小であることを確かめてみそ。
その結果を考察に書いておけば＋１０点の花丸ハンバーグだ！

みなさんほんと詳しく教えていただきありがとうございます。
ほんと助かりました。

ところで、Fortranにオブジェクトデータ構造はいつ入るんですか．

>>113
FORTRAN2003。コンパイラーの実装待ち。

>>113
オブジェクトデータ構造ってなんですか？
オブイェークト

ERK=ER(K)*DCOS(GANMA2(K))-EI(K)*DSIN(GANMA2(K))
EIK=EI(K)*DCOS(GANMA2(K))+ER(K)*DSIN(GANMA2(K))

実数 ERK=ERK*DCOS(GANMA3(K))-EIK*DSIN(GANMA3(K))
虚数 EIK=EIK*DCOS(GANMA3(K))+ERK*DSIN(GANMA3(K))

ERK=ERK*DCOS(GANMA4(K))-EIK*DSIN(GANMA4(K))
EIK=EIK*DCOS(GANMA4(K))+ERK*DSIN(GANMA4(K))

この式をEXPを用いた複素数型に変換できなくこまってます。
どなたか優しい方おしえていただけませんか（＞＜）

>>116
やぁ！やさしいお兄さんだよ！援助してあげるから、今日からパパと呼ぶんだよ！
定義：　
ze = er + i ei
zg = COS(Γ) + i SIN(Γ) = EXP(iΓ)

積：
ze * zg = ( er + i ei ) * ( COS(Γ) + i SIN(Γ) )
　　　　　= ( er * COS(Γ) - ei * SIN(Γ) ) + i ( ei * COS(Γ) + er * SIN(Γ) )

これより
COMPLEX*8 :: ze, zek
ze = (er, ei)
zek = ze * EXP( CMPLX(0.0d0, gamma2) )

これでＯＫさ！

ただ細かいことを言うと、倍精度複素数はＦ７７標準規格に存在しないので、
完全に規格に乗っ取ろうとすると複素変数は使えないのさ！
でも、大概のコンパイラは倍精度複素数があるのでＯＫさ。
ＥＸＰとＣＭＰＬＸは総称名なので、精度は引数に合わせてくれるはず。

でもマニュアルを調べて総称名が、ちゃんと倍精度複素数になってるかどうか
確認したほうがいいかもね。

おっと
ze = CMPLX(er, ei)
だったね。

さぁ、>>116　おっぱいうｐうｐ
　　　 _ 　∩
　　(　ﾟ∀ﾟ)彡　おっぱい!おっぱい!
　　(　 ⊂彡
　 　|　　　|　
　 　し ⌒Ｊ

なんだかしりませんが
Fortran2003って普通にかな〜りすごそうなんですけど
ほんとに実装されるんですか？

>>118
つ ttp://kahouweb.hp.infoseek.co.jp/foods/food_230.html

>>119
FORTRANの場合スーパーコンピュータやCPUのキラーアプリとして
単独での採算を度外視されて開発されているので、多分実装されるでしょう。

INTELが旧DECのFORTRANコンパイラーチームをCompaqから買い取ったのは
そういう方向からだし。

次の宿題もってこい！

ｆ９０を使い初めてプログラムの複素数化をしているのですがうまくいきません、どなたか助言をおねがいします。
PROGRAM TTT
IMPLICIT NONE
INTEGER :: I,N,M
REAL(8) :: EI,ER(0:2047),EI(0:2047),A,B
REAL(8) :: FACTR,H,PI,RN1,RN2,EP0,C,WL,TINTY,TH,S(0:2047),TT(0:2047),GEFF,ERI,EII
!
PI=4.0*DATAN(DBLE(1.0))
C=3.0E+8
WL=1.55E-6
DO I=0,N-1
FACTR=-H*PI*RN2*RN1*EP0*C/WL
TINTY=(ABS(S(I)*S(I))+ABS(TT(I)*TT(I)))/DBLE(2.0)
TH=FACTR*TINTY
ERI=EXP(-GEFF*H)*(1*DCOS(TH)-1*DSIN(TH))
EII=EXP(-GEFF*H)*(1*DCOS(TH)+1*DSIN(TH))
ER(I)=ERI
EI(I)=EII
OPEN(20,FILE='TH.dat')
WRITE(20,*)ER(I),EI(I)
END DO
END PROGRAM TTT
上のプログラムからSINやCOSをけしてEXP等で表現
し複素数でのプログラムにするにはどうしたらいいでしょうか？

>>123
まずその前に・・・・・
C=3.0E+8
WL=1.55E-6
これだと単制度で定数が入るので、ここで数値が腐ってしまう。E→Dにせよ。

吸収係数かなんかの計算か？
発想が逆さまだ。もともとの公式を素直にプログラムするほうがよっぽど楽だぞｗ

ほいで
(1 + i) * EXP(iθ) = (1 + i) * (cosθ+ isinθ)
　　　　　　　　　　= (cosθ+ isinθ) + (icosθ - sinθ)
　　　　　　　　　　= (cosθ - sinθ) + i(cosθ +sinθ)
ゆえに
ZE = EXP(-geff * h) * (1.0d0, 1.0d0) * EXP(CMPLX(0.0d0, th))
こうじゃねーか?

間違っていてもシラネｗ

>>124
サンクスです、再度考えてみます

補足
(1+i) = SQRT(i) * SQRT(2) = EXP( iπ / 2 )

(i + i) * EXP(iθ) = EXP( i(θ + π/2) )

ごめｗ　今借り物でFEPが違くて勝手が違うｗ　間違い膜リング
(1+i) = SQRT(i) * SQRT(2) = EXP( iπ / 2 )　* SQRT(2)
(1 + i) * EXP(iθ) = SQRT(2) * EXP( i(θ + π/2) )

つーかやっぱ、一度元の公式に戻れｗ

糞ごめｗ　π/2 →　π/4 うっはーｗｗｗｗｗ　信用度０ｗｗｗｗ

ZE = EXP(-geff * h) * (1.0d0, 1.0d0) * EXP(CMPLX(0.0d0, th))
　 = EXP(CMPLX(-geff * h, th + pi / 4.0d0)
こうかな？ｗ

スレ汚し　すまんこﾉｼ！

>>124
できましたありがとうございますm(__)m

さぁ、>>123　おっぱいうｐうｐ
　　　 _ 　∩
　　(　ﾟ∀ﾟ)彡　おっぱい!おっぱい!
　　(　 ⊂彡
　 　|　　　|　
　 　し ⌒Ｊ

>>131
つ ttp://members.at.infoseek.co.jp/rimssecret/oppai2.htm

実数部分 ERK=ER(K)*DCOS(GANMA2(K))-EI(K)*DSIN(GANMA2(K))
虚数部分 EIK=EI(K)*DCOS(GANMA2(K))+ER(K)*DSIN(GANMA2(K))

ERK=ERK*DCOS(GANMA3(K))-EIK*DSIN(GANMA3(K))
EIK=EIK*DCOS(GANMA3(K))+ERK*DSIN(GANMA3(K))

ERK=ERK*DCOS(GANMA4(K))-EIK*DSIN(GANMA4(K))
EIK=EIK*DCOS(GANMA4(K))+ERK*DSIN(GANMA4(K))
を同じ考えで
EK=E(K)*EXP(CMPLX(0.0d0,(GANMA2(K)+GANMA3(K)+GANMA4(K))))
としたら同結果にならないのですがご教授おねがいします
本当に連投すみませんm(__)m

>>133
EK = E(K) * (1.0d0, 1.0d0)**3 * EXP( CMPLX(0.0d0, GAMMA2 + GAMMA3 + GAMMA4) )
= E(K) * SQRT(8.0d0) * EXP(CMPLX(0.0d0, GAMMA2 + GAMMA3 + GAMMA4 + 3.0d0 * pi / 4.0d0) )

でね？

>>133はもう一度>>124を読み直しなさい！！ｗ

誰かネタを投下しろです！

前半の距離h/2にわたる線形空間伝搬
[exp( j (r2 + r3 + r4 ) h/2 }F{Φ(z,t)}]
　
元の式がこれなのですが複素数でやろうとするとうまくいかないです・・・
(1.0d0, 1.0d0)の部分に変数をいれればいいのですかね？

>>137
その式のままプログラムすればいい。
phi = phi_func(z, t)
zf = F(phi)
zres = EXP( CMPLX（0.0d0, 0.5d0 * h * (r2 + r3 + r4) ) ) * zf

3行に分けたのはdebug時に、１行ごとにエラーを知るため。
うまく行くのがはっきりしたら、１行にまとめればいい。

>>138
コンパイラによってはｆｕｎｃｔｉｏｎを別のfunctionやサブルーチンの引数に取れないかもな。

1 名前：ぶつわよ！φ ★[] 投稿日：2006/02/08(水) 18:26:54 ID:???0
毎日新聞「開かれた新聞」委員会の１２、１月度見解を報告します。
今回は、昨年１２月２５日に山形県庄内町で起きたＪＲ羽越線の特急脱線転覆事故を扱った
毎日新聞の社説「安全管理で浮ついてないか」（１２月２７日朝刊）が
「風の息づかいを感じていれば、事前に気配を感じていたはずだ」と指摘したことなどに対し、
読者から多数の批判意見が寄せられた問題を中心に取り上げます。

※とっても長いので全文はソースで
http://www.mainichi-msn.co.jp/shakai/wadai/news/20060207ddm012070089000c.html

プログラミングの宿題です。助けて下さい。

2種の生物x(t)とy(t)が同一場所に生息している。
x(t)の出生率はa>0であり、x(t)はy(t)に食べられる、その死亡率はby(t)<0であるとする。
一方x(t)を食べるy(t)の死亡率はc<0であり、食料となるx(t)に比例して出生率がdx(t)>0であるとする。
このとき次の連立微分方程式系のシミュレーションのプログラムを作れ。
dx/dt=(a+by)x
dy/dt=(c+dx)y
ただし、t=0、x=3、y=1、h=0.2、a=1、b=-1、c=-1、d=1(hは時間tの刻み幅)
　
Fortran77の形式でお願いします。よろしくお願いします。

>>141 ほれ！ティンティンうｐねｗ
　微分方程式の積分を最高簡単にしておいた。誤差が溜まって結果怪しい。
　ルンゲ＝クッタとか真面目なのに自分で治しとけｗ

PROGRAM oppai
PARAMETER (h = 0.2, a = 1.0, b = -1.0, c = -1.0, d = 1.0)
t = 0.0
x = 3.0
y = 1.0
DO 10 i = 1, 100
t = i * h
dx = (a + b * y) * x * h
dy = (c + d * x) * y * h
x = MAX(0.0, x + dx)
y = MAX(0.0, y + dy)
WRITE(6, *) x, y, t, i
10 CONTINUE
END

>>142　インデント

結果はというと、ｔ＝１０くらいで、生命体ｘが滅亡し、その後餌のなくなった
生命体ｙもすぐに滅亡したｗ

>>141
メッシュを細かくして計算したら、ｘとｙは交互に増えたり減ったりして振動したｗ
この結果がもっともなので、>>142のままではやはりﾔｳﾞｧｽｗｗｗｗｗｗｗｗｗｗ

fortran(パソコン？）に初めて触ると言うくらいの初心者です。
http://www.coastal-env.k.u-tokyo.ac.jp/koibuchi/fortran95/fortran77.htmにそって、".Net"が無い為、fortran77をダウンロードしました。
書いてあるとおり、"F77.exe"の位置を指定して、プログラムをコンパイル(実行）を行った結果、画面には…


「コンパイルできません。
ファイル" …(ファイル名).exe"は、存在しません。」

と表示され、何も反応してくれません。
どうしたらよいでしょうか？教えてください…。

ちなみに、f77.exeは、"C:\Program Files\Fortran77\ftn77.exe"に保存してあります。
保存先は、"…\My Documents\Fortran"に指定しました。

F95でも、同様な事をすれば"fcpad231.lzh"でうまく行くと書いてありますが、こちらも教えて頂けましたら、幸い名のですが…。
(今、自分は研究者なのですが…遅まきながらの数値計算の導入で、かなり四苦八苦しています。よろしくお願いします。）

>>144
それは
「あなたの作ったプログラムに間違いがあるのでコンパイルできませんでした。
（コンパイルできなかったので）実行ファイルが存在せず、実行できません」
という意味です。
「メッセージ」の欄にエラーの起こった行とその理由が表示されるので
それを参考にもう一度よくプログラムを見直して訂正して下さい。

>>145
FTN95 コンパイル時のオプション指定の方法が FTN77 と違うため、
そのままでは CPad for Salford FTN77 では使えません。
FTN95付属のIDEを使うか（但し、日本語は使えません）、
適当なエディタを使って下さい。

>>146
ありがとうございます。そうだったんですね。

メッセージには、
■C:\Documents and Settings…\デスクトップ> ftn77 /link ex.for
[CPad: コンパイル終了]
とあります。

どうしたらよいのでしょうか…？

あと、適当なエディタを用いるとありますが、どのようなものを使用すればよいのでしょうか？
本当に初歩的な質問ですいません…。

もう一つ質問させてください。
f77よりf90でパンコイルしたほうが高速ですか？

>>148
初心者レベルで問題になるような著しい差はない。

エディターに関しては、初心者だというなら

何がしたいのか分からないが、新たに言語を習得したいというならF90を
学ぶほうがいい。

エディターに関しては、初心者だというなら予約語（？）が色づけされる
タイプのものがいいはずなので、統合環境のものを使うのが吉だろう。

MSの開発環境は複雑化しすぎて初心者には荷が重いだろう。

無料のものといえば、旧ソ連の金のない連中が、g77や著作権切れの
本のPDFなどをセットにしたFORTRAN統合環境を出していたがｗｗｗ
ま、FTNでいいんでないの？

>>147
>どうしたらよいのでしょうか…？

あなたが使っているPCの環境やプログラムの内容がわかるエスパーが現れるまで
しばらくお待ち下さい。

>>141　今日は2次のルンゲ=クッタにしてみたｗ
PROGRAM oppai2
PARAMETER (h = 0.2, a = 1.0, b = -1.0, c = -1.0, d = 1.0)
t = 0.0
x = 3.0
y = 1.0
hh = h / 2.0
DO 10 i = 1, 200
t = i * h
dx0 = (a + b * y) * x * hh
dy0 = (c + d * x) * y * hh
dx = (a + b * (y + dy0) ) * (x + dx0) * hh
dy = (c + d * (x + dx0) ) * (y + dy0) * hh
x = MAX(0.0, x + dx)
y = MAX(0.0, y + dy)
WRITE(6, *) x, y, t, i
10 CONTINUE
END

>>152　インデント

とりあえず、振動が続くようになった。捕食のモデルから予想されるものに
なっていると思う。
ま、普通は４次のルンゲ=クッタなんだがｗ

>>152　4次ルンゲ食った
PROGRAM oppai4
PARAMETER (h = 0.2, a = 1.0, b = -1.0, c = -1.0, d = 1.0)
t = 0.0
x = 3.0
y = 1.0
hh = h / 2.0
DO 10 i = 1, 2000
t = i * h
dx1 = (a + b * (y ) ) * (x ) * hh
dy1 = (c + d * (x ) ) * (y ) * hh
dx2 = (a + b * (y + dy1) ) * (x + dx1) * hh
dy2 = (c + d * (x + dx1) ) * (y + dy1) * hh
dx3 = (a + b * (y + dy2) ) * (x + dx2) * hh
dy3 = (c + d * (x + dx2) ) * (y + dy2) * hh
dx4 = (a + b * (y + dy3) ) * (x + dx3) * hh
dy4 = (c + d * (x + dx3) ) * (y + dy3) * hh
x = MAX(0.0, x + dx1 / 6.0 + dx2 / 3.0 + dx3 / 3.0 + dx4 / 6.0)
y = MAX(0.0, y + dy1 / 6.0 + dy2 / 3.0 + dy3 / 3.0 + dy4 / 6.0)
WRITE(9, *) x, y, t, i
10 CONTINUE
END

>>153　インデント
振動しつつ、ｘ、ｙとも１に収束してゆく。　ｔ＝３５０　くらいでほぼ収束か。

ポトペタ可能な開発環境はありませんか？

ポトペタ　ってなんですか？

>>154
こんなんがあるな。使ったこと無いけどｗ
ttp://www.hulinks.co.jp/software/ginomenustudio/index.html

つーか、GUI は他の言語で作った方が楽だけどな。

>>147
・・・いっぽうロシア人は鉛筆を使った。・・・・

金なしロシアのｆ７７開発環境だｗ
http://www.imamod.ru/~vab/vfort/

多変数多項式のGCD計算がなかなか作れないのですが、どなたか見本みたいな感じでソースを見せて頂けませんか？
あと、グレブナ基底を求めるプログラムもあれば頂きたいです。

>>158
　【FORTRAN７７】　ﾌﾟﾛｸﾞﾗﾑを起動させて！！
http://science4.2ch.net/test/read.cgi/sim/1057990786/68

はい、次の方どうぞー

>>158　いろいろあった気がするが。　　岩波から本が出てなかったか？　買え！

To search the Netlib repository, enter your search query in the form below and
then press the Go button.

Look for : any of the wordsall the words

Query "groebner" found 1 matches.
1. toms/628
keywords: groebner basis, polynomial ideals, rational integers
gams: C3b
title: GROEB
for: canonical (or Groebner) bases of polynomial ideals
by: F. Winkler et al.
ref: ACM TOMS 11 (1985) 66-78
Score: 100%

http://www.netlib.org/cgi-bin/search.pl

すごい資料が出たらしいぞｗｗｗｗ

http://bbs.enjoykorea.jp/tbbs/read.php?board_id=thistory&nid=1578859

これをみれ！　そして１００回saveしろ！

>>159
あ〜、すみません。それ俺です…
向こうは人が少ないし、偶々ここを見つけて、こっちの方が本職かな？と思いまして。

>>160
岩波からですね。漁ってみます。
自分は脳ミソの出来がよろしくないので、わからないときはまた質問させて頂きます m(_ _)m

連続ですみません。

>>160
この検索すごいですね。grobnerもGCDも一発で見つかりました。
まだ中身を細かく見てないですが、
岩波の本を買って参照しながら解読してみます。
助かりました。

>>162
岩波でなかったかもｗ　間違ってたらすまんｗ
プログラムはよくネットに落ちているので検索すればもっと出てくると思う。

>>163
NetLibはインターネット以前からある数学系のサブルーチン集だ。
古いのがごみのように詰まっているｗ　

>>164
それらしい本で
岩波から計算代数と計算幾何(佐々木,今井,他,共著)
共立から計算機数学(野崎著)
を本屋の検索で発見したので、明日図書館で調べてみようかなと思います
groebnerのプログラムが長すぎて脳がスーパーノヴァしてるので、今日は寝ることにします(_ _;)
数値解析の本ばかりで、GCD計算は代数系の入門書からだったので、本当に助かりました

>>165
古くても重宝しますね、これは。
数にビビリました。
ルンゲも、ガウスも、DE公式も、SVDも、なんでもありますね！
無償で頂けるなんて親切な人達がいるもんですね(^-^)

>>166
昔の大型計算センターはどこでもライブラリが大量に公開されているのが普通だった。

ＦＯＲＴＲＡＮの世界は歴史が古いので、ソフト囲い込み運動と
その反動のオープンソース運動とは無関係にソース公開共有の
伝統がまだ残っている。田舎で鍵をかけないで出かける世界？ｗ

さぁ次の宿題もってこい！

配列を返す関数はどのように定義するのですか？

>>169
program hoge
　implicit none
　real, dimension(3) :: x, y, z
　real, dimension(5) :: u, v, w
　x = (/1, 2, 3/); y = (/2, 3, 4/)
　u = (/4, 5, 6, 7, 8/); v = (/10, 20, 30, 40, 50/)
　z = adda(x, y)
　w = adda(u, v)
　print *, z
　print *, w
　stop
contains
　function adda(a, b)
　　real, intent(in), dimension(:) :: a, b
　　real, dimension(size(a)) :: adda
　　integer :: na, nb
　　na = size(a); nb = size(b)
　　if (na > nb) then
　　　adda(1:nb) = a(1:nb) + b(1:nb); adda(nb+1:na) = 0.0
　　else if (na < nb) then
　　　adda = a(1:na) + b(1:na)
　　else
　　　adda = a + b
　　end if
　end function adda
end program hoge

>>170
レスありがとうございます！わかりました！

pathとファイル名指定して、そこでファイルオープンするような処理を
書きたいのですけど、char型に突っ込んだ後の、うしろの余ったスペースが
どうにもならないっす。
みなさんはうしろの余ったスペースどうやって消してますか？

age忘れたので、ageさせてください。すいません。

>>172
! TRIM 関数, LEN_TRIM 関数の使い方
program huge
implicit none
character(len=30) :: path, fname
character(len=70) :: r1, r2

print *, "path:"; read *, path
print *, "file name:"; read *, fname
r1 = path // fname ! そのまま連結
r2 = trim(path) // fname ! 末尾の空白を除いて連結

print *, r1
print *, r2
print *, len(r2) ! 文字列の長さを返す
print *, len_trim(r2) ! 末尾の空白を除いた文字列の長さを返す
stop
end program huge

>>174
わざわざプログラムまで書いていただいて、ありがとうございます！
おかげでわかりました！

1 FORMAT(1H ,I5)
2 FORMAT(1H ,E15.7)
3 FORMAT(1H ,I5,1H ,I5,1H ,7E15.6)
4 FORMAT(1x,I5,1X,I5,1X,E15.6,E15.6)
5 FORMAT(B20)
6 FORMAT(I5,1H ,I5)
7 FORMAT(1H ,I5,1H ,3E15.6)
11 FORMAT(7X,I5,1X,E15.6,E15.6)

fortran90からはじめたのですが上のFORMATｎ意味がよくわかりません
どなたかご教授おねがいします

HはホレリスのH。廃止勧告が出ていて、いまではHな人しか使わないｗ
基本的には文字列を囲む引用符”　”の祖先で、以降が文字列ということを示している。
文字数はHの前の数字で表している。
FORTRAN66までは文字型が無かったという事情がある。

１文字目が１Hになっていることが多いのは、かつてのFORTRANでは出力の１文字目は
プリンター制御文字になっているという事情がある。改行、改ページ、改行せずができる。
（ビデオコンソールの普及は７０年代以降なのでコンソール出力を前提とした発想は捨てよ）。

それゆえ、この１文字目を空白で飛ばさないと、頭の１文字が消えてしまうという現象がおきる。
いまでもコンパイラーごとに１文字目の扱いは異なる。VisualFOTRANなどではコンパイラ
オプションで１文字目の扱いを変えることができる。

最近のコンパイラーは画面出力では１文字目をプリンター制御文字と解釈しない傾向がある。
そういうコンパイラーで実行すると、常に行頭に１文字空白が付いたり、時々０，１が現れることになる。

プリンター制御文字　４通り
０　２行改行
１　改ページ
＋　改行せず
その他　１行改行（普通の出力）


その他
Xは入力では頭の数だけカラムをスキップ、出力では頭の数だけ空白出力。
Bは２進法で出力。B２０なら２０桁の幅で2進法出力。
I,ESは基本なのでマニュアル嫁。

１H　が冒頭に連発するのは初心者には理解できない躓きポイントではあるが、
これは上代の古語が出てきたと思って、ももしきや〜と古きをしのんで喜ぶべき
ポイントでもある。歴史と伝統のFORTRANの醍醐味でもある。
あなとうと、あなうるわし。

FORMATは昔からFORTRANの難所の一つとされてきたが、
最近はFORMATのほうも進化してきたので学び甲斐もあるべ。

配列出力で、：や（）による反復指定を使って思い通りに出力できると、
他の言語では得られない快感が得られる。

ありがとうございました＾＾
がんばって理解していこうと思います

COMMON /M/!,N,LZ
COMMON /P/P,GAMMA,GPH

とゆうCOMMON文があるのですがCOMMON文自体を理解していなく
何の処理をしているのかいまいちわかりません
ｆ９０ではMODULE文にできるらしいのですがどうのように
書き直せばいいのかどなたかご教授お願いします。

>>181
COMMONブロック P の先頭から順に実数型変数 P, GAMMA, GPH が割り当てられる。
（同じプログラム単位中でP, GAMMA, GPHの定義がない場合）

COMMONブロック M は変数名の並びがコメントアウトされているので無効？
（エラーになる？）

モジュールにするとこんな感じ
module hoge
implicit none
private
!integer, public, save :: N, LZ
real, public, save :: P, GAMMA, GPH
end module hoge

但し、他のプログラム単位中でCOMMONブロックM, Pが他の用途に使われていたり
P, GAMMA, GPHの並びが違っていたりすると正常に動作しないことがある。

COMMON文はグローバル変数としての機能もあるが、その概念を外延すると
理解しにくいこともある。むしろ同じメモリー領域の手動割付と考えると分かりやすい。

昔メモリーが少なくて貴重だった時代には、手動で半動的にメモリー割付を
してメモリーを有効活用していた。こういう場合同じCOMMONブロックに全然
対応の付かない変数が書かれていたりする。

COMMONブロックでは、書いた変数の名前の順に先頭から順番にメモリーが
割り振られてゆく。

COMMON /P/P,GAMMA,GPH

この場合Pという名前のブロックに、（P,GAMMA,GPHを単精度４バイトとして）
P,GAMMA,GPHに対応する４バイトの領域が３つとられる。

このコモンブロックを別のサブルーチンでは、COMMON/P/A,B,C　と受けることも
出来る。この場合P,GAMMA,GPHの内容はA,B,Cと一致する。（変数の番地が
一緒になる）

またこのコモンブロックを別の形で受けてもいい。
COMMON/P/A(3)
で受けた場合、P,GAMMA,GPHの内容はA(1),A(2),A(3)に対応する。

重要なのは順番で変数名ではない。これはMODULEの場合と異なる。
COMMON /P/GAMMA,GPH,P
で受けたとすると、元のP,GAMMA,GPHはそれぞれGAMMA,GPH,Pに対応する。

このように違う名前でCOMMMONを受けていなければ、単純にMODULEに直せる。
そうでない場合は慎重に作業する必要がある。

この式のプログラムできませんか？さっぱりわからんのです。。。
F(P)=(9/(128*PAI**2))*(1/P)|EXP(-P*(1/TAN(X/2)))*(1+COSX**4)dX
[|]はインテグラルで、範囲は０からπ/２です。シンプソンのプログラムはあるんですが、もうさっぱりです！
どなたかお願いします。。。

>>184　
tan(X/2)=1/Yとおくと積分範囲は∞〜１

dX=-2/(1+Y^2)dY

cos(X)=(Y^2-1)/(Y^2+1)

1/(1+cos(X)^4) = (Y^2+1)^4/(Y^8+6Y^4+1)

代入して
(1/p)∫exp(-p/tan(X/2)) / (1+cos(X)^4) dX = -(2/p)∫(Y^2+1)^3 exp(-pY) / (Y^8+6Y^4+1) dY
積分範囲を１〜∞に反転させれば頭のマイナスが消える。

とりあえず原点付近の１／ｔａｎの発散を避けるようにしてみた積分範囲が無限大にまでになるが、
被積分項は急速に小さくなるので適当なところで打ち切れるだろう。

適当なので間違ってたらスマソ。

fortranは、コマンドライン引数使って実行できないのですか？

>>186
iargc()
getarg

>>186
規格に入ったのはFortran2003から、Fortran95までの規格には存在しない。
しかしたいていはベンダー拡張ルーチンでできることが多い。
マニュアルを調べてみるといい。POSIX互換の形式をとることが多いので、>>187が
書いたようなコマンド名になることが多い。

>>187, 188
なるほど、勉強になります。ありがとうございます。

>>182 183
ありがとございます m(__)m

VAXVMS上で動く，Fortran IVで書かれたプログラムを
Fortran77以降のものに自動変換してくれるプログラムとか
無いですか？

>>191
見たこと無い。

宿題を頼んだ人は事後報告をするように。

>>193
宿題に困るレベルの人が総括（事後報告）できると思う？
「ありがとうございます」で納得してあげてちょ。

解答もらったらそれっきり音沙汰ない奴もおるからのぉ。

さて、遠慮せず次の宿題を持ってきたまへ。

宿題たまってるけど眠いのでもう寝ます。
お休みなさい。

Fortranで書かれた数値計算のプログラムをC#に移植しようとしています。
至る所で引数に配列の要素数をもつ関数が定義されているのですが、
C#のGetLengthやLengthをつかってその都度判定するような仕様にしては問題あるのでしょうか。

つまり、GetLengthなどの関数を実行する時間を節約するために要素数を引数で与えているのでしょうか。
それともFortranでは配列数を取得するような関数が無いためなのでしょうか。

>>198
FORTRNでは、FORTRAN９０になるまで配列の大きさを知る命令は存在しない。
FORTRAN９０で関数SIZEが導入された。

FORTRAN７７以前では、配列の大きさは副プログラム（サブルーチン、関数）に、
引数として明示的に渡すか、COMMON等で渡すか、あるいは”知っている”もの
としてサブルーチン側で勝手に使うかした。

FORTRANでは変数はアドレス渡しであり、FORTRAN７７以前の書式を使うと
配列であっても渡しているのは先頭要素のアドレスだけで、配列の大きさもの情報も
次元も渡していない。（FORTRAN９０で導入された新しい書式を使うと配列の
サイズや次元の情報も渡される。）

ここで気をつけなければならないのは、FORTRAN６６時代のプログラムでは、
サブルーチン側の引数の仮配列の大きさを実際の長さとは無関係に１０などにして
内部では好き勝手な長さとして使用する書き方が標準的に採用されていたことだ。
（FORTRAN77では、サイズを＊で定義するやり方に相当するものかな？）
こういう書法の伝統を引いている場合、>>198の移植方法では大混乱をきたすだろう。

参考　（ｔabが消えるので７カラム目までの空白は脳内で補完してください）
＊FORTRAN66時代にありがちなプログラム
REAL A(20, 20), B(200)
CALL ZERO(200, B)
CALL ZERO(400, A)
CALL ZERO(20, A(1, 2))
CALL ZERO(20, A(1, 3))
CALL ZERO(20, A(1, 4))
END
SUBROUTINE ZERO(N, X)
REAL X(10)
DO 10 I = 1, N
X(I) = 0.0
10 CONTINUE
END

これは昔よく使われたやり方。サブルーチンに渡されるのは配列の先頭要素の
アドレスだけなので、次元が整合していなくても問題はない。また配列の途中
だけを渡して部分的に利用しても問題はない。（FORTRANの配列はCとは
逆順で列方向に対して連続的）

CALL ZERO(200, B)　これは素直な呼び方。
CALL ZERO(400, A)　これは２次元配列を１次元的に還元して渡している。
CALL ZERO(20, A(1, 2))　これは２次元配列の２列目だけを渡している。
CALL ZERO(20, A(1, 3)) この呼び方は列ベクトルの使用法として良く使われる。
CALL ZERO(20, A(1, 4))　同様に３，４列目を渡している。

＊FORTRAN77時代
SUBROUTINE ZERO(N, X)
REAL X(N)
または
REAL X(*)

１９９，２００で書いたこととは無関係に一般論的な結論を言わせてもらうが、
fortranでは宣言した配列の大きさとは無関係な大きさで配列を使うことが
ままあるので、引数として渡している大きさを配列の宣言長に直すのは、
とても危険である。

FORTRANのライブラリなどでやたらと配列の寸法を引数として求められるのも、
配列は先頭要素のアドレスしか渡しておらず、そのほかの情報を一切渡していない
からである。

またFORTRANでは、この性質を大いにプラスに利用しているので
（というかこれが無いとまともな実用的なプログラムが書けない）これを非難するのは
お門違いである。

FORTRAN９０では、配列の次元や寸法も渡しているがその性で、>>200で書いた
ような列ベクトルを渡す時に、余分な処理が必要になって速度が低下する。
こういう時にFORTRAN77式の先頭アドレス渡しの簡潔さが心地よく感じられる。

>>199

丁寧なご解説、痛み入ります。
なるほど、配列に関する色々な情報ではなく、配列のアドレスを渡しているだけなのですね。
そういえば他の箇所で多次元配列が1次元の配列として処理されていたり、
首をかしげたおぼえがあります。
このあたり、移植の際には気をつけないと駄目ですね。
ありがとうございました。ご指摘の点を踏まえて頑張ってみます。

10年以上前の話だけど、
なんとなく入学した大学で、入学直後からFORTRAN77の宿題を大量に課せられた。
もちろん僕にはFORTRANが一体何なのかもよくわからなかった。
とりあえず、テキストをマネしながらキーボードを叩いてはみたものの、
なんでこんな勉強しなきゃいけないのかわからなかった。
正直言って、楽しくなかった。
同学部の友達が居なくて、先生も冷たくって、僕は１人だった。
結局、大学には通わなくなり、僕はその場所を去った。
その後、僕はプログラムとは無縁の人生を送ることになった。

このスレにFORTRANの宿題を持ってくる人・・・
それに対して丁寧に解説をされている方々・・・
ログを上から読んでみたけど、なんだかスゲー楽しそうじゃん？
なんだか羨ましいよ。マジで。

あの頃、この2chのFORTRAN（超初心者用）スレがあったら・・・・。
なんて、思ってしまったけど、仮定法過去完了で人生を語るのはやめようと思う。

皆さんお元気で。
ひっそりと心の奥底で応援いたします。

>>203
昨日の涙を今日の笑顔に変えて、でっかい夕日に向かって走るんだ！

それはそうとFORTRANの教科書や教育課程はろくなものが無い。
面白くなくて訳分からなくなるのが標準だ。

やさしく分かりやすく書くということを、内容を整理して書くことではなく、
本文と無関係な下手糞な挿し絵をちりばめる事だと思っている世界だｗ

さぁ宿題をもってきたまえ

>>203
僕はＦＯＲＴＲＡＮスレで助けてもらったおかげで
大学院を無事修了できました。
ここの住人の方々には本当に感謝しています。
僕もがんばって、知恵で人を助けられるような人間になりたい

>>204
うむうむ
「作法」が執筆された時に 77 があったら良かったのに。
ratfor で見た目がちょっと違うというだけで、
初学者には敬遠されてしまうから。

これも仮定法過去完了だぬ。

>>206
おめ！
これからも山崎渉を・・じゃなかったFORTRANスレを応援してくださいｗ

>>204
『FORTRAN 入門』浦昭二　が名著だと思う。
うちの大学ではどこでも転がっている標準本。
Fortran90 への対応を希望する。

>>209
オレもそれのFORTRAN77版で習ったんだけどさ、
あれが名著と言われるようでは日本もお終いだと考えるべきだと思う罠ｗｗｗ

前半例題で後半が文法という形式以外は例題もいまいちだしいいとこなしだと思うよ。

日本語で書かれたFortran90以降の機能が詳しく解説されてる本ってない？
近所で探してみたけど、どれも全くの入門書止まりだったり
FORTRAN77に毛が生えたような本しか見つからなくて。

>>211
これが結構いい。
コードも比較的Fortran90の機能だけで書こうとしている。

入門 Fortran90実践プログラミング SOFTBANK BOOKS
東田 幸樹 (著), 山本 芳人 (著), 熊沢 友信 (著)
ttp://www.amazon.co.jp/exec/obidos/ASIN/4890524932/249-3009705-2437946


しかし絶版だったｗｗｗｗ

>>210
そうかな？
文法事項を段階的にかつ網羅的に適用した、
適切な入門例題だと思うが。
科学技術計算を意識した例題が多いというのも、
FORTRAN という言語の性質上、好ましい。

>>211
そうだよな。
Fortran90 以降の名著が少ない（あっても絶版）というのは、
Fortran 市場を反映しているようで、残念でならない。
『FORTRAN 入門』浦昭二を90向けに改訂するなら、喜んで手伝うが。

>>209-210
あの本はFORTRAN77版だけじゃないんですか？
詳細ｷﾎﾞﾝﾇ．（"・ω・`）ﾉ

>>214
あの本はFORTRAN６６版がある。
アマゾンとかで検索してFORTRANというだけの７７のついていないのがそれ。
それを元にFORTRAN77版を書いているので、
わりとすっきりしていたFORTRAN６６版に比べるとFORTRN７７版は
いまいちな気がしないでもない。


FORTRAN90の教科書だったら、FとかLahey ELFのようにFORT７７以前との
互換性を捨て去ったバージョンで書くのがすっきりしていいと思う。Fなんかは
色々入門書が出ていたと思う。もっとも肝心のFがNAGのフロントエンドを入れた
せいで、当初の過去互換性を切り捨てた簡潔な言語になるという理想から
逸脱して迷走しているのがずっこける。
http://www.fortran.com/F/
http://www.fortran.com/

FORTRAN９０で導入された新しい機能を用いれば、昔PASCALなんかで
やっていた正則な構造化（モジュール化）されたプログラム書法が身につくと思うのだが、
いかんせん現実問題としてはFORTRAN77やFORTRAN６６のプログラムをポンと渡されて、
動かせと言われたりするから、まあどうしたものか分からん。

PROGRAM GGGG
IMPLICIT NONE
INTEGER*8 :: N,K,NH,M,I
REAL*8 :: GANMA2(0:10000),GANMA3(0:10000),GANMA4(0:10000)
REAL*8 :: BETA2,BETA3,BETA4,SECHYP,EAMP0,ER(0:2047),EI(0:2047)
REAL*8 :: COS,SIN,ERK,EIK,T0,T,PI,FRQ,FQ,H,C,CD,CDD,WL,ROU,TAU0
COMPLEX*16 :: E(0:4094),EK
M=11
N=2**M
NH=N/2
PI=4.0D0*DATAN(DBLE(1.0D0))
T=5.12D-12
TAU0=0.22D-12
C=3.0D+8
CD=-42.0D-1
CDD=0.0D0
WL=1.55D-6
H=5.0D-5
BETA2=-CD*WL*WL/(2.0D0*PI*C)
BETA3=WL*WL*WL*(WL*CDD+2.0D0*CD)/(4.0D0*PI*PI*C*C)
BETA4=0.0D0
T0=TAU0/(2.0D0*DLOG(DSQRT(DBLE(2.0D0))+DBLE(1.0D0)))
EAMP0=10.0D0

DO K=0,N-1

SECHYP=DBLE(1.0D0)/DCOSH((T/T0)*(DBLE(K)/DBLE(N)-DBLE(0.5D0)))

E(K)=CMPLX(EAMP0*SECHYP,0.0D0)

END DO
DO K=0,NH-1
ROU=DBLE(K)*T0/T

FRQ=PI*ROU/T0
GANMA2(K)=-H*BETA2*PI*PI*ROU*ROU/(T0*T0)
GANMA3(K)=-H*BETA3*PI*PI*PI*ROU*ROU*ROU/(T0*T0*T0)*2.0D0/3.0D0
GANMA4(K)=-H*BETA4*FRQ*FRQ*FRQ*FRQ*1.0D0/3.0D0
FQ=ROU/T0

END DO

DO K=0,20,1
EK=E(K)*CEXP(CMPLX(0.0d0,(GANMA2(K))))

EK=EK*EXP(CMPLX(0.0d0,(GANMA3(K))))

EK=EK*CEXP(CMPLX(0.0d0,(GANMA4(K))))

E(K)=EK
END DO

END PROGRAM GGGG

DO K=0,20,1
ERK=ER(K)*DCOS(GANMA2(K))-EI(K)*DSIN(GANMA2(K))
EIK=EI(K)*DCOS(GANMA2(K))+ER(K)*DSIN(GANMA2(K))

ERK=ERK*DCOS(GANMA3(K))-EIK*DSIN(GANMA3(K))
EIK=EIK*DCOS(GANMA3(K))+ERK*DSIN(GANMA3(K))

ERK=ERK*DCOS(GANMA4(K))-EIK*DSIN(GANMA4(K))
EIK=EIK*DCOS(GANMA4(K))+ERK*DSIN(GANMA4(K))

ER(K)=ERK
EI(K)=EIK

END DO

END PROGRAM GGGG

連投すみません、216,217と計算が実数か複素数かだけの違いなのですが
217ではなぜか範囲を広げるほど虚数の値がかけ離れていきます
桁落ちでもないのですが、どなたか考えられる原因を教えてくださいm(__)m

>>218
範囲を広げるという意味が良くわからんので、適当に書いてみるｗ

EK=E(K)*CEXP(CMPLX(0.0d0,(GANMA2(K))))
EK=EK*EXP(CMPLX(0.0d0,(GANMA3(K))))
EK=EK*CEXP(CMPLX(0.0d0,(GANMA4(K))))

この式で、EK(K)は最後に掛けたほうがいい。なぜなら、それ以外の部分は
全部位相因子なので幾ら掛けても、絶対値は１になるはずの量だから。
まずはそこでチェックしてみたらどうか？

つまり
EK = EXP(CMPLX(0.0d0,(GANMA2(K))))
EK = EK * EXP(CMPLX(0.0d0,(GANMA3(K))))
EK = EK * EXP(CMPLX(0.0d0,(GANMA4(K))))

ここでABS(EK)＝１．０ｄ０のはず（誤差を抜かして）　
もう一つの式でも同様にして、その段階で虚部を比較すれば、
より限定的な情報が得られて問題解決に近づくだろう。

EK = EK(N) * EK　
（ここで後ろのEKはpahase_factorとかの変数名にしたほうがいいと思うが
まぁ好みの問題だ）


あとCEXPは総称名EXPに統一したほうがいいだろう。
精度の規定値が気になるなるが、KIND=8をつければ強制できるはず。

> ERK=ERK*DCOS(GANMA3(K))-EIK*DSIN(GANMA3(K))
> EIK=EIK*DCOS(GANMA3(K))+ERK*DSIN(GANMA3(K))

とあるけど、一行目で、ERKを書き換えてしまってはダメじゃないの？

>>220
実数、虚数にわけてる計算はあっているので、ただしいはずです・・・

>>219
ありがとうございます、試してみます
EK=EXP(CMPLX(0.0d0,(GANMA2(K)),8))
といった感じで精度をあわせてしまえばいいんですよね？

>>221
いや、１行目でERKが書き換わっているので、正しいとは思えない。
その2行で１つの計算のはずだから、虚部を求める時に新しくなった実部を
使ってはいけない。

これで正しいとすると、二つが一致しないのは当然。
ただこれが正しいことはまずありえないと思う。

>>221
いや、式はあってるんだけど、1行目でreal partを計算する際に、
元の変数のreal partを書き換えてしまっているのでは？
2行目の右辺に出てくるERKはもはや元の変数のreal partではないないと思うけれど。

かぶったw

式があっているってのは、実部虚部への分け方はあっている、ってことね。
そのプログラムは変だと思うよ。

>>223~226
気づかなかった；；；サンクスです。

もりあげようぜ！ｗ

盛り上げようにも大学が春休みで宿題がない予感。

欠測のある平均値の求め方を教えてください。
10行ごとの各列の（欠測を除外した）平均値と欠測数を求める問題です。

時間,A,B,C
1,23.5,23.2,22.7
2,24.5,25.2,22.1
3,25.4,24,27.5
4,26.4,,22.4
5,23.5,22.4,
と2000行、続きます。

>>230 適当なのでもっとうまい方法もあるやも知れぬ。参考までに。
INTEL　FORTRAN　Ver.9

PROGRAM vipp
IMPLICIT NONE
INTEGER :: i, j, k
REAL :: x(10)
i = 0
DO
i = i + 1
x = 0.0
READ(9, '(i5, 10F15.5:)', ADVANCE = 'NO', EOR = 98,END = 99) k, (x(j), j = 1, 10)
98 PRINT *, j, ":", SUM(x(1:j)) / REAL(j)
END DO
99 PRINT *, "EOF"
STOP
END PROGRAM vipp

実行結果
C:\Documents and Settings\All Users\Documents\FORTRAN\vip\Debug>vip
3 : 23.13333
3 : 23.93333
3 : 17.63341
3 : 16.26667
3 : 15.30000
EOF

>>230
ごめｗ　列ごとの平均かｗ　

風呂から上がったらまた考えるｗ

>>232
トライ、ありがとうございます。
温度なので、23度あたりになるはずです。
5列分の平均ですが、
A 24.66 , B 23.7 , C 23.675

内部ファイルとLEN 関数をつかって、
うまくできないかとチャレンジ中です。

くたばれFORTRAN

がんばれFORTRAN

>>233
単純に思いつく方法は
１．手抜き方式、０．０が読み取られたときを欠損値とするというやり方
２．内部ファイルに読み込んで、INDEX関数で”，”を区切りとして切り取るやり方

汎用性からみると２かな？

第３の方法を求めて
FORMAT文がF95で進化したので、空データの読み取りが検知できないか今調べているｗ

>>233
使用コンパイラ　Compaq VFでは、
空白を欠測ではなく0と認識してしまったので、ダメでした。
文字列として読み込んで、^@ の処理を考えます。

こんな感じッスかね？

program hoge
implicit none
integer, parameter :: N = 10, NT = 2000
real, parameter :: NG = -1.0E10
real :: A(N), B(N), C(N), t, aa, ab, ac
integer :: ka, kb, kc, i, j, ios
character(len=80) :: tmp
open(unit=11,file="huge.dat",status="OLD",blank="NULL")
loop_j: do j = 1, NT/N
A = NG; B = NG; C = NG
loop_i: do i = 1, N
read(unit=11,fmt="(A)",iostat=ios) tmp
if (ios < 0) then
! エラー処理は省略...
exit loop_j
end if
read(unit=tmp,fmt=*,iostat=ios) t, A(i), B(i), C(i)
end do loop_i
ka = count(A <= NG); aa = sum(A, mask=(A > NG)) / (N-ka)
kb = count(B <= NG); ab = sum(B, mask=(B > NG)) / (N-kb)
kc = count(C <= NG); ac = sum(C, mask=(C > NG)) / (N-kc)
print *, "欠測 = ", ka, kb, kc
print *, "平均 = ", aa, ab, ac
end do loop_j
close(11)
end program hoge

test

>>238
ありがとうございます！
無事に解決できました。
10個とも欠測の時間帯の平均値が NaN と表記されて、
いたれりつくせりです。

count(A <= NG)　と　mask=(A > NG)
がミソですね。
0割になりそうでならないのが少し不思議です。

出遅れたかｗｗｗｗｗｗｗｗｗｗｗ

PROGRAM vipp
IMPLICIT NONE
INTEGER :: i, j, k
CHARACTER(LEN = 136) :: buff
REAL :: a, x(10) = 0.0
INTEGER :: n(10) = 0
DO
READ(9, '(a)', END = 99) buff
k = INDEX(buff, ',')
READ(buff(:k), '(i10)') j
buff = buff(k + 1:)
DO i = 1, 10
k = INDEX(buff, ',')
IF (k == 0 .AND. LEN_TRIM(buff) == 0) EXIT ! end by ","
IF (k /= 1) THEN ! k == 1 is missing data : ",,"
READ(buff, '(F15.0)') a
x(i) = x(i) + a
n(i) = n(i) + 1
END IF
buff = buff(k + 1:)
IF (k == 0) EXIT ! end of record
END DO
END DO
99 PRINT *, "EOF"
PRINT *, n
PRINT *, x / REAL(n)
STOP
END PROGRAM vipp

>>241

C:\Documents and Settings\All Users\Documents\FORTRAN\vip\Debug>vip
EOF
5 4 4 0 0 0
0 0 0 0
24.66000 23.70000 23.67500 NaN NaN
NaN NaN NaN NaN NaN

>>241
わかりやすい注釈までいれていただいて、ありがとうございます。
k = INDEX(buff, ',')
buff = buff(k + 1:)
LEN_TRIM(buff) == 0
など、文字列処理の強力さはさすがですね。

配列手続といい、文字列処理といい、
f90 の強力さに、あらためて惚れなおしました。

>>241
今、調べなおしてみると、
LEN_TRIM 以外は、f77 の機能のようですね。
f77 の文字列処理の先見の明に驚いてしまいます。

>>243
場合わけはあまりエレガントでなかった。もう少し冗長でも分かりやすくすべきかも。


buff = buff(k + 1:)
これで長さが足りない分に空白が入るかどうか確かめておいた。
ちゃんと空白が入っているようなので、ごみが尻のほうに残ることは無い。

文法書にもそう書いてあった気がするのだが、記憶がアイマイミー。

>>245
INDEX を用いた切り分けや、テクニカルな文字列処理の手本
を示していただいただけでも、感謝感激です。
本当にありがとうございました。

>>246
喜んでもらえてうれしいお（＾ω＾）ﾉｼ

>>238も勉強になった。漏れは固定FORMATの愛好者なのでCVS形式はほとんど使ったこと無い。
自由書式の時にCVS形式の空データは読み飛ばすとは知らなかった。０が入ると思っていた。



F77はINDEX関数がとても便利だったが、TRIM関数が無かったせいで初心者時代に苦労した。
ある程度習熟してからは、不定長の文字引数で尻から先頭に向かって空白以外が出るまで
数える関数などを自分で用意できるようになったが、そこに到達するまではかなり時間がかかった。


内部ファイルの使い方も共同研究の外国人の書いたプログラムを読んでいて初めて知った。
それまでこんな便利なものがあるとはまったく知らなかったｗ
後になって改めて文法書を読むと、文字変数を入力装置にできるとはさりげなく書いてあるのだが、
それがとても便利なものだということはまず分からない。


この辺のノウハウを徒弟制度や口伝に頼っているFORTRAN界の未来は暗いなと思ったりしてｗ

ところでFortranに分割コンパイルっていつ入るんですか？

>>248
単純に複数のプログラム単位をコンパイル＆リンク出来るという意味でなら
FORTRAN IIで副プログラム呼び出しとBinary Symbolic Subroutine Loaderが
使えるようになった1958年から。

複数のプログラム単位のコンパイルを自動化出来るという意味でなら
FORTRAN IIモニタの出来た1959年から。

近代的なモジュールの機能が使えるようになったという意味でなら
Fortran 90の規格が出来た1991年から。

ということでよろしいかな？

非線形連立代数方程式をとくFORTRANのコードをC++に移植しています。
解の一部がどうしても両プログラムで一致せず、
丸一日かけて色々と試行錯誤したところ、Fortranでreal型を使っていたところ、
c++でdouble型にしてしまっていたことが原因でした。

このような場合、c++への移植の際にもfloat型を使うなどしてもとのコードと精度を合わせた方がよいのでしょうか。
なんとなくより大きな範囲を扱える変数型の方がよさそうな気がするのですが、
移植元が倍精度を使っていないので迷っています。
移植元のコードは1980年頃に書かれたもので、計算時間との兼ね合いから単精度を選択したと捉えたほうがよいでしょうか。

また、計算誤差の発生しないdecimal型があると知りました。
一般に数値計算をする場合はこのような型を使うのでしょうか。
精度が高い代わりに計算速度が遅いと聞きましたので、
一般にはどのような基準で使い分けて折り合いをつけていくのか興味があります。

すみませんが、どなたかご存知の方がいらっしゃいましたらお教えいただけますか。
よろしくお願いいたします。

>>250
精度が出てるのを確かめられたならdouble型にすべし。
floatにする必要はないと思う。
C++の組み込み型にdecimalがあるとはきいたことないが、
だた10進数の内部表現ができるってだけで(0.1とかが正確に表せるのかな？)
数値計算自体に誤差が出ないわけではない。はず・・・たぶん

ところでFORTRANにリアルタイムイベントハンドラはいつ入るんですか？

>>250
1980年代のコードである場合、浮動小数点には今とは違う事情がある。
現在の浮動小数点はIEEEの形式にほぼ統一されているが、かつてはメーカーごとに
実数の内部形式は異なっていた。そのために、同じワード長でも仮数部により多くの
ビットを割り当てているものもあれば、指数部により多くのビットを割り当てているもの
もあって、機種によって有効桁の範囲が違っていた。
（IBMは仮数部に、CRAYは指数部により多くのビットを割り当てていた。
さらに言えば、ＮＥＣのＡＣＯＳのあるシリーズのように１ワード＝３６bitとか、
CDCのいくつかのシリーズのように１ワード５６？bitとかでワード長自体が
違うものもあった。）

というわけで、書かれた対象機種が分かれば必要精度に関するヒントになる。

ただ数値演算プロセッサが倍精度で作られている以上
基本的に倍精度にしておくのがよろしいと思う。

まぁSSEで高速化を優先するならまた話は別だが・・・

>>251, 253

ありがとうございます。
頂いたアドバイスに従って
とりあえず、倍精度で移植を進めてみようと思います。

>現在の浮動小数点はIEEEの形式にほぼ統一されているが、かつてはメーカーごとに
>実数の内部形式は異なっていた。

そういえば、コンピュータの精度に応じてコメントアウトすべしというような注釈のついた関数がありました。
そういうことだったんですね。勉強になります。

>>250
DECIMAL型について言えば、あれは商業用の機能で、数値計算的にはあまり
利益がないと思う。特に計算スピードはがくっと落ちるはずである。

専門ではないから正確ではないが、商業の世界では利息で生じた小数点以下の
端数をあるところで切り捨てる（または四捨五入？）する慣例がある。浮動小数
点では２進法を使っている関係で丸め誤差から小数点以下の端数の切り捨てで
１０進法の時と微妙なズレが生じ、大量の金額に対して計算すると慣例とのズレ
が無視できなくなる。これに対する対策がdecimal型すなわち１０進型である。

一般的には、１桁ごとに１０進法の数字を１個割り当てて計算するので極めて
計算効率が悪い。１桁につき最低４ビット必要。同じビット数で表される数の
範囲が狭くなってしまう。
（確か８０８７コプロにもBCD演算機能があったはずだが、固定小数点だったような
気がする。）

数値計算的には、倍精度、倍々精度という方向に行くことはあっても、
１０進演算の方向に行くことは特殊な場合を除いてないと思う。
実際FORTRANで１０進演算はできないが、利用者からの不満は聞かない。
一方COBOLやPL/Iなどでは１０進演算が出来るようになっている。

>>254
＞そういえば、コンピュータの精度に応じてコメントアウトすべしというような
＞注釈のついた関数がありました。

昔は精度に関する機種依存性が結構あって移植の問題になっていた。

Fortran90で数値の精度に関する関数が大量に導入されたのは、
この辺でみんな苦労していたからだと思う。ただ９０年代に入ってから
IEEE形式に集約されて行ったので、ややありがたみが減った。


最近は０割りが発生しても強行突破して計算を続けてくれるので便利だと思う。
投機的実行の為にわざわざIF文で排除したはずの０割りを実行して糸冬了とか
脳の血管切れそうな事態が起きなくなった。

Fortran2003（一部はFortran95の拡張規格）でIEEE754に関する様々な関数が
導入されているので、将来はよりマニアックな要請に対応できるはずである。

PROGRAM TEST
IMPLICIT REAL*8 (A-H,O-Z)
IMPLICIT INTEGER (I-N)
*
DIMENSION G(50)
*
WRITE(*,*)'個数(入力数+1)'
READ(*,*)N
*
WRITE(*,*)'整数'
DO 10 I = N, 0, -1
READ(*,*)G(I)
10 CONTINUE
*
DO 500 I = N, 0, -1
WRITE(*,*)I,G(I)
500 CONTINUE
*
DO 510 J = N, 0, -1
WRITE(*,*)J,G(J)
510 CONTINUE
*
END

これ実行するとＧ(０)の値がちゃんと出力されないんですが、何故ですか？
ＤＯ５００とＤＯ５１０で、出力が違う値で出てきてしまいます。
あまりにも簡単すぎる質問で申し訳ないですが、誰か教えてください。

frtで実行するとこんなのが出ます。
↓

jwe0019i-u The program was terminated abnormally with Exception Code EXCEPTION_A
CCESS_VIOLATION.
error occurs at _MAIN__ line 19 loc 0040108c offset 0000008c
_MAIN__ at loc 00401000 called from o.s.
error summary (Fortran)
error number error level error count
jwe0019i u 1
total error count = 1

>>257
配列の使い方を復習しましょう。
ttp://www.g.dendai.ac.jp/lecture/fortran/fort07.html

>>259
なるほど・・・初歩の初歩ですね。
とてもよくわかりました。
助かりました。ありがとう。

>>257
数学科や計算機学科での自然な数は０から始まるのだが、
FORTRANの自然数は１から始まる古典的な自然数体系を取っている。

したがってG(50)と宣言したときはG(1)〜G(50)までの領域が取られており、
G(0)はout of rangeになる。

エラーメッセージにACCESS_VIOLATIONと出ているが、これはは大抵
配列のはみ出しと思っていい。

G(0)がどうしても使いたければG(0:50)という形で宣言すればいいが、
まあ素直に１から始まる配列になれたほうがいいだろう。

エラーメッセージから見る限り富士通のコンパイラを使っていると思うが、
デバッグ用のオプションに配列はみだしをチェックするサブスクリプトチェックの項目が
あるだろうから、デバッグ段階ではそれをオンにしておけばはみだしが発生したとき
そのソース行を教えてくれる。エラーが取れたら、チェックを外せばよい。

>>261
めちゃめちゃ親切に説明して頂いて助かります。
富士通のコンパイラまでわかるなんて、凄いです。
エラーメッセージの意味が理解できなかったので、まだまだでした…精進します。

自分、257でした。
259さん、すみません。

64bitに対応したトリップ検索プログラムありますか？

TextSS　のWindowsXP(Professional)64bit化おながいします

もしくは64bitにネイティブ対応したテキスト置換ソフトありますか？

この板でも爆撃喰らうんだな。ＶＩＰでは日常茶飯事だがｗ

xxxx01.txt
xxxx02.txt
....
みたいな連番になっているファイルを順番に読み込むようにしたいのですが、
読み込むファイル名の部分をどのような書式にすれば良いのでしょうか？
（（文字列+整数）をまとめて文字列として認識する方法が分かりません）

>>266
内部ファイルを使うのが一番手っ取り早いだろう。

CHARACTER*10 FNAME
INTEGER I

DO I = 1, 20
　WRITE(FNAME,'(A, I2.2, ''.txt'')') 'xxxx', I
　OPEN(10, FILE=FNAME, STATUS='OLD')
　：
　ファイルの読み込み
　：
　CLOSE(10)
END DO

>>267
どうもありがとうございます。助かります。

申し訳ないですが、もうひとつ質問、
連番になっていなくても、ディレクトリにあるファイル全部を
読み込むような手段はあるでしょうか？

>>268
Fortran95 までの標準仕様ではディレクトリの操作や検索は出来ない。

しかし、大抵の処理系にはファイル操作関係の拡張手続きが付いているので
それらを使えば可能。
以下は Compaq Visual Fortran での例：

　use DFLIB
　implicit none
　character(len=100) :: dir, files, fname
　integer(kind=4) :: handle, fnlen
　type(FILE$INFO) :: info
　：
　dir = "D:\DataDir\"
　files = trim(dir) // "*.txt"
　handle = FILE$FIRST
　do
　　fnlen = GETFILEINFOQQ(files, info, handle)
　　if ((handle == FILE$LAST) .or. (handle == FILE$ERROR)) exit
　　fname = trim(dir) // info%NAME
　　open(unit=10, file=fname, status="OLD")
　　：
　　ファイルの読み込み
　　：
　　close(10)
　end do
　：

詳しくは処理系のマニュアル等を参照のこと。

>>269
どうもありがとうございました。
連番に直してから読み込むことにします。

format 文について、教えてください。
配列Aのサイズ変化にあわせて、
format を可変にしたい場合どうすればよいのでしょうか？

program main
integer, parameter :: NUM=4
real*8 A(NUM)
integer n

do n=1,NUM
A(n)=0.1*n
end do

! write (6,'(1x, NUMf5.1)') (A(n),n=1,NUM) ! これはエラー
write (6,'(1x, 4f5.1)') (A(n),n=1,NUM) ! これなら通る
end

write文で DO型並びを使わず愚直に DO文をかく。

write(6, '(1x, $)')
do n=1, NUM
write(6, '(f5.1, $)') A(n)
end do
write(6, '(/)')


FMT記述子を文字列変数にするという手もあるが、
こちらの方がシンプル。

>>272
ありがとうございます。無事にできました。
改行制御編集記述子と斜線編集記述子を組み合わせれば、
できるのですね！

仕様上はFORMAT書式に対応する項目がなくなった時点で書式制御が終了するので
出力する項目に対して多めに指定してやっても良い（ちょっと強引だが）
write(6, '(1X, 100F5.1)') (A(n),n=1,NUM)

処理系依存だが、<> でくくった中に整数式を書ける処理系もある。
write(6, '(1X, <NUM>F5.1)') (A(n),n=1,NUM)

ついでに、コロン編集記述子を使えば入出力項目がなくなった時点で
以降の書式を無視して終了させられる。
write(6, '(1X, 100(F5.1, :, ", "))') (A(n),n=1,NUM)

>>272
＄　指定子ってFortranの規格外でないかい？

>>274-275
情報、ありがとうございます。
過大定数の仮代入やコロン編集記述子は、アリですね！

<> でくくった中に整数式を書ける処理系
Compaq Visual Fortran ではできました！すごい！
処理系依存はできるだけ避けたいと思っているので、
考えどころですね。

改行制御編集記述子は規格外でしたか！
なかなか使い勝手のよさそうな機能なのに。

$の代わりにこっちでどうかな？
write(*,'(f5.1)',advance='no') A(n)

>>277
ありがとうございます。
こちらでもできました！
停留入出力という技もありましたね！

F90 の書式仕様の理解は、手を抜いていましたが、
まじめに取り組むべきだと思いました。

教えてください！！
初めて書き込みさせていただきます。
Ｆｏｒｔｒａｎのプログラムの課題が出ているのですがわからなくて・・・
簡単な内容なんだと思うのですが、
みなさんのお力貸していただけないでしょうか？？

内容は
配列データａ(i) (i=0,Num-1)　が与えられたとき、
このデータの平均値および標準偏差を計算するプログラムを書け。

というものです。
お願いします！！

>>279
平均値と標準偏差を求める式は知っているな？
それをそのまま書けばよい。

avg = sum(a(0:Num-1)) / Num
v = sum((avg - a(0:Num-1))**2) / Num
s = sqrt(v)

>>280
ありがとうございます！！
書いてみます！

学校の課題のこのプログラムが、全体として何が行われているか、どなたか教えてください！
お願いします。

real::x,x1,x2,func
real::dy,delt=1.e-6

read(*,*) x1,x2
if(func(x1)*func(x2) > 0.0) then
write(*,*) 'boo boo'
stop
end if
dy = sqrt(-func(x1)*func(x2))

do while (dy > delt)
x=(x1+x2)/2.0
if(func(x1)*func(x) > 0.0) then
x1=x
　else
x2=x
end if
dy = sqrt(-func(x1)*func(x2))
end do

x=x1
if(abs(func(x1)) > abs(func(x2))) x=x2
write(*,*) x
end

>>282
方程式の解、二分法
ttp://www.arch.kumamoto-u.ac.jp/hagane/yamanari/joho3/bisec.html

>>283
ありがとうございます！
ＵＲＬとっても参考になりました。
これでやっとレポート書けそうです。
ありがとうございます！！

しばらくぶりに・・・・・・宿題もっと持って来い！

こんにちわ。

02.dat,05,dat,08,datなどのように、
飛び飛びで01.datから10.datまでの間の番号のついたファイルがあります。
ただし、どの番号のファイルがあって、どの番号が無いのかは
時々刻々と変わります。
このファイルを、01から、10まで、順番にOPENしてREADしたいのですが、
このとき、存在しない番号のファイルをスキップするにはどうしたらよいでしょうか。
status='old'を試しましたが存在しない番号のところで止まってしまいます。
ヒントだけでもよいのでよろしくお願いします。

追記
status無しのただのOPEN文だと、存在しない番号ファイルを
空ファイルで作ってしまいます。これを作らないようにしたいのです。

>>286

logical :: exist
character(len=80) :: filename
do i=1,10
write(filename,'(i2.2,".dat")') i
inquire(file=filename,exist=exist)
if(exist) then
open(io,file=filename,status='old')
!
close(io)
endif
enddo

以下の複素多項式p(z)（zは複素変数）のゼロ点を全て、ニュートン法で求めよ。
(1)p(z)=z**2+4*i
(2)p(z)=z**4-2*z**3+2*z**2-2*z+1

お願いします。

PROGRAM newton
IMPLICIT NONE
INTEGER :: i, n
REAL :: pi
COMPLEX :: z
pi = 4.0 * ATAN(1.0)
n = 4 !2
DO i = 1, n
z = EXP( 2.0 * pi * CMPLX(0.0, REAL(i) / REAL(n)) )
DO
IF ( ABS(f(z))**2 < EPSILON(REAL(f(z)))**2 + EPSILON(IMAG(f(z)))**2 ) EXIT
z = z - f(z) / df(z)
END DO
PRINT *, 'finished', z, f(z)
END DO
STOP
CONTAINS
COMPLEX FUNCTION f(z)
IMPLICIT NONE
COMPLEX, INTENT(IN) :: z
f = z**4 - 2.0 * z**3 + 2.0 * z**2 - 2.0 * z + (1.0, 0.0) ! f = z**2 + (0.0, 4.0)
RETURN
END FUNCTION f
!----------------------------------
COMPLEX FUNCTION df(z)
IMPLICIT NONE
COMPLEX, INTENT(IN) :: z
df = 4.0 * z**3 - 6.0 * z**2 + 4.0 * z - (2.0, 0.0) ! df = 2.0 * z
RETURN
END FUNCTION df
!----------------------------------
END PROGRAM newton

解
finished (-1.414214,1.414214) (0.0000000E+00,1.3691417E-07)
finished (1.414214,-1.414214) (0.0000000E+00,1.3691417E-07)
Press any key to continue

finished (-7.1054274E-15,1.000000) (-2.8421709E-14,0.0000000E+00)
finished (0.9998764,-2.0247151E-11) (0.0000000E+00,1.0012824E-14)
finished (1.1924881E-08,-1.000000) (4.7699519E-08,0.0000000E+00)
finished (1.000000,1.7484555E-07) (0.0000000E+00,5.6843419E-14)
Press any key to continue


ふつう、ｎ次代数方程式の初期値は、１のｎ乗根から出発するのだが、
２番目の問題はそれがすでに解だったりして、newton法を使う前にLOOPから抜けてしまうｗ
問題設定が悪いｗ

Newton法では微係数が必要なのだが、ここでは手動で微分して関数として与えた。
発展形としては数値微分するもよし、多項式なら機械微分するもよし。

収束判定条件にはEPSILON関数を使ってみた。　あまり真面目に考えてない。

どうもありがとうございます。
学校でやってみます。

>>288様
ありがとうございます。
教えていただいたコードをそのままやってみたのですが、
existの中身がファイルのある無しに関わらず常に「F」になって
if文の中に入っていきません。原因として何が考えられるでしょうか。
使っているコンパイラはFTN95というやつです。
INQUIREというヒントをいただきましたので、自分でも考えてみます。
ありがとうございました。

できました。
ファイルの拡張子指定ミスでした。
大変失礼いたしました。orz

あげ

そろそろ　宿題持ってこないか？ｗ

質問です。
私のやりたいことは、
file_001.txt
file_002.txt
file_003.txt ・・・・
とナンバーの振られたファイルを順次読み込んでいくプログラムを作りたいのです。

do文などで何とかならないものでしょうか？

>>297
character*12 fname
・・・

do i=1,20
　write(fname, '(''file_'', I3.3, ''.txt'')') i
　open(unit=11, file=fname, status='OLD')
　・・・

　close(11)
end do

こんな感じでどうでしょう。

できました！ありがとうございます。

ここでwrite文はなにをしてるんでしょうか？
fnameの定義ですか？

>>299
内部ファイルを使って文字列変数 FNAME に
　　'file_' + 3桁の整数 + '.txt'
を書き込んでいる。
つまり、ループカウンタ I に応じて
ファイル名 'file_001.txt' 〜 'file_020.txt' を生成している。

内部ファイルについての詳細は↓やマニュアル等を参照のこと。
ttp://docs.sun.com/source/806-4843/Cp2_io.html

Fortran90で直行座標(x，y，z)を極座標(r，θ，Φ)に変換するプログラムを作れ。
って課題が出たのですがわかりません。教えて下さいm(._.)m

>>301
「直交座標 極座標」でググればすぐ出て来るだろう、というのは置いといて
なんのひねりもなく作ると、

program hoge
　implicit none
　real, parameter :: PI = 3.14159265
　real :: x, y, z, r, th, ph
　
　print *, "x, y, z:"
　read *, x, y, z
　r = sqrt(x**2 + y**2 + z**2)
　th = acos(z / r)
　ph = atan2(y, x)
　print *, "r, θ, φ"
　print "(3F12.3)", r, th*180.0/PI, ph*180.0/PI
end program hoge

はじめて書き込みます。
大学の先輩にFortranのソースプログラム(？)をもらって、それをインストールしたCOMPAQ VISUAL FORTRANで動かそうとしてます
でもコンパイルすると、
Can't open indirect file C:DOCUME~1
みたいなエラー文が出ます。これはどういう事なんでしょうか？

>>303
ソースのパスが空白や日本語を含まないようにして（例えば"D:\TEST\HOGE.FOR"）
再トライ

303です
304さんありがとうございました
エラーなくコンパイルし、実行できました

実行ファイルで条件を入力し、Enterを押すとすぐにウィンドウが消えてしまうのですが、どうしたらいいでしょうか？違うパソコンではそういう状態になりませんでした。

>>306
いくら何でも情報少なすぎ。
少なくともどんな環境で実行したか位は書いとかんと、まともの答えてもらえんぞ。

まあとりあえず、プログラムの最後に
read *
とか
pause
とか書いとけや。

学校のソフトを使ったので詳しい名前はわかりませんが、fujitsu fortran90
とかいう名前だったと思います。学校のパソコンでファイルを実行したとき
は、条件の値を入力したら答えが出て、答えが出た状態でキーボードに
触れるとプログラムウィンドウが閉じていたんですが、自分のパソコン
（NEC VALUESTAR Pentium4）で実行したときは条件の値を入力した瞬間
にプログラムウィンドウが閉じてしまいます。307さんのアドバイス通りに
最後にpauseと書いたらウィンドウが勝手に閉じることはなくなりましたが、
それでは根本的に解決していないのでまだ困ってます。

>>308
何がしたいのか今一よくわからないんですが…
使っているOSと自分のパソコンに入ってるコンパイラは何ですか？

学校のFujitsu Fortranの開発環境で実行して終了後一時停止していたプログラムを
自分のパソコンでexeファイルをダブルクリックして起動したときでもソースを変更せずに
同じように終了後一時停止するようにしたい。

というのなら>>307の様にして終了直前にキー入力待ちにして止めるか、
コマンドプロンプトを開いてその上で実行しないと無理ですよ。

初心者ですがフリーのをDLしてやってみたんですが、あるプログラム
を実行したら入力させる画面が出て、入力したら答えが出ないで、プログラム
が終了してしまいます。どうしたらいいでしょうか？ちなみに正常にコンパイル
されています。

>>310
>>309

>>309
＞コマンドプロンプトを開いてその上で実行

ありがとうございます、解決しました。

m組の変数xと関数f(x)の値を読み込み、x0に対するf(x0)の値をLagrange内挿する
関数副プログラムを作りなさい。また主プログラムを付けて、
次のデーターに対してf(0.525)を内挿してみなさい。

x f(x)
0.0 1.0
0.1 0.90484
0.2 0.81873
0.3 0.74082
0.4 0.67032
0.5 0.60653
0.6 0.54881

1次〜6次の内挿による解の精度も調べなさい。

参考

dimension XI(N),FI(N)
RLAGRN=0.0
do 20 I=1,N
F=FI(I)
do 10 J=1,N
if(J.ne.I) P=P*(X-XI(J))/(XI(I)-XI(J))
10 continue
20 RLAGRN=RLAGRN+P
return
end

御願いします

放物運動で時刻ｔ、角度θの任意の物体の位置のソースファイル？はあるのですが角度を３０度、４５度、６０度の３通りに
ついて計算し、すべての軌跡を出力せよ。との問題に全く手がでません。
ヒントとして配列を用いて３つの角度をプログラム内で定義する。とあるのですがよくわかりません。
どなたか教えて下さい。
お願いします。

>>314
元のソースがどんなのか知らんが、「配列を用いて」というところからすると
角度と配列を渡すとその配列に結果を格納するサブルーチンを作って
それを３回CALLすればいいんじゃないの？

すいません初心者なんでCALLとかわかんないです・・・
テキストには
dimension angle(3)
data angle/30.,45.,60./
・・・・・・・
do 20 i=1,3
theta=angle(i)*PI/180.
t=0.
10 continue
・・・・・・・
if(y.lt.0.) then
go to 20
else
go to 10
end if
20 continue
close(1)
end
と書いてあるんですが
・・・・・・
を考える課題なんですけど本当にわかりません。
教えて下さい。

>>316
配列って角度の方を配列にしてるのね・・・

まあ、プログラムを見る限り、最初の「・・・・・・・」の部分で
・データファイルをOPEN文で開く
次の「・・・・・・・」の部分で
・角度theta、時間t、などから位置x、yを求める

ということなのだろうけど
・初速度や初期位置は既知なのか、自分で決めるのか、ファイルから読み込むのか
・データファイルの書式は？
など不明な点があるのでこれ以上は自分でおやりなさい。

いずれにせよ情報を小出しにされたのではｴｽﾊﾟｰでもない限り「本当にわかりません」。

1からＮまでの整数について、次のそれぞれの分類に属する整数の数を数えて表示するプログラムを作成せよ。
・３の倍数であり、５の倍数でもあるもの
・３の倍数であるが、５の倍数ではないもの
・３の倍数でも５の倍数でもないもの



誰か教えてください。

プログラムというか数字の扱いの勉強したら？
高校数学あたり。

>>319
ここはどんな質問にもマジレスする擦れです。
宿題も教えまくりますｗ

>>318　FORT77で書いてみたｗ　英語はうんこｗ

PROGRAM vip
WRITE(*, *) ' INPUT N'
READ (*, *) N
N1 = 0
N2 = 0
N3 = 0
DO I = 1, N
M3 = MOD(I, 3)
M5 = MOD(I, 5)
IF (M3 .EQ. 0 .AND. M5 .EQ. 0) N1 = N1 + 1
IF (M3 .EQ. 0 .AND. M5 .NE. 0) N2 = N2 + 1
IF (M3 .NE. 0 .AND. M5 .NE. 0) N3 = N3 + 1
END DO
WRITE (*, *) ' The number of multples of both 3 and of 5', N1
WRITE (*, *) ' The number of multples of 3 but not of 5 ', N2
WRITE (*, *) ' The number of multples not of 3 nor of 5 ', N3
STOP
END

>>320　専ブラ用インデント

>>314,316
出された情報だけでは解けない。
初速度が必要だが、それがファイルに入っているのだろうか？
重力常数も必要だがどうするか？（まぁ知ってはいるが、問題が地球とは限らないかもｗｗ）

もうちょっと情報出せ。
テキストプログラムを見る限り出題者はセンス悪いｗ　
このスレで素直な良い手筋を学べ

>>313　F90でエラーチェックなしウンコプログラムｗ
PROGRAM Lagrange
IMPLICIT NONE
INTEGER, PARAMETER :: n = 7
INTEGER :: i
REAL :: xin(n), yin(n), x, y
OPEN (10, FILE = 'input.dat')
DO i = 1, n
READ(10, *) xin(i), yin(i)
END DO
x = 0.525
y = alagran(x, 6, xin, yin)
WRITE(*, *) ' x, f(x) =', x, y
STOP
CONTAINS
REAL FUNCTION alagran(x, n, xin, yin)
IMPLICIT NONE
REAL , INTENT(IN) :: x, xin(:), yin(:)
INTEGER, INTENT(IN) :: n
REAL :: tmp
INTEGER :: i, j
alagran = 0.0
DO i = 1, n
tmp = yin(i)
DO j = 1, n
IF (j /= i) tmp = tmp * (x - xin(j)) / (xin(i) - xin(j))
END DO
alagran = alagran + tmp
END DO
END FUNCTION alagran
END PROGRAM Lagrange

>>323　ごめw　まちがい　このままだと５次ｗ

y = alagran(x, 6, xin, yin)
↓
y = alagran(x, n, xin, yin)

でｗ

とりあえず、答えは
x, f(x) = 0.5250000 0.591555
のはずｗ

間違ってても知らんｗ
ほぼ直線に乗っているのでまぁこんな感じ？

ラグランジュ法はデータの差分の漸化式で計算するほうが一般的なだが
とりあえず参照として出されたサブルーチンと同じく公式どおりに書いた。

１〜６次までの誤差の見積もりのほうは、課題の雰囲気からすれば、まじめな誤差解析は要求して
いないだろうから、０．６と０．５・・・・０．０まで逆進する方向でデータを増やしてサブルーチンを呼んだ
結果を書いて何事かムジャムジャ言えばよかろうｗ

>>325

DO i = 1, 6
! y = alagran(x, n, xin, yin)
y = alagran(x, i + 1, xin(7 - i:7), yin(7 - i:7))
WRITE(*, *) ' i, x, f(x) =', i, x, y
END DO

i, x, f(x) = 1 0.5250000 0.5921000
i, x, f(x) = 2 0.5250000 0.5915309
i, x, f(x) = 3 0.5250000 0.5915560
i, x, f(x) = 4 0.5250000 0.5915547
i, x, f(x) = 5 0.5250000 0.5915551
i, x, f(x) = 6 0.5250000 0.5915558

こんな感じか？　
誤差が０．１＾Nに比例しているとか書けばいいのだろうか？ｗ

ttp://www.library.cornell.edu/nr/bookfpdf/f3-1.pdf
この辺でも読みやがれですぅ

n個の実数データを配列datに読み込み、読み込んだデータを降順に出力する。
ただし、データを読み込んだ配列の中身は変更せず、別に整数型の配列indexを
用意して、大きい順の並びが
dat(index(1)),dat(index(2)),･･･,dat(index(n))
であたえられるようにする。つまり次のような出力命令
do 30 i=1,n
write(*,*)dat(index(i))
30 continue
によって、データが降順に出力されるようにせよ。

という課題がでたのですが、fortran77の形で教えていただけませんか？
よろしくお願いします<m(__)m>

ＮＡＧのFortranBuilderの試用版が無料セミナーで配布されるらしい
以下参照
【日時】　5月31日（水）13時30分〜（13時開場）
【場所】　英国大使館　（東京都千代田区一番町1）
【地図】　http://www.nag-j.co.jp/fbreleasemap.htm
【内容】　
１：Fortran Builderのご紹介(Fortran Builder開発者　中野 剛）
Fortran Builderの機能、製品コンセプト、事例ご紹介.デモンストレーションなど
２：プログラミング学習とe-Learning
　（法政大学情報メディア教育研究センター所長　竹内　則雄 教授）
プログラミング言語学習用電子ブックと開発環境が連動した体験的なFortraｎ
プログラミング学習環境の概要と、対面授業やe-Learningでのその効用について
３： NAGWare f95 and reliable, portable programming
　（英国ＮＡＧ NAGWare Fortranコンパイラ開発者、マルコム・コーエン）
How to detect errors using NAGWare f95,and how to write portable,
reliable programs.Support for features from the latest Fortran standard
and plans for future releases. 　　
４： Fortran Builder を使った数値計算 (筑波大学 長谷川 秀彦 助教授）
行列計算パッケージ LAPACK の紹介と Fortran Builder を利用したLAPACK の手軽な使い方について
５：Fortran Builder 体験コーナー・30日トライアル試用版配布・
ノベリティ配布など もちろん無料です。
英国大使館の警備都合上事前にご登録をお願いしております。
下記ＵＲＬフォームにて参加登録ができます。

>>326
特に指定がなかったので解りやすく(?)バブルソートでやってみますた。
　　　parameter (n=10)
　　　real dat(n)
　　　integer index(n)
C
　　　do 10 i=1,n
　　　　read(*,*) dat(i)
　　　　index(i) = i
　 10 continue
C
　　　do 22 i=1,n-1
　　　　do 21 j=1,n-i
　　　　　if (dat(index(j)) .lt. dat(index(j+1))) then
　　　　　　ntmp = index(j)
　　　　　　index(j) = index(j+1)
　　　　　　index(j+1) = ntmp
　　　　　end if
　 21　 continue
　 22 continue
C
　　　do 30 i=1,n
　　　　write(*,*) dat(index(i))
　 30 continue
　　　stop
　　　end

処理系によっては配列の宣言部分でエラーになるかも知れないので
その時は適時修正のこと。
#F77形式で、って言われるとここが一番難しいかもｗ

>>328選択法でプログラミングしてもらえませんか？

>>329
オレが口をだすべきではないと思うが、苦言を一言。
モノの頼み方を知ってるのか？

このスレはフォートランの超初心者向けではあるが、
人生の超初心者（不躾な厨房）向けではないぞ。

質問をするなら、回答者が答えやすいように、できるだけ詳しく書く。
最も簡単な方法は、（要約などせず）出題内容をすべて書く。
直接選択法の解が必要なら、最初から書け。

労を厭わず、回答案を出してもらえたら、
最初の一言は「ありがとうございます」だ。
次に、再要望であったり、回答内容に対する質問だ。

>>330
申し訳ありません。
バブルソートでの解に関しましては、御礼申し上げます。
もし、気分を害されてないようでしたら、直接選択法の解をお願いできませんか？

>>330
まぁそうなんだけどさ，ここは２ｃｈだし，いいんでないのｗ

Ｆｏｒｔｒａｎ人口を増やすためにも甘い水で誘うべしｗ

>>331
オレじゃなくて、>>328 に　感謝して、非礼を詫びろ。

>>332 たしかにそうとも思うが、
そういう教育はだれかがしなければ、
会社に入って恥をかくまで治らないと思う。
社会人になっても治らないものもいるが。orz

Fortran 人口を増やすには、説教より回答が一番だな。
328 えらい。オレばか。orz

>>331
ほれ、>>328の真ん中のループの所を入れ替えれ。

do 22 i = 1, n-1
imax = i
do 21 j = i+1, n
if (dat(index(imax)) .lt. dat(index(j))) imax = j
21 continue
ntmp = index(i)
index(i) = index(imax)
index(imax) = ntmp
22 continue

君の行動を称えて以下の文章を贈ろう。
コレをよく読んで今後も精進に励め。（どう励むかはお前次第だ）
http://www.redout.net/data/osietekun.html

>>333
まあそう気にしなさるな。
ここに宿題を丸投げして後で苦労するのも、社会に出てからちんぷんかんぷんな
質問をして恥をかくのも、全て本人の自己責任ってやつだ。

　 WRITE(*,'(A5$)')name(i)
のような書式で、表示位置を決める定数"5"の部分をREAD文で読み取って
任意の値に変更できるようにしたいのですが､
どうすればいいでしょうか？

ちなみに
　 num=5
　1000 FORMAT(Anum$)
　 WRITE(*,1000)name(i)
と書いたのですがうまくいきませんでした。

>>335
1. 上の方でも出てきているが、<>の中に変数名（整数式）を書く（ただし、一部の処理系の拡張形式）
　　　write(*,'(A<num>$)') name(i)

2. 文字列の書式を書き込む
　　　character*10 str
　　　・・・
　　　write(str, '(''(A'',I3,''$)'')') num
　　　write(*,str) name(i)

>>336
内部ファイルは便利なのに知られていないよう棚。もっと啓蒙が必要なようだ。
FAQにでも入れておきたいところだｗ

>>335
■必要な第1の知識
FORTMATは
１．FORMAT文で与える。
２．文字常数として与える。
３．文字変数として与える。
という、３通りの方法がある。

たとえば
.....1.FORMAT(A13)
.......WRITE(6, 1) 'WE ARE VIPPER'
は
.......WRITE(6, '(A13)') 'WE ARE VIPPER'
と等価だし
......CHARACTER*10 fmt
......fmt = '(A13)'
......WRITE(6, fmt) 'WE ARE VIPPER'
とも等価である。

■必要な第２の知識
内部ファイルというものがある。これは文字列がファイルのように読み書きできるということを意味する。

WRITE(9, '（A)' ) ' chin chin man man ' ! これはファイルの９番に文字列を書き出している。
WRITE(fmt, '(A)') ' okki ' 　　　　　　　　　! これは文字列fmtに文字列を書き出している。

ちなみにREAD文でも同様のことができる。
これを利用すれば、好みのｆｍｔ文を実行時に動的に生成できる。
文字列操作を使ってもいいが、この方が数値を文字列として書き込むのが楽勝でできる。

無料のWindows用のコンパイラはありませんか？

>>338
g77（gccに入ってる）
Salford f77

できました！短時間に二人もありがとうございます！
>337さんの
WRITE(9, '（A)' ) ' chin chin man man '
が顔に見えてちょっとうれしいです。

>338
g95
ttp://www.g95.org/

>>334
ありがとうございます。
よい勉強になりました。

fortranの遺産があるんですが・・・
なぜか配列の出力が20行ごとｗ
むかしの名残なんでしょうか？

ってか計算だけfortranでやらせてGUI？はほかの言語でやるにはどうしたらいいの？
ｵｼｴﾃ

あとゎぁ書式の指定で
format(500i4)
ここの500のところを定数とかで指定できないの？

>>344
最初から読めとまでは言わないけれど、
10前のレスくらい読もうよ。

>>335-336
答えそのものではないけれど、ほぼ回答。

>>337
少し難しいけれど、極めて有益な情報。

ｻﾝｸｽｗｗｗ
くそふるいマニュアル本みててもわからないわけだ・・・ｗ

>>343
大体はＦＯＲＴＲＡＮ側をＤＬＬにしてＶＢやＤｅｌｐｈｉなんかで処理しているみたいだな．

まぁ，ＧＵＩもＦＯＲＴＲＡＮでやってる連中もそれなりにいるようだ．

Ｉｎｔｅｌの質問掲示板に行けばうまくいかねーとウンウンうなっている人たちに出会えるｗ
ttp://softwareforums.intel.com/ids/board?board.id=5

FTN95でコンパイルしたんですが，
オブジェクトファイルができただけで，
この後どうやって実行すればいいのか教えてください．

>>348
1. コンパイル時に「/LINK」オプションを付ける。
例： ftn95 hoge.f90 /link

2. コンパイルが正常終了したら、おもむろにプログラム名を入力してEnterを押す。
例： hoge

3. スプラッシュが消えるまで全裸正座で待機。
例：
　+　　　+
　　 ∧＿∧ 　+
　　（0ﾟ・∀・）　　　ﾜｸﾜｸﾃｶﾃｶ
　　（0ﾟ∪ ∪ + 　　　　　　　
　　と＿_）__）　+

>>349
おお！できました．ありがとうございます。
ちょうど風呂からでたところなので助かりましたｗ

　　　　　　（＿＿＿＿＿＿＿＿＿＿＿＿＿＿　J
　そ　君　 （__　　　　　　　　 　┌―‐―┐　　　　）　き
　ん　！　　(_　　　　 　 　 　 　 |`l TT了|　　　 　｝　 ゃ
　な　　　　 (_　　　　　　　　　　 j .| .|:| .l |　　　　 /　 あ
　か　待　 （_　　　　　　　　　 　 | | .|j .j |　　　　 ｲ　　|
.　っ　ち　　(,_　　　　　　 　 　 　 |.| .|l .|:|　 ,. -‐ /_.　 っ
　こ　 た　 （　　　　 　 　 　 　 　 |.|. |! |／　　 　　/　！
　う　 ま　　(｀　　　　　　　　　　　｀ｰ /..:::::＼≧,,,､:::７___
　で　え　　（―――――――――（:::::::＞'´ ==　＼::⌒l＾⌒
_　　　　　(⌒￣￣￣￣￣￣￣￣￣ ﾉ く彡///　∪,ﾉ　　　;|
.ﾚ⌒Y^'⌒｀＼＿＿＿＿＿＿＿＿ く:::::∧　'_,. -､ く／::::::::/
　　 |:::|　＼xく　　　　_,,,...,_　　　　　 　＼:::::l､ヽ ,ノ　 ＼,,∠,,__
　＼|:::| _,....!,,_ ＼　 iれ__,.､ヽ　 　 　 ｌF〒｀ヾ.＼,,..イ　 　 |::::::::,
　　 ｀7´ ＿,,.ｨ　　ヽ{|iｭ　ｪｯﾘ　 　 　 |　||　　_,..-／7ﾞｈ　_|:::::://
　＼.｛ｎ|.ｨｧ it}　　　 ',　_'_　ｊ)　　 　r'"三¨7´＼|　　　 |´.|:::://
　　 |:::ﾄl、 rｭｊ　.　 　 ト ﾆ ｲl、　　/ ﾟ`.|ｎ./　 .イl　 　,∧ |:://
　　 |::,|　 'ーｹﾄr'TTlｲ　　/＿｀ヽﾞtっ　r'lﾞ　 　 /⌒`lくミV /
　,r1´|｀'六´ //｀￣´ ｀Y´ 　 　 |└┬ｼｊ　 .／ 7ヽ〈 　/ヾ)＜
./　|　∨|::|∨ !　｛ r　　,､　＿,ｼ ／ﾞ丁〈　/　 　 　 }　{　｛ ＼
　　| 　 ',|::|/　!　 ,ゝ-＜ (　　　/ 　　.|　|/　　　　 ∧　＼|
　　 l　　.Ｙ。 .|　 |`　 〃￣￣⌒　　/　〈　　　　 /!　', __,,....::-‐
　　.∧.　 |。 ｛　　ﾞ爪｀ ' ‐- ､..,,,...イ　　　'、 　 / .|　 ｀|::::::::::::::::
＼/　 ｌ　 |。./　　,l | l,　　.|　　.　 ||　　　　｀'ｰ' i　|　　j:::::::::::::::::
ヽ､｀'::､L.∧/　　/　|.｛　　u　　　〈.|　 　 　 　 イ 〈　 /::::/:::::::::::
::::::::｀ヽ､ ∨　　/￣|　| ､　　　/　 l:l.　　　　 　 |　j　/::::/:::::::::::::
:::::＼::::..... ￣｀|　　 ｌ　|　 鬱鬱　　.l:l/／　　　 l|　∨:::/::::::::::::::::

いいかお前らにためになるお言葉を言い渡しておこう。
「問題はレベルの高いプログラム言語を使うか、レベルの低いプログラム言語をつかうか
　といったことではない、そうではなくて、
　レベルの高いプログラム言語を使うか、もしくはレベルの高いプログラム言語処理系を低レベル言語で作成するか、
　ということだ。」

つまりフォートランなんかいまさら使っている時点で負け組み。
それで半ばコンピュータの専門か気取ってる痛いやつが多すぎなんだよｗ

日本語でおｋ

さて、諸君。次の問題を持ってきたまへ。

RoR じゃなくて FoR はありませんか

駐車料金を計算するプログラム
料金は30分単位で切り上げて計算される
最初の一時間まで３００円/３０分
一時間を越えると２００円/３０分
入力は実数型　2,35（2時間35分）
この実数を引数で関数の値が料金のFunctionをつくり使用する

↑のような宿題が出たのですがFortran77の型で誰かお願いします

>>356
program parking
read(*,*) hour,fmin
write(*,*) ifee(hour,fmin)
stop
end

function ifee(hour,fmin)
totalmin=hour*60.+fmin
if(totalmin.le.60.) then
ifee=300*int(totalmin/30)
if(totalmin-30.*int(totalmin)/30.gt.0.) ifee=ifee+300
else
ifee=600+200*int(totalmin-60.)/30
if(totalmin-30.*int(totalmin-60.)/30.gt.0.) ifee=ifee+200
endif
return
end

>>357
ありがとうございます

あ、忘れてたのですが最高課税額が4000円になるようにしてもらえませんか
本当に申し訳ない＞＜

>>359
357 じゃないけど。
せっかくプログラムがあるんだから、
いろいろいじってみようよ。
ヒント：if 文　をどこかに一行いれる。

慌てて作ったら間違っていたみたい。
それを修正するのも良い勉強になるでしょう。
あと、変数に7文字以上使っているのも問題かも。
実は厳密なFORTRAN77を知らない。

>>361
それ以前に小文字を使っているぞいｗ

そういやそうだったねorz。字下げもしない方が良かったかな。

1/√π×exp（-ｘ^2）は中央値0、偏差1の正規分布曲線を与える。
〔-1、1〕および〔-2、2〕の範囲に存在する確率はいくらか。台形法で求めよ。

学校の課題なのですがわかる方いたらお願いします。
先生いきなり難しいよ・・・

小学生の算数じゃなかったか？＞台形面積

>>363
いや字下げは文法にかなっているからおｋ
FORTRAN７７までは大文字のみが文法で規定されていた。
小文字は見にくくてｷﾓｽｗ

>>365
おまえは>>1を読み返せ！

このスレは糞既出FAQにも丁寧にマジレスし、宿題も文句も言わずに解いてやり、
FORTRAN依存症を作り出して初心者をFORTRAN Userに引き込む秘密結社の陰謀だ。

まぁ本スレの方ならなじった上で、こっちのスレに誘導してくれればおｋｗｗｗ

>>364　Fortran90で難しく書いてみた。が、少しイモくさいｗ

MODULE mod_integ
IMPLICIT NONE
REAL :: pi
CONTAINS
REAL FUNCTION trapez(n, xmin, xmax)
IMPLICIT NONE
INTEGER, INTENT(IN) :: n
REAL , INTENT(IN) :: xmin, xmax
INTEGER :: i
REAL :: h, x
h = (xmax - xmin) / n
trapez = 0.0
DO i = 1, n
x = xmin + (i - 1) * h
trapez = trapez + sub_trapez(x, h)
END DO
RETURN
END FUNCTION trapez

REAL FUNCTION sub_trapez(x, h)
IMPLICIT NONE
REAL, INTENT(IN) :: x, h
sub_trapez = 0.5 * h * ( func(x) + func(x + h) )
RETURN
CONTAINS
REAL FUNCTION func(x)
IMPLICIT NONE
REAL, INTENT(IN) :: x
func = EXP( -x**2 ) / SQRT(pi)
RETURN
END FUNCTION func
END FUNCTION sub_trapez
END MODULE mod_integ
!
PROGRAM oppai
USE mod_integ
IMPLICIT NONE
INTEGER :: i
REAL :: s, s1, s2
pi = 4.0 * ATAN(1.0) ! arctan(1.0) = pi / 4; tan(pi/4) = 1.0
! [-1..1]
s1 = trapez(100, -1.0, 0.0)
! [-2..2]
s2 = trapez(200, -2.0, 2.0)
! [-inf..+inf]
s = trapez(1000, -10.0, 10.0)
PRINT *, s, s1, s2
STOP
END PROGRAM oppai
>>367-368

実行結果
0.9999998 0.4213469 0.9953194
Press any key to continue

微妙に問題の偏差の定義などが一般的な統計の場合とずれているが、
ここでの定義を信じてゆくことにするｗ

正規分布の−∞〜＋∞の積分は１になるように規格化してあるはず。
（√π）で割っている部分。それを確かめたのが変数ｓ。
ただし積分範囲は−１０〜＋１０までにしている。
無限大じゃないじゃんか！というかもしれないが、これでも有効数字の範囲では
もう無限大と変わりない。というか、関数は０になっている。
誤差は台形公式の誤差だろう。分割数ｎ（分割の間隔）を変えると調べられる。

-1..1の積分は約０．４２１、-2..2までの積分は約０．９９５になっている。
つまり、−∞から＋∞まで積分すると1.0になるはずものだが、
−２〜＋２までの範囲で約９９．５％が尽くされることが分かる。

プログラムでは、メインルーチンは積分範囲を指定して積分ルーチンを呼んでいる。
サブルーチンはモジュールの中に入れてある。
台形（ｔｒａｐｅｚｏｉｄ）積分は、等間隔にｘ軸を切ってそれを台形の高さとし、それぞれの
小さな台形の面積の和を返している。小さな台形の面先はさらに台形面積を返す
ルーチンを呼び出している。小さな台形の面積を返すルーチンは正規分布関数の
値を返すルーチンを呼んでいる。

・プログラム
　・積分ルーチン
　　　・台形面積ルーチン
　　　　　・正規関数ルーチン

FORTRAN77で書いてほしければ、はじめから遠慮せずに言うようにｗ

なお、間違っていても謝罪も賠償もしないので４６４９！

ありがとうございます！77とか90とかよくわかりません。学校のパソコンはどっちだろう・・・
これをテラパッドに入力してシグインで実行するわけですな

>>371
! [-1..1]
s1 = trapez(100, -1.0, 0.0)
これまちがいｗ　ごめｗ　でも賠償は（ｒｙ
正しくは、s1 = trapez(100, -1.0, 1.0)　面積は０．８４２６７２９

積分範囲がまちがって-1.0〜0.0になってたｗｗｗｗｗｗｗ
なんか値がおかしいと思ったのだがｗ　標準偏差の定義が標準と異なるので
表で値を確かめるときに換算を間違えておｋと思ってしまったｗｗｗ
ま、宴会から帰ってきた夜中だったので許せｗ

＊やや無関係な話題ｗ
πの値はパラメータで与えてもいいが、ここでは４．０＊ＡＴＡＮ（１．０）を使った。
これは昔からよく使われている定石だが、どこかで１回計算しておかないと
いけない。それがめんどい。だがこうしておくと精度を機械任せにできる。
今のプログラムではやや不自然で、結果醜くなっている。

Ｆｏｒｔｒａｎ２００３ならパラメータ文の中で関数が使えるようになるのだが・・・・
（パラメータは実行時ではなくコンパイル時に評価されるので、
今のところ四則演算とベキぐらいしか使えない）

＊正規分布の−∞から＋∞への積分
∫exp(-x^2)dxは√πになるのだが、不定積分があらわせないので、
必死に積分しようとしないようにｗ

多重積分のテクニックで直交座標を円座標にして出すのが伝統芸。
∫exp(-x^2)dx∫exp(-y^2)dy=∬exp(-(x^2+y^2))dxdy=∬exp(-r^2)ｒdrdθ
=-1/2[exp(-r^2)]^∞_0 * 2π　=　1/2*2π = π
つまり∫exp(-x^2)dx=√πだお。　

数値積分では無限大の積分は、普通変数変換して有限範囲にしてから実行する。
ここでの-10~10の積分は、あくまで確認のためのお遊びということで。

>>372
ありがとう。使わせてもらいます。

DAIKEIHOU
gakusekibangou NAME
sekibun.f conded at 29 May
DIMENSION NJUN(J)
DATA NJUN/10,50,100/
A=-1,0
B=1,0
WRITE(6,100)A,B
FORMAT(5x,'SEKIBUN FOR 1/SQR(PI)*EXP(-X*X) A'
,F5,2,3X,'B',F5,2)
DO 10 K=1,3
N=NJUN(K)
H=(B-A)/N
S=FX(A)+FX(B)/2,0
DO 20 J* J=1,N-1
X=A+BJ
S*S+FX(X)

ちなみにこれが学校の先生が見本としてやってたのを急いで写したものですが（途中まで）
あってますか？

微妙に写し間違えているｗ
DIMENSION NJUN(3)
...
S=( FX(A) + FX(B) ) / 2.0
DO 20 J = 1, N - 1
X = A + H * J
S = S + FX(X)
20 CONTINUE
S = S * H
WRITE(6, *) 'BUNKATSU = ', NJUN(I), 'MENSEKI =', S
10 CONTINUE
STOP
END
こんな感じになるのではないかな．
台形公式を>>364の様に素直に小さな台形の和を足して求めると，
関数の値を２回同じ値に対して求めなければならない．

そこを改良すると，ここでの例題のように，両端だけを半分にして，
残りを一回計算して足し上げてゆけば良いことになる．
絵を描いて考えればすぐ分かると思う．

NJUNは分割数だろう．何個の台形に分けるかを指定している．
数学的には分割が細かいほど正確な値になるが（積分は分割→０の極限で定義
されている）現実には数値誤差や計算時間の関係からほどほどの分割数がいい．
いくつがいいかは時と場合によるので一概には言えない．

台形公式の改良版としてロンバーグ加速ってのがある。
台形公式の幅を半分ずつにしていって結果がどのように変化するかを調べて、
幅無限小の時の値を補外するというもの。
関数の形にもよるけど、台形公式とほぼ同じ処理量で凄い精度が出る。

>>50

ポインタが無いぶん最適化はしやすいと聞くが、
正直微々たる差だろうな。

GCCの場合、どちらも一旦中間コードに直されてから最適化されるので、同じ構造のプログラムを書いた場合
基本的に最適化に差はほとんどない。

>>378
ＧＣＣに限らず最近のコンパイラは、中間言語を出力する字句解析のフロントエンド部と、
機械語変換のバックエンド部に分離しているから、その理屈ではどんな言語でも
似たようなスピードになる。実際そういう傾向はある。

問題は今後のハードの発展が、ＣＰＵの並列化や、非連続な分散メモリー、
ＳＳＥのようなベクトル処理の方向に向かっていることで、言語自体が向いていないと
最適化がかかりにくくなってくる。

ＦＯＲＴＲＡＮの規格はハードウェアの進化のフィードバックを改定を受けているが、
Ｃ言語などは最近のネットセキュリティーの重視を受けてセキュアーさのフィードバックを
強く受けているようなので長期的にはどうなるか・・・

最近のマルチコア・マルチスレッド・計算用専用プロセッサの利用など（ＧＰＵ流用や
ＣＥＬＬとか）を見ていると、コンパイラによる自動並列化というより、ＭＰＩのような
ＡＰＩ呼び出しで解決しようとしているようだし、どうなることかｗ

すみませんフォートランの始め方からしてわかりません・・・

http://www.g.dendai.ac.jp/lecture/fortran/fort01.html
このサイトにはFTN77 Personal Editionをダウンロードしろとあるのですが、
77でなく95しか置いてないみたいなんです。

Cpadが77向けなので95と77では合わずに不具合が生じてしまう気がするのですが・・・

どなたかご回答願います・・・

>>380
FORTRAN III
http://pc8.2ch.net/test/read.cgi/tech/1104724162/

861 名前：デフォルトの名無しさん[] 投稿日：2006/05/27(土) 21:37:28
FTN77 Personal Edition ってもうダウンロードできないの？
SalfordのHPやSilverfrostのHPからもリンクが消えてるんだけど…

862 名前：デフォルトの名無しさん[] 投稿日：2006/05/29(月) 01:52:14
>>861
FTN95がＦｒｅｅになったから消したんでないのかな？

863 名前：デフォルトの名無しさん[] 投稿日：2006/06/03(土) 02:31:53
もうどこでもダウンロードできないんですかね?

864 名前：デフォルトの名無しさん[] 投稿日：2006/06/04(日) 00:35:52
>>863
ロシアか中国あたりにうｐられていそうな気もするがｗ
漏れは使ったこと無いけど、FTN77はそんなにいいものなのかい？
フリーのFTN95はバナーがでてウザかった気はするが使えたよ。
F77だったらOPEN WATFORやｇ７７などがあるので代替手段が無いわけではない。
WATFORもDOS EXTENDER時代に本来の開発が止まったものなので、昔懐かしい
感じを保っているはずだが・・・ｗ

865 名前：デフォルトの名無しさん[sage] 投稿日：2006/06/04(日) 01:38:40
「FTN77PE.EXE」でぐぐるとそれらしきものがいくつか見つかるが
どれもなんか怪しい雰囲気がｗ
>>864
単体で使う分には取り立ててよいというほどではないと思うけど
CPad for Salford FTN77と一緒に使うとお手軽にプログラムの
作成・実行が出来るので重宝しています。
このセットで紹介しているHPも結構あるみたいだし。

およよ・・・ありがとうございますです。
がんばってみます。

すでにめげそう；；

>>380
文法的にはFTN95はFTN77を包括している。
（完全に規約に準拠していると、７７の機能の一部が削られていることになるが）
Cpadの設定を自分で書き直せば何とかなるだろう。

まぁ、そろそろ９０、９５に移行することをお勧めする。
といっても色々事情もあるだろうから、好きになさいｗｗｗ

安易な方法としては、違法くさいのがうｐされているので、具具ってそれを落としてくれば
とりあえずは足りるべｗ　　
いつまであるかはわからんがｗ

みなさまありがとうございます。
がんばりますです。

95と77の違いってどんなところ？？

>>385
９５はモジュール言語の特徴を大きく取り入れた。
ほとんど温泉旅館の新館と旧館のようになかば違う言語がくっついた感じになっている。
Fortran90の時点ではFORTRAN77を完全に包括していた。
Fortran95はFortran90のBugFix的な位置づけだが、一部Fortran90で廃止勧告が
出ていた機能を削っている。ただし、ほとんどのベンダーは過去互換性を重視して
そういう機能を残している。

古い人なら、FORTRAN77がpascalと合体したと思えばいい。

どもっす

行列がわかんねーorz

>>388
それは FORTRAN の問題じゃなくて数学の問題ではなかろうか。

DAIKEIHOU
gakusekibangou NAME
sekibun.f conded at 29 May
DIMENSION NJUN(J)
DATA NJUN/10,50,100/
A=-1,0
B=1,0
WRITE(6,100)A,B
FORMAT(5x,'SEKIBUN FOR 1/SQR(PI)*EXP(-X*X) A'
,F5,2,3X,'B',F5,2)
DO 10 K=1,3
N=NJUN(K)
H=(B-A)/N
S=FX(A)+FX(B)/2,0
DO 20 J* J=1,N-1
X=A+BJ
S*S+FX(X)
DIMENSION NJUN(3)
...
S=( FX(A) + FX(B) ) / 2.0
DO 20 J = 1, N - 1
X = A + H * J
S = S + FX(X)
20 CONTINUE
S = S * H
WRITE(6, *) 'BUNKATSU = ', NJUN(I), 'MENSEKI =', S
10 CONTINUE
STOP
END

これ実行できないです。どこか間違ってますか？

>>390
コメントは行の先頭に"C"を付けろとか
小数点(".")とカンマ(",")を混同してるとか
文番号のないFORMATなんて何に使うんだとか
勝手に行を折り返すなとか
実行文の途中に配列宣言入れるなとか
ツッコミどころ多数でどこもかしこも間違っている。

>>390　もう少し自力でがんばれｗ
PROGRAM DAIKEIHOU
C gakusekibangou NAME
C sekibun.f conded at 29 May
DIMENSION NJUN(3)
DATA NJUN/10,50,100/
A = -1.0
B = 1.0
WRITE(6,100) A,B
100 FORMAT(5X,'SEKIBUN FOR 1/SQR(PI)*EXP(-X*X) A', F5.2, 3X, 'B',F5.2)
DO 10 K = 1, 3
N = NJUN(K)
H = (B - A) / N
S = ( FX(A) + FX(B) ) / 2.0
DO 20 J = 1, N - 1
X = A + H * J
S = S + FX(X)
20 CONTINUE
S = S * H
WRITE(6, *) 'BUNKATSU = ', NJUN(K), 'MENSEKI =', S
10 CONTINUE
STOP
END
C
FUNCTION FX(X)
PARAMETER (PI = 3.14159265)
FX = EXP(- X * X) / SQRT(PI)
RETURN
END
>>392

小数点をカンマにするのは、ドイツ・フランス。中点にするのがイギリス。
あと日本の縦書きの新聞とかも中点にする。

しか〜しｗ　FORTRANやアメリカ・日本では小数点はピリオドだｗｗｗ

FORMAT文はWRITE文の中に書き込んだほうがいいのだが、
学校教育では、まぁ文番号つきのFORMAT文でやっとけｗ

関数副プログラムを使ったが、今の場合FX(X)は文関数でも書ける。

質問があれば答えるから、まぁなんちゅーの、もう少し理解するようにしろｗｗｗｗ

>>389
行列の表示の仕方がわからない・・・

>>394
real A(2,3)
data A/1.,2., 3.,4., 5., 6./
do i = 1, 2
write(*,*) (A(i,j), j=1,3)
end do
end
とか
program hoge
implicit none
real :: A(2,3) = reshape((/1.,2., 3.,4., 5.,6./), (/2,3/))
integer :: i
do i = 1, 2
print *, A(i,:)
end do
end program hoge
とか

関係ないがFortran 2003では (/.../) のかわりに [...] と書けるので
可読性が少し良くなる（のかな？）

>>395
解決しました！ありがとうございます。

>>394
もっと具体的に質問すれば、華麗なるFORMATを教えて進ぜようゾナｗ

>>392
ありがとうございます。家のパソコンにはシグインないので学校行ったらチャレンジしてみます。

>>398
シグインってcygwinのことかな？ｗ

念のため言っておくけどｗ　変数のAとBが積分範囲の下限と上限なので、
−１，１というのを−２，２にすれば当初の課題が達成されるｗ

シグインなんてDLすればいいジャンジャン

Evca(Ω,x)を求めるプログラムを作りたいのですが教えてください。
Evca(Ω,x)=Ei(Ω,x)+E0(Ω,x)+E1(Ω,x)+E2(Ω,x)
Ei(Ω,x)=-1.79186*Z**(5/3)/r
E0(Ω,x)=Z(2.21/r**2-0.916/r+Ecorr)
Ecorr=-(0.115-0.031logr)
E1(Ω,x)=2πZe**2((1-X)Z**Cs*(Rm**Cs)**2(1+2/3*u**Cs)+xZ**K(Rm**k)**2(1+2/3u**K))
E2(Ω,x)=1/2��(4πe**2/ΩG**2)**-1(V(G))**2*1/(1-f(G))*(1/ε(G)-1)
�狽ﾍG/=0
Gは逆格子ベクトルG(n1,n2,n3)
V(q)=(1-x)V**Cs+xV**k(q)
Vi(q)=-4πZe**2/Ωq**2((1+u)cos(qRm**i)-usin(qRm**i)/qRm) (i=Cs,K)
u=qRmcos(qRm)/sin(qRm)-qRmcos(qRm)
わからない値は後から自分で入れるので教えてください
状態方程式の計算をするらしいです

E2(Ω,x)=1/2��(4πε**2/ΩG**2)**-1(V(G))**2(1/(1-f(G)))(1/ε(G)-1)
�狽ﾍG/=0です
Ｇは逆格子ベクトルです。G(n1,n2,n3)
これを求めるプログラムはどうしたらいいか教えてください
わからない値は自分でいれます。お願いします。

>>401
質問が曖昧すぎるぞい！ｗ
もう少しうまく質問しないと助けたくとも助けられないｗ

見た感じ結晶の電子相関込みの凝集エネルギーかなんかを計算しているようだが，
どの辺のなにがわからんのか，求めたいものは何なのか，もうちょっと他人が分かる
質問をしろｗｗｗｗ

これだけで答えられる奴はいねーよｗ

丸投げしたいなら持ってる情報を全部出せ！ｗ　
誰かが回答するかどうかは別としてｗ

>>401-402
セシウムとカリウムのアルカリ金属合金の凝集エネルギーを求めて
なにかしたいようであることは分かった．


まずは欲張らずに単体のアルカリ金属の凝集エネルギーは求めろよｗ
そしたら，ｘ＝１．０または０．０の時の極限にも対応するので
あとからチェックにも使える．

一体どの辺が，どのくらいわかんないのよｗ　

すいません。あほなんで基本的には全然わかりません。
状態方程式のＰを求めるためにやってると思うんですけど。
圧力Ｐは結晶エネルギーＥvca(Ω,x)の体積Ωに関係する微係数で与えられるらしいのです。
a=5.225A A=10**-8
Ω=(a**3)/2 rs=(9π/4)**(1/3)/Kf
Ef=h**2*Kf**2/2m c=2(1+0.0155π/Kf)**-1
f(q)=q**2/2(Q**2+cKf**2)
ε(q)=1+(1-f(q))*(4πZe**2/Ωq**2)*(3/2Ef)*(1/2+((4Kf**2-q**2)/8Kfq)*log(2Kf+q/2Kf-q))
CsのＨ近似の場合Rmは4.610 uは-0.9605
KのＨ近似の場合Rmは3.019 uは-0.5682

せめて�濡/=0 のＧベクトルが０以外の時を全部足す方法が知りたいです。
お願いします。

k=5.225
a=k/(0.529177)
o=a**3/2
z=1
rs=(3*o/z*4*3.14)**(1/3)
Kf=(9*3.14/4)**(1/3)/rs
Ef=Kf**2
c=2(1+(0.0155*3.14)/Kf)**-1
Ei=(1.79186*z**(5/3))/3*o*rs
E0=z*(2.21/rs-0.916/rs+0.031/3*o)
Rmk=3.019
uk=-0.5682
E1=-2*3.14*z*2/o*(z*(Rmk)**2*(1+2/3*uk))
k=5.225
a=k/(0.529177)
o=a**3/2
z=1
rs=(3*o/z*4*3.14)**(1/3)
Kf=(9*3.14/4)**(1/3)/rs
Ef=Kf**2
c=2(1+(0.0155*3.14)/Kf)**-1
Ei=(1.79186*z**(5/3))/3*o*rs
E0=z*(2.21/rs-0.916/rs+0.031/3*o)
Rmk=3.019
uk=-0.5682
E1=-2*3.14*z*2/o*(z*(Rmk)**2*(1+2/3*uk))

あとはE2(Ω,x)=1/2��(4πe**2/oG**2)**-1(V(G))**2*1/(1-f(G))*(1/ε(G)-1)
�狽ﾍG/=0 を数値微分して全部足すと出るんですけどここがわからないです。
G=2π/a(n1*n1+n2*n2+n3*n3) 収束すると思うのである程度までの結果を出していきたいです。
n1=n2=n3/=0の時以外を計算させる方法がわかりません。数値微分もわかりません。
V(G)=-4*Z*2/oG**2((1+uk)cos(G*Rmk)-uk*sin(G*Rm)/G*Rm)
ε(G)=1+(1-f(G))*(4*3.14*z*2/o*G**2)*(3/2*Ef)*(1/2+(4*Kf**2-G**2)/(8*Kf*G)*log((2*Kf+G)/2*Kf-G)
f(G)=G**2/4*(G**2+c*Kf**2)+G**2/4*c*Kf**2
です。

>>405
ＯＫ
エネルギーを体積で微分すれば圧力になるのは熱力学の帰結だ．
しかしこれはFORTRANというより物理の質問だなｗ

＞せめて�濡/=0 のＧベクトルが０以外の時を全部足す方法が知りたいです。
ＦＯＲＴＲＡＮスレ的には，これは３次元逆格子ベクトルに関する和なので，
もっとも素朴には３重のＤＯ　ＬＯＯＰを回せばよいとい答えになる．
DO n1 = -n, n
　DO n2 = -n, n
　　DO n2 = -n, n
IF(n1==0 .and. n2==0 .and. n3==0) CYCLE　
云々
END DO
END DO
END DO
実際は結晶の対象性があるので，このループの１／８〜１／２４位回せばいいはず．


実際の計算には逆格子ベクトルを求める必要がある．これは結晶構造に依存する．
もし今アルカリ金属が対象ならbcc構造ということになるだろう．
その逆格子ベクトルはfccと同じものになる．
これは固体物理の教科書に載っているから，一般式をつくるか，
手で求めて表を配列に入れることになる．

DO　LOOPの中にn1,n2,n3の組み合わせが逆格子ベクトルにふさわしいかどうかを
判別する式を，G=0を飛ばすのと同じように，いれてやればいい．

>>407
Gは逆格子ベクトルそのものではなくて、その長さだから（ルート忘れてるみたいだが）
たんに逆格子点に相当するn1,n2,n3が分かっていればいい。

後は式の通りに書いて、逆格子点について足せばいい。

FORTRANの疑問ではないので、ここで聞く前に周囲の人に聞いた方が早いと思う。


数値微分は単に、格子間距離aを、特定のaとａ＋ｄａに取ったときに計算して
差分をｄａでわればいい。ｄａは小さい数だ。

fortran90で、数式が二行に渡るとエラーがでるのですが、
どうしたらいいですか？

>>410
一行目のしりに＆をつけろ。

C言語やPASCALと違って、FORTRANは終端文字（；）が不要なかわり行センシティブ。
ゆえに継続行の記号を明示する必要がある。

漏れはFortran方式が好き。

Ruby や Haskell のような方式もあるぞい

>>408
>>409
わざわざありがとうございます。物理の質問でしたねｗ
ＣＹＣＬＥとはどういう意味ですか？
結果にＮａＮとでるのはどういうことですか？
お願いします。

>>413
CYCLEというのはFortran90で導入された命令で，DO文の中で使う．
それ以降の文を実行すると無くDO行へ戻る．
つまり，次のループ計算を実行する．
これと類するものにEXIT命令があって，これはEND　DOの次の行に飛んで
DO　LOOPを脱する．

NaNとは数値に非ず．Not A　Numberの略号だ．配列がはみ出してとんでもない値を
書いているか，オーバーフロー等を起こしている．

SUBCHKなどのデバッグオプションをつけてみそ．

オーバーフローだと普通は ±inf じゃない？
NaN は実数での負値の平方根とか、未定義の演算結果かな。

＞オーバーフローだと普通は ±inf じゃない？
おっしゃるとおり．漏れの間違いだ．
負の平方根とか関数の定義域外呼び出しでもでるね．

漏れが出会うNaNは，ゼロ割か，配列はみ出しでとんでもないメモリー領域の
アクセスしたときが多いかな．個人的にはゼロ割で実行停止しなくなったのでとても助かる．



NaNを厳密に追求したければ，IEEE745の浮動小数点の定義を調べればいいが気力那須ｗｗ

最近、大学でFortranの授業が始まりました。
今日で授業が三回目だったのですが次のような課題がだされました。
おそらく基礎中の基礎だと思うのですが、授業はプリントに書かれた
プログラムを説明もなくそのままうつして実行させるような授業ですので
この課題を読んだ時日本語かどうかすら疑いました。
　いろいろと調べているのですがなかなかよくわかりません。よろしければ
どなたかこの課題について解答していただけませんか？
　あと、本当に自分勝手で申し訳ないのですが、少しでいいので解説もして
いただけたらとても嬉しいです。よろしくお願いします。
（１）subroutineを呼ぶときの、変数の引き渡し方について説明せよ。
特に、配列変数を引き渡す方法についてプログラム名をあげて説明せよ。
（２）Fortranで文字列を取り扱うにはどうすればいいか調べて、プログ
ラム例をあげて説明せよ。
（３）ＲＥＡＤ文やＷＲＩＴＥ文のＩ，Ｅ，Ｆ，Ａ形式について調べて
解説せよ。またプログラム例と印字例を示せ。

すいません。（１）のなかのプログラム名は間違いでプログラム例でした。

調べてって書いてるんだから、調べれば？
プリントか教科書あるんじゃないの？

プリントは説明がなくて・・・教科書もないんです・・・やっぱりちょっと
調べればわかることみたいですね。

FORTRAN で検索したら色々出てくるよ。

ありがとうございます！いろいろ調べて自分なりの解答作ってまた書き込み
にくるんでよかったらみてください。

本物のプログラマって現存しますか?

あれは最初からジョークだからw

色々示唆には富んでるけどね。

>>419
教えてやれよｗ

（１）なんて、参照渡しと値渡しの説明を求められていたら初心者には大変だぜ！
しかし、配列の先頭要素のアドレスしか渡していないことを知らないと７７では
実戦的なプログラムは書けない。

（２）はFORTRAN７７以降ならCHARACTERがあるが、FORTRAN６６だと整数型に
ENCODE／DECODE文で文字列突っ込んだりしてるし。

（３）FORMAT書式は初心者がつまづく、難所中の難所だろう。



と、>>417をびびらせることを書いてみたりしてｗ　うふｗ

data init/'12345'X,'23456'X,'34567'X,'45678'X/
↑
'整数'X　が何を意味するのか分かりません。
よろしくお願いいたします。

'整数'Xは16進数を表すものでした。
お手数をかけました。

>>417
大阪大学工学部応用理工学科二回生乙。

これが噂に聞いた「ソシャる」というやつか

>>422
臆せずに書き込むように。

みんなの勉強になる。

次の宿題もってこいや！！

>>392を使わせてもらったんですけど、実行できませんでした。
cygwinを起動してコンパイルするとこまでは行けたんですけど(g77 report.f)、
./a.exewoをしてもcommand not found になります。どこかエラーがあるのでしょうか？

./a.exewo→./a.exeの間違いです。

>>432
command not found
ということは、a,exeが見つからないということなので、コンパイルに失敗している。

それ以前に何かエラーメッセージが出ていないか見て味噌。
まずはｌｓしてa.exeができているか見て味噌。

どなたか、教えてください;:il:il|;l|;il:i(-ω-｀；)ll|l|il|;:il|!
マジでなやんでます…。
1/√(π)exp(-x^2) 中央値０　偏差１　の正規分布表をあたえる。
[−１．１]および[−２．２]の範囲に存在する確立はいくらか？
で台形法で解けと…。

プログラムがわかんないです…。
たすけてください…

>>435
>>364

一緒だ…。

なんとなく、わかりました。ありがとーございます<(_ _)>

>>434
ありがと。もっかいやってみます。Invalid first character at (^)とか出るんですよね・・・。

>>438
同じ授業を受けてる気がする・・・。特定はやめてね・・・

＞＞439
それ、思った…。

>>440
ここでもその名前でお願い
http://ex11.2ch.net/test/read.cgi/campus/1149336703/

>>439
もしかしてプログラムを「１桁目から」書いてないか？

2chの表示の都合上、ここではプログラムが行の先頭から書かれているように見えるが
実際のFORTRAN77のソースは行の各桁に何を書くかも文法上で決まっている。
そこら辺は初歩の初歩でどんな入門書にも載っているようなことなので
少しは教科書なりWeb上の入門記事なりを読んで勉強してくれ。

>>442
Σ(ﾟДﾟ)！すいません、そのままコピペしてました。

>>443
君は２ｃｈ専用ブラウザを使ってるかい？
使っているなら>>392の引用符を見れば一応ある程度インデントが効いて
本来の姿に近くなる。

>>392の中で自分を引用しているのは、そうやって見るためだ。

もうちょっと最小限のところは勉強しろｗｗｗｗｗｗ
他のスレならぶちキレられてるぞｗｗｗｗｗｗ
このスレはVIPレベルなのでOKだがな（＾ω＾；）

>>435のケロも遠慮せずに書き込めよｗ

>>414
わかりました。ありがとうございます。
一応作ってみたのですがエラーがでてコンパイルできません。
初めて作るのでめちゃくちゃかもしれませんがおかしいところを教えてください。
お願いします。
k=5.225
a=k/(0.529177)
o=a**3/2
z=1
rs=(3*o/z*4*3.14)**(1/3)
Kf=(9*3.14/4)**(1/3)/rs
Ef=Kf**2
c=2*(1+(0.0155*3.14)/Kf)**-1
Ei=(1.79186*z**(5/3))/3*o*rs
E0=z*(2.21/rs-0.916/rs+0.031/3*o)
　Rmk=3.019
uk=-0.5682
E1=-2*3.14*z*2/o*(z*(Rmk)**2*(1+2/3*uk))

E2 = 0
E4 = 0
da = 0.001

WRITE(*,*)'(n4)'
READ(*,*) n4
DO 10 n1= -n4 , n4
DO 10 n2= -n4 , n4
DO 10 n3= -n4 , n4

　IF(n1==0.and.n2==0.and.n3==0) CYCLE
G=2*3.14/a*(n1*n1+n2*n2+n3*n3)
f=G**2/4*(G**2+c*Kf**2)+G**2/4*c*Kf**2
e=1+(1-f)*(4*3.14*z*2/o*G**2)*(3/2*Ef)*(1/2+(4*Kf**2-G**2)/
&(8*Kf*G))*log(2*Kf+G/2*Kf-G)
V=(-4*Z*2/oG**2)*((1+uk)*cos(G*Rmk)-uk*sin(G*Rmk)/G*Rmk)
E3=(4*3.14*2/o*G**2)**-1*V**2*(1/(1-f))*(1/e-1)
E2 = E2 + E3
V1=(-4*Z*2/(o+da)G**2)*((1+uk)*cos(G*Rmk)-uk*sin(G*Rmk)/G*Rmk)
e1=1+(1-f)*(4*3.14*z*2/(o+da)*G**2)*(3/2*Ef)*(1/2+(4*Kf**2-G**2)/
&(8*Kf*G))*log(2*Kf+G/2*Kf-G)
E5 =(4*3.14*2/(o+da)*G**2)**-1*V1**2*(1/(1-f))*(1/e1-1)
E4 = E4 + E5　
10 CONTINUE
E6 = E2/2
E7 = E4/2
E8 = (E7-E6)/da
P=Ei+E0+E1+E8

WRITE(6,601) P
601 FORMAT(1H , F15.8)
STOP
END
エラー
fortcom: Error: Illegal character in statement label field [
fortcom: Error: Illegal character in statement label field [・
fortcom: Error: keita.f, line 55: Syntax error, found END-OF-FILE when expecting one of: <LABEL> <END-OF-STATEMENT> ; BLOCK BLOCKDATA PROGRAM TYPE COMPLEX BYTE CHARACTER ...
.
^
compilation aborted for keita.f (code 1)

あほすぎて迷惑かけます。

-1以上1以下の値ｘを入力してe^xの近似式を出力するプログラム。e^xの近似式はe^x=1+x+x^2/2!+x^3/3!+..+x^n/n!を利用。
(最後の項に絶対値が0.00001より大きい間、項を増やして足し続ける)また、適正なxの値が入力される限り、何度もこの計算ができるようにすること。って課題が出されたんですけどまったく分かりません。誰か助けてください！

>>445
我慢して聞いてくれｗ

基礎が全くなっていないので，ちょっと見ただけでもまともに動くとは思われない．
周りに聞く人がいるなら聞いた方がいい．掲示板で説明出来る範囲を超えているｗｗｗｗ

１．まず変数の型がむちゃくちゃ．第一行目から，整数型のＫに実数を入れている．
式の途中に（１／３）みたいなのがいっぱいあるが，整数の場合はこれは０になる．
ようするに型に関する認識が全くなってないｗ　整数と実数は違うものだと思った方がいい．
実数に整数を代入する方は，一応型の変換が自動でなされるが，明示的に実数にしておいた
ほうがいい．

２．有効桁に関する認識もなってない．円周率を３．１４で入れているが，これでは有効桁が
小数点以下２桁になってしまう．そもそも常数は，ＰＩとかちゃんとした名前の付いた変数
ないしパラメータ変数にすべきもの．

３．式も所々演算子が抜けているし，演算子の優先度に関して間違っているとしか
思えないところがある（割り算の分母）．

４．微分を取りにゆくときのセンスも悪い．まず特定の格子間距離に対してエネルギーを
求めるサブルーチンをつくって，それを２回異なる格子間距離について呼ぶべき．

とにかく，なんというか，プログラミングの基本的センスがなっていないし，
問題自体に対する理解もなってないようなので，
今のままでは他者からのかなりの助けがないと完成はおぼつかない．

周りに聞く人がいれば，頭を下げて一から教えてもらえ．
とにかく道はティッペラリーより遙かかなただｗ

>>448
ほれ、７７形式で昔のテキスト風に芋くさく書いておいたぞい。

PROGRAM VIPPER
1 WRITE(*, *) 'INPUT x (-1.0<= x <= 1.0)'
READ (*, *) x
IF (ABS(x) > 1.0) STOP
e = 1.0
f = x
i = 1
10 CONTINUE
e = e + f
IF (ABS(f) < 0.00001) GOTO 999
i = i + 1
f = f * x / REAL(i)
GOTO 10
999 CONTINUE
WRITE(*, '( 2(a, f8.5) )') 'EXP(', x, ')', e
WRITE(*, *) EXP(x)
GOTO 1
END

すまんｗ　ついうっかりFortran90の癖を出してしまった。
７７としては以下の部分を置き換えておいてくれｗ

IF (ABS(x) .GT. 1.0) STOP

IF (ABS(f) .LT. 0.00001) GOTO 999

大学はいってから、別にプログラム関係の学部でもないのに
fortranをやらされ、まったくわからない…

なんか時分秒の引き算のプログラムを作らなければいけない。

　開始時刻と終了時刻を指定すると、経過時間
を計算し端末に表示するプログラムをつくりなさい。
時刻はhh:mm:ssの形式で指定されるものとする。
終了時刻が開始時刻より早い場合には、その旨を画面に
表示して、再入力させること。

だそうだ。

PRPGRAM TIME
REAL A,B,C,D,E,F,G,H,I
READ(*,*)A,B,C,D,E,F
X=A*60*60+B*60+C
Y=D*60*60+E*60+F
IF(X.GT.Y)THEN
Z=Y-X
G=Z/60/60
M=Z-G*60*60
H=M/60
I=M-H*60
WRITE(*,*)G,H,I
ELSE
WRITE(*,*)'ERROR'
END IF
STOP
END

こんなのが限界で、hh:mm:ssと入力できないし、出力もできない。
２４時間以降もERRORにしなければならないが、どうやるかわからない。
;←を打つには書式を変えろって言われたけど習ってないし、
教えてもくれなかった＿|￣|○

助けてください。こんな簡単なことを聞いてすいません…

>>453
ほれ、出来るだけ藻舞のオリジナルを尊重しつつ77で作った。
変数は整数型にしないといかんべ。
MODは余りを求める関数だ。マニュアルみそ。

入力をhh:mm:ssの形式にしたければ、入力エラーは無いことを仮定すれば、
ＷＲＩＴＥ文のＦＯＲＭＡＴをちょっと直してＲＥＡＤ文の方にも書いてやればいい。



PROGRAM TIME
INTEGER IH0, IH2, IM0, IM1, IS0, IS1, IT0, IT1
INTEGER IH, IM, IS, IT
1 READ(*, *) IH0, IM0, IS0, IH1, IM1, IS1
IT0 = IH0 * 60 * 60 + IM0 * 60 + IS0
IT1 = IH1 * 60 * 60 + IM1 * 60 + IS1
IF(IT1 .GT. IT0)THEN
IT = IT1 - IT0
IH = IT / (60 * 60)
IT = MOD(IT, 60 * 60)
IM = IT / 60
IS = MOD(IT, 60)
WRITE(*, '( I2.2, A1, I2.2, A1, I2.2 )' ) IH, ':', IM, ':', IS
ELSE
WRITE(*,*)'ERROR'
GOTO 1
END IF
STOP
END


>>454

２４時間がどうのこうの謂う所は自分でやってちょｗ
ＩＦ文で飛ばすだけだし、>>452を見る限り自力で出来る感じがするｗ

>>411
遅レスですいません！たすかります、どうもっす！

>>452
＞大学はいってから、別にプログラム関係の学部でもないのに
＞fortranをやらされ、まったくわからない…

来るべき２１世紀の脱工業化社会には、電子計算機の知識が欠かせないのだ。
計算機言語のサンスクリット語たるＦＯＲＴＲＡＮを学べることを誇りに思って励め！

>>457
電子計算機の基礎ですよね…がんばらなきゃ（；゜〇゜)
>>454
すいませんが、( I2.2, A1, I2.2, A1, I2.2 )の.2って何ですか？
I2は２行表示で、A1は…？教科書はめちゃくちゃ高いのに
ぜんぜん詳しく書いてない��(￣Д￣)

教わってる教授の本だからしょうがないか…＿|￣|○

>>458
I2 と　I2.2 の違いは，桁が余ったときの処理の違い．　後者は余ったところを０で埋める．

時計で言えば，１時８分８秒を前者で出すと　１：　８：　８となる，
後者で出すと０１：０８：０８となる．

リターンしちまったｗ

Ａ１の方は文字１桁ということ，今の場合は１文字．　
数字を省略すると，文字列長に自動的に合わさる．

そうなんですか
わかりました。
入力をhh:mm:ssの形式にするには
READ(*,*)を
READ(*,10)
10 FORMAT(I2,':',I2,':',I2)
って変ですか？わからない…
なんかだめな気がするので、助言お願いします。

>>461
1 WRITE(*, *) 'INPUT START TIME hh:mm:ss '
READ(*, '(I2, 1X, I2, 1X, I2)') IH0, IM0, IS0
WRITE(*, *) 'INPUT END TIME hh:mm:ss '
READ(*, '(I2, 1X, I2, 1X, I2)') IH1, IM1, IS1

たとえばこんな感じにできる。１Xというのは１文字分飛ばせということなので、
実際は　：　の代わりに空白でも何でもいい。とにかく無視される。

ちゃんと　：　を認識して処理したければ、CHRACATER型の変数に読んでチェック
すればいいのだが、その場合時分秒がちゃんと０〜２４、０〜６０に入っているかを
チェックするほうが常識的に先になる。

(0,1)の擬似乱数をMT19937で発生させて、
その平均値と標準偏差が0.5と0.25になるのを確認しようとしました。
平均値は0.5にかなり近づくんですが、標準偏差が0.28程度になります。
乱数の個数を増やしていくと、0.28程度で収束しているように見えます。
octaveでもやってみたんですが、0.25程度に収束しません。
適合度検定にパスするとは思うのですが。

擬似乱数はばらつき具合はうまく表現できないのでしょうか？
よろしくお願いします。

>>463
漏れはシミュレーションとかやらないので乱数知らないのだが、
擬似乱数って用途によってアルゴリズムを選べと言われていたはず。


乱数マニアのサイトを覗くと粘着した情報が得られると思うｗｗ
ものによって最初の数万個は捨てろとか色々個別の事情があるようだ。

返事ありがとうございます。
乱数マニアのサイトをもう少し見てみます。

１文字分飛ばせばいいんですね！
ちゃんと　：　を認識して処理したければ、CHRACATER型の変数に読んでチェック
って言うのは…？
もっと勉強したいと思います。
時分秒がちゃんと０〜２４、０〜６０に入っているかを
チェックするのは、IF文で０〜２４、０〜６０じゃなければERRORにすればいいですよね。

>>466
ういうい。

CHARACTER DUMMY1*1, DUMMY2*1
READ(*, '(I2, A1, I2, A1, I2)') IH, DUMMY1, IM, DUMMY2, IS
IF (DUMMY1 .NE. ':') STOP 'ERROR'
・・・・・・・・
IF (IH .GT. 24 .OR. IH .LT. 0) STOP 'ERROR'
云々

まあ、いきなりストップせず、WRITE 文で入力値を出力して、
どこがおかしいのか書いたほうが親切ではある。

alpha 10.0
beta 0.00
....
gamma 4E-3

↑という感じのファイルを読み込みたいのですけど、
一行に文字と色々なフォーマットの数字が並んでるんですけど、
できますか？

文字は20文字固定で、余りは半角で埋めようと思っています。

>>468
できる。
もし文字と数値の間に必ず空白が空いているなら簡単。

CHARACTER NAME*20, TEXT*100
READ(9, '(a)') TEXT
IPOS = INDEX(TEXT,' ')
NAME = TEXT(1:IPOS - 1)
READ(TEXT(IPOS + 1:IPOS + 20), '(F20.0)') VALUE

こんな感じで、一度１行全体を文字列として読み、名前と値の境目を見つけたら
そこより前は文字変数に突っ込み、境目から後ろは文字列を内部ファイルとして
もちいて、数値として数値変数に読み出せばいい。
（ここでは２０桁の数値を仮定した。F20.0にしたので小数点がないと
とんでもない値になるが、この辺もINDEXで尻を見つけるとか、
FORTRAN90でいいならTRIMや何やらを使えば回避できる。）

[0,1]の一様乱数の標準偏差は0.28...で、0.25ではありませんでした。
お手数をかけました。

-1以上1以下の値ｘを入力してe^xの近似式を出力するプログラム。e^xの近似式はe^x=1+x+x^2/2!+x^3/3!+..+x^n/n!を利用。
(最後の項に絶対値が0.00001より大きい間、項を増やして足し続ける)また、適正なxの値が入力される限り、何度もこの計算ができるようにすること。って課題が出されたんですけどまったく分かりません。誰か助けてください！
って昨日打ったんですけどc++でお願いします！！

>>471
ちょｗｗｗ　おまｗｗｗｗ
Ｆｏｒｔｒａｎスレに聞くなよｗ
>>450を自分で直せｗｗｗｗ

テキストファイルから読み込んだ値を変数に代入したいのですが、
テキストの値の書式が
1,22.489,445,67888,23
のようにカンマ区切りの場合、どのように読み込むのが良いでしょうか？

>>473
書式無しの自由書式。この場合コンマで区切って読んでくれる。

それで問題がある場合は、文字列として読んで、自分でコンマのところで切って
文字列を内部ファイルとして数値に変換する。

>>474
integer i,j,k,l,m
real r
open(10,file='hoge.txt')
read(10,*)i,r,j,k,l,m
のような感じで良いのでしょうか？？

>>475
いいはず

>>476
ありがとう

c言語にはfortranでいうcycle機能ってあるの？

continue;

すいません。昨日からはじめた超初心者です。
初心者用のサイトを見ながらやってますが、
書いてある通りに行かないので質問させてください。
FTN95でコンパイルする際、
「ftn95.exe プログラムのファイル名.for95」と入れれば.exeファイルができる
と説明してあるのですが、
上のように入力しても、「.exe」ではなく「.obj」ファイルになってしまいます。
どうすれば.exeになるのでしょうか？

>>480
コンパイル時に/linkオプション。つまり
ftn95.exe プログラムのファイル名.for95 /link

>>481
ありがとうございます。
出来ました。

>>467
ありがとうございました。

real X
read(*,*)X
として、小数点がつかない
5
などを読み込んだとき、Xにとんでもない値が入ってしまうのを
回避する方法はあるでしょうか？

行き詰ってしまっているので、失礼します

1行に
(整数):(整数)(3文字の文字列)(8ケタの整数)(17文字の文字列)(空白1)(14文字の文字列)
というデータが並び、これが何百行と続くデータを読みこみ、
(8ケタの整数)というもののみを、他ファイルに出力したいと思います。
(整数)というのはケタ数が決まっていません。これがエラーでコンパイルできません。ご指導お願いします

program read
implicit none
integer::i,j,k,l
open(unit=11,file='input.dat')
open(unit=22,file='output.txt')
　do l=1,100
　read(11,'(I8,1X,I15,3X,I8,32X)i,j,k
　write(22,*)k
　end do
close(11)
close(22)
end program read

l=1,100　となっているのは、100行読み取りたいとき、ということにしてます

修正です
　read(11,'(I8,1X,I15,3X,I8,32X)')i,j,k
でお願いします

また失礼します。コンパイルはできるのですが、実行ができないのでした・・

fmt: read unexpected character
apparent state: unit 11 named input.dat
last format: (I8,1X,I15,3X,I8,32X)
lately reading sequential formatted external IO
aborted (core dumped)

と表示されます。環境はcygwinで、g77でコンパイルをしています。ぐだぐだと失礼しました

事故解決しました。すみません

>>484おながいいたします
あと、"を文字として書くときってどうすれば良いでしょうか

>>484
普通の処理系なら
read(*,*) X
に対して
5
と入力した場合、小数部を省略したとみなして X は 5. となるはずだが？

>>489
>"を文字として書くとき
print *, '"'
print *, """"
print *, 1H"

4:345:23:346543
346:344:4:455

のような、ケタ数が一致しない複数の数値が:で区切られた値を読み取るにはどうしたら良いのでしょうか

>>491
各行に含まれる数値の数が一定なら、以下のように一旦文字列として
読み込んで加工し、内部ファイルを使えば良い。
program hoge
implicit none
integer i, k,l,m,n, ios
character*80 str
open(11,file='data.dat',status='OLD')
100 continue
read(11, fmt='(A)', iostat=ios) str
if (ios .ne. 0) goto 999
do 200 i=1, len(str)
if (str(i:i) .eq. ':') str(i:i) = ','
200 continue
read(str,*) k,l,m,n
write(*,*) k,l,m,n
goto 100
999 continue
close(11)
stop
end

>>492
ありがとうございます。内部ファイル・・？
これを参考にして頑張ります

>>493
>read(*,*)X
>として、小数点がつかない
>5
>などを読み込んだとき、Xにとんでもない値が入ってしまうのを
>回避する方法はあるでしょうか？
自由書式の場合は問題ない。
FORMATを与えていた場合に問題が生ずることがある。
小数点があれば与えたFORMATよりも、入力の小数点が優先されるが、
小数点がないと与えたFORMATの桁数にしたがって小数点位置が決められる。

>>491の安易な解決法は：を、に置換して自由書式で読み出すことだｗ

>>491の問題は>>492の方法で解決できるが、Fortran90以降を使うと
文字列関数が豊富な分やりやすい。ただし77の範囲でも実現は可能。
だらだらと書くもよし、：までを切り出した後、文字変数に右詰に
して返す関数を作るもよし。

内部ファイルと外部ファイルの違う点は、ファイルの読み込み位置が
毎回先頭に戻ってしまう事で、これがちょっとうざい。
外部ファイルとなんでも一緒のつもりで行って昔ちょっとはまったｗ

文字列の配列ってできますか？
character*50 c(10)
みたいな感じで。

あっ、できました。
また質問いたします。

integer i
character*50 c
do i=1,4
write(c,'(i4)')i
end

みたいな感じにして、
cに「1234」のような文字列が入るようにしたいのですが、
上手くいきません。
どうすればよろしいのでしょうか？

>>496
integer i
character*50 c
write(c,'(4(I1))') (i, i=1,4)
write(*,*) c
end

>>496
別解
DO i = 1, 4
WRITE(c(i:i), '(i1)') i
END DO

>>497
>>498
どうもサンクス

もう一個、
読み込んだファイルの全行数を知ることはできませんか？

>>499

多分1回１行づつ数えるしかないと思うのだが。ほかに方法あるのかな？


iline = 0
DO
READ(9, iostat = io, *)
IF (io < 0) EXIT
iline = iline + 1
END DO

こうかな？

巻き戻し
REWIND(9)

>>500
あ、できました。どうもありがとう
READ(9,*, iostat = io)
ですね

今一番新しいfortranの規格ってfortran2003ってやつですか？
これ普及すんのかなぁ

>>502
する！

まだ実装したコンパイラないけどなｗ

ﾏｼﾞｽｶ
まだ77でさえしぶとく残ってるのにね…

>>504
６６だってまだまだ健在だぜ！

９０で書いてるって言ってる奴も、珍妙な７７や６６の入り混じったハイブリットだったりするし。

そこまでいくとｶｺｲｲなｗ

fortran90で
ｎ個の配列を準備してＲＡＮＤＯＭ＿NUMBER関数を用いて
ランダムな値を配列に代入する。
ばらばらの順になっている配列の実数を、小さな順に並び替えるサブルーチンを作成し、
そのサブルーチンをcallする。
結果を出力する。


でばらばらの順になっている配列の実数を、小さな順に並び替えるサブルーチンを作成し、
そのサブルーチンをcallする。
というところがさっぱりわかりません。
どうすればいいでしょうか？？

>>507
とりあえず、「ソート アルゴリズム」でググれば
いろいろとソート（並び替え）の方法が出てくるよ。
あとはそのうち簡単なものを選んで実装するのみ。

>>507
作ってみた
program rndm
real,dimension(:),allocatable :: x; integer :: i,n,stat
read(*,*) n
allocate(x(n),stat=stat)
if(stat/=0) stop "could not allocate x"
do i=1,n
call random_number(x(i))
enddo
print *,"before sort"; call printx(x)
call srt(x)
print *,"after sort"; call printx(x)
deallocate(x)
contains
subroutine printx(x)
real,dimension(:) :: x; integer :: i
do i=1,size(x); print *,x(i); enddo
end subroutine printx
recursive subroutine srt(x)
integer :: i,iflag; real :: tmp
real,dimension(:) :: x
iflag=0
do i=1,size(x)-1
if(x(i)>x(i+1)) then
iflag=1; tmp=x(i); x(i)=x(i+1); x(i+1)=tmp
endif
enddo
if(iflag/=0) call srt(x)
end subroutine srt
end program rndm

>>509
おまいFortranプログラマーでないだろうｗ

行内複数命令、暗黙の型からの逸脱、再帰関数・・・・・

他の言語だと推奨されるところだが、Fortranでは嫌われるものばかりｗ

育ちが分かるぜｗ　うっしっしｗ

>>510
一応77育ちで最近90に手を出したところ。
行内複数命令はやりたくなかったけど、「改行が多過ぎます」の
書き込みエラーが出たのを回避するため。行数を減らす時に
「implicit none」を消してしまった。
>>507がfortran90と書いていたので、90で取り込まれた仕様を
使ってみた。call random_number(x(i))のところはdoループに
しないでcall random_number(x)でも良かったみたい。

行内複数命令を使わずに出来るだけ短くしてみた
program rndm
real,dimension(:),allocatable :: x
read(*,*) n
allocate(x(n))
call random_number(x)
print *,"before sort"
print *,x
call srt(x)
print *,"after sort"
print *,x
deallocate(x)
contains
recursive subroutine srt(x)
real,dimension(:) :: x
iflag=0
do i=1,size(x)-1
if(x(i)>x(i+1)) then
iflag=1
x(i:i+1)=(/x(i+1),x(i)/)
endif
enddo
if(iflag/=0) call srt(x)
end subroutine srt
end program rndm

>>512
x(i:i+1)=(/x(i+1),x(i)/)
うわ！　こういうswap tips があったのか。

>>513
隣接で良ければ x(i:i+1) = x(i+1:i:-1) という手もあるな．

>>514
Fortran 90 以降の並列型配列処理の強力さが全開！！

>>515
強力なんだが、暗黙の内にテンポラリ配列に展開するので、
超巨大配列に適用すると、すゲー時間がかかったりメモリーを食ったりする。

サブルーチンに渡すときに７７の時のつもりで先頭要素のアドレスだけ渡している気になると
実はテンポラリにコピー、引き渡し、受け取り、テンポラリから元の配列へコピーとやってたりするｗ

すいません、くそ初心者なんですが文の途中に改行するときに
入れる記号を教えてください！

>>517
すいません、くそ初心者なんですが文の途中に改行するときに &
入れる記号を教えてください！

>>517
Fortran90なら>>518にあるように　＆　を尻に付ける．

FORTRAN７７なら　６　カラム目に何か文字を入れる．＋とか数字がよく使われるかな．

>519
> FORTRAN７７なら　６　カラム目に何か文字を入れる．＋とか数字がよく使われるかな．
カーニハン、プローガー、プログラム書法、第2版、木村泉訳、共立出版、p.21では、
継続行を表すのに$が使われています。
その理由が次のとおり書かれています。
（ここで継続行をあらわすのに「$」を使っているのは、「$」がFortranの規格上
公式に認められている文字のうち、構文上特別の意味を持たないただ一つの
ものであることによる。こうすれば混乱が生ずる機会を最小にすることができる
とともに、欄をまちがったときにそれに気づく可能性が大きい。）

なるほど…
90でも$にしたほうがいいの

主プログラムの中で配列の番号指定して
それをfunctionに渡してやりたいんだけど・・・うまく渡せないんだけど･･･
どうしたらいいの？（＞＜）

program UNKO
real chinko(10)　←ファイルから読み込ませる
i=5
なんちゃらかんちゃら

function sex(ぱらめーた)
sex=chinko(i)*manko



chinko(i)がゼロになるｗ

>>520
そういえばプログラム書法にはそうあったような気がするな。
漏れはＡＣＯＳっ子だったので＆を６カラム目に置かされていたｗｗｗｗ

数字を使う人たちは１，２，３，４・・と行が継続行の数だけ数字を増やして行く流派の人たちだった。
＋は昔のＦＯＲＴＲＡＮのテキストに多かった気がする。

>>521
９０では前の行の尻に＆だ。＆以外は許されない。
続いた行頭においてもよかったはずだが。

９０で固定カラムというのならまぁ＄でもいいがｗｗｗｗ

>>522
いくつかやり方がある。

い）必要な要素だけを渡す。
program UNKO
real chinko(10)　←ファイルから読み込ませる
anal = sex( chinko(5)、ぱらめーた　 )

function sex(chinpo, ぱらめーた)
sex=chinpo*manko

ろ）配列全部渡す。
program UNKO
real chinko(10)　←ファイルから読み込ませる
anal = sex( chinko、ぱらめーた　 )

■７７式　甲
function sex(chinpo, ぱらめーた)
real chinpo(10)
i = 5
sex=chinpo(i)*manko
■７７式　乙
function sex(chinpo, ぱらめーた)
real chinpo(*)
i = 5
sex=chinpo(i)*manko
■９０式
function sex(chinpo, ぱらめーた)
real, intent(in) :: chinpo(:)
i = 5
sex=chinpo(i)*manko

下品すぎます

下品ですいません（＞＜）
一瞬で答えが出ましたねすごいですｗ
必要な要素だけを渡す場合にchinko(5)のところはchinko(i)でもいいんですか？
具体的には主プログラム中の繰り返し回数iによってfunction中の値も変えたいんです（＞＜）

>>526
OK

DO impo = 1, 69
anal = sex( chinpo(i) , ....... )
.......
END DO

やっぱりわかりません泣　どこが違うのでしょうか？

　　H1 = 0.0639
H2 = 0.0001
PAI = 3.14159
T = 0.0164
OMG = 2.0 * PAI / T
D = 0.5
V = 20.0


WRITE(12,10)
10 FORMAT(' T X1(t) X2(T)')
DO 100 N=0,300,1
X = FLOAT(N) / 200.0
Y1 = EXP(-H1 * OMG * X) * (D * COS(SQRT(1.0 - H1**2) * OMG * X)
$ + (D * H1 + V / OMG) / SQRT(1.0 - H1**2)
$ * SIN(SQRT(1.0 - H1**2) * OMG * X) * EXP(-H1 * OMG * X)
Y2 = EXP(-H2 * OMG * X) * COS(SQRT(1.0 - H2**2) * OMG * X)

WRITE(12,20) X,Y1,Y2
20 FORMAT(3F7.3)

100 CONTINUE

>>528

Y1 = EXP(-H1 * OMG * X) * (D * COS(SQRT(1.0 - H1**2) * OMG * X)

この行のＤの前の括弧　（　に対応する　閉じ括弧がないから叱られている．

継続行の問題ではない．もしくは，継続行にしているうちに壊してしまったのだろう．


>>524　への追加だが，　
ＦＯＲＴＲＡＮ７７での配列の渡し方としては，配列の大きさを渡すのが一番一般的だと思う．
数学ライブラリのたぐいはみんなこうだろう．

■７７式　丙
anal = sex( n, chinko, ..... )

function sex( n, chinko, ..... )
real chinko(n)
.......

ついでながら，
PAI = 3.14159
これだと，単精度の有効桁にちょっと足りない．６〜７桁はあるのでもう１，２数字つけておくと安心だ．
まぁ，本来的には他の数字の有効桁にもよるのだが，常数だからギリギリまで正確に与えた方がいい．

PAI = 4.0 * ATAN(1.0)

倍精度なら

PAI = 4.0d0 * ATAN(1.0d0)

がよく使われる．これだと，機械によらずに有効な桁まで数値が与えられる．
倍精度の時にπを与えるのがかったるいので便利である．

すいません、教科書見ながらやってみたのですが全く手に負えません。
この課題お願いします。

ｍ行ｎ列の行列（実数）Ａ，Ｂをキーボードから読み込んで、
行列の積Ｃ＝ＡＢを計算し、結果を画面に表示させなさい。

序盤からなかなかうまいことがいかず・・・お願いします。

>>531
m≠n のときはどうするんだｗ

>>531
Fortran90ならＭＡＴＭＵＬ　という組み込み関数があるぞｗ

>>532
m≠nのときも考えるのだと思います。

一応はいりだけですが・・
DIMENSION A(M,N)
WRITE(*,*) 'M,N='
READ(*,*) M,N

DO 1 I=1,M
DO 2 J=1,N
WRITE(*,*) I,J, '要素を入力'
READ(*,*) A(I,J)
2CONTINUE
1CONTINUE

(つД`)たすけて〜

>>534
助けてやってもよいが、行列の積を考えると　ｍ＊ｎ　と　ｍ＊ｎ　の積は取れないと
>>532は言っている。

ｍ’　＊　ｎ　と　ｎ　＊　ｍ　のように　積をとる部分では次元がそろってないといけない。

もう少し問題文を正確に教えてくれ。

>>534
一応ＦＯＲＴＲＡＮ７７で作ってみた。若干中級レベルの技法を使っている。
まぁ、素直に多重ＤＯ　ＬＯＯＰをまわすのもよし。
ＦＯＲＴＲＡＮ７７では行列の大きさを実行時に決められない。Ｆｏｒｔｒａｎ９０ならできる。

PROGRAM matrix
PARAMETER (m = 2, n = 3)
REAL a(n, m), b(m, n), c(n, n)
C input A
WRITE(*, *) ' input Matrix A '
WRITE(*, *) ' size =', n, ' * ', m
READ(*, *) ((a(i, j), j = 1, m), i = 1, n)
C input B
WRITE(*, *) ' input Matrix B '
WRITE(*, *) ' size =', m, ' * ', n
READ(*, *) ((b(i, j), j = 1, n), i = 1, m)
C C = A * B
DO 10 i = 1, n
DO 20 j = 1, n
c(i, j) = 0.0
DO 30 k = 1, m
c(i, j) = c(i, j) + a(i, k) * b(k, j)
30 CONTINUE
20 CONTINUE
10 CONTINUE
C print C
DO 40 i = 1, n
WRITE(*, '(9F8.4)') ( c(i, j), j = 1, n )
40 END DO
STOP
END
>>536　

>>535
ｍ行ｎ列のＡとｎ行ｍ列のＢということだと思われます。
>>536
ありがとうございます。
理解してみようと思います。

宿題ではないのですが、
WINDOWSのスタートアップに入れておいて、
PC起動してから約1時間後に勝手に計算を始めるFORTRANのプログラムをつくりたいのです。

そこで、以下のような概略を考えました。

・スタートアップでプログラムスタート
・現在時刻を取得して格納・DOループスタート。
このループ中で時刻取得を繰り返し、上で取得した最初の時刻の1時間後過ぎになったらループ終了
・やりたい計算を実行
・プログラム終了

と、こんな感じで考えたのですが、PC起動後ずっと時刻取得のループをまわし続けるのが
どうも無駄のような気がします。
こういう場合、どういうやり方が上手なんでしょうか

>>538
・案１
Windows のタスク機能を使う。
[スタート] ボタンをクリックし、[プログラム]、[アクセサリ]、[システム ツール] の順に
ポイントし、[タスク] をクリックします。

・案２　ＭＳ／ＤＥＣ／Ｃｏｍｐａｑ／ＨＰ／ＩＮＴＥＬ　ＦＯＲＴＲＡＮ　の場合
ＳＬＥＥＰ　関数を呼んで眠っていてもらう。

行列ﾃﾗｳｻﾞｽ

>>540
Ｆｏｒｔｒａｎ９０なら行列気持ちよくいじれる．オナニーできるくらい．

いいなー。ところで、
Gauss-Jordanのピボット選択ってどうやるんですか
絶対値が最大のと０になるところを入れ換える。。
バカなんでさっぱりわからん。。

>>542
最低限度で良ければ，０が出た時に，そこから未処理分の列の係数をみていって
絶対値最大の行と入れ替えればいいだけだろう。

Ａｘ＝ｂ
のｂの対応行も入れ替え。

入れ替えた行を覚えておいて元に戻すのがめんどいかな。

課題が出たのですが、まったくわかりません。
助けてください。

m行n列の２つの実行れい、A、Bをキーボードから読み込んで
行列の席C＝ABを計算し、結果を画面に出力しなさい。

留意点
　汎用性のあるプログラムを設計すること。
　出力設計については結果が見やすいことが肝心である。
　問題の分析、再定義が必要（特に配列の扱いについて）
　行列A、Bの各要素は**.*のように整数２桁、小数１桁の実数である。

と書いてありますが、まったくわかりません。
お願いします。

>>544　昨日の>>531とは別人か？
ＦＯＲＴＲＡＮ７７とＦｏｒｔｒａｎ９０でどっちがいいね？

その問題は設定が悪いｗ　　Ｍ＊Ｎのサイズがあいまい。
>>536を見て、もっと詳しい要望を出せ。出力ＦｏｒｍａｔをＦ４．１に変えればおｋ。

＞行列A、Bの各要素は**.*のように整数２桁、小数１桁の実数である。
これは整数と小数の言葉を間違って使っているｗ

あぁー　ホント行列わかんねぇよ
考えてんだけど・・　Doループで頭がクルクルパー・・

>>543の入れ替えた行を覚えておいて元に戻すのがめんどいかな
って具体的にどうするの？

すいません>>536についてですが
c(i, j) = 0.0
というのはＣの初期値ということでよろしいのでしょうか？

>>546
ググルと色々サンプルプログラムが出てくるので参照してみたらいかがかな。
ピボッティングだったら教科書の類にもよく出てきているので、他人のプログラムを
読んでみるといいと思うんこ。


>>547
このループでｃ（ｉ，ｊ）に次々足してゆくわけだが、一番最初に０．０になっていないと困る。
DO 30 k = 1, m
c(i, j) = c(i, j) + a(i, k) * b(k, j)
30 CONTINUE

Ｆｏｒｔｒａｎでは変数の初期値は不定なのでごみが入っていることがある。
最近のコンパイラ（というか正確にはローダーだが）変数をゼロに初期化してくれることが
多いが、これを期待してはいかんざき。変数はみんな初期化する癖を付けるのがよろしいアル
カイダ。

>>545ＦＯＲＴＲＡＮ７７です。
L＊N行列のものと、N＊M行列のものの積で汎用性を持たせろということなので
10*6と6*9や3*3と3*3もできなければならないものと思いました。
>>548ググってみましたがサンプルプログラムが汎用性があるのがわかりません＿|￣|○
汎用性…

>>549
まず『汎用性』というものを具体的に定義しなければいけない。
ＦＯＲＴＲＡＮ７７での汎用性であるから『行列の大きさを自由に変えられる
サブルーチン形式』と解釈しておけばいい。

以下例・・・

PROGRAM matrix
PARAMETER (nmax = 40)
REAL a(nmax, nmax), b(nmax, nmax), c(nmax, nmax)
CALL inpmat(nmax, nra, nca, a)
CALL inpmat(nmax, nrb, ncb, b)
IF (nca .NE. nrb) STOP 'matrix size mismatch'
CALL matmlt(nmax, nra, nca, ncb, a, b, c)
CALL primat(nmax, nra, ncb, c)
STOP
END
C==== C = A * B ==============================================================
SUBROUTINE matmlt(nmax, n, m, l, a, b, c)
REAL a(nmax, m), b(nmax, l), c(nmax, l)
DO 10 i = 1, n
DO 20 j = 1, l
c(i, j) = 0.0
DO 30 k = 1, m
c(i, j) = c(i, j) + a(i, k) * b(k, j)
30 CONTINUE
20 CONTINUE
10 CONTINUE
RETURN
END
>>551

C=====================================================================
SUBROUTINE inpmat(nmax, nrow, ncol, a)
REAL a(nmax, nmax)
WRITE(*, *) 'input Matrix Size ( row & col ).'
READ(*, *) nrow, ncol
IF (nrow .GT. nmax .OR. ncol .GT. nmax) THEN
WRITE(*, *) 'matrix size must be smaller than', nmax, '.'
STOP 'Increase nmax!'
END IF
DO 10 irow = 1, nrow
WRITE(*, *) 'input Matrix element of ', irow, ' th row.'
READ(*, *) (a(irow, icol), icol = 1, ncol)
10 CONTINUE
RETURN
END
C=====================================================================
SUBROUTINE primat(nmax, nrow, ncol, a)
REAL a(nmax, nmax)
WRITE(*, *)
DO 10 irow = 1, nrow
WRITE(*, '(1x, 40F5.1)') ( a(irow, icol), icol = 1, ncol )
10 CONTINUE
RETURN
END
>>552

ホントは整合配列をきっちりしてやんなきゃないが、面倒なのでいんちきしたｗ
すぐ直せるから調べて直してくれｗ

色々エラーチェックとかも適当ｗ
じゃあ

かなり難しいですね。。
わからない部分がほとんどですが、１つずつわからないところを
つぶしていきたいと思います。
わからなくなったら聞きますのでよろしくお願いします。
ありがとうございました。

サイズも渡せばいいだけじゃね？

>>555
うみゅ。
渡す配列毎に、使用した縦・横方向と、定義した縦方向の大きさが必要。

FORTRAN77流にお行儀悪くやってみたお

PROGRAM HOGE
DIMENSION A(100), B(100), C(100)
WRITE(*,*) 'INPUT L, N, M:'
READ(*,*) L, N, M
WRITE(*,*) 'INPUT A(',L, '*', N, '):'
CALL MTINP(A, L,N)
WRITE(*,*) 'INPUT B(', N, '*', M, '):'
CALL MTINP(B, N,M)
CALL MTMUL(A, B, C, L, N, M)
WRITE(*,*) 'C ='
CALL MTPRT(C, L,M)
STOP
END

>>557の続きだお

SUBROUTINE MTINP(A, M, N)
DIMENSION A(M,N)
READ(*,*) ((A(I,J), J=1,N), I=1,M)
END
C
SUBROUTINE MTPRT(A, M, N)
DIMENSION A(M,N)
CHARACTER F*20
WRITE(F,'(''(1X,'',I3,''F7.1)'')') N
WRITE(*,F) ((A(I,J), J=1,N), I=1,M)
END
C
SUBROUTINE MTMUL(A, B, C, L, N, M)
DIMENSION A(L,N), B(N,M), C(L,M)
DO 10 J=1, M
DO 20 I=1, L
C(I,J) = 0.
DO 30 K=1, N
C(I,J) = C(I,J) + A(I,K) * B(K,J)
30 CONTINUE
20 CONTINUE
10 CONTINUE
END

>>554は、たぶん配列の添え字と縦横の関係も分かっていないと思われｗ

Ｆｏｒｔｒａｎでは、第一添え字が縦方向、第二添え字が横方向を指定する。

たとえば３＊３の場合

Ａ（１，１）　Ａ（１，２）　Ａ（１，３）
Ａ（２，１）　Ａ（２，２）　Ａ（２，３）
Ａ（３，１）　Ａ（３，２）　Ａ（３，３）

という並びになっている。これはＣ言語とは逆でＦｏｒｔｒａｎ型のほうが数学記号と合う。

また『行　（row　ロウ）』というのは、こういう横方向の並びをいう。
Ａ（１，１）　Ａ（１，２）　Ａ（１，３）

こういう縦方向の並びは『列　（column　カラム)』という。
Ａ（１，１）
Ａ（２，１）
Ａ(３，１）

■やや高度な話題
Fortranでは、２次元の配列は、メモリー上ででは、１列づつ縦書きで日本語とは逆に
左から右へ読む形で１次元に並んでいる。つまり、

Ａ（１，１）　Ａ（２，１）　Ａ(３，１）　Ａ（１，２）　Ａ（２，２）　Ａ(３，２）　Ａ（１，３）　Ａ（２，３）　Ａ(３，３）

これはサブルーチンに配列を渡すときに知っていると役に立つ。

ただしFortran90以降では、分散メモリー型の並列機に配慮したせいか、
規格上ではこの種の並びが保障されなくなった。とはいえ、まず信じて大丈夫ｗｗ

どなたか、VB.NETでビリヤード程度のゲームのソースを
お持ちの方、またサイトを知っておられる方
教えていただけませんか??

どなたか↓をC言語にしてもらえませんか？F2C: FORTRAN to C Translator
を使いこなせなくて・・・。よろしくお願いします。

C*** 　report2
C*** 　gakusekibangou
DIMENSION YA(30),VA(30),YS(30),VS(30),YR(30),VR(30)
OPEN(UNIT=16,FILE="dy.d")
OPEN(UNIT=17,FILE="dv.d")
G=9.8
R=0.3
CANALYTIC SOL.
DO 10 K=1,30
T=K
EP=EXP(-R*T)
VA(K)=G/R*(1.0-EP)
YA(K)=G/R*T-G/R/R*(1.0-EP)
10CONTINUE
YS0=0.0
YS1=0.0
YS0=0.0
VS0=0.0
YR0=0.0
VR0=0.0
H=0.01
DO 20 K=1,3000

CSABUN
T=0.01*K
YS2=2.0*YS1-YS0+G*H*H-R*(YS1-YS0)*H
VS2=(YS2-YS1)/H
CRUNGE-KUTTA
BK1=VR0
BM1=G-R*VR0
BK2=VR0+H*BM1/2.0
BM2=G-R*(VR0+H*BM1/2.0)
BK3=VR0+H*BM2/2.0
BM3=G-R*(VR0+H*BM2/2.0)
BK4=VR0+H*BM3
BM4=G-R*(VR0+H*BM3)
YR1=YR0+H*(BK1+2.0*BK2+2.0*BK3+BK4)/6.0
YR1=YR0+H*(BM1+2.0*BM2+2.0*BM3+BM4)/6.0
IF(MOD(K,100).NE.0) GO TO 30
KK=K/100
YS(KK)=YS2
VS(KK)=VS2
YR(KK)=YR1
VR(KK)=VR1
30 YS0=YS1
YS1=YS2
VS0=VS0
VS1=VS2
YR0=YR1
VR0=VR1

20 CONTINUE
WRITE(17,100)
DO 40 K=1,30
WRITE(17,150) K,VA(K),YS(K),VR(K)
40 CONTINUE
WRITE(16,200)
DO 50 K=1,30
WRITE(16,150) K,YA(K),YS(K),YR(K)
50 CONTINUE
100 FORMAT(5X,"T",5X,"VA",11X,"VS",11X,"VR")
150 FORMAT(2X,T3,3X,F10.5,3X,F10.5,3X,F10.5)
200 FORMAT(5X,"T",5X,"YA",11X,"YS",11X,"YR")

STOP
END

>>560-563 UZEEEEEEEEEEEEeeeeeeeeeeeeeeeeeewwwwwwwwwww

やっとこさコンパイルできたと思ったら実効できない・・・。

run-time error
error format/data mismatch
00401000 MAIN# [+03e4]

これはどういうエラー？？

>>565
FORMAT文のエラーでないかな。

オプションをつけるとソース行が出ることがある。
トレースバック　というデバッグオプションをオンにしてみそ。

デバッグオプションがわかんない場合は、
１）近くのわかっている人に聞くか、
２）あきらめて、プログラムにＰＲＩＮＴ　＊、’ｕｎｋｏ　１’　等のＰＲＩＮＴ文をバラまいてエラー行を特定せよｗ

すいません、real*8の変数がNaN(Not a Number)になったとき、プログラムを
終了させるにはどうしたらいいですか・

Gauss Jordan法を用いて逆行列を求めなさい
ピボット選択：i行以降の行j=i〜nでa(j,i)の絶対値が最大である行を捜し
　　　　　　　その行とi行の入れ替えを行う。

　よろしくお願いします(TT)

>>567
Fortran2003にならないと標準命令ではできない。
ベンダー拡張でできる場合もある。
IEEEのFormatを知っていれば自分でBITを調べる関数は作れる。

またランタイムのデバッグオプションで何とかなる処理系もある。


一般的にはNaNが出るのは配列のはみ出し等なので、まずはバグをなくせ！ｗ

>>569
そうなんすか、どうもっす！

NaN は、負数の sqrt でよく出てくるな。
数値誤差でギリギリ負になる場合ってのもあるから、
式上は非負でも安心できない事が。
負なら 0 とする関数を作っておく必要がある。

functionの引数にfunctionの値入れるにはどうしたら良いの？
f(x)=3x+g(x)+a+b
g(x)=x**2+c

program test
real x,a,b,c
external g

S=f(3,g,4,4,4)

function f(x,g,a,b,c)
xeternal g
f=3x+g(x,c)+a+b
return
end

function g(x,c)
g=x**2+c
return
end


いちおこれでいけたんだけど
g(x)=(x,h(x),d,e)

とかって3連コンボにするとうまくいかないんです（＞＜）

function hの引数をgやｆの引数いれたらいいんだよね？
つまりfの引数が一番多くなるってことでいいんだよね？
でもできないんです（＞＜）

やりかたおしえて♪

>>572
いまいち質問がはっきりしないので外しているかもしれないが・・・

この種の問題は他の言語から来た人などがよく出会う問題だ。
FORTRANでfunctionを引数とすることには、規格上のあいまいさの問題が
あることを知っておくとよい。

FORTRANでは副プログラムに変数を引き渡すとき、参照渡し（アドレス渡し、CallByReference)
になっている。一方、関数自体は値を値として持っているのかどうか規格で定められていない。
C言語やPASCALなんかでは値を返していると思う。（たしか、そうだったと、思うｗ）

だから、引数に関数を置いたときの挙動は処理系に依存している。
大昔では、HITAC暴走、FACOMコンパイラエラー、ACOS、VAXはOKだったか？

また、値が返ってくるような引数に（F90のINTENT（INOUT）やINTENT(OUT））
のような関数を置くのは明らかにおかしい。
======================================================================
さて、>>572のプログラムは、つづり等の間違いを除けば一応正しい。
ただ、引数の型をきっちり合わせていないとうまくいかない可能性がある。
つまり　S=f(3.0,g,4.0,4.0,4.0)　

また　g(x,h(x),d,e)　こういう形にしたい場合は、第二引数は普通の変数として
定義しておけばいいはず。　F90以降ならINTENT(IN)属性を付けておけば
コンパイラも多分おこらないだろう。

ま、一番確かなのはコンパイラのマニュアルを読むことだｗ

質問が、関数ポインターで引数リストも渡したい、というものならFortranではできない。
Fortran2003を待つしかない。

>>572のプログラム自体は、externalで関数を渡しているので問題ない。
ただ、EXTERNALは関数ポインタではないので、引数の数（形式）までは渡せない。

固定の引数の形式（引数リスト？）をもったものを渡せるだけ。

難しい解説ありがとうございますｗ
でもちんぷんかんぷんです＾＾；
f(x)=3x+g(x)+a+b
g(x)=x**2+ｈ(x)+c
h(x)=4x+d
こういうのをやりたいってことなんです（＞＜）

解説がよくわからないんですがCOMMONでa,b,c,dを定義しないとだめってことですか？

S=f(3,g,4,4,4,)

function f(x,g,a,b,c,d)
xeternal g
f=3x+g(x,h,c,d)+a+b

function g(x,h,c,d)
external h
g=x**2+h(x,d)+c

function h(x,d)
h=4x+d
これじゃぁﾀﾞﾒだったんです（＞＜）

>>575
たとえば、これならうまくいくはず。

program test
external g

S = f(3.0, g, 4.0, 4.0, 4.0)
PRINT *, S
END

function f(x,g,a,b,c)
external g
f=3+x+g(x,c)+a+b
print *, x, g(x, c), a,b,c
return
end

function g(x,c)
g=x**2+c
return
end

>>575
質問がうんこだと、回答もウンコしかもらえないよｗ

とりあえず、これがやりたいことか？動くよ。
program test
external g, h
S = f(3.0, g, h, 4.0, 4.0, 4.0, 1.0)
PRINT *, S
END

function f(x, g, h, a, b, c, d)
external g, h
f = 3.0 * x + g(x, h, c, d) + a + b
return
end

function g(x, h, c, d)
external h
g = x**2 + h(x, d) + c
end

function h(x, d)
h = 4.0 * x + d
end

ｳﾝｺでｺﾞﾒﾝｗ
でもうごかなかったんだよ・・・
も一度確認してみるね(・∀・)
普通3連コンボとかこういう風にやるもんなの？
もっとスマートな方法とかあるの？
でも関数の中に関数が定義できないってのがしょぼくね？fortran

もちろん、しょぼい。

>>579
関数の中に関数は定義できる。
引数の形が任意の関数を渡すことはできない。
数学でいう汎関数のようなものはできない。

できないのは、ほとんど必要とされていないから。

普通は、以下のようにサブルーチンは固定のままで話が済む。
program test
S = f(3.0, 4.0, 4.0, 4.0, 1.0)
END

function f(x, a, b, c, d)
f = 3.0 * x + g(x, c, d) + a + b
return
end

function g(x, c, d)
g = x**2 + h(x, d) + c
end

function h(x, d)
h = 4.0 * x + d
end

どもできましたｗ
function　fの引数にhを入れるのを忘れてましたｗ
でもひとつうんこなしつもんなんだけどぉ
functionの中に配列渡すにはどうしたら良いの？

real Un(100)
integer Ko
Ko=36
X=UNK(Un,Ko)
end

function UNK(Un,Ko)
real Un
integer Ko
UNK=Un(Ko)
end

>>582
ネタでなかったらそんな聞き方はしない筈

>>582
少しは前レス嫁ｗｗ　いいけどなｗｗ

受け取る側で配列であることを示せばいい。
function UNK(Un,Ko)
real Un(*)
UNK=Un(Ko)
end
たとえば、こういうの。配列の最後の次元については、＊で示してもいい。

ただ、普通はもう少し上品に配列の定義長を渡す。
多次元の場合は最後の次元以外でどうしても必要になる。
real Un(100)
Ko=36
X=UNK(100,Un,Ko)
end

function UNK(Mn, Un, Ko)
real Un(Mn)
UNK=Un(Ko)
end

FORTRAN77では配列の先頭要素のアドレスしか渡しておらず、配列のサイズに関する情報は
別個に明示的に渡さないといけない。いろいろ便利だが、バグの原因にもなるので、やれるなら
F90で書いてきっちりやったほうがいい。

けど、好きにやればいいさｗ

ども
access violationとかでる（＞＜）
変数宣言の順番変えたらうまくいったりする場合もあったり・・
これってバグってやつ？
f90でやったほうがいいの？それともCってやつのほうがいいの？
fortran以外だとdelphi pascalしかしらない・・・

つまり配列を渡す場合にはいったん外部で配列の値を入力した変数を渡してやれば良いってことジャン？？
めんどくせーーｗ
こんなところで10時間近く悩んでたなんてあほだなｗｗｗｗｗｗｗ
fortranってしょぼくね？ｗ

real Un(100)
integer Ko
Ko=36
Un=Un(ko)
X=UNK(Un)
end

function UNK(Un)
real Un
UNK=Un

あ　これじゃだめかU=Un(Ko)とかにしないとだめなんだな

fortranまじしょぼいｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗｗ
引数の順序変えたらまともに動かなくなったり動くようになったりｗｗｗ
なんでこんなに不安定なの？
使い物にならないジャン？

>>588
しょぼいのはFORTRANではない、お前だｗｗｗ

単に使い方を誤解しているだけだ。

引数の順番を変えてうまく行ったり行かなかったりするというのは、
じつはどっちにしろうまく行ってないということだｗ

PASCALができるなら、F90の方が向いている。Module言語的になっているから。

しかし、書き込みを見ている限りPASCALで出来るとも思われぬがｗうひょｗｗｗ

ごめんやっぱしょぼいのおれだったｗｗｗｗ
ってか72文字しかかけないんだって？
だったら７２文字以上書いた時点でエラー出せよｗ
コンパイルスンナｗ
しかも変数の宣言間違ってたらエラー出せよｗ

まあ、F77 がしょぼいのは確かだけどな。

>>590
FORTRAN７７には予約語とかはないし、空白とかをトークンとして区切るわけでもない。
コンパイラはまず空白を削除して頭から解釈してゆく。
したがって７２文字より先の部分が切り捨てられても、解釈が通ればエラーにはならない。

大体、コンパイラが厳しければ厳しいで文句を言うんだから、ゆるいのを喜べ！
FORTRAN６６なんかDO　LOOPから飛び出して、なんかやった後またそのLOOPの中に
戻っていけんるんだぞ！こんなことできる言語はなかなかないぜ！

それは利点じゃねーw

>>590
下品でない変数名にしようよ。
貴婦人が見ているかもしれないんだから。

>>594
KISS MY ASS！　F*CK　U！

FORTRANの変数名が６文字の制限下では、Four Letter Words がぴったりだぜ！