function f= fSSR(a, Sm, v0m)
v0p= a(1)*Sm.^3./(a(2)+Sm.^3);
f= sum((v0m-v0p).^2);