%q1
r=rand(9)
b=r>mean(mean(r))
q1=sum(sum(b))

%q2
a=linspace(-1,1);
b1=sin(a)./(a.^2+1);

plot(a,b1)

fplot('sin(a)./(a.^2+1)',[-1 1])

ezplot('sin(10*x)-sqrt(2*x)')

% q3 : second func
figure(2)
fplot('log(p+2)',[-1 1])

% q4: third and fourth funcs
fplot('[sin(p.^2) cos(p.^2)]',[0 2*pi])
% using hold
fplot('sin(p.^2) ',[0 pi])
hold on
fplot(' cos(p.^2)',[0 pi])
hold off

% q5
subplot(2,2,1)
ezplot('sin(a)./(a.^2+1)',[-1 1])
subplot(2,2,2)
ezplot('log(p+2)',[-1 1])
subplot(2,2,3)
ezplot('sin(p.^2) ',[0 pi])
subplot(2,2,4)
ezplot(' cos(p.^2)',[0 pi])


% 6
diag(diag(((magic(4)+rand(4)+12)\[4 0 2 -1]')*ones(1,2)+rand(4,2)))
A_logic = logical( diag(diag(((magic(4)+rand(4)+12)\[4 0 2 -1]')*ones(1,2)+rand(4,2))) );
n = sum(sum( A_logic ))
 



% q1 extra
% num of cells > avg of their col
a=pascal(5)
a(3,1)=12
a(4,3)=3
a(3,4)=3
b=mean(a)
m=repmat(b,5,1)
sum(sum(a>m))