روی این دستور :
v1(S)=v(S)-tol;

این خطا رو میده :
??? Error using ==> minus
Matrix dimensions must agree.


مهندس فکر کنم همینطوری رفع اشکال کنیم بتونیم جواب بگیریم!