clc; clear all; close all
A=[1 2 3; 4 5 6; 7 8 9]
%Y=rot90(rot90(rot90(rot90(A))))
%Y1=fliplr(A)
%Y2=flipud(A)
% Y3=reshape(A',1,9)
% diag(Y3)
% triu(A,1)
% tril(A,-1)
% A=[1 2 3]
% B=[1 2 3]
% dot(A,B)
% A*B'

% x=0:0.1:1
% y=x
% subplot(2,2,1)
% plot(x,y,'color',[1 0 0])
% subplot(2,2,2)
% plot(x,y,'color',[1 0 0])
% clf
% 
% subplot(2,2,4)
% plot(x,y,'color',[1 0 0])
% [x,y]=ginput 
%axis equal
% comet(x,y)

[th,fi]=meshgrid(0:0.01:pi,0:0.02:2*pi);
x=1*sin(th).*cos(fi);
y=1*sin(th).*sin(fi);
z=1*cos(th)
surf(x,y,z,x)
shading interp
axis equal