新歓中止

ロード新歓は懸念とされていた雨で中止になっていまったし、さらに朝走ろうって言ってたのに起きられなかったのが悔やまれる。そして次の日、朝練があった。

6：30に起きてご飯を高速で食べ、ギターだけ持って自転車に乗る。スタジオの前まで行ってみると誰もいなかった。なんだーみんな起きられなかったのかと思い、コンビニで立ち読み開始。「京大M1物語」を読んだあたりで、もしかしたら誰か来てるかも？と思って外に出てみると、丁度タカフが到着していた。中に入ってちょっと待ったけど誰も来ないのでシールドとか借りて入った。

用意していたらてっちゃんが来た。バイトがありえなかったと盛んに言っていた。WASA合わせたりしてたけど相変わらず俺が大変になるとギターを弾かなくなるので、3ピースはきつい気がした。てっちゃんがギターを弾いてベースなしでやったときが一番まともだった気がする。その後ひろきが来て合わせたりして、9：00になって終わる。なんでいきなり歌詞忘れるんだろうなあ。覚えたと思ったのに。あと他の曲はどうなったんだろう。

終わってから家に帰って仮眠したけど相変わらず眠いままだった。時間になって学校へいく。昼ごはんはまたもや高速でご飯だけ。

3・4限は実験。今日は材料のやつで、最初に顕微鏡で組織観察をして（イスが変なやつだった）、それから部屋を移動してビッカース硬さ試験をやった（イスが狭い）。俺の誤ったイメージでは硬さの試験ってもっとガツン！ってぶつけるのかと思ってたけど、実際は結構まったりだった。最後に回転曲げ疲労のやつを見て終わる。ぬるい実験だった。

帰ってきて眠くなってしまって寝て、日が変わるころに起きた。よくない傾向だけどしょうがない。起きてからひたすらFortranで断面積形状からマッハ数、静温、静圧を求めるプログラムを書いていた。何回か頑張ったら成功したけどスロート部でマッハ数が1に達しないのでどう考えてもおかしいことが分かったけれど、一応結果が出たので満足して寝た。下にコードを載せるけど、もし誰か読めたらどこがおかしいのか教えて。最初はニュートン法の初期値が問題なのかと思ったけどそうじゃないっぽい。

parameter(eps=1.0d-12,shoki=1.0)
dimension a(2370)
f(x)=x**6+15*(x**4)+75*(x**2)-0.2162*a(j)*x+125
fd(x)=6*(x**5)+60*(x**3)+150*x-0.2162*a(j)
g(x)=x-f(x)/fd(x)
t(x)=(1500+273.15)/(1+0.2*(x**2))
p(x)=0.6/((1+0.2*(x**2))**(3.5))
open(10,file='areaonly.d',status='old')
do i=1,2370
read(10,*)a(i)
end do
close(10)
open(11,file='MachNumber.d',status='unknown')
open(12,file='Temp.d',status='unknown')
open(13,file='Pres.d',status='unknown')
do j=1,2370
x2=shoki
times=0
100   times=times+1
x1=x2
x2=g(x1)
if(abs(x1-x2)>eps)then
goto 100
else
y=j*0.1
write(11,*)x2,y
write(12,*)t(x2),y
write(13,*)p(x2),y
end if
end do
close(11)
close(12)
close(13)
end