天津科技大でOctaveを教えてほしいとのことなので,調べました。早速,サンプルを実行。ロジスティック写像がなんと16行でできるんです。なおOctaveは音楽ソフトでもMIDI生成ソフトでもありません。
% ロジスティック写像: osim18.m
% x(n+1) = a*x(n)*(1-x(n))
% a: コントロールパラメタ
% x0: 初期値, N: 繰り返し数
function osim18(a,x0,N)
x(1) = x0;
printf("n \t x(n)\n");
for n = 1:N
x(n+1) = a*x(n)*(1-x(n));
printf("%d \t %f\n",n,x(n+1))
end
n = 1:N;
plot(n,x(n+1));
axis([1 N 0 1]);
xlabel(‘n’);
ylabel(‘X_n’);
これはきっとcastor先生がはまりますよ。
【Octave】
Octave
GNU Octave: 日本語マニュアル for 2.1
図 Octaveはgnuplotを呼び出すんです。