%
% source.m
%
function source(c,label,orient)
t = 0:.01:2*pi;
if strcmp(orient,'v')
   plot([c(1) c(1)],[c(2) c(2)-1],'k','linewidth',1.5)
   hold on
   plot(0.8*cos(t)+c(1),0.8*sin(t)+c(2)-1.8,'k','linewidth',1.5)
   plot([c(1) c(1)],[c(2)-2.6 c(2)-3.6],'k','linewidth',1.5)
   text(c(1),c(2)-1.8,label)%,'fontsize',24)
   text(c(1)+.25,c(2)-.75,'+')
else
   plot([c(1) c(1)+1],[c(2) c(2)],'k','linewidth',1.5)
   hold on
   plot(0.8*cos(t)+c(1)+1.8,0.8*sin(t)+c(2),'k','linewidth',1.5)
   plot([c(1)+2.6 c(1)+3.6],[c(2) c(2)],'k','linewidth',1.5)
   text(c(1)+1.6-length(label)/10,c(2),label)%,'fontsize',24)
   text(c(1)+2.5,c(2)+1,'+')
   %text(c(1)+1,c(2)+1,'-')
end
axis equal
return