function [H,Xp]=houghradon(i,theta)

[H,Xp]=radon(i,theta);  