Histogram computation from gradient image.

Hello everyone

I want to compute the histogram using image gradient.let say that I have an image

Img=imread(myimage);

[gx gy]=edge(Img,'sobel');

Mag=sqrt(gx,gy);

Dir=atan2(gy,gx);

now I want to compute the orientation bins using 12 orientations each covering 30degree.

how can I do it.After that I need to sum(Img,Dir) from 1 to number of bins 12.

Any help is appreciated.

[Walter Roberson]

sqrt() cannot take two numeric inputs.

[Algorithms Analyst]

I will correct it. The question is about histogram computatiom