\documentstyle[11pt,epsf]{article}
\input{psfig}
\textwidth 14.0cm
\hoffset -2cm
\headsep 1.5cm
\parindent 1.2em
%\renewcommand{\baselinestretch}{2.0}
\catcode`@=11
\def\seceqaa{\@addtoreset{equation}{section}
\def\theequation{A\arabic{equation}}}
\def\seceqbb{\@addtoreset{equation}{section}
\def\theequation{B\arabic{equation}}}
\def\seceqcc{\@addtoreset{equation}{section}
\def\theequation{C\arabic{equation}}}
\catcode`@=12
\begin{document}
\begin{center}
%\hfill hep-th/0108038\\
%\hfill IP/BBSR/2001-20\\
{\Large \bf 
The Peculiarity of a  Negative Coordinate Axis in Dyonic Solutions 
of Noncommutative ${\cal N}=4$ Super Yang-Mills}
\vskip 0.1in
{Aalok Misra\\ 
Harish Chandra Research Institute,\\
Jhusi, Allahabad - 211 019, India\\email: aalok@mri.ernet.in}
\vskip 0.5 true in
\end{center}

\begin{abstract}

We show that in the neighborhood of a negative coordinate axis,  the $U(1)$
sector of the static dyonic solutions to the noncommutative $U(4)$ ${\cal N}=4$
Super Yang-Mills (SYM) can be consistently decoupled from the $SU(4)$ 
to {\it all orders in the noncommutativity parameter}. 
We show the above decoupling in two ways. First, we show 
the noncommutative dyon being the same as the commutative dyon, is a 
consistent solution to noncommutative equations of motion in the abovementioned
region of noncommutative space.  Second, as an example of 
decoupling of a non-null $U(1)$  sector, 
we also obtain a family of solutions with nontrivial 
$U(1)$ components for all components of the gauge 
field in the same region of noncommutative space. 
In both cases, the $SU(4)$ and $U(1)$ components separately 
satisfy the equations of motion.
\end{abstract} 

\section{Introduction}

Nonperturbative solutions in Super Yang-Mills (SYM) theories have been an 
important area of work. More recently, static dyonic solutions 
in $SU(N), N\geq3,4$ ${\cal N}=4$ SYM corresponding to planar and 
non-planar (respectively) string junctions have been
obtained in \cite{SI1,SI2} (See also references therein). Also, solitonic
solutions in  noncommutative SYM have been considered by others before 
(See \cite{hata1,hata2} and references therein). Now, the gauge group
$SU(N)$, in general is not allowed in noncommutative gauge
theories \cite{ncggps} 
(See \cite{wess} though, where the gauge parameter in
noncommutative gauge tranformations is allowed to depend on the usual $SU(N)$
gauge paramater and gauge fields. The explicit form of the map to all orders
in the noncommutativity parameter [which is what we are interested in,
in this paper] is quite involved and has not been worked out.). 
One of the questions that we address in this work is whether it is possible
to construct solutions to the equations of motion of a noncommutative
supersymmetric gauge theory in which one can decouple the $U(1)$ components
from the $SU$ components of the fields at least in some region of the
noncommutative space, {\it to all orders in the noncommutativity paramater}. 
There are the following two additional motivations 
for this work. The equality of the angular momentum and (null) 
quadrupole moment of the commutative and noncommutative monopoles of (S)U(2) 
$N=2$ SYM  up to $O(\theta^2)$ was shown in \cite{noncomcomeq} 
and in all of space. This makes it natural to pose the following question:
is it possible to construct solutions of a noncommutative supersymmetric
gauge theory that are characterized by quantities that receive no
noncommutativity corrections {\it to all orders in the noncommutativity
paramater}? Finally, it is of obvious interest to construct solutions
that are valid up to all orders in the noncommutativity paramater.

In this paper, we show that for noncommutative $U(4)$ ${\cal N}=4$ SYM, 
the static dyonic solutions in the neighborhood of a negative coordinate axis
of the noncommutative space, are such that one can decouple 
the $U(1)$ from the $SU(4)$, and construct
a 6-parameter family of noncommutative dyonic solutions (for a non-null 
$U(1)$ decoupling). Even though, the decoupling has been shown for a very
limited region of the noncommutative space, the reason the solutions are of
interest is because they have the same electric and magnetic charges as  and 
energy differing infinitesimally from
the commutative dyonic solutions of \cite{SI1}, {\it to all orders in the
noncommutativity parameter}. This points to the peculiar nature of the
the negative $x^3$ axis when constructing dyonic solutions to noncommutative
${\cal N}=4$ SYM, which requires further investigation.

The paper is organized as follows. In section 2, we discuss the basics
needed for getting static dyonic solutions in commutative ${\cal N}=4$ 
Super Yang-Mills (SYM). In section 3, 
we discuss the modifications in the techniques of section 2 to
get static dyonic solutions of noncommutative 
${\cal N}$=4 SYM. It has two subsections: (a short) 
3.1 on decoupling of a null $U(1)$ from the $SU(4)$, and 3.2 on
decoupling of a non-null $U(1)$ from the $SU(4)$. In section 4, we 
give a summary of and a discussion on the results obtained in the paper, 
and indicate possible future directions for work.

\section{Basics}

In this section, we give a brief review of the techniques given
in \cite{SI1} to obtain dyonic solutions in commutative $SU(4)$ ${\cal N}=4$
SYM.

It is advisable to write the equations of motion as 
if one were working in curved space-time.  The reason is that for 
solutions of the equations of motion, it is convenient to go to 
the stereographic coordinate system (See Fig 1: $z=e^{i\phi}tan{\theta\over2}$
) in which the metric is not
constant. The equations of motion for commutative $SU(4)$ ${\cal N}=4$ SYM
are given by:
\begin{eqnarray}
\label{eq:comeoms}
& &
{D_i}(\sqrt{g}g^{ii_1}g^{jj_1}{F}_{i_1j_1})=i\sqrt{g}\sum_J
[{D^j}{\Phi^J},{\Phi^J}];\nonumber\\
& &
{D_i}(\sqrt{g}g^{ij}{D}_j
{\Phi^J})=\sqrt{g}[{\Phi}^I,[{\Phi}^I,
{\Phi}^J]].
\end{eqnarray}
The component-form of the above equations is given in \cite{SI1}
which we won't be repeating here. Following the close parallel
with Skyrme  models\cite{projskr}, one uses the harmonic-map anzatz to
solve the equations of motion - the Higgs and gauge fields are 
assumed to have the following form:
\begin{eqnarray}
\label{eq:AzPhidefs}
& & {A_z}= {A_z}^a T^a 
=-(1+c_1-c_2-c_0)(r)(\partial P_0) P_0 + (c_2-c_1)(r)[P_1,\partial P_1] 
+ (1-c_2(r))[P_2,\partial P_2],\nonumber\\
& & {A_r}=0;\nonumber\\
& & {\Phi^I}={\Phi}^{I,a} T^a =
(b^I_2+b^I_0+b^I_1)(r)(P_0
-{1\over4}{\bf 1}_4) + (b^I_2+b^I_1)(r)(P_1-{1\over4}{\bf 1}_4)
+b^I_2(r)(P_2-{1\over4}{\bf 1}_4),\nonumber\\
& & 
\end{eqnarray}
$T^a$ being the generators of $SU(4)$, $I$ indexing the number of scalars
that are turned on, and for this paper taking values 1,2 and 3,
$c_{0,1,2}(r)$ and $b^I_{0,1,2}(r)$ being the `profile functions' of the 
gauge field and Higgs respectively, and $P_i$ the orthogonal set of 
${\bf CP}^3$ sigma model projectors defined as:
\begin{eqnarray}
\label{eq:defsprojs}
& & P_0={ff^\dagger\over f^\dagger f};\nonumber\\
& & P_j={(\bigtriangleup^j f)(\bigtriangleup^j f)^\dagger\over
|\bigtriangleup^j f|^2},\ {j=1,2},
\end{eqnarray}
where 
\begin{eqnarray}
\label{eq:Nablafdef}
& & \bigtriangleup f\equiv \partial f - f{f^\dagger\partial f
\over|f|^2};
\nonumber\\
& & f\equiv(1,\sqrt{3} z, \sqrt{3} z^2, z^3)^t
,  
\end{eqnarray}
and $\bigtriangleup^2f=\bigtriangleup(\bigtriangleup f)$.
They can be thought of as hermitian maps from $S^2$ to ${\bf CP}^3$. 
The explicit forms of the projectors $P_i$s can be worked out
using Mathematica or Maple.

The substitution of the harmonic map anzatz into the equations of motion
(written in components) either results in an equation of motion being
identically satisfied, or a system of 
non-linear second-order coupled differential equations
involving the Higgs and gauge field profile functions
that can be solved using techniques referred to in \cite{SI1}.
The solutions corresponding to (non-)planar string networks  
(for non-planar, assuming an infinitesimal non-planarity corresponding to
an infinitesimal deformation from the BPS nature of the planar solution)
were then obtained in \cite{SI1}.

\section{Noncommutative dyonic solutions}

In this section we derive a class of static noncommutative dyonic solutions
assuming  $r-z$ and $r-{\bar z}$ noncommutativity only. Thus, the
star product of two functions $f(r,z,{\bar z})$ and $h(r,z,{\bar z})$ will
be:
\begin{equation}
\label{eq:rzrzbarncdieuf}
f(r,z,{\bar z})*h(r,z,{\bar z})
=e^{{i\over2}[\Theta^{rz}(\partial_r^1\partial^2 - \partial_r^2\partial^1)
+\Theta^{r{\bar z}}(\partial_r^1{\bar\partial}^2 
- \partial_r^2{\bar\partial}^1)]}f(r_1,z_1,{\bar z}_1)
h(r_2,z_2,{\bar z}_2)|_{1=2}.
\end{equation}
We show that using stereogrpahic coordinates, in the neighborhood
of the South Pole, or equivalently the neighborhood of the negative $x^3$ axis,
one can consistently decouple in the equations of motion,
the $U(1)$ from the $SU(4)$ components of the Higgs and the
gauge fields.

The noncommutative equations of motion are given by:
\begin{eqnarray}
\label{eq:noncomeoms}
& &
\hat{D_i}*(\sqrt{g}*g^{ii_1}*g^{jj_1}*\hat{F}_{i_1j_1})=i\sqrt{g}*\sum_J
[\hat{D^j}*\hat{\Phi^J},\hat{\Phi^J}]_*;\nonumber\\
& &
\hat{D_i}*(\sqrt{g}*g^{ij}*\hat{D}_j*
\hat{\Phi^J})=\sqrt{g}*[\hat{\Phi}^I,[\hat{\Phi}^I,
\hat{\Phi}^J]_*]_*.
\end{eqnarray}
As we are considering only $r-z, r-{\bar z}$-noncommutativity,
the projectors of \cite{SI1} are still valid here. 
We assume the following anzatz for the gauge field $\hat{A_z}$, 
$\hat{A_r}$ and the Higgs field $\hat{\Phi^I}$:
\begin{eqnarray}
\label{eq:AzPhihatdefs}
& & \hat{A_z}= \hat{A_z}^a T^a +i\epsilon F_z(\Theta)
g_z(r)\hat{{\cal A}^0_z}{\bf 1}_4;\nonumber\\
& & \hat{A_r}=i\epsilon 
F_r(\Theta)g_r(r)\hat{{\cal A}^0_r}{\bf 1}_4;\nonumber\\
& & \hat{\Phi^I}=\hat{\Phi}^{I,a} T^a+i\epsilon F_\Phi(\Theta)g_\Phi^I(r)
\hat{{\cal\phi}^0}{\bf 1}_4,\nonumber\\
\end{eqnarray}
where $F_{z,r,\Phi}(\Theta)$ are functions of the noncommutativity
parameter $\Theta$ which vanish as $\Theta$ is set to zero and are bounded
as $\Theta\rightarrow\infty$, and
\begin{eqnarray} 
\label{eq:anztzSU}
& & \hat{A_z}^a
T^a=-(1+c_1-c_2-c_0)(r)*(\partial P_0) P_0 + (c_2-c_1)(r)*[P_1,\partial P_1] 
+ (1-c_2(r))*[P_2,\partial P_2],\nonumber\\
& & \hat{\Phi}^{I,a}T^a=
(b^I_2+b^I_0+b^I_1)(r)*(P_0
-{1\over4}{\bf 1}_4) + (b^I_2+b^I_1)(r)*(P_1-{1\over4}{\bf 1}_4)
+b^I_2(r)*(P_2-{1\over4}{\bf 1}_4)
,\nonumber\\
& & 
\end{eqnarray}
$c_{0,1,2}(r)$ and $b^I_{0,1,2}(r)$ being the `profile functions' of the 
gauge field and Higgs respectively as in (\ref{eq:AzPhidefs}).  
We have thus taken an anzatz in which the $U(4)$ 
fields have infinitesimal $U(1)$ components. We will
work up to $O(\epsilon)$. Hence, the noncommutative $U(1)$ field
strengths are the same as the commutative $U(1)$ field strengths.
There is no star product in the anzatz (\ref{eq:AzPhihatdefs}) in the $U(1)$ 
components of the fields because as argued
below, the star products can be dropped in the neighborhood of the South 
Pole in the $SU(4)$-components of the fields, and we assume
similar behavior for the $U(1)$ and $SU(4)$ components. The anzatz
for the $U(1)$ components that we choose in (\ref{eq:U1anztz}) below,
is such that if one were to put in 
star product between the radial and angular functions, then
it (the star product) can not be dropped in the neighborhood of the
South Pole.

