function [sol, val] = example1Fn(sol,options)
% Demonstration evaluation function used in the example.
% f(x)= 20+x + 10*sin(4*x)+8*cos(3*x)
%
% function [val,sol] = example1Fn(sol,options)
% 
% val - the fittness of this individual
% sol - the individual, returned to allow for Lamarckian evolution
% options - [current_generation]


%x=sol(1);

x1=sol(1);
x2=sol(2);
x3=sol(3);
x4=sol(4);
x5=sol(5);

gen2
%pointv(:)=pointvr(:)+pointvi(:)*i;      %變成復數
%pointv2(:)=pointv2r(:)+pointv2i(:)*i;    %變成復數
%dfrs(:,:)=dfrsr(:,:)+dfrsi(:,:)*i;        %變成復數
numVar=size(sol,2)-1;
nmsr=22;
for k=1:nmsr
  pointv(k)=pointvr(k)+pointvi(k)*i;      %變成復數
  pointv2(k)=pointv2r(k)+pointv2i(k)*i;   %變成復數
  for m=1:numVar
     dfrs(k,m)=dfrsr(k,m)+dfrsi(k,m)*i;    %變成復數
  end
end


for k=1:nmsr
ess(k)=0.0;
for j=1:numVar              %   !1 TO NXX
ess(k)=dfrs(k,j)*sol(j)+ess(k);		%!K步 散射場	  J+1
end
end
for k=1:2:nmsr	 %!計算測量相位差，K步計算相位差
   j=(k-1)/2+1;
	dcpointv(j)=atan2(imag(pointv(k)),real(pointv(k)))...
 	  -atan2(imag(pointv(k+1)),real(pointv(k+1)));
	djpointv(j)=atan2(imag(ess(k)+pointv2(k)),...
	  real(ess(k)+pointv2(k)))...
 	  -atan2(imag(ess(k+1)+pointv2(k+1)),...
       real(ess(k+1)+pointv2(k+1)));			 %!pointv2(i)
end
%pause



%for i=1:numVar
%  vlr(i)=0.0;
%end
%for i=1:numVar
%   for j=1:numVar
%      vlr(i)=vlr(i)+vrs(i,j)*sol(j);
%  end
%   vlr(i)=(vlr(i)-rsvol(i))^2;
%end
h1=0.0;
h2=0.0;
%for i=1:numVar
for k=1:nmsr/2
   h1=h1+(dcpointv(k)-djpointv(k))^2;
   h2=h2+dcpointv(k)^2;
end

%ag=sqrt(vlr/vr);
ag=sqrt(h1/h2);
val=1.0/ag;

%if ag<6.0
%val=6.0-ag;
%else
%val=0.0;
%end