mupad2casadi.m

function [casadi_expr] = mupad2casadi(mupad_expr,mupad_var,casadi_var)
  s = size(mupad_expr);
  assert(numel(s)<=2, 'Tensors not supported. Make a vector or matrix from the mupad expression first.');
  mupad_expr = mupad_expr(:);
  
  matlabFunction(symfun(mupad_expr,mupad_var),'File','temp_fun.m');
  clear temp_fun
  args = vertsplit(casadi_var);
  
  casadi_expr = temp_fun(args{:});
  %casadi_expr = reshape(casadi_expr,s);
end