plot_kitagawa    a simple calculation and a plot
