clc
clear
x=[70;115;150;110;90;180];
y=[140;115;150;170;190;210];
z=[115.4;123.1;113.8;110.5;107.2;131.78];
X=110;
Y=150;
i=length(x);
A=0.0;
for j=1:i
xi=x(j);
yi=y(j);
zi=z(j);
distance=sqrt((X-xi).^2+(Y-yi).^2);
A=sum(1/distance+A);
end
weight=0.0;
distance=0.0;
Z=0;
for j=1:i
xi=x(j);
yi=y(j);
zi=z(j);
distance=sqrt((X-xi).^2+(Y-yi).^2);
weight=(1/distance)/A;
Z=sum(weight.^zi+Z)
end