def METHOD_NAME(ax, theta, f=None):
  dim = theta.shape[1]
  for i in range(dim):
    for j in range(i + 1, dim):
      if f:
        ax[i, j].scatter(
            theta[:, i], theta[:, j], marker='o', s=3, alpha=0.5, c=f)
      else:
        ax[i, j].plot(theta[:, i], theta[:, j], '.', markersize=3)
      ax[i, j].grid(b=True, which='both')
      ax[i, j].set_xlabel("F"+str(i))
      ax[i, j].set_ylabel("F"+str(j))