clc
close all
clear all
z = importdata('C:\Users\amin\Desktop\New Text Document.txt');
[e s]=size(z);
Minx=min(z(:,2));
Maxx=max(z(:,2));
Miny=min(z(:,3));
Maxy=max(z(:,3));
R=input('please input your resolution number:');
m=floor((Maxx-Minx)/R)+1;
n=floor((Maxy-Miny)/R)+1;
b=1;
for i=0:m
    for j=0:n
      S(b,:)=[(i*R+Minx) (Maxy-j*R)];
    b=b+1;
    end 
end
[r t]=size(S);
p=1;
c=1;
for p=1:e
for i=1:r
        W1(c,p)=(S(i,1)-z(p,2))^2+(S(i,2)-z(p,3))^2;;
      c=c+1;
end 
c=1;
end
P=0;
Q=0;
b1=1;
[u l]=size(W1);
for j=1:u
   for i=1:e
    if (W1(j,i))^0.5<200
   P=(z(i,4)/W1(j,i))+P;
   Q=Q+1/W1(j,i);
   H(b1,1)=P/Q;
   b1=b1+1;
   end
   end
end
disp(H)