Setwindow([0,21], [0,23]);
Setscaling(0.8);  // 縦横比を縱0.8:橫1に

data=[[ 0,   3.9*10^6], [ 1,   5.3*10^6], [ 2,   7.2*10^6], [ 3,   9.6*10^6],..
      [ 4,  12.9*10^6], [ 5,  17.1*10^6], [ 6,  23.2*10^6], [ 7,  31.4*10^6],..
      [ 8,  38.6*10^6], [ 9,  50.2*10^6], [10,  62.9*10^6], [11,  76.0*10^6],..
      [12,  92.0*10^6], [13, 106.5*10^6], [14, 123.2*10^6], [15, 132.1*10^6],..
      [16, 152.3*10^6], [17, 180.7*10^6], [18, 205.1*10^6], [19, 226.5*10^6]];
data2=list()
for i=1:length(data)/2
  data2($+1)=[data(2*i-1), data(2*i)/10^7];
end
//seq=Listplot(data2);

//Windisp(seq);

N0=3.9*10^6;
mu=197*10^6;
r=0.3134;

function N=N(t);
  N=mu/(1+(mu-N0)/N0*exp(-r*t))/10^7;
endfunction

G=Paramplot('[t, N(t)]', 't=[0,21]', 'N=300');
L=Plotdata('mu/10^7', 'x=[Xmin(), Xmax()]');  // Gの漸近線
L1=Plotdata('5', 'x=[Xmin(), Xmax()]');   // 補助線
L2=Plotdata('10', 'x=[Xmin(), Xmax()]');  //
L3=Plotdata('15', 'x=[Xmin(), Xmax()]');  //

//Windisp(G);

Openfile('tex/tick_mark_sample.tex');
Beginpicture('0.5cm');
// 座標軸のラベルの設定
Setax('', '\text{year}', 's4',..
      '\text{population}\,(\times 10^6)', 'n2', '$1790$', 's2w');

Drwline(G, 1.5);
Dashline(L, 1, 1);
Dottedline(L1, L2, L3, 1, 1);  // 補助線
Setpen(5);
Drwpt(data2);
Setpen(1);

// 目盛りを付ける
for i=1:22
  Vtickmark(i, string(10*i));  // 縦軸の目盛り
end

for i=0:19
  Htickmark(i, ' ');  // 横軸の目盛り
end
for i=0:9  // ひとつ飛ばしで数値を入れる
  Expr([2*i+1,0], 's2', string(1800+i*20));
end

Expr([4,mu/10^7], 's2', '\mu=197\times 10^6');
Expr([16,8], 'c', '\displaystyle N=\frac{\mu}{1+e^{-(\gamma t+C)}}');
Endpicture(1);
Closefile();

Setscaling(1);
