How does the SAS code shown below look when translated into python?
/* Plot the K-L div versus lambda for a sequence of Poisson(lambda) models */
lambda = do(4, 16, 0.1);
KL = j(1, ncol(lambda), .);
do i = 1 to ncol(lambda);
poisPDF = pdf("Poisson", Levels, lambda[i]);
KL[i] = KLDiv(Proportion, poisPDF);
end;
title "K-L Divergence from Poisson(lambda)";
call series(lambda, KL) grid={x y} xvalues=4:16 label={'x' 'K-L Divergence'};