clear;clco=0.7;v=0.1;r=0.5;c=0.01;[pe,w]=meshgrid(linspace(0,1));dr=(1-r)*(1-o-w+pe)/(2*(1-o));de=(1-r)*(o*(1-o)+o*w-(2-o)*pe)/(2*o*(1-o))+r*(1-pe-w);z=(w-c).*dr+(pe-c).*de;z(~(w>=(2-o)/o*pe-1+o & w>pe+2*v-1+o & pe>=w))=nan;mesh(pe,w,z);