x=1:6;
y=4:5:100;
Z=[31171,29724,27356,75532,143344,147435,113016,79544,43183,30086,18252,11323,10698,6973,3652,1870,952,460,181,50,3;
37569,35805,32825,80779,163233,175792,141087,97463,49917,34216,19457,13862,13442,7566,3538,1628,779,395,164,49,5;
29047,25166,18817,80665,154653,146284,111750,66396,29970,19994,13177,8948,7899,4910,2432,1114,509,244,96,21,3;
76842,64177,46193,480551,800303,563641,333781,165926,69940,48645,28671,16290,13773,10484,7362,4588,2382,1122,348,102,8;
52710,44361,31808,294039,452028,353441,223506,115277,49576,34785,20761,12056,10272,8029,5541,3468,1893,937,319,83,18;
4319,4141,3298,19908,38428,31350,21961,12897,6052,3910,2074,1277,1051,759,440,234,135,58,25,6,0];
Z=Z(:,1:20);
[X Y]=meshgrid(x,y);
surf(X,Y,Z')