Pr(B<Btarget)のみを返す単純なやつ

get.kobemat(fout, N = fout$input$N, nyear = fout$input$nyear,
  Btarget = 0, fmulti = seq(from = 0.3, to = 1, by = 0.1))