00001 function y=ms(x,s,v,n) 00002 00003 myEps=sqrt(eps); 00004 00005 00006 00007 y=(fs(s+myEps*x,v,n)-fs(s,v,n))/myEps;