Как построить допустимую область задачи нелинейной оптимизации в Matlab

У меня есть нелинейная задача оптимизации с двумя переменными с ограничением неравенства.

min  f(x,y)
s.t. g(x,y) >= 0
          x >= 0    % y free

с

a = 0.9;
c = 20; 
h = 5;
j = 1;
l = 1;
m = 10;
s = 3;

f = c.* integral(@(z) (z-x).*normpdf(z,m,s),x,inf,'ArrayValued',true) + h*(y+x-j*m-l*x*(1- normcdf(x,m,s))-l.*integral(@(z) z.*normpdf(z,m,s),-Inf,x));
g = normpdf(x,m,s).*normcdf(y,j*m,sqrt(j)*s) + integral( @(z) normpdf(x-z,m,s).*normcdf(y+z,j.*m,sqrt(j).*s),0,x,'ArrayValued',true) - a;

Я хочу построить возможную область на трехмерном графике, чтобы лучше понять пространство решений. Каков наилучший или наиболее разумный способ сделать это в Matlab?

Я попытался следовать документации в https://de.mathworks.com/help/optim/ug/solve-constrained-nonlinear-optimization-problem-based.html, но моя проблема плохо представлена ​​на контурном графике.

Есть ли лучший способ сделать это или с адаптацией к 3D-графике?