BursErr(FrameNo,Burst)

function [Err_Seq Burst_out] = BursErr(FrameNo,Burst)
%len=FrameNo*8*8*204*2;   % assigning the lenght of sequence double size of
%due to mother coderate 1/2
A0 = Burst(1);% 0.2;
A1 = Burst(2);%0.8;
TargetAv = Burst(3);% 0.1;
p01plusp10 = Burst(4);%0.1;
emis = [A0 (1-A0); A1 (1-A1)];
p0 = (TargetAv - A1)/(A0 - A1);
p1 = 1 - p0;
p10 = p0*p01plusp10;
p01 = p1*p01plusp10;
trans = [(1-p01),p01
         p10,(1-p10)];

 %FrameNo*8*8*204  =>  Length of bits in the data
[Err_Seq,states] = hmmgenerate(FrameNo*8*8*204*2+100,trans,emis,'Symbols',[1,0]);

Burst_out(1)=p0;
Burst_out(2)=p1;