epoch=100; for i=1:epoch [m2,T2,P2]=PTNK(m1,T1,P1); .............. .............. .............. [m1,T1,P1]=HLR(m4,T4,P4) end