As {\it expanding the star product in (\ref{eq:anztzSU}) will involve taking
derivatives of ($\partial) P_i$}, using the explicit 
forms of the projectors that can be worked out using Mathematica or
Maple,  one sees that 
{\it one can drop $(\partial)\partial P_i$ relative to $(\partial) P_i$ 
as one approaches the neighborhood of the South Pole}, as the former
is suppressed relative to the latter by terms of $O({1\over|z|^m}), m>0$ 
that vanish as $|z|\rightarrow\infty$ (i.e. in the 
neighborhood of the south pole).  {\it One thus has for the $SU$
part of $\hat{A_z},\hat{\Phi}$ the same expression as in
\cite{SI1}.} 
One can argue similarly that the * products in (\ref{eq:noncomeoms}) can also
be dropped (in the neighborhood of the South Pole).
We now discuss the decoupling from $SU(4)$, 
of a null $U(1)$ in {\bf 3.1} and non-null $U(1)$ in {\bf 3.2}.

\subsection{Decoupling of null $U(1)$}

Using (\ref{eq:AzPhihatdefs}), for null $U(1)$ components,
one can require consistently that the $SU(4)$ and $U(1)$ 
components, separately, satisfy the equations of motion, 
and one can show there is no mixing between the $SU(4)$ and $U(1)$ sectors.  
Hence, one can consistently set the $U(1)$ components to zero.
The entire analysis of \cite{SI1} to show the connection between (non-)planar 
stable (non-)BPS string junctions and dyon solutions, then goes through 
even in the neighborhood of the South Pole of the abovementioned 
noncommutative space if one notices and uses the following 
transformation that takes one from the North Pole to 
the South Pole (and is also a symmetry of the  equations of motion of the 
$SU(4)$ part of the solutions, as given in reference two of \cite{SI1}):
\begin{equation}
\label{eq:sym}
c_0\leftrightarrow c_2;\ b^I_0\leftrightarrow b^I_2.
\end{equation}
One also uses the fact that if $\hat{\Phi}$ is a solution then
$-\hat{\Phi}$ is also a solution to the equations of motion.

\subsection{Decoupling of non-null $U(1)$}

Here we discuss in detail, the decoupling of non-null $U(1)$ components 
of the Higgs and the gauge fields from their $SU(4)$ components in the 
neighborhood of the South Pole of the noncommutative space.
To show the decoupling of a non-null $U(1)$ from $SU(4)$,
we choose the following anzatz for the $U(1)$ components (of the 
Higgs and gauge fields):
\begin{eqnarray}
\label{eq:U1anztz}
& & \hat{{\cal\phi}^0},\hat{{\cal A}_{z}^0}=
Det[(P_0-a_{\Phi,z} {\bf 1}_4)(P_1-M_{\Phi,z})(P_2-M_{\Phi,z})]+
{\rm cc}
\nonumber\\
& & \hat{A_r^0}= z Det[(P_0-a_{r} {\bf 1}_4)(P_1-M_{r})(P_2-M_{r})]+{\rm
    cc}.\nonumber\\
\end{eqnarray}
where
\begin{equation}
\label{eq:M_idef}
M_i=\left(\begin{array}{cccc}\\
0&b_i&0&0\\
c_i&0&0&0\\
0&0&0&e_i\\
0&0&f_i&0\\
\end{array}\right),
\end{equation}
where $i\equiv r,z,\Phi$, and the parameters $a_i, b_i, c_i,  e_i$ and $f_i$
have been introduced to get non-zero determinants using the projectors
and some constant matrices; the parameters  get constrained later.
We now assume that one approaches the South pole along any path very close to
$\phi(\equiv$ azimuthal angle)$={\pi\over2}$. This implies that $Re z$ is
finite, but $Im z$ approaches infinity. Hence, one sees that 
$\partial{\bar\partial}\hat{{\cal A}_r^0}$ (as shown in  
(\ref{eq:laphatArzinf}))  and $(\partial-{\bar\partial}) 
(\hat{\phi}^0,\hat{{\cal A}_z^0})$ 
(as shown in  (\ref{eq:diffzinf})) are real. 
If 
\begin{equation}
\label{eq:cond1}
b_i=f_i,\ c_i=e_i,
\end{equation}
then  $\partial\biggl((1+|z|^2)^2(\partial+{\bar\partial})(\hat{\phi}^0,
\hat{{\cal A}_z^0})\biggr)$ 
(as shown in (\ref{eq:derz^4}))
will be finite and real as $Im\ z\rightarrow\infty$. These considerations
become useful when one solves the equations of motion for the $U(1)$ components
as done in (\ref{eq:U1eqns1}) and (\ref{eq:U1eqns2}) in the neighborhood of
the South Pole.

From power (of $|z|$ or $Im z$) counting 
arguments using expressions given in the  appendix, one 
can require consistently that the $SU(4)$ and $U(1)$ 
components, separately, satisfy the equations of
motion, and one can show there is no mixing between the $SU(4)$ and 
$U(1)$ sectors. As the decoupling for non-null $U(1)$ is non-trivial,
to illustrate the point, we consider the example of $\hat{F_{rz}}$.
We see that up to $O(\epsilon)$:
\begin{equation}
\label{eq:decplngod1ex}
\hat{F_{rz}}^{U(4)}=\hat{F_{rz}}^{SU(4)}
+\hat{F_{rz}}^{U(1)}-i[\hat{A_r}^{U(1)},\hat{A_z}^{SU(4)}]_*,
\end{equation}
Now,
\begin{eqnarray}
\label{eq:decplgod2ex}
& & [\hat{A_r}^{U(1)},\hat{A_z}^{SU(4)}]_*=i\Theta^{rz}\Biggl(
\partial_r\hat{A_r}^{U(1)}\partial\hat{A_z}^{SU(4)}-
\partial\hat{A_r}^{U(1)}\partial_r\hat{A_z}^{SU(4)}\Biggr)
\nonumber\\
& & 
+i\Theta^{r{\bar z}}\Biggl(
\partial_r\hat{A_r}^{U(1)}{\bar\partial}\hat{A_z}^{SU(4)}-
{\bar\partial}\hat{A_r}^{U(1)}
\partial_r\hat{A_z}^{SU(4)}\Biggr)
+O(\Theta^2).
\end{eqnarray}
Using expressions in the appendix, one sees that:
\begin{eqnarray}
\label{eq:decplgod3ex}
& & \partial_r\hat{A_r}^{U(1)}\equiv O(1),\ \partial\hat{A_r}^{U(1)},
{\bar\partial}\hat{A_r}^{U(1)}\equiv O(1);\nonumber\\
& & \partial\hat{A_z}^{SU(4)}\equiv O({1\over|z|^3})\times M,\
{\bar\partial}\hat{A_z}^{SU(4)}\equiv O({1\over|z|^4})\times
{\bf 1}_4,\
\partial_r\hat{A_z}^{SU(4)}\equiv O({1\over|z|^2})\times 
M;\nonumber\\
& & \partial_r\hat{A_z}^{U(1)}\equiv O(1),\
\partial\hat{A_z}^{U(1)},{\bar\partial}\hat{A_z}^{U(1)}\equiv
O({1\over|z|^2})({\rm relevant\ 
for\ }[\hat{A_r}^{U(1)},\hat{A_z}^{U(1)}]_*),\nonumber\\
& & 
\end{eqnarray}
where
\begin{equation}
\label{eq:defM}
M\equiv\left(\begin{array}{cccc}\\
0&1&0&0\\
0&0&1&0\\
0&0&0&1\\
0&0&0&0\\
\end{array}\right).
\end{equation}
Thus
\begin{equation}
\label{eq:decplgod4ex}
\hat{F_{rz}}^{U(1)}\equiv O(1);\hat{F_{rz}}^{SU(4)}
\equiv O({1\over|z|^2})\times M,
\end{equation}
and in Table 1, we show which terms can be dropped relative to
which terms as $Im z\rightarrow\infty$.
\begin{table} [htbp]
\centering
\caption{Consequences of power counting done in (\ref{eq:decplgod3ex})}
\begin{tabular}{|c|c|} \hline
Term that can be dropped & relative to \\ \hline
$\partial_r\hat{A_r}^{U(1)}\partial
\hat{A_z}^{SU(4)}$ & $\hat{F_{rz}}^{SU(4)}$ \\
$\partial\hat{A_r}^{U(1)}\partial_r\hat{A_z}^{SU(4)}$ &
$\hat{F_{rz}}^{SU(4)}$\\
$\partial_r\hat{A_r}^{U(1)}{\bar\partial}\hat{A_z}^{SU(4)}$
& $\hat{F_{rz}}^{U(1)}$\\
${\bar\partial}\hat{A_r}^{U(1)}\partial_r\hat{A_z}^{SU(4)}$&
$\hat{F_{rz}}^{SU(4)}$\\ \hline
\end{tabular}
\end{table}
Similar reasoning can be used to argue that terms of higher order
in $\Theta$, can also be dropped relative to the commutative
counterparts. Hence, as $Im z\rightarrow\infty$,
\begin{equation}
\label{eq:decplgod5}
\hat{F_{rz}}^{U(4)}\sim \hat{F_{rz}}^{SU(4)}+\hat{F_{rz}}^{U(1)},
\end{equation}
showing the stated decoupling.  

