function F=equation(landa)
global N m t x0  q alpha GAMMA shape B
     o = ones(1,length(t));
      rx(:,:,1) = t*o - (t*o)';  % signed distance matrix 
      r(:,:,1) = abs( rx(:,:,1));  %  distance matrix
       B(:,:,1)= mq(r(:,:,1),shape);
     for k=2:m+1
      rx(:,:,k)= (t*q^(k-1))*o - (t*o)';
      r(:,:,k)= abs(rx(:,:,k));
       B(:,:,k)= mq(r(:,:,k),120*N);
     end
for i=2:N  
    s1=0;
    for j=1:N
    s1=s1+landa(j)*B(i,j,1);  
    end
V=0;
R=0;
for n=1:m
    clear D
    D=(1-q^(n))/(1-q^(n-1+alpha));
    for k=1:m
        D=D*(1-q^(n+k))/(1-q^(n-1+alpha+k));
    end
    s2=0;
   
    for j=1:N
        s2=s2+landa(j)*B(i,j,n);
    end
    s2=s2^2;   
    V=V+q^(n-(3/2)*n+(1/2))*D*s2;
end
R=((t(i)^(alpha-(3/2)))*(1-q)/GAMMA)*V;
F(i)=s1-x0-R;
end
    i=1;
    s1=0;
    for j=1:N
s1=s1+landa(j)*B(i,j,1);
end
  F(1)=s1-x0;
end

