%2009年人口分布
       %1-10歲人口總數
       clear;
       clc;
    i=[0.13730000000000   0.13720000000000   0.14200000000000   0.14480000000000];
    j=[ 
        0.99959890780000
   0.99997102180000
   0.99997135740000
   0.99997448160000

  



];
  for  k=2:4
        l(1,1)= 0.99959890780000;
      l(1,k)=l(1,k-1)*j(k,1);
  end
  l
 a1=sum(l.*i);
  %10-90歲人口總數
  
  qq1=[
       0.99959890780000
   0.99997102180000
   0.99997135740000
   0.99997448160000
   0.99998170060000
   0.99998226840000
   0.99998260300000
   0.99998212380000
   0.99998447160000
   0.99998922520000
   0.99998009820000
   0.99998673980000
   0.99998428860000
   0.99997356700000
   0.99997243400000
   0.99997387780000
   0.99997593280000
   0.99996979080000
   0.99996828320000
   0.99997615940000
   0.99997291260000
   0.99996753000000
   0.99996793340000
   0.99996369780000
   0.99997257140000
   0.99996786980000
   0.99996776760000
   0.99991722540000
   0.99995410820000
   0.99994150420000
   0.99993725960000
   0.99994318080000
   0.99993362860000
   0.99993231660000
   0.99993418120000
   0.99992510140000
   0.99991386480000
   0.99992392120000
   0.99991113320000
   0.99990958880000
   0.99990408060000
   0.99990129540000
   0.99992022320000
   0.99991731540000
   0.99989850020000
   0.99989549060000
   0.99988952980000
   0.99988415020000
   0.99985624660000
   0.99984350100000
   0.99983634080000
   0.99983597420000
   0.99981427820000
   0.99983449820000
   0.99982848560000
   0.99980826480000
   0.99981171560000
   0.99983210000000
   0.99978389340000
   0.99978713440000
   0.99975925860000
   0.99976886940000
   0.99973009780000
   0.99970850820000
   0.99969022980000
   0.99965806860000
   0.99966612880000
   0.99962281660000
   0.99959544260000
   0.99955510720000
   0.99949930320000
   0.99951838440000
   0.99949040340000
   0.99951042140000
   0.99951945300000
   0.99948134840000
   0.99955297420000
   0.99951216400000
   0.99954049680000
   0.99956356140000
   0.99952490460000
   0.99957870480000
   0.99954648400000
   0.99964882240000
   0.99966972380000
   0.99970738640000
   0.99974535840000
   0.99980299960000
   0.99979688060000
   0.99984942880000
   0.99934253520000
  ];

 uu1=zeros(1,87);
  for  nn1=1:87
   uu1(1,nn1)=qq1(nn1+1,1)*qq1(nn1+2,1)*qq1(nn1+3,1)*qq1(nn1+4,1);
  end

  p=[
      0.14035372189598
   0.13805039410699
   0.13308175801683
   0.13768820387250
   0.14901047873717
   0.15405451882324
   0.15053928788807
   0.16704175142020
   0.16951941339914
   0.17502569554059
   0.19653063295855
   0.18737977488504
   0.20842977859113
   0.21555058525905
   0.23529639188826
   0.26605386569084
   0.24241262331688
   0.21407599362998
   0.21479497444874
   0.17259954544947
   0.15140776921230
   0.14785231604185
   0.15326163191736
   0.18186773412613
   0.16357569761409
   0.16198718298575
   0.17208666855821
   0.17100200458666
   0.16490607161051
   0.18411974402934
   0.19466763997089
   0.21356465850849
   0.22498230085943
   0.23544386455759
   0.24590604512357
   0.27061583111449
   0.25031178516027
   0.27646956177907
   0.22318456553506
   0.25164107982666
   0.25805563142396
   0.24684861473294
   0.29165536871310
   0.21896854896901
   0.12145575814297
   0.15503330337872
   0.14130915037006
   0.18488851542135
   0.21009411353305
   0.19168747819822
   0.20976900197092
   0.20387127882014
   0.18924432889072
   0.18968956809073
   0.16091185220288
   0.16104752363142
   0.15509363091662
   0.13263009124613
   0.13104495720945
   0.11913554259837
   0.11579132017412
   0.10618980837349
   0.09892648429136
   0.09617837884843
   0.09631681400081
   0.09820861237015
   0.08003844298868
   0.08883347966051
   0.08476972933919
   0.08229177232353
   0.08359696610728
   0.07204974173070
   0.07502536268646
   0.06347571202392
   0.05634753981631
   0.05860172614307
   0.04583273678731
   0.04560878398228
   0.04001384062527
   0.03275085398104
   0.03333423329672  
   0.02715591062403
   0.02309610284956
   0.02025577348327
   0.01641991851384
   0.01357958914755
   0.00987902576760
   
  ];
  size(uu1);
  size(p')
  q=uu1.*p'
  sum(q)+a1
  ss=[(l.*i)';q']