function val = ThirdOrderCorrection(s,v,u,m)

n   = length(u);
val = s^2;
k = 1:n;
for j=n+1:m
    a = 1./((j+k-1).*(j+k)/2-(k-1));
    val = val + (a*v)^2;
end
for j=n+1:m
    a = 1./((k+j-1).*(k+j)/2-(j-1));
    val = val + (a*u)^2;
end
val = sqrt(val);

return