2013-10-27 22:24:23 +00:00
|
|
|
def euler(y, a, b, h)
|
|
|
|
|
a.step(b,h) do |t|
|
|
|
|
|
puts "%7.3f %7.3f" % [t,y]
|
|
|
|
|
y += h * yield(t,y)
|
2013-04-10 16:57:12 -07:00
|
|
|
end
|
|
|
|
|
end
|
|
|
|
|
|
2013-10-27 22:24:23 +00:00
|
|
|
[10, 5, 2].each do |step|
|
|
|
|
|
puts "Step = #{step}"
|
|
|
|
|
euler(100,0,100,step) {|time, temp| -0.07 * (temp - 20) }
|
|
|
|
|
puts
|
|
|
|
|
end
|
