abs = [0:0.1:30];
 
a=A(abs);
b=B(abs);
cb = comp(b);
 
hold on
plot(abs,a);
plot(abs,b);
plot(abs,union(a,b));
hold off
 
figure
hold on
plot(abs,a);
plot(abs,b);
plot(abs,inter(a,b));
hold off
 
figure 
hold on
plot(abs,b);
plot(abs,cb);
hold off
 
figure 
hold on
plot(abs,a);
plot(abs,cb);
plot(abs,inter(a,cb));
hold off
 
function ua = A(x)
    ua=[1:300];
    for i= 1:length(x)
        n=x(i);
        if n<7
            un = 0;
        else
            un = (n-7)/(n-5);
        end
        ua(i)=un;
    end
end
 
function ub = B(x)
    ub = 1./(1+(x-12).^2);
end
 
function AuB = union(A,B)
    if length(A) ~= length(B)
        ME = MException("Union:wrongArraySize","The sizes of the two arrays do not match");
        throw(ME);
    end
    AuB = max(A,B);
end
 
function AnB = inter(A,B)
    if length(A) ~= length(B)
        ME = MException("Union:wrongArraySize","The sizes of the two arrays do not match");
        throw(ME);
    end
  
    AnB=min(A,B);
end
 
function cA = comp(A)
    cA=1-A;
end