GS = {A -> 3.41, B -> 0.54, C -> 2.27, D -> 1.07};GS = {A -> 14.47, C -> 9.57, B -> 2.28, D -> 4.51}; PhaP[a_, b_] :=   1/2((A + B) b^2 + (C + B) a^2 +         Sqrt[((A - B) b^2 - (C - B) a^2)^2 + 4 (D + B)^2 a^2 b^2]);Muir[a_, b_] := (a^4 v^4 + (1 + q)a^2 b^2 v^2 h^2 + b^4 h^4)/(a^2 v^2 +         b^2 h^2);MPP = {v -> Sqrt[C], h -> Sqrt[A], q -> ((B + D)^2 + B(C - B))/(A(C - B))};WeP[a_, b_] := vp^2 (1 + 2 e b^4 + 2 d a^2 b^2);Th = {vp -> Sqrt[C], vs -> Sqrt[B], e -> (A - C)/(2 C),     d -> ((B + D)^2 - (C - B)^2)/(2 C (C - B))};Fow[a_, b_, x_] := (v^2 a^2 + h^2 b^2)(1 - x) +     x Sqrt[(v^2 a^2 + h^2 b^2)^2 + 2(q - 1) a^2 b^2 v^2 h^2/x];ParametricPlot[{a 180/Pi, 100        Abs[Sqrt[WeP[Cos[a], Sin[a]]/PhaP[Cos[a], Sin[a]]] - 1]} /. Th /.     GS, {a, 0, Pi/2}, PlotStyle -> AbsoluteDashing[{2}]];ParametricPlot[{a 180/Pi, 100          Abs[Sqrt[Muir[Cos[a], Sin[a]]]/Sqrt[PhaP[Cos[a], Sin[a]]] -               1]} /. MPP /. GS, {a, 0, Pi/2}, 	      PlotStyle->AbsoluteDashing[{5}]];ParametricPlot[{a 180/Pi, 100          Abs[Sqrt[Fow[Cos[a], Sin[a],1/2]]/Sqrt[PhaP[Cos[a], Sin[a]]] -               1]} /. MPP /. GS, {a, 0, Pi/2}];Show[{%,%%,%%%},Frame->True,FrameLabel->{"Phase Angle (degrees)","Relative Error (%)",None,None},PlotLabel->"Phase Velocity Error"];Display["junk_ma.eps",%,"EPS"];