clc
clear
A=zeros(4,36);

for i=1:36
    A(1,i) = (i-16)*0.1;
end


for i=1:36
    x=A(1,i)
    A(2,i) = sin(x);
    A(3,i) = cos(x);
    A(4,i) = (sin(x))^2;
end

subplot(3,1,1)
plot(A(1,:), A(2,:), 'r')


mtlb_hold on;
subplot(3,1,2)
plot(A(1,:), A(3,:), 'g')



subplot(3,1,3)
plot(A(1,:), A(4,:), 'b')


subplot(3,1,4)
plot(A(1,:), A(4,:), 'b')