Yeechs...try something like ( (1+RR) . P2 --> (1+RR) * P2 ) xprime : (1+RR) * (1-RR) * P1 + (1+RR) * P2 + 2 * (1-RR) * SQ * matrix([x[1]],[x[2]],[x[3]]); Barton