      subroutine solve(u1,u2,t,c1,c2)
*
*  called from within stitch
*
      implicit double precision(a-h,o-z)
      dimension t(2,2)
 
      d=t(1,1)*t(2,2)-t(2,1)*t(1,2)
      u1=(t(1,2)*c2-t(2,2)*c1)/d
      u2=(t(2,1)*c1-t(1,1)*c2)/d

      return
      end 

************************************************************************

