6.3  Subroutines LSEM3 and LSEMA