radial.m

to find stastical moments

function [xr,yr]=radial(x0,y0,x,y)
[xr,yr] =intline(x0,x,y0,y);