BeginPackage["YMat`"]
Unprotect@@Names["YMat`*"];
ClearAll@@Names["YMat`*"];
YMos::usage="YMos[n_Integer,list_?VectorQ, i_]"
YBjt::usage="YBjt[n_Integer,list_?VectorQ, i_]"
myPseudoInverse::usage="myPseudoInverse[m_List?MatrixQ/;Precision[m]===Infinity]"
Begin["`Private`"]
YMos[n_Integer,list_?VectorQ, i_]:=
Module[{M,idx,idx1, Y,gm,ro, RG},idx=Pick[Range[Length[list]],UnitStep[list-1] UnitStep[n-list],1];
M=ConstantArray[0,{n,n}];
idx1=list[[idx]];
gm=ToExpression["gm"<>IntegerString[{i}]];