%Forming Phi_i for each component 
function [tmpPhi]=tPhi(xy,LSgridx,LSgridy,p)
    at_tht = atan2(xy(4)-xy(2),xy(3)-xy(1));
    L = ((xy(3)-xy(1))^2+(xy(4)-xy(2))^2)^0.5;
    xcrd = LSgridx - xy(1)/2 - xy(3)/2;
    ycrd = LSgridy - xy(2)/2 - xy(4)/2;
    tmpPhi = 1 - (2*L^-1*(xcrd*cos(at_tht) + ycrd*sin(at_tht))).^p - (2*(-xcrd*sin(at_tht) + ycrd*cos(at_tht))/xy(5)).^p;
end