00001 function jac=numJacV(s,v,n)
00002
00003 N=n*n;
00004
00005 x=ones(N,1);
00006
00007 jac=zeros(N,N);
00008
00009
00010 for i=1:N
00011
00012 x=zeros(N,1);
00013 x(i)=1;
00014
00015 jac(:,i)=mv(x,s,v,n);
00016
00017 end
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029