By requiring that the $SU$ components of the Higgs 
and gauge fields still satisfy the commutative
equations of motion (having dropped the star product at the 
South pole), the $U(1)$ components should satisfy these equations:
\begin{eqnarray}
\label{eq:U1eqns1}
& & {dg_z\over dr}(\partial-{\bar\partial})\hat{{\cal A}^0_z}-
2g_r\partial{\bar\partial}\hat{{\cal A}^0_r}
=0;\nonumber\\
& & {d^2g_z\over dr^2}
\hat{{\cal A}^0_z} - g_r\partial\hat{{\cal A}^0_r} 
+{g_z\over2r^2}
\partial\biggl((1+|z|^2)^2(\partial+{\bar\partial})\hat{A_z}
\biggr)=0
\nonumber\\
& & {1\over r^2}{d\over dr}(r^2{dg_\Phi^I\over dr})\hat{{\cal\phi}^0}
+{(1+|z|^2)^2\over r^2}g_\Phi^I
\partial{\bar\partial}\hat{{\cal\phi}^0}=0.
\end{eqnarray}
The * product has been dropped as it will generate $O({1\over{|z|^m}})$-
corrections (relative to $O(1)$),
which can be dropped in  the neighborhood of the 
south pole. We now discuss  a family of dyonic
solutions to (\ref{eq:U1eqns1}).

One gets (setting $c_r=e_r$ and $b_r=f_r$ for convenience):
\begin{eqnarray}
\label{eq:U1eqns2}
& & {dg_z\over dr}
+2g_r{[(b_rc_r)^3  
(2c_r - b_r)]\over(b_zc_z)^3(c_z-b_z)}
=0
\nonumber\\
& &  2{d^2g_z\over dr^2} (-1+a_{z})a^3_{z}
(b_{z}c_{z})^4 -
g_r(-1+a_r)a_r^3(b_rc_r)^4
+{
g_z\over r^2}(3(-1+a_z)a_z^3b_z^2c_z^2(b_z-c_z)^2)
=0\nonumber\\
& &  2{1\over r^2}{d\over dr}(r^2{dg_\Phi^I\over dr})(-1+a_{\Phi})a^3_{\Phi}
(b_{\Phi}c_{\Phi})^4
\nonumber\\
& & 
-{6\over r^2}|z|^2g_\Phi^I(-1+a_{\Phi,z})a_{\Phi}^3
(c_{\Phi}^4b_{\Phi}^2+b_{\Phi}^4c_{\Phi}^2)=0.
\end{eqnarray}
By eliminating $\hat{{\cal A}^0_r}$ from the first two equations, one gets the
following differential equation for $g_z$:
\begin{eqnarray}
\label{eq:diffeqgz}
& & 
{d^2g_z\over dr^2} 
+{dg_r\over dr}{(-1+a_r)a_r^3(b_rc_r)^4\sqrt{3}
(-1+a_z)a_z^3[(b_zc_z)^3(c_z-b_z)]\over
4(-1+a_z)^2a^6_z(b_zc_z)^4[b_r^3 c_r^3 
\sqrt{3} (2c_r - b_r)]}
\nonumber\\
& & 
+
{g_z\over r^2}(3(-1+a_z)a_z^3{(b_z-c_z)^2\over2 b_z^2c_z^2})=0
\end{eqnarray}
The equation (\ref{eq:diffeqgz}) is of the form
\begin{equation}
\label{eq:bessleqn}
r^2 y''(r) + \alpha r^2 y'(r) + \beta  y(r)=0,
\end{equation}
the solution to which (using maple) is given by:
\begin{eqnarray}
\label{eq:besselsolnnotZ}
& & 
y(r)=
{\it C_1}\sqrt{r}
{\it I}_{{1\over2}\sqrt{1-4\beta }}({1\over2}\alpha r)
{e^{-{1\over2}\,\alpha r}}\nonumber\\
& & +
{\it C_2}
\sqrt {r}{\it I}_{-{1\over2}\sqrt{1-4\beta }}
({1\over2}\alpha r){e^{-{{1\over2}\,\alpha r}}}
\end{eqnarray}
for ${1\over2}\sqrt{1-4\beta }\not\in{\bf Z}$and
\begin{eqnarray}
\label{eq:besselsolnZ}
& & y(r)=
{\it C_1^\prime}\sqrt{r}
{\it I}_{{1\over2}\sqrt{1-4\beta }}({1\over2}\alpha r)
{e^{-{1\over2}\,\alpha r}}\nonumber\\
& & +
{\it C_2^\prime}
\sqrt {r}{\it K}_{{1\over2}\sqrt{1-4\beta }}
({1\over2}\alpha r){e^{-{{1\over2}\,\alpha r}}}
\end{eqnarray}
for ${1\over2}\sqrt{1-4\beta }\in{\bf Z}$. Lets consider the two cases 
separately. 

For the former, to make sure that the solution goes to zero as
$r\rightarrow
\infty$, one sets:
\begin{equation}
\label{eq:condCs}
C_1=-C_2.
\end{equation}
Even though this is a singular solution (at $r=0$), as we shall see,
it is a finite energy solution.

The derivative of the solution involving $I_{\pm\nu}(\alpha r)$ is given by:
\begin{eqnarray}
\label{eq:derBesselI}
& & {d\over dr}
\Biggl(e^{-{\alpha r\over2}}\sqrt{r}\biggl(
I_{{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})-I_{-{1\over2}
\sqrt{1-4\beta}}({\alpha r\over2})\biggr)
\Biggr)\nonumber\\
& & 
={e^{{-\alpha r\over2}}\over4\sqrt{r}}\biggl(\alpha 
r I_{-1+{1\over2}\sqrt{1-4\beta}}(
{\alpha r\over2}) +\alpha r I_{1+{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})
+2(-1+\alpha r)I_{-{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})\nonumber\\
& & -\alpha r I_{-1-{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})
-\alpha r I_{1-{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})
-2(-1+\alpha r)I_{{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})\biggr).
\end{eqnarray}
For the purpose of getting finite energy, as we will very shortly
see, one requires to impose that $r^2({dy(r)\over dr})^2\rightarrow0$
as  $r\rightarrow\infty$, 
and that it is finite at $r=0$.
Using the series expansions around $r=0$ and the asymptotic expansions
for $I_\nu(r)$, the former is satisfied identically, 
and the latter is satisfied if one imposes:
\begin{equation}
\label{eq:constbeta}
\sqrt{1-4\beta}<1 \Leftrightarrow\ 0<\beta<{1\over4}.
\end{equation}

The  solution for $g_\Phi^I(r)$ that one gets is: 
\begin{equation}
\label{eq:gPhisoln}
g_\Phi^I(r)=0.
\end{equation}
Hence, the final result for the $U(1)$ components of the Higgs and
gauge fields, in the neighborhood of the South Pole (approaching it
along $\phi$ close to ${\pi\over2}$) are:
\begin{eqnarray}
\label{eq:final1}
& & \hat{A}^0_z(r,z)=i \epsilon
f_z(\Theta)\sqrt{r}{e^{-{1\over2}\,\alpha r}}\biggl(
{\it I}_{{1\over2}\sqrt{1-4\beta }}({1\over2}\alpha r)-
{\it I}_{-{1\over2}\sqrt{1-4\beta }}({1\over2}\alpha r)\biggr)
\nonumber\\
& & \times
Det[(P_0-a_z{\bf 1}_4)(P_1-M_r)(P_2-M_r)] +
{\rm cc});\nonumber\\
& & \hat{A^0_r}=i \epsilon f_r(\Theta)\nonumber\\
& & \times{e^{{-\alpha r\over2}}\over4\sqrt{r}}\biggl(\alpha r
I_{-1+{1\over2}\sqrt{1-4\beta}}(
{\alpha r\over2})
+\alpha r I_{1+{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})
+2(-1+\alpha r)I_{-{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})\nonumber\\
& & -\alpha r I_{-1-{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})
-\alpha r I_{1-{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})
-2(-1+\alpha r)I_{{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})\biggr)
\nonumber\\
& & \times(z 
Det[(P_0-a_z{\bf 1}_4)(P_1-M_r)(P_2-M_r)] +
{\rm cc});\nonumber\\
& &
\hat{\Phi^0}=0,
\end{eqnarray}
where 
\begin{eqnarray}
\label{eq:defs}
& & \beta\equiv(3(-1+a_z)a_z^3{(b_z-c_z)^2\over2
  b_z^2c_z^2})<0;\nonumber\\
& & \alpha\equiv{(-1+a_r)a_r^3(b_rc_r)^4\sqrt{3}
(-1+a_z)a_z^3[(b_ze_z)^3(c_z-b_z)]\over
4(-1+a_z)^2a^6_z(b_zc_z)^4[b_r^3 c_r^3 
\sqrt{3} (2c_r - b_r)]};\nonumber\\
& & 0<a_z<1;\ c_{z,r}=e_{z,r};\
b_{z,r}=f_{z,r}.
\end{eqnarray}

