meeting 2014/06/17
implement lsor4c.f90 but with bugs