def plotpairs(features, k, m):
    TD_ind=np.where(diagnosis==0)[0]
    ASD_ind=np.where(diagnosis==1)[0]

    betw1_TD = features[TD_ind, k]
    betw1_ASD = features[ASD_ind, k]
    betw2_TD = features[TD_ind, m]
    betw2_ASD = features[ASD_ind, m]

    plt.scatter(betw1_TD, betw2_TD, color='red', marker='o')
    plt.scatter(betw1_ASD, betw2_ASD, color='blue', marker='o')

    plt.show()