% Data couples
tj = [  0  0.07 0.14 0.21 0.28 0.35 0.42 0.49 0.56 0.63 0.70];
Vj = [ 194 184  177  156  142  160  168  166  170  178  187];
% Trigonometric interpolation
n = length( tj ) - 1; % n=10, n+1 = 11
m = 100;
t_values = linspace( 0, 0.77, m*(n+1) );
V_values = interpft( Vj, m*(n+1) );
plot( tj, Vj,'o', t_values, V_values, 'r--' );
% Ejection fraction 
Ef = ( max( V_values ) - min( V_values ) ) / max( V_values )
%  Ef =
%      0.2714

% Correction of data V_7 
Vjc = Vj;
Vjc(7+1) = ( Vj(6+1) + Vj(8+1) ) / 2; %(Note: MATLAB index starts from 1)
% Trigonometric interpolation
Vc_values = interpft( Vjc, m*(n+1) );
hold on; plot( t_values, Vc_values, 'k-' );
legend('Data','Uncorrected','Corrected','Location','SouthEast')
% Ejection fraction (post-correction)
Efc = ( max( Vc_values ) - min( Vc_values ) ) / max( Vc_values )		
%  Efc =
%      0.2720