diff --git a/ShoulderCase/@Muscles/plot.m b/ShoulderCase/@Muscles/plot.m
new file mode 100644
index 0000000..61b8737
--- /dev/null
+++ b/ShoulderCase/@Muscles/plot.m
@@ -0,0 +1,37 @@
+function plot(obj)
+
+
+  SCase = obj.shoulder.SCase;
+  dataPath = fullfile(SCase.dataCTPath,'muscles');
+
+  rotatorCuffCrossSection = imread(fullfile(dataPath,'rotatorCuffCrossSection.png'));
+  IS = imread(fullfile(dataPath,'IS','mask.png'));
+  SC = imread(fullfile(dataPath,'SC','mask.png'));
+  SS = imread(fullfile(dataPath,'SS','mask.png'));
+  TM = imread(fullfile(dataPath,'TM','mask.png'));
+
+  leg = '';
+  imshow(rotatorCuffCrossSection);
+  hold on
+  if (max(IS,[],'all') > 0)
+      leg = [leg "IS"];
+      contour(IS,'color','g');
+  end
+  if (max(SC,[],'all') > 0)
+      leg = [leg "SC"];
+      contour(SC,'color','r');
+  end
+  if (max(SS,[],'all') > 0)
+      leg = [leg "SS"];
+      contour(SS,'color','b');
+  end
+  if (max(TM,[],'all') > 0)
+      leg = [leg "TM"];
+      contour(TM,'color','y');
+  end
+  if length(leg)
+      leg = leg(2:end);
+  end
+  legend(leg);
+  % saveas(fig,['./output/' ctFolder(i).name]);
+end