English: plot generated using the following maxima commands:

array (check, 99, 99); for m:1 thru 100 do for n:1 thru 100 do check[m-1,n-1] : 0 $

array (result, 99, 99); for m:1 thru 100 do ( for n:1 thru 100 do ( if (check[m-1,n-1] = 0) then ( tmp : 1-1/365**(m+n)*sum(sum(stirling2(m,i)*stirling2(n,j)*product(365-k,k,0,i+j-1),j,1,n),i,1,m), result[m-1,n-1] : tmp, result[n-1,m-1] : tmp, check[m-1,n-1] : 1, check[n-1,m-1] : 1, if (mod(n, 10) = 0 or mod(m, 10) = 0) then display(n,m) )))$

plot3d (result[floor(x)-1,floor(y)-1], [x,1,100],[y,1,100], [zlabel,""], [mesh_lines_color,false], [elevation,0], [azimuth,0], color_bar, [grid,99,99], [ztics,false], [color_bar_tics,1])$