【東6東1】F科 part6【ﾚｰｻﾞｰ】

403 ：電通太郎：2008/05/06(火) 20:36:19

c
c filename: f3.f
c contents: f(x) = exp(-(x**2))
c
program main
integer i
real*8 x, y, z1, z2, d
d = 0.1
x = -5.0
do 10 i=1, 100
x = x + d
y = exp(-(x**2))

c z1 = (6 - 2*x) / (6 + 4*x + x*x)
c Z2 = (6 - 4*x + x*x) / (6 + 2*x)

c write(6,20) x, y, z1, z2
write(6,20) x, y
10 continue
c 20 format(F6.2,3(1X,1PE11.4))
20 format(F6.2,1(1X,1PE11.4))
stop
end
c

>g77 f3.f -o f3

>./f3 > f3.dat

>cat f3.plt
set xzeroaxis
set yzeroaxis
set xrange[-3:3]
set yrange[-1:1]
set xtics axis
set ytics axis
plot "f3.dat" with lines

>gnuplot -persist f3.plt

	■ したらば のおすすめアイテム ■ ビューティフルコード - Brian Kernighan 2ch閉鎖騒動時にこんなことを言った人がいるそうです。 「なぜコードを最適化するのですか？」 「・・・そこにコードがあるから」彼らは静かに答えた。 開発者コメント：したらばのコードも買収当事は相当ひどかったぜｗｗ
	この欄のアイテムは掲示板管理メニューから自由に変更可能です。