I found another way. First I created a subplot with a legend, then placed the legend on the appropriate place on the plot, I created a M-code from the plot and got the position of the legend, then used the position for all subplots. . . .
end
leg1=legend(legend);
set(leg1,'Position',[p1 p2 p3 p4]);
j=j+5;
end