if def_nr in [3, 4]:
                            if (Nx == 1):
                                # Joint density (TMd,TMm) given the Max and min
                                # Note the density is not scaled to unity
                                pdf[0, ts, tn] = fxind[0]  # *CC