Matlab program for pda algorithm