That is not an error, it is advice. You can improve your code speed by using
s1 = f(si);
Note: I suspect that your code is not correct for this case, but I do not understand it well enough to be sure.
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!