probability.m

safsdrghfydjg kukugilihg b nmvhjtdyewrtwert

function [P]=Probability(k,ps,ns,PS,NS,OP,La)

N=length(PS);
for j=1:N
    if(ps==PS(j))
        i=j;
        break;
    end
end

if(ns==NS(i,1))
    P=0.5*(1 + OP(i,1,k)*tanh(0.5*La));
elseif(ns==NS(i,2))
    P=0.5*(1 + OP(i,2,k)*tanh(0.5*La));
end