EASolverで結果を出しても、数字だけだとピンときませんね。 そこで周期とその確かさを表すようにしてみました。 plot(d[m]/(N1+(1:(k+1))-1),1/dev[1:(k+1)],pch=16,col=3,grid(col=4),xlab="P") こういう行を、繰り返しの後にでも入れます。 複雑に見えま…