function area = calcarea(rad)
diam = rad * 2;
area = pi * (diam/2)^2;
end

