plot Rmax and c50 in a loop in plotfits.m