maximize
  logexp_sum(spline_sum(matrix_parameters_vars, matrix_data))
constraint: <= 10
  linear(matrix_linear)

calculate
Point: point_problem_1
  logexp_sum(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))
  logistic(spline_sum(matrix_parameters_vars, matrix_data, matrix_data_knots))