Самое близкое, что я могу придумать, это:
function q55322965
% Evaluate the equation (half domain!)
t = linspace(-pi,0,50);
c = exp(-cos(t));
% Turn vectors into a mesh:
[TT,CC] = meshgrid(cos(t),c);
% Clear all points that are above the curve:
CC(CC > c) = NaN;
% Fill in the rectangle between the chart and zero:
CC(end+1,:) = 0;
TT(end+1,:) = TT(end,:);
% Plot:
figure(); mesh(TT,CC,CC,'FaceColor','interp','EdgeColor','interp'); view([0,90]);
Что дает:
Если вы хотите получить менее зазубренный вид при построении этого метода, вы можете увеличить разрешение в t
. Например, если мы используем 500
вместо 50
в linspace
, мы получим: