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 ************************************************************************