cop calculation and plot