Lets consider the energy of the above solution. It is given by (we
follow the conventions of \cite{SI1}):
\begin{eqnarray}
\label{eq:energy1}
& & E=-{1\over4\pi}\int r^2 dr {dz d{\bar z}\over(1+|z|^2)^2}\
 Tr\biggl[(\hat{D_i}*\hat{\Phi})^2
+{1\over2}(\hat{F_{ij}})^2\ _*\biggr]\nonumber\\
& & =E_{SU(4)}+E_{U(1)}.
\end{eqnarray}
But, in the neighborhood of the South Pole, where one can drop the * product,
\begin{eqnarray}
\label{eq:energy2}
& & E_{U(1)}
=-\int r^2 dr {dz d{\bar
    z}\over{(1+|z|^2)^2}}\biggl[(\partial_r
\hat{{\cal A}^0}_z -\partial\hat{{\cal A}_r^0})^2+
(\partial_r
\hat{{\cal A}^0}_z +{\bar\partial}\hat{{\cal A}_r^0})^2\nonumber\\
& & +((\partial+{\bar\partial})\hat{{\cal A}_z^0})^2\biggr]\geq0.
\end{eqnarray}
Hence, in the neighborhood of the South Pole, 
even though one is dealing with non-BPS solutions, one gets the
following inequality:
\begin{equation}
\label{eq:energy3}
E_{U(4)}\geq E_{SU(4)},
\end{equation}
with saturation of the inequality occuring for null $U(1)$ component.
Substituting (\ref{eq:final1}) in the $Im\ z\rightarrow\infty$ limit into
(\ref{eq:energy2}), and integrating in the neighborhood of the South
Pole as described earlier, i.e., around $\phi={\pi\over2}$ and
$\theta\in[\pi-\epsilon,\pi], \lim_{\epsilon\rightarrow0}$, 
one gets an infinitesimal result, provided one gets  a finite result 
from the $r$ integration. The $r$-integral should 
converge to get a finite-energy solution.  One sees that one has to evaluate:
\begin{eqnarray}
\label{eq:rintfinite}
& & E_{U(1)}\sim
\int_0^\infty dr r^2 {e^{-\alpha r}\over16r}\biggl(\alpha r
I_{-1+{1\over2}\sqrt{1-4\beta}}(
{\alpha r\over2})
+\alpha r I_{1+{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})
+2(-1+\alpha r)I_{-{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})\nonumber\\
& & -\alpha r I_{-1-{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})
-\alpha r I_{1-{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})
-2(-1+\alpha r)I_{{1\over2}\sqrt{1-4\beta}}({\alpha
  r\over2})\biggr)^2.\nonumber\\
\end{eqnarray}
We have not been able to evaluate the above integral. But the integrand
for $0<\beta<{1\over4}$ is analytic for all values of
$r\in[0,\infty)$, and it vanishes at $r=0$ and $r\rightarrow\infty$. We
hence assume that it is convergent corresponding to a finite energy
solution. We have verified the same numerically using mathematica for
several values of $\beta<{1\over4}$.

Let us now consider the case ${1\over2}\sqrt{1-4\beta}\in{\bf Z}$. In
that case in order to get a solution that vanishes as
$r\rightarrow\infty$, one sets $C_1^\prime=0$ in
(\ref{eq:besselsolnZ}).

One then requires:
\begin{eqnarray}
\label{eq:finenergZ}
& & E_{U(1)}\sim \int_0^\infty 
dr r^2
{e^{-\alpha r}\over16r}
\biggl(-\alpha r 
K_{-1+{1\over2}\sqrt{1-4\beta}}({\alpha r\over2}) 
- 2 (-1 + \alpha r) K_{{1\over2}\sqrt{1-4\beta}}({\alpha r\over2})  
\nonumber\\
& & - 
   \alpha r K_{1 + {1\over2}\sqrt{1-4\beta}}({\alpha r\over2})\biggr)^2<\infty.
\nonumber\\
& & 
\end{eqnarray}
By using the series expansion of $K_\nu(x)$ about $x=0$ (See
\cite{Bateman}), one sees that one can not find 
a $\beta: {1\over2}\sqrt{1-4\beta}\in{\bf Z}$ that will satisfy 
(\ref{eq:finenergZ}) because of the singular nature of the integrand
at $r=0$.

Now, lets turn to calculating the commutative Higgs 
and gauge field: $\Phi^I,A_{r,z,{\bar z}}$. As
we have thus far continued working to all orders in the
noncommutativity parameter $\Theta$, one hence requires an all-order (in
$\Theta$) Seiberg-Witten map\cite{SW,Liu}. There is a conjecture for that as
given in \cite{Liu}. We will not explicitly construct the SW map, but
conjecture that in the neighborhood of the negative $x^3$ axis (or the
neighborhood of the South Pole of the Riemann Sphere, approaching it
the way described in the paper), the noncommutative hatted fields are the
same as the commutative unhatted fields.  Noting that as 
$Im z\rightarrow\infty$, $\hat{A_{z,{\bar z}}}^{SU(4)},
\partial,{\bar\partial}\hat{\Phi}^{SU(4)},
\partial\hat{A_z}^{U(1)}$, $\partial^2\hat{A_r}^{U(1)}\rightarrow0$,
and dropping terms of $O(\epsilon^2)$ in $\hat{A_{r,z}}^{U(1)}$ in
equations (3.3), (3.4)  and (3.7) for the Higgs and equations  (A3), (A6) and
(A8) for the gauge field of  \cite{hatath2}, the conjecture can be  
explicitly verified to $O(\Theta^2)$.  Hence, in the neighborhood of 
the South Pole, approaching it along an azimuthal angle close to ${\pi\over2}$
\begin{eqnarray}
\label{eq:noncomcomequiv}
& & \hat{\Phi}=\hat{\Phi}^a T^a=\Phi^a T^a;\nonumber\\
& & \hat{A_i}=\hat{A_i^a} T^a + \hat{A_i^0}{\bf 1}_4=A_i^a T^a +
A_i^0{\bf 1}_4.
\end{eqnarray}
The electric charge is then the same as that given in \cite{SI1}. The
magnetic charge will receive a $U(1)$ contribution proportional to 
$(\partial+{\bar\partial}){A_z}(z,{\bar z},r\rightarrow\infty)$, which using
(\ref{eq:final1}), vanishes. Hence, the noncommutative dyon in the
neighborhood of the South Pole, carries the same electric and magnetic
charges as the commutative dyon of \cite{SI1}. Their energy will be
more than that of the commutative dyons in the neighborhood of the 
South Pole, though by an infinitesimal amount. As the equality of the
commutative and noncommutative electric and magnetic charges is valid
only in the neighborhood of the South Pole, these noncommutative dyonic
solutions are likely to be stable.

\section{Summary and Discussion}

To summarize, we have shown that noncommutative (with $r-z$ and
$r-{\bar z}$ noncommutativity) $U(4)$ SYM possesses static dyonic 
solutions in the neighborhood of the South Pole of the 
Riemann sphere or equivalently the negative $x^3$ axis, 
that either are the same as the static commutative $SU(4)$ dyonic
solutions of \cite{SI1} (corresponding to {\it decoupling of a null $U(1)$
from the $SU(4)$}) or, have the same electric and magnetic charges as 
the commutative dyons and differ infinitesimally in their energy (corresponding 
to {\it decoupling of a non-null $U(1)$ from the $SU(4)$}), 
to {\it all orders in the noncommutativity parameter}. 
We have thus shown the decoupling of the $U(1)$ sector, both for null and 
infinitesimal $U(1)$ components in the neighborhood of the negative $x^3$ axis. 
This decoupling to $O(\Theta)$ for the whole of space, was 
also seen in the context of noncommutative monopole solutions in 
SYM in \cite{hata1}.  The noncommutative hatted fields are
the same as the commutative unhatted fields (we conjecture to all orders in the
noncommutativity parameter and show explicitly up to $O(\Theta^2)$).  
Thus, we see that the neighborhood of the 
negative $x^3$ axis, acts as ``denoncommutatifier". 
It will be interesting to study the same problem using noncommutative ADHMN
construction (See \cite{NCADHMN} and references therein), and 
to be able to establish the existence of such
solutions using techniques of \cite{Taubes} using Morse theory. 
Also, finding similar solutions with $z-{\bar z}$ noncommutativity will be nice.

\section*{Acknowledgements} 

We thank A.Kumar for useful discussions. We
also thank K.L.Panigrahi, R.R.Nayak, S.S.Pal and G.Bhattacharya 
for bringing to our notice references \cite{SI1}, \cite{Liu} and \cite{ncggps}.

\appendix
\section*{Appendix A}
\setcounter{equation}{0}
\seceqaa

In this appendix, we summarize the asymptotic (w.r.t to $z$)
expressions for  $A^{U(1)}_{r,z}$
and some  expressions involving them necessary to obtain the
$U(1)$ components of the gauge field and Higgs. We have used
mathematica for arriving at the expressions given in this appendix.

Using (\ref{eq:U1anztz}), one can explicity write out the expressions for
$\hat{{\cal\phi}^0},\hat{A_{z}^0}=
Det[(P_0-a_{\Phi,z}
{\bf 1}_4)(P_1-M_{\Phi,z}).(P_2-M_{\Phi,z})] +{\rm cc}$ and 
$\hat{A_r^0}=z Det[(P_0-a_{r}
 {\bf 1}_4)(P_1-M_{r}).(P_2-M_{r})]+{\rm cc}$, using mathematica,
 but because of the length of the expressions, we will give only
 their asymptotic forms.

One gets as $Im z\rightarrow\infty$:
\begin{equation}
\label{eq:hatPhiAzzinf}
\hat{\cal\phi}^0,\hat{A_z^0}\sim 2(-1+a_{\Phi,z})a^3_{\Phi,z}
(b_{\Phi,z}c_{\Phi,z}e_{\Phi,z}f_{\Phi,z})^2.
\end{equation}

Simlarly, as $z\rightarrow\infty$, one gets:
\begin{equation}
\label{eq:hatArzinf}
\hat{A_r^0}\sim2(-1+a_r)a_r^3(b_rc_re_rf_r)^2Re z.
\end{equation}
Hence, we see that if $z\rightarrow\infty$ as $z=(Re z)_0+i\infty$,
where $(Re z)_0$ is finite, then $\hat{A_r}$ is non-singular, else it
is singular.

One sees that as $z\rightarrow\infty$,
\begin{equation}
\label{eq:dhatArzinf}
\partial\hat{{\cal A}^0_r}\sim(-1+a_r)a_r^3(b_re_rc_rf_r)^2.
\end{equation}

