% plot1.m 96/07/14
%
% Program to plot simple graph of integral of surface
% 

    [X,Y]=meshgrid(1:.1:3);
    Z=2-.5*(X-2).^2-.5*(Y-2).^2;

    hold off;
    colormap('white');
    mesh(X,Y,Z);
    axis([0 4 0 4 0 4]);
    hold on
%    plot3([1 1],[1 3],[0 0],'w--'); 
%    plot3([3 3],[1 3],[0 0],'w--'); 
%    plot3([1 3],[1 1],[0 0],'w--'); 
%    plot3([1 3],[3 3],[0 0],'w--'); 
%    plot3([1 1],[1 1],[0 3],'w--');
%    plot3([1 1],[3 3],[0 3],'w--');
%    plot3([3 3],[1 1],[0 3],'w--');
%    plot3([3 3],[3 3],[0 3],'w--');
%    text(1.1,.5,0,'a');
%    text(3.0,.5,0,'b');
%    text(.4,1.2,0,'c');
%    text(.4,3.2,0,'d');
%    text(1.3,2.5,0,'R')
    text(3.5,2,1.5,'S');
    xlabel('x');
    ylabel('y');
    zlabel('z');
    arrow3([2;2;2],[2;2;3],.2);
    arrow3([1;1;1],[0;1;2],.2);
    arrow3([3;1;1],[3.3;1;2],.2);
    arrow3([1;3;1],[1;3;2],.2);
    hold off;
