function [hf] = PlotScapHumRhythm(HumAng,ScapAng)
%PlotScapHumRhythm - plots scapular rotation with respect to humeral
%position to display scapulohumeral rhythm
% Dr Martin Warner, University of Southampton, 2014

%% Ensuring data has been trime normlised

if length(HumAng)~= 101;
    error('Error: Humeral kinematic data are not time normalised');
end

if length(ScapAng)~= 101;
    error('Error: Scapula kinematic data are not time normalised');
end


%% Plotting scapular rotation with respect to humeral evelation angle

hf = figure; hold on
plot(HumAng,ScapAng)
xlabel('Time (percentage)')
ylabel('Degrees')

