def METHOD_NAME(self, n_elem_curve, defor=False):
    curve = np.zeros((n_elem_curve, 3))
    t_vec = np.linspace(0, 2, n_elem_curve)
    for i in range(n_elem_curve):
        t = t_vec[i]
        for idim in range(3):
            if defor:
                polyfit, _, _ = algebra.get_polyfit(self.coordinates_def, self.ordering)
            else:
                polyfit, _, _ = algebra.get_polyfit(self.coordinates_ini, self.ordering)
            polyf = np.poly1d(polyfit[idim])
            curve[i, idim] = (polyf(t))
    return curve