Как ограничить решения системы линейных уравнений? (linalg.solve)

Я решаю систему линейных уравнений с помощью функции linalg.solve из numpy. Хотелось бы поставить на эти решения определенные ограничения, например, чтобы все найденные x лежали в интервале от 0 до 1. Вот что я получаю без попытки поставить это ограничение (я не знаю как это сделать):

А хотелось бы иметь что-то типа: 0.25, 0.56, 0.82 ... или 0, 0.55, 0.79 и тд...