S_dis:=ln(N!/(N-M)!/M!);

plot3d(S_dis,N=0..100,M=0..N);

S_con:=N*ln(N)-(N-M)*ln(N-M)-M*ln(M);

plot3d(S_con,N=1..100, M=0..N,grid=[100,100]);

S_con_big:=10000*ln(10000)-(10000-M)*ln(10000-M)-M*ln(M);

plot(S_con_big,M=0..10000);