A..4 課題4

prog4.jl

# prog4.jl --- 漸化式 a_{n+1}=3a_n+1 で定義される数列
# Float64 は C の double と同じで 10 進法16桁,  10^{-308}〜10^{308}
# 実は BigInt を使うと正確に計算できる
function prog4(m)
  a=zeros(Float64,m) # a=zeros(BigInt,m)
  a[1]=1
  println("a_1=$(a[1])")
  for n=2:m
    a[n] = 3*a[n-1]+2
    println("a_$n=$(a[n])")
  end
end

prog4.jl

# prog4.jl --- 漸化式 a_{n+1}=3a_n+1 で定義される数列
# Float64 は C の double と同じで 10 進法16桁,  10^{-308}〜10^{308}
# 実は BigInt を使うと正確に計算できる
function prog4(m)
  a=1.0 # a=BigInt(1)
  println("a_1=$a")
  for n=2:m
    a = 3*a+2
    println("a_$n=$a")
  end
end