clear all
close all

r=ones(20,1)*linspace(0,1,20);
p=r';
[r,p]=rotate2d(r,p,45,0.5,0.5);
rp=r+.5*p-1.5*r.*p;

pa=surf2patch(r,p,rp);

h=patch(pa);
alpha(.9)

set(h,'facecolor',[.5,.5,.5],'edgecolor',[.5,.5,.5])

grid on
daspect([1 1 .5]);
view(-66,18);
axis tight

camlight right
%camlight(-24,12);
lighting phong
xlabel('r')
ylabel('p')
zlabel('payoff')
print -djpeg rp

