function outM = Myobjective(x,y,title1,n)
    
    parfor i=1:n
%           M = csvread(title1(i,1:25),0);
%           rv1 = M(:,1);
%           zv1 = M(:,2);
%           thetaGDv = M(:,3);
% 
%           nz1 = 20;
%           nr1 = 10;
%           rm1 = reshape(rv1,nr1,nz1);
%           zm1 = reshape(zv1,nr1,nz1);
%           thetaGD = reshape(thetaGDv,nr1,nz1);
%           r1 = rm1(:,1);
%           z1 = zm1(1,:);
          
%           poolobj = gcp;
%           addAttachedFiles(poolobj,{title1(1,1:25),title1(2,1:25),title1(3,1:25),title1(4,1:25),title1(5,1:25),title1(6,1:25)})
          
          M = csvread(title1(i,1:25),0);
          rv1 = M(:,1);
          zv1 = M(:,2);
          thetaGDv = M(:,3);

          nz1 = 20;
          nr1 = 10;
          rm1 = reshape(rv1,nr1,nz1);
          zm1 = reshape(zv1,nr1,nz1);
          thetaGD = reshape(thetaGDv,nr1,nz1);
          r1 = rm1(:,1);
          z1 = zm1(1,:);
          
         [thetaG,CL,CG,r,z,nr,nz] = modelThetaG(x,y(i,:));
          thetaGMt = griddata(r,-z,thetaG',r1,z1);
          thetaGM = thetaGMt';


          out(i) = norm(thetaGM -thetaGD);
    end
    
    outM = mean(out);