t = (0:0.01:2*pi);
a = log(2+t+t.^2);
b = exp(t).*(1+cos(3*t));
c = cos(t).^2+sin(t).^2;
d = tan(t).^(-1);
e = cot(t);
f = sec(t).^2+cot(t)-1;

subplot(2,2,1); % pod-plotovi, arg: br. redova, br. kolona, 
%pozicija sledeceg plota (popunjavanje sa leva na desno, odozgo na dole
plot(t,a);
subplot(2,2,2);
plot(t,b);
subplot(2,2,3);
plot(t,c);
subplot(2,2,4);
plot(t,d);