% Prop9InvMultOpRun: Inversion of Frame multiplier according to Prop. 9 
% Run of the function Prop9InvMultOp(c,r,TPhi,TG,m,e) 
% with given input-parameters
%function y= Prop8InvMultOpRun;
ltfatstart;


% c=3; % the number c of the frame vectors
% r=2; % the number r of the coordinates of the frame vectors
% TPhi=[1 2 6; 2 3 5]; %the synthesis matrix (rxc) of the frame Phi
% TG= [1 2 -6; 2 -3 5]; %the synthesis matrix (rxc) of the frame G (Psi-Phi)
% m=[1 2 3]; %the elements of the positive symbol m (c numbers in a row): 
% e=1/100000; %the aimed error bound

c=4; % the number c of the frame vectors
r=3; % the number r of the coordinates of the frame vectors
TPhi=[4 -2 6 3; 2 3 5 -1; 3 7 2 -9]; %the synthesis matrix (rxc) of the frame Phi
TG= [3 2 -5 1; 6 -5 -2 -2; 8 -3 4 6]; %the synthesis matrix (rxc) of the frame G (Psi-Phi)
m=[2 -5 3 -4]; %the elements of the positive symbol m (c numbers in a row): 
e=1/100000; %the aimed error bound


[m,TPsi,M0,M1,M2,M0inv,n0,M1inv,M2inv,n] = Prop9InvMultOp(c,r,TPhi,TG,m,e);




M0inv % the inverse matrix of M0 via the iterative algorithm in Prop9InvMultOp
disp('Compare M1inv to the inverted matrix via Matlab:')
M0invMatlab=inv(M0)



M1inv % the inverse matrix of M1 via the iterative algorithm in Prop9InvMultOp
disp('Compare M1inv to the inverted matrix via Matlab:')
M1invMatlab=inv(M1)

M2inv % the inverse matrix of M2 via the iterative algorithm in Prop9InvMultOp
disp('Compare M2inv to the inverted matrix via Matlab:')
M2invMatlab=inv(M2)

n0
n


%end