10/05/23 17:10:04
>>146
毎回最初からと言っているのは、READではなくて、補間位置を決める、2つのデータに挟まれている位置を
求める部分の話。
>>107での2重DO-LOOP構造では、DO内でK=1にしているので毎回X(0)から探し始めている
DO I=0,INT(X(MP-1)-X(0))
XX=X(0)+real(I)
K=1
do while ((X(K-1)-XX)*(X(K)-XX) > 0.0d0)
K=K+1
end do
..YY=SUB2(G,H,X,Y,XX,K,MP)
end do
..
end do
一応>>132では(Xの単調増加を仮定して)外に出してある。
>>146
読み込んだデータを延々と配列にため込んでいるので、毎回現在のポインタ位置のような意味でNを渡す
必要があるのではないかと思う。