% Directorio: PFCz_2mp % Archivo: YB2mp_1 %[Y1]=YB2mp_1(Minc,z1) function [Y1]=YB2mp_1(Minc,z1) y1=inv(z1); Y1=Minc.'*y1*Minc;