%Task 4 - Determine frequencies at which amplitude is 1/2 maximum
    disp(' '),disp(''),disp('Task 4: Determine frequencies at which amplitude is 1/2 maximum')
    poly = [1/Hz^4,0,(c^2-2)/Hz^2,0,1-1/(maxA/2)^2];
    disp(poly)
    format short e
    disp('The frequencies in rpm:')
    rt = roots(poly)*60;
    disp(rt)