clc;clearrho = 0:0.1:10;theta = [0:0.1:2*pi,0];[rho1,theta1] = meshgrid(rho,theta);x=rho1.*cos(theta1);y=rho1.*sin(theta1);z = rho1.^1.*(exp(sin(theta1))+exp(cos(theta1)))/2;contourf(x,y,z)colormap(jet(16))colorbaraxis offaxis equal