syms x y
f(x,y) = sqrt(x*y);
g(x,y) = -sqrt(x*y);
fsurf(f, [0 4 0 4])
hold on
fsurf(g, [0 4 0 4])
hold off
pbaspect([1 1 1])
camorbit(0, 0, 'data', [0 0 1]);
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 0 1]); 
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 1 0]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 0 1]); 
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [1 0 0]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [1 1 0]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 1 1]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [1 0 1]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 0 1]); 
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 1 0]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 0 1]); 
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [1 0 0]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [1 1 0]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 1 1]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [1 0 1]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 0 1]); 
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 1 0]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 0 1]); 
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [1 0 0]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [1 1 0]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [0 1 1]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end
for angle = linspace(0, 360, 36)
       camorbit(1, 0, 'data', [1 0 1]);
       drawnow; % Update the plot
       frames(fix((angle/360)*36) + 1) = getframe(gcf); % Capture the frame
end