function [res] = FF(M,V,s1)
% This function is used by Steepest Descent (Algorithm 10.3)
% to define the function g given F.
syms('W','I','r');
W = zeros(1,M);
if M == 2
   r = 0;
   for I = 1:M
   s = s1{I};   
   g = inline(s,'y1','y2');
   r = r + g(V(1),V(2))^2;
   end;
elseif M == 3
   r = 0;
   for I = 1:M
   s = s1{I};   
   g = inline(s,'y1','y2','y3');
   r = r + g(V(1),V(2),V(3))^2;
   end;
elseif M == 4
   r = 0;
   for I = 1:M
   s = s1{I};   
   g = inline(s,'y1','y2','y3','y4');
   r = r + g(V(1),V(2),V(3),V(4))^2;
   end;
elseif M == 5
   r = 0;
   for I = 1:M
   s = s1{I};   
   g = inline(s,'y1','y2','y3','y4','y5');
   r = r + g(V(1),V(2),V(3),V(4),V(5))^2;
   end;
elseif M == 6
   r = 0;
   for I = 1:M
   s = s1{I};   
   g = inline(s,'y1','y2','y3','y4','y5','y6');
   r = r + g(V(1),V(2),V(3),V(4),V(5),V(6))^2;
   end;
else
   r = 0;
   for I = 1:M
   s = s1{I};   
   g = inline(s,'y1','y2','y3','y4','y5','y6','y7');
   r = r + g(V(1),V(2),V(3),V(4),V(5),V(6),V(7))^2;
   end;
end;
res = r;