For the purpose of calculation of energy and  imposing it to be
finite,
one requires to calculate
${\bar\partial}\hat{\cal A}^0_r$, and its asymptotic expression is given by:
\begin{equation}
\label{eq:dervbarhatArzinf}
{\bar\partial}\hat{\cal A}_r^0\sim(-1+a_r)a^3_r b_r^2 c_r^2 e_r^2 f_r^2 
\end{equation}

One also needs to evaluate $\partial{\bar\partial}\hat{{\cal A}^0_r}$
as $z\rightarrow\infty$, and is given by:
\begin{equation}
\label{eq:laphatArzinf}
\partial{\bar\partial}\hat{{\cal A}^0_r}
\sim(-1+a_r)a_r^3
{b_r^2 c_r e_r^2 f_r 
(\sqrt{3} c_r - \sqrt{3} f_r){\bar z}^2+\sqrt{3}(b_re_rc_r)^2f_rz^2\over|z|^4}
\end{equation}
Also, one one sees that:
\begin{equation}
\label{eq:laphatPhiAzzinf}
\partial{\bar\partial}(\hat{\cal\phi}^0,\hat{A_z^0})
\sim-6(-1+a_{\Phi,z})a_{\Phi,z}^3
{(((c_{\Phi,z}e_{\Phi,z})^2b_{\Phi,z}f_{\Phi,z}+(b_{\Phi,z}f_{\Phi,z})^2c_{\Phi,z}
e_{\Phi,z})\over|z|^2}.
\end{equation}

For the purpose of solving for the $U(1)$ components of the gauge
field, and for getting the magnetic charge of the dyon, one needs to
evaluate $(\partial+{\bar\partial}){\cal A}^0_z$, and 
as $z\rightarrow\infty$, one gets:
\begin{eqnarray}
\label{eq:dersumzinf}
& & (\partial+{\bar\partial})\hat{{\cal A}^0_z}\sim
{b_z c_z^2 e_z f_z^2(\sqrt{3} b_z  - \sqrt{3} e_z){\bar z}^2
+(c_z f_z b_z^2(-e_z f_z\sqrt{3}e_z+ c_z e_z\sqrt{3} e_z)
-
3c_zf_zb_ze_zc_ze_z)z^2
\over|z|^4}
\nonumber\\
& & 
\end{eqnarray}

Two other quantities that one
needs to calculate are:
$(\partial-{\bar\partial})(\hat{\Phi},\hat{A_z})$
and
$ \partial\biggl((1+|z|^2)^2(\partial+{\bar\partial})
(\hat{\Phi},\hat{A_z})\biggr)$.
The former as $z\rightarrow\infty$, is given by:
\begin{equation}
\label{eq:diffzinf}
(\partial-{\bar\partial})(\hat{{\cal\phi}^0},\hat{{\cal A}_z^0})\sim
{\sqrt{3}(-1+a_i)a_i^3\biggl[(b_ie_i)^2c_if_i(c_i-f_i)
{\bar z}
+(-b_ic_i^2e_i^2f_i(f_i-b_i)+ e_ib_i^2c_if_i^2(c_i-e_i)
)z\biggr]
(z-{\bar z})\over|z|^4},
\end{equation}
The latter as $z\rightarrow\infty$ is given by:
\begin{eqnarray}
\label{eq:derz^4}
& & \partial\biggl((1+|z|^2)^2(\partial
+{\bar\partial})(\hat{{\cal\phi}^0},\hat{{\cal A}_z^0})\biggr)\sim
\nonumber\\
& & -2(-1+a_i)a_i^3\biggl(
\sqrt{3}[c_i^2e_i^2b_if_i(-f_i+b_i)+b_i^2f_i^2c_ie_i(c_i-e_i)]z
-(3b_if_i(c_ie_i)^2+3c_ie_i(b_if_i)^2)\nonumber\\
& & +{3b_ic_ie_if_i(b_ie_i+c_if_i){\bar z}^4\over|z|^4}\biggr)
\end{eqnarray}


\begin{thebibliography}{99}
\bibitem{SI1}T.~Ioannidou and P.~M.~Sutcliffe,
{\it Non-BPS string junctions and dyons in N = 4 super-Yang-Mills},
Phys.\ Lett.\ B {\bf 467}, 54 (1999)
[hep-th/9907157]; {\it Non-Bogomolny SU(N) BPS
  monopoles}, Phys.\ Rev.\ D {\bf 60}, 105009 (1999) [hep-th/9905169].
\bibitem{SI2}T.~Ioannidou and P.~M.~Sutcliffe,
{\it Monopoles and harmonic maps,}
J.\ Math.\ Phys.\  {\bf 40}, 5440 (1999)
[hep-th/9903183].
\bibitem{hata1}K.~Hashimoto, H.~Hata and S.~Moriyama,
{\it Brane configuration from monopole solution in non-commutative super  
Yang-Mills theory},JHEP {\bf 9912}, 021 (1999)
[hep-th/9910196].
\bibitem{hata2}H.~Hata and S.~Moriyama,
{\it String junction from non-commutative super Yang-Mills theory},
JHEP {\bf 0002}, 011 (2000)
[hep-th/0001135].
\bibitem{ncggps}K.~Matsubara,
{\it Restrictions on gauge groups in noncommutative gauge theory},
Phys.\ Lett.\ B {\bf 482}, 417 (2000)
[hep-th/0003294].
\bibitem{wess}
B.~Jurco, L.~Moller, S.~Schraml, P.~Schupp and J.~Wess,
{\it Construction of non-Abelian gauge theories on noncommutative spaces,}
hep-th/0104153.
\bibitem{projskr}
T.~Ioannidou, B.~Piette and W.~J.~Zakrzewski,
{\it Low energy states in the SU(N) Skyrme models},
hep-th/9811071.
\bibitem{Bateman} {\it Bateman Manuscript Project, Higher
    Transcendenta Functions, Vol II}, A.Erde$^\prime$lyi (editor),
  McGraw-Hill Book Company (NY), 1953.
\bibitem{SW}N.~Seiberg and E.~Witten,
{\it String theory and noncommutative geometry,}
%x
JHEP {\bf 9909}, 032 (1999)
[hep-th/9908142].
\bibitem{Liu} H.~Liu, 
{\it *-Trek II: *n operations, open Wilson lines and the Seiberg-Witten  map,}
hep-th/0011125.
\bibitem{ambSW}T.~Asakawa and I.~Kishimoto,
{\it Comments on gauge equivalence in noncommutative geometry,}
JHEP {\bf 9911}, 024 (1999)
[hep-th/9909139].
\bibitem{hatath2}
S.~Goto and H.~Hata,
{\it Noncommutative monopole at the second order in theta,}
Phys.\ Rev.\ D {\bf 62}, 085022 (2000)
[hep-th/0005101].
\bibitem{noncomcomeq}E.~S.~Na, 
{\it Superpartner Solutions of a BPS Monopole in Noncommutative
  Space}, hep-th/0106278.
\bibitem{NCADHMN}
D.~Bak,
{\it Deformed Nahm equation and a noncommutative BPS monopole,}
Phys.\ Lett.\ B {\bf 471}, 149 (1999)
[hep-th/9910135].
\bibitem{Taubes}
C.~H.~Taubes,
{\it Monopoles And Maps From S(2) To S(2): The Topology Of The Configuration Space,}
Commun.\ Math.\ Phys.\  {\bf 95} (1984) 345; 
{\it 
The Existence Of A Nonminimal Solution To The SU(2) Yang-Mills Higgs
Equations On $R^3$}, Commun.\ Math.\ Phys.\  {\bf 86}, 257 (1982).
\end{thebibliography}
\newpage
\begin{figure}[htbp]
\vskip -1in
\centerline{\mbox{\psfig{file=Riemann.ps,width=0.8\textwidth}}}
\caption{The stereographic coordinate system}
  \end{figure}
\end{document}
% Psfig/TeX 
\def\PsfigVersion{1.9}
% dvips version
%
% All psfig/tex software, documentation, and related files
% in this distribution of psfig/tex are 
% Copyright 1987, 1988, 1991 Trevor J. Darrell
%
% Permission is granted for use and non-profit distribution of psfig/tex 
% providing that this notice is clearly maintained. The right to
% distribute any portion of psfig/tex for profit or as part of any commercial
% product is specifically reserved for the author(s) of that portion.
%
% *** Feel free to make local modifications of psfig as you wish,
% *** but DO NOT post any changed or modified versions of ``psfig''
% *** directly to the net. Send them to me and I'll try to incorporate
% *** them into future versions. If you want to take the psfig code 
% *** and make a new program (subject to the copyright above), distribute it, 
% *** (and maintain it) that's fine, just don't call it psfig.
%
% Bugs and improvements to trevor@media.mit.edu.
%
% Thanks to Greg Hager (GDH) and Ned Batchelder for their contributions
% to the original version of this project.
%
% Modified by J. Daniel Smith on 9 October 1990 to accept the
% %%BoundingBox: comment with or without a space after the colon.  Stole
% file reading code from Tom Rokicki's EPSF.TEX file (see below).
%
% More modifications by J. Daniel Smith on 29 March 1991 to allow the
% the included PostScript figure to be rotated.  The amount of
% rotation is specified by the "angle=" parameter of the \psfig command.
%
% Modified by Robert Russell on June 25, 1991 to allow users to specify
% .ps filenames which don't yet exist, provided they explicitly provide
% boundingbox information via the \psfig command. Note: This will only work
% if the "file=" parameter follows all four "bb???=" parameters in the
% command. This is due to the order in which psfig interprets these params.
%
%  3 Jul 1991	JDS	check if file already read in once
%  4 Sep 1991	JDS	fixed incorrect computation of rotated
%			bounding box
% 25 Sep 1991	GVR	expanded synopsis of \psfig
% 14 Oct 1991	JDS	\fbox code from LaTeX so \psdraft works with TeX
%			changed \typeout to \ps@typeout
% 17 Oct 1991	JDS	added \psscalefirst and \psrotatefirst
%

% From: gvr@cs.brown.edu (George V. Reilly)
%
% \psdraft	draws an outline box, but doesn't include the figure
%		in the DVI file.  Useful for previewing.
%
% \psfull	includes the figure in the DVI file (default).
%
% \psscalefirst width= or height= specifies the size of the figure
% 		before rotation.
% \psrotatefirst (default) width= or height= specifies the size of the
% 		 figure after rotation.  Asymetric figures will
% 		 appear to shrink.
%
% \psfigurepath#1	sets the path to search for the figure
%
% \psfig
% usage: \psfig{file=, figure=, height=, width=,
%			bbllx=, bblly=, bburx=, bbury=,
%			rheight=, rwidth=, clip=, angle=, silent=}
%
%	"file" is the filename.  If no path name is specified and the
%		file is not found in the current directory,
%		it will be looked for in directory \psfigurepath.
%	"figure" is a synonym for "file".
%	By default, the width and height of the figure are taken from
%		the BoundingBox of the figure.
%	If "width" is specified, the figure is scaled so that it has
%		the specified width.  Its height changes proportionately.
%	If "height" is specified, the figure is scaled so that it has
%		the specified height.  Its width changes proportionately.
%	If both "width" and "height" are specified, the figure is scaled
%		anamorphically.
%	"bbllx", "bblly", "bburx", and "bbury" control the PostScript
%		BoundingBox.  If these four values are specified
%               *before* the "file" option, the PSFIG will not try to
%               open the PostScript file.
%	"rheight" and "rwidth" are the reserved height and width
%		of the figure, i.e., how big TeX actually thinks
%		the figure is.  They default to "width" and "height".
%	The "clip" option ensures that no portion of the figure will
%		appear outside its BoundingBox.  "clip=" is a switch and
%		takes no value, but the `=' must be present.
%	The "angle" option specifies the angle of rotation (degrees, ccw).
%	The "silent" option makes \psfig work silently.
%

% check to see if macros already loaded in (maybe some other file says
% "\input psfig") ...
\ifx\undefined\psfig\else\endinput\fi

%
% from a suggestion by eijkhout@csrd.uiuc.edu to allow
% loading as a style file. Changed to avoid problems
% with amstex per suggestion by jbence@math.ucla.edu

\let\LaTeXAtSign=\@
\let\@=\relax
\edef\psfigRestoreAt{\catcode`\@=\number\catcode`@\relax}
%\edef\psfigRestoreAt{\catcode`@=\number\catcode`@\relax}
\catcode`\@=11\relax
\newwrite\@unused
\def\ps@typeout#1{{\let\protect\string\immediate\write\@unused{#1}}}
\ps@typeout{psfig/tex \PsfigVersion}

%% Here's how you define your figure path.  Should be set up with null
%% default and a user useable definition.

\def\figurepath{./}
\def\psfigurepath#1{\edef\figurepath{#1}}

%
% @psdo control structure -- similar to Latex @for.
% I redefined these with different names so that psfig can
% be used with TeX as well as LaTeX, and so that it will not 
% be vunerable to future changes in LaTeX's internal
% control structure,
%
\def\@nnil{\@nil}
\def\@empty{}
\def\@psdonoop#1\@@#2#3{}
\def\@psdo#1:=#2\do#3{\edef\@psdotmp{#2}\ifx\@psdotmp\@empty \else
    \expandafter\@psdoloop#2,\@nil,\@nil\@@#1{#3}\fi}
\def\@psdoloop#1,#2,#3\@@#4#5{\def#4{#1}\ifx #4\@nnil \else
       #5\def#4{#2}\ifx #4\@nnil \else#5\@ipsdoloop #3\@@#4{#5}\fi\fi}
\def\@ipsdoloop#1,#2\@@#3#4{\def#3{#1}\ifx #3\@nnil 
       \let\@nextwhile=\@psdonoop \else
      #4\relax\let\@nextwhile=\@ipsdoloop\fi\@nextwhile#2\@@#3{#4}}
\def\@tpsdo#1:=#2\do#3{\xdef\@psdotmp{#2}\ifx\@psdotmp\@empty \else
    \@tpsdoloop#2\@nil\@nil\@@#1{#3}\fi}
\def\@tpsdoloop#1#2\@@#3#4{\def#3{#1}\ifx #3\@nnil 
       \let\@nextwhile=\@psdonoop \else
      #4\relax\let\@nextwhile=\@tpsdoloop\fi\@nextwhile#2\@@#3{#4}}
% 
% \fbox is defined in latex.tex; so if \fbox is undefined, assume that
% we are not in LaTeX.
% Perhaps this could be done better???
\ifx\undefined\fbox
% \fbox code from modified slightly from LaTeX
\newdimen\fboxrule
\newdimen\fboxsep
\newdimen\ps@tempdima
\newbox\ps@tempboxa
\fboxsep = 3pt
\fboxrule = .4pt
\long\def\fbox#1{\leavevmode\setbox\ps@tempboxa\hbox{#1}\ps@tempdima\fboxrule
    \advance\ps@tempdima \fboxsep \advance\ps@tempdima \dp\ps@tempboxa
   \hbox{\lower \ps@tempdima\hbox
  {\vbox{\hrule height \fboxrule
          \hbox{\vrule width \fboxrule \hskip\fboxsep
          \vbox{\vskip\fboxsep \box\ps@tempboxa\vskip\fboxsep}\hskip 
                 \fboxsep\vrule width \fboxrule}
                 \hrule height \fboxrule}}}}
