using new version of matrix form of conic sections, matlab script still not consistent with convention used in ltaem fortran