\fi
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% file reading stuff from epsf.tex
%   EPSF.TEX macro file:
%   Written by Tomas Rokicki of Radical Eye Software, 29 Mar 1989.
%   Revised by Don Knuth, 3 Jan 1990.
%   Revised by Tomas Rokicki to accept bounding boxes with no
%      space after the colon, 18 Jul 1990.
%   Portions modified/removed for use in PSFIG package by
%      J. Daniel Smith, 9 October 1990.
%
\newread\ps@stream
\newif\ifnot@eof       % continue looking for the bounding box?
\newif\if@noisy        % report what you're making?
\newif\if@atend        % %%BoundingBox: has (at end) specification
\newif\if@psfile       % does this look like a PostScript file?
%
% PostScript files should start with `%!'
%
{\catcode`\%=12\global\gdef\epsf@start{%!}}
\def\epsf@PS{PS}
%
\def\epsf@getbb#1{%
%
%   The first thing we need to do is to open the
%   PostScript file, if possible.
%
\openin\ps@stream=#1
\ifeof\ps@stream\ps@typeout{Error, File #1 not found}\else
%
%   Okay, we got it. Now we'll scan lines until we find one that doesn't
%   start with %. We're looking for the bounding box comment.
%
   {\not@eoftrue \chardef\other=12
    \def\do##1{\catcode`##1=\other}\dospecials \catcode`\ =10
    \loop
       \if@psfile
	  \read\ps@stream to \epsf@fileline
       \else{
	  \obeyspaces
          \read\ps@stream to \epsf@tmp\global\let\epsf@fileline\epsf@tmp}
       \fi
       \ifeof\ps@stream\not@eoffalse\else
%
%   Check the first line for `%!'.  Issue a warning message if its not
%   there, since the file might not be a PostScript file.
%
       \if@psfile\else
       \expandafter\epsf@test\epsf@fileline:. \\%
       \fi
%
%   We check to see if the first character is a % sign;
%   if so, we look further and stop only if the line begins with
%   `%%BoundingBox:' and the `(atend)' specification was not found.
%   That is, the only way to stop is when the end of file is reached,
%   or a `%%BoundingBox: llx lly urx ury' line is found.
%
          \expandafter\epsf@aux\epsf@fileline:. \\%
       \fi
   \ifnot@eof\repeat
   }\closein\ps@stream\fi}%
%
% This tests if the file we are reading looks like a PostScript file.
%
\long\def\epsf@test#1#2#3:#4\\{\def\epsf@testit{#1#2}
			\ifx\epsf@testit\epsf@start\else
\ps@typeout{Warning! File does not start with `\epsf@start'.  It may not be a PostScript file.}
			\fi
			\@psfiletrue} % don't test after 1st line
%
%   We still need to define the tricky \epsf@aux macro. This requires
%   a couple of magic constants for comparison purposes.
%
{\catcode`\%=12\global\let\epsf@percent=%\global\def\epsf@bblit{%BoundingBox}}
%
%
%   So we're ready to check for `%BoundingBox:' and to grab the
%   values if they are found.  We continue searching if `(at end)'
%   was found after the `%BoundingBox:'.
%
\long\def\epsf@aux#1#2:#3\\{\ifx#1\epsf@percent
   \def\epsf@testit{#2}\ifx\epsf@testit\epsf@bblit
	\@atendfalse
        \epsf@atend #3 . \\%
	\if@atend	
	   \if@verbose{
		\ps@typeout{psfig: found `(atend)'; continuing search}
	   }\fi
        \else
        \epsf@grab #3 . . . \\%
        \not@eoffalse
        \global\no@bbfalse
        \fi
   \fi\fi}%
%
%   Here we grab the values and stuff them in the appropriate definitions.
%
\def\epsf@grab #1 #2 #3 #4 #5\\{%
   \global\def\epsf@llx{#1}\ifx\epsf@llx\empty
      \epsf@grab #2 #3 #4 #5 .\\\else
   \global\def\epsf@lly{#2}%
   \global\def\epsf@urx{#3}\global\def\epsf@ury{#4}\fi}%
%
% Determine if the stuff following the %%BoundingBox is `(atend)'
% J. Daniel Smith.  Copied from \epsf@grab above.
%
\def\epsf@atendlit{(atend)} 
\def\epsf@atend #1 #2 #3\\{%
   \def\epsf@tmp{#1}\ifx\epsf@tmp\empty
      \epsf@atend #2 #3 .\\\else
   \ifx\epsf@tmp\epsf@atendlit\@atendtrue\fi\fi}


% End of file reading stuff from epsf.tex
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% trigonometry stuff from "trig.tex"
\chardef\psletter = 11 % won't conflict with \begin{letter} now...
\chardef\other = 12

\newif \ifdebug %%% turn me on to see TeX hard at work ...
\newif\ifc@mpute %%% don't need to compute some values
\c@mputetrue % but assume that we do

\let\then = \relax
\def\r@dian{pt }
\let\r@dians = \r@dian
\let\dimensionless@nit = \r@dian
\let\dimensionless@nits = \dimensionless@nit
\def\internal@nit{sp }
\let\internal@nits = \internal@nit
\newif\ifstillc@nverging
\def \Mess@ge #1{\ifdebug \then \message {#1} \fi}

{ %%% Things that need abnormal catcodes %%%
	\catcode `\@ = \psletter
	\gdef \nodimen {\expandafter \n@dimen \the \dimen}
	\gdef \term #1 #2 #3%
	       {\edef \t@ {\the #1}%%% freeze parameter 1 (count, by value)
		\edef \t@@ {\expandafter \n@dimen \the #2\r@dian}%
				   %%% freeze parameter 2 (dimen, by value)
		\t@rm {\t@} {\t@@} {#3}%
	       }
	\gdef \t@rm #1 #2 #3%
	       {{%
		\count 0 = 0
		\dimen 0 = 1 \dimensionless@nit
		\dimen 2 = #2\relax
		\Mess@ge {Calculating term #1 of \nodimen 2}%
		\loop
		\ifnum	\count 0 < #1
		\then	\advance \count 0 by 1
			\Mess@ge {Iteration \the \count 0 \space}%
			\Multiply \dimen 0 by {\dimen 2}%
			\Mess@ge {After multiplication, term = \nodimen 0}%
			\Divide \dimen 0 by {\count 0}%
			\Mess@ge {After division, term = \nodimen 0}%
		\repeat
		\Mess@ge {Final value for term #1 of 
				\nodimen 2 \space is \nodimen 0}%
		\xdef \Term {#3 = \nodimen 0 \r@dians}%
		\aftergroup \Term
	       }}
	\catcode `\p = \other
	\catcode `\t = \other
	\gdef \n@dimen #1pt{#1} %%% throw away the ``pt''
}

\def \Divide #1by #2{\divide #1 by #2} %%% just a synonym

\def \Multiply #1by #2%%% allows division of a dimen by a dimen
       {{%%% should really freeze parameter 2 (dimen, passed by value)
	\count 0 = #1\relax
	\count 2 = #2\relax
	\count 4 = 65536
	\Mess@ge {Before scaling, count 0 = \the \count 0 \space and
			count 2 = \the \count 2}%
	\ifnum	\count 0 > 32767 %%% do our best to avoid overflow
	\then	\divide \count 0 by 4
		\divide \count 4 by 4
	\else	\ifnum	\count 0 < -32767
		\then	\divide \count 0 by 4
			\divide \count 4 by 4
		\else
		\fi
	\fi
	\ifnum	\count 2 > 32767 %%% while retaining reasonable accuracy
	\then	\divide \count 2 by 4
		\divide \count 4 by 4
	\else	\ifnum	\count 2 < -32767
		\then	\divide \count 2 by 4
			\divide \count 4 by 4
		\else
		\fi
	\fi
	\multiply \count 0 by \count 2
	\divide \count 0 by \count 4
	\xdef \product {#1 = \the \count 0 \internal@nits}%
	\aftergroup \product
       }}

\def\r@duce{\ifdim\dimen0 > 90\r@dian \then   % sin(x+90) = sin(180-x)
		\multiply\dimen0 by -1
		\advance\dimen0 by 180\r@dian
		\r@duce
	    \else \ifdim\dimen0 < -90\r@dian \then  % sin(-x) = sin(360+x)
		\advance\dimen0 by 360\r@dian
		\r@duce
		\fi
	    \fi}

\def\Sine#1%
       {{%
	\dimen 0 = #1 \r@dian
	\r@duce
	\ifdim\dimen0 = -90\r@dian \then
	   \dimen4 = -1\r@dian
	   \c@mputefalse
	\fi
	\ifdim\dimen0 = 90\r@dian \then
	   \dimen4 = 1\r@dian
	   \c@mputefalse
	\fi
	\ifdim\dimen0 = 0\r@dian \then
	   \dimen4 = 0\r@dian
	   \c@mputefalse
	\fi
%
	\ifc@mpute \then
        	% convert degrees to radians
		\divide\dimen0 by 180
		\dimen0=3.141592654\dimen0
%
		\dimen 2 = 3.1415926535897963\r@dian %%% a well-known constant
		\divide\dimen 2 by 2 %%% we only deal with -pi/2 : pi/2
		\Mess@ge {Sin: calculating Sin of \nodimen 0}%
		\count 0 = 1 %%% see power-series expansion for sine
		\dimen 2 = 1 \r@dian %%% ditto
		\dimen 4 = 0 \r@dian %%% ditto
		\loop
			\ifnum	\dimen 2 = 0 %%% then we've done
			\then	\stillc@nvergingfalse 
			\else	\stillc@nvergingtrue
			\fi
			\ifstillc@nverging %%% then calculate next term
			\then	\term {\count 0} {\dimen 0} {\dimen 2}%
				\advance \count 0 by 2
				\count 2 = \count 0
				\divide \count 2 by 2
				\ifodd	\count 2 %%% signs alternate
				\then	\advance \dimen 4 by \dimen 2
				\else	\advance \dimen 4 by -\dimen 2
				\fi
		\repeat
	\fi		
			\xdef \sine {\nodimen 4}%
       }}

% Now the Cosine can be calculated easily by calling \Sine
\def\Cosine#1{\ifx\sine\UnDefined\edef\Savesine{\relax}\else
		             \edef\Savesine{\sine}\fi
	{\dimen0=#1\r@dian\advance\dimen0 by 90\r@dian
	 \Sine{\nodimen 0}
	 \xdef\cosine{\sine}
	 \xdef\sine{\Savesine}}}	      
% end of trig stuff
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\def\psdraft{
	\def\@psdraft{0}
	%\ps@typeout{draft level now is \@psdraft \space . }
}
\def\psfull{
	\def\@psdraft{100}
	%\ps@typeout{draft level now is \@psdraft \space . }
}

\psfull

\newif\if@scalefirst
\def\psscalefirst{\@scalefirsttrue}
\def\psrotatefirst{\@scalefirstfalse}
\psrotatefirst

\newif\if@draftbox
\def\psnodraftbox{
	\@draftboxfalse
}
\def\psdraftbox{
	\@draftboxtrue
}
\@draftboxtrue

\newif\if@prologfile
\newif\if@postlogfile
\def\pssilent{
	\@noisyfalse
}
\def\psnoisy{
	\@noisytrue
}
\psnoisy
%%% These are for the option list.
%%% A specification of the form a = b maps to calling \@p@@sa{b}
\newif\if@bbllx
\newif\if@bblly
\newif\if@bburx
\newif\if@bbury
\newif\if@height
\newif\if@width
\newif\if@rheight
\newif\if@rwidth
\newif\if@angle
\newif\if@clip
\newif\if@verbose
\def\@p@@sclip#1{\@cliptrue}


\newif\if@decmpr

%%% GDH 7/26/87 -- changed so that it first looks in the local directory,
%%% then in a specified global directory for the ps file.
%%% RPR 6/25/91 -- changed so that it defaults to user-supplied name if
%%% boundingbox info is specified, assuming graphic will be created by
%%% print time.
%%% TJD 10/19/91 -- added bbfile vs. file distinction, and @decmpr flag

\def\@p@@sfigure#1{\def\@p@sfile{null}\def\@p@sbbfile{null}
	        \openin1=#1.bb
		\ifeof1\closein1
	        	\openin1=\figurepath#1.bb
			\ifeof1\closein1
			        \openin1=#1
				\ifeof1\closein1%
				       \openin1=\figurepath#1
					\ifeof1
					   \ps@typeout{Error, File #1 not found}
						\if@bbllx\if@bblly
				   		\if@bburx\if@bbury
			      				\def\@p@sfile{#1}%
			      				\def\@p@sbbfile{#1}%
							\@decmprfalse
				  	   	\fi\fi\fi\fi
					\else\closein1
				    		\def\@p@sfile{\figurepath#1}%
				    		\def\@p@sbbfile{\figurepath#1}%
						\@decmprfalse
	                       		\fi%
			 	\else\closein1%
					\def\@p@sfile{#1}
					\def\@p@sbbfile{#1}
					\@decmprfalse
			 	\fi
			\else
				\def\@p@sfile{\figurepath#1}
				\def\@p@sbbfile{\figurepath#1.bb}
				\@decmprtrue
			\fi
		\else
			\def\@p@sfile{#1}
			\def\@p@sbbfile{#1.bb}
			\@decmprtrue
		\fi}

\def\@p@@sfile#1{\@p@@sfigure{#1}}

\def\@p@@sbbllx#1{
		%\ps@typeout{bbllx is #1}
		\@bbllxtrue
		\dimen100=#1
		\edef\@p@sbbllx{\number\dimen100}
}
\def\@p@@sbblly#1{
		%\ps@typeout{bblly is #1}
		\@bbllytrue
		\dimen100=#1
		\edef\@p@sbblly{\number\dimen100}
}
\def\@p@@sbburx#1{
		%\ps@typeout{bburx is #1}
		\@bburxtrue
		\dimen100=#1
		\edef\@p@sbburx{\number\dimen100}
}
\def\@p@@sbbury#1{
		%\ps@typeout{bbury is #1}
		\@bburytrue
		\dimen100=#1
		\edef\@p@sbbury{\number\dimen100}
}
\def\@p@@sheight#1{
		\@heighttrue
		\dimen100=#1
   		\edef\@p@sheight{\number\dimen100}
		%\ps@typeout{Height is \@p@sheight}
}
\def\@p@@swidth#1{
		%\ps@typeout{Width is #1}
		\@widthtrue
		\dimen100=#1
		\edef\@p@swidth{\number\dimen100}
}
\def\@p@@srheight#1{
		%\ps@typeout{Reserved height is #1}
		\@rheighttrue
		\dimen100=#1
		\edef\@p@srheight{\number\dimen100}
}
\def\@p@@srwidth#1{
		%\ps@typeout{Reserved width is #1}
		\@rwidthtrue
		\dimen100=#1
		\edef\@p@srwidth{\number\dimen100}
}
\def\@p@@sangle#1{
		%\ps@typeout{Rotation is #1}
		\@angletrue
%		\dimen100=#1
		\edef\@p@sangle{#1} %\number\dimen100}
}
\def\@p@@ssilent#1{ 
		\@verbosefalse
}
\def\@p@@sprolog#1{\@prologfiletrue\def\@prologfileval{#1}}
\def\@p@@spostlog#1{\@postlogfiletrue\def\@postlogfileval{#1}}
\def\@cs@name#1{\csname #1\endcsname}
\def\@setparms#1=#2,{\@cs@name{@p@@s#1}{#2}}
%
% initialize the defaults (size the size of the figure)
%
\def\ps@init@parms{
		\@bbllxfalse \@bbllyfalse
		\@bburxfalse \@bburyfalse
		\@heightfalse \@widthfalse
		\@rheightfalse \@rwidthfalse
		\def\@p@sbbllx{}\def\@p@sbblly{}
		\def\@p@sbburx{}\def\@p@sbbury{}
		\def\@p@sheight{}\def\@p@swidth{}
		\def\@p@srheight{}\def\@p@srwidth{}
		\def\@p@sangle{0}
		\def\@p@sfile{} \def\@p@sbbfile{}
		\def\@p@scost{10}
		\def\@sc{}
		\@prologfilefalse
		\@postlogfilefalse
		\@clipfalse
		\if@noisy
			\@verbosetrue
		\else
			\@verbosefalse
		\fi
}
%
% Go through the options setting things up.
%
\def\parse@ps@parms#1{
	 	\@psdo\@psfiga:=#1\do
		   {\expandafter\@setparms\@psfiga,}}
%
% Compute bb height and width
%
\newif\ifno@bb
\def\bb@missing{
	\if@verbose{
		\ps@typeout{psfig: searching \@p@sbbfile \space  for bounding box}
	}\fi
	\no@bbtrue
	\epsf@getbb{\@p@sbbfile}
        \ifno@bb \else \bb@cull\epsf@llx\epsf@lly\epsf@urx\epsf@ury\fi
}	
\def\bb@cull#1#2#3#4{
	\dimen100=#1 bp\edef\@p@sbbllx{\number\dimen100}
	\dimen100=#2 bp\edef\@p@sbblly{\number\dimen100}
	\dimen100=#3 bp\edef\@p@sbburx{\number\dimen100}
	\dimen100=#4 bp\edef\@p@sbbury{\number\dimen100}
	\no@bbfalse
}
% rotate point (#1,#2) about (0,0).
% The sine and cosine of the angle are already stored in \sine and
% \cosine.  The result is placed in (\p@intvaluex, \p@intvaluey).
\newdimen\p@intvaluex
\newdimen\p@intvaluey
\def\rotate@#1#2{{\dimen0=#1 sp\dimen1=#2 sp
%            	calculate x' = x \cos\theta - y \sin\theta
		  \global\p@intvaluex=\cosine\dimen0
		  \dimen3=\sine\dimen1
		  \global\advance\p@intvaluex by -\dimen3
% 		calculate y' = x \sin\theta + y \cos\theta
		  \global\p@intvaluey=\sine\dimen0
		  \dimen3=\cosine\dimen1
		  \global\advance\p@intvaluey by \dimen3
		  }}
\def\compute@bb{
		\no@bbfalse
		\if@bbllx \else \no@bbtrue \fi
		\if@bblly \else \no@bbtrue \fi
		\if@bburx \else \no@bbtrue \fi
		\if@bbury \else \no@bbtrue \fi
		\ifno@bb \bb@missing \fi
		\ifno@bb \ps@typeout{FATAL ERROR: no bb supplied or found}
			\no-bb-error
		\fi
		%
%\ps@typeout{BB: \@p@sbbllx, \@p@sbblly, \@p@sbburx, \@p@sbbury} 
%
% store height/width of original (unrotated) bounding box
		\count203=\@p@sbburx
		\count204=\@p@sbbury
		\advance\count203 by -\@p@sbbllx
		\advance\count204 by -\@p@sbblly
		\edef\ps@bbw{\number\count203}
		\edef\ps@bbh{\number\count204}
		%\ps@typeout{ psbbh = \ps@bbh, psbbw = \ps@bbw }
		\if@angle 
			\Sine{\@p@sangle}\Cosine{\@p@sangle}
	        	{\dimen100=\maxdimen\xdef\r@p@sbbllx{\number\dimen100}
					    \xdef\r@p@sbblly{\number\dimen100}
			                    \xdef\r@p@sbburx{-\number\dimen100}
					    \xdef\r@p@sbbury{-\number\dimen100}}
%
% Need to rotate all four points and take the X-Y extremes of the new
% points as the new bounding box.
                        \def\minmaxtest{
			   \ifnum\number\p@intvaluex<\r@p@sbbllx
			      \xdef\r@p@sbbllx{\number\p@intvaluex}\fi
			   \ifnum\number\p@intvaluex>\r@p@sbburx
			      \xdef\r@p@sbburx{\number\p@intvaluex}\fi
			   \ifnum\number\p@intvaluey<\r@p@sbblly
			      \xdef\r@p@sbblly{\number\p@intvaluey}\fi
			   \ifnum\number\p@intvaluey>\r@p@sbbury
			      \xdef\r@p@sbbury{\number\p@intvaluey}\fi
			   }
%			lower left
			\rotate@{\@p@sbbllx}{\@p@sbblly}
			\minmaxtest
%			upper left
			\rotate@{\@p@sbbllx}{\@p@sbbury}
			\minmaxtest
%			lower right
			\rotate@{\@p@sbburx}{\@p@sbblly}
			\minmaxtest
%			upper right
			\rotate@{\@p@sbburx}{\@p@sbbury}
			\minmaxtest
			\edef\@p@sbbllx{\r@p@sbbllx}\edef\@p@sbblly{\r@p@sbblly}
			\edef\@p@sbburx{\r@p@sbburx}\edef\@p@sbbury{\r@p@sbbury}
%\ps@typeout{rotated BB: \r@p@sbbllx, \r@p@sbblly, \r@p@sbburx, \r@p@sbbury}
		\fi
		\count203=\@p@sbburx
		\count204=\@p@sbbury
		\advance\count203 by -\@p@sbbllx
		\advance\count204 by -\@p@sbblly
		\edef\@bbw{\number\count203}
		\edef\@bbh{\number\count204}
		%\ps@typeout{ bbh = \@bbh, bbw = \@bbw }
}
%
% \in@hundreds performs #1 * (#2 / #3) correct to the hundreds,
%	then leaves the result in @result
%
\def\in@hundreds#1#2#3{\count240=#2 \count241=#3
		     \count100=\count240	% 100 is first digit #2/#3
		     \divide\count100 by \count241
		     \count101=\count100
		     \multiply\count101 by \count241
		     \advance\count240 by -\count101
		     \multiply\count240 by 10
		     \count101=\count240	%101 is second digit of #2/#3
		     \divide\count101 by \count241
		     \count102=\count101
		     \multiply\count102 by \count241
		     \advance\count240 by -\count102
		     \multiply\count240 by 10
		     \count102=\count240	% 102 is the third digit
		     \divide\count102 by \count241
		     \count200=#1\count205=0
		     \count201=\count200
			\multiply\count201 by \count100
		 	\advance\count205 by \count201
		     \count201=\count200
			\divide\count201 by 10
			\multiply\count201 by \count101
			\advance\count205 by \count201
			%
		     \count201=\count200
			\divide\count201 by 100
			\multiply\count201 by \count102
			\advance\count205 by \count201
			%
		     \edef\@result{\number\count205}
}
\def\compute@wfromh{
		% computing : width = height * (bbw / bbh)
		\in@hundreds{\@p@sheight}{\@bbw}{\@bbh}
		%\ps@typeout{ \@p@sheight * \@bbw / \@bbh, = \@result }
		\edef\@p@swidth{\@result}
		%\ps@typeout{w from h: width is \@p@swidth}
}
\def\compute@hfromw{
		% computing : height = width * (bbh / bbw)
	        \in@hundreds{\@p@swidth}{\@bbh}{\@bbw}
		%\ps@typeout{ \@p@swidth * \@bbh / \@bbw = \@result }
		\edef\@p@sheight{\@result}
		%\ps@typeout{h from w : height is \@p@sheight}
}
\def\compute@handw{
		\if@height 
			\if@width
			\else
				\compute@wfromh
			\fi
		\else 
			\if@width
				\compute@hfromw
			\else
				\edef\@p@sheight{\@bbh}
				\edef\@p@swidth{\@bbw}
			\fi
		\fi
}
\def\compute@resv{
		\if@rheight \else \edef\@p@srheight{\@p@sheight} \fi
		\if@rwidth \else \edef\@p@srwidth{\@p@swidth} \fi
		%\ps@typeout{rheight = \@p@srheight, rwidth = \@p@srwidth}
}
%		
% Compute any missing values
\def\compute@sizes{
	\compute@bb
	\if@scalefirst\if@angle
% at this point the bounding box has been adjsuted correctly for
% rotation.  PSFIG does all of its scaling using \@bbh and \@bbw.  If
% a width= or height= was specified along with \psscalefirst, then the
% width=/height= value needs to be adjusted to match the new (rotated)
% bounding box size (specifed in \@bbw and \@bbh).
%    \ps@bbw       width=
%    -------  =  ---------- 
%    \@bbw       new width=
% so `new width=' = (width= * \@bbw) / \ps@bbw; where \ps@bbw is the
% width of the original (unrotated) bounding box.
	\if@width
	   \in@hundreds{\@p@swidth}{\@bbw}{\ps@bbw}
	   \edef\@p@swidth{\@result}
	\fi
	\if@height
	   \in@hundreds{\@p@sheight}{\@bbh}{\ps@bbh}
	   \edef\@p@sheight{\@result}
	\fi
	\fi\fi
	\compute@handw
	\compute@resv}

%
% \psfig
% usage : \psfig{file=, height=, width=, bbllx=, bblly=, bburx=, bbury=,
%			rheight=, rwidth=, clip=}
%
% "clip=" is a switch and takes no value, but the `=' must be present.
\def\psfig#1{\vbox {
	% do a zero width hard space so that a single
	% \psfig in a centering enviornment will behave nicely
	%{\setbox0=\hbox{\ }\ \hskip-\wd0}
	%
	\ps@init@parms
	\parse@ps@parms{#1}
	\compute@sizes
	%
	\ifnum\@p@scost<\@psdraft{
		%
		\special{ps::[begin] 	\@p@swidth \space \@p@sheight \space
				\@p@sbbllx \space \@p@sbblly \space
				\@p@sbburx \space \@p@sbbury \space
				startTexFig \space }
		\if@angle
			\special {ps:: \@p@sangle \space rotate \space} 
		\fi
		\if@clip{
			\if@verbose{
				\ps@typeout{(clip)}
			}\fi
			\special{ps:: doclip \space }
		}\fi
		\if@prologfile
		    \special{ps: plotfile \@prologfileval \space } \fi
		\if@decmpr{
			\if@verbose{
				\ps@typeout{psfig: including \@p@sfile.Z \space }
			}\fi
			\special{ps: plotfile "`zcat \@p@sfile.Z" \space }
		}\else{
			\if@verbose{
				\ps@typeout{psfig: including \@p@sfile \space }
			}\fi
			\special{ps: plotfile \@p@sfile \space }
		}\fi
		\if@postlogfile
		    \special{ps: plotfile \@postlogfileval \space } \fi
		\special{ps::[end] endTexFig \space }
		% Create the vbox to reserve the space for the figure.
		\vbox to \@p@srheight sp{
		% 1/92 TJD Changed from "true sp" to "sp" for magnification.
			\hbox to \@p@srwidth sp{
				\hss
			}
		\vss
		}
	}\else{
		% draft figure, just reserve the space and print the
		% path name.
		\if@draftbox{		
			% Verbose draft: print file name in box
			\hbox{\frame{\vbox to \@p@srheight sp{
			\vss
			\hbox to \@p@srwidth sp{ \hss \@p@sfile \hss }
			\vss
			}}}
		}\else{
			% Non-verbose draft
			\vbox to \@p@srheight sp{
			\vss
			\hbox to \@p@srwidth sp{\hss}
			\vss
			}
		}\fi	



	}\fi
}}
\psfigRestoreAt
\let\@=\LaTeXAtSign

