function [Pr]=JProbability(t,w,ps,ns,PS,NS,OP,P)

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

if(ns==NS(i,1))
    f=1;
    k=1;
elseif(ns==NS(i,2))
    f=1;
    k=2;
else
    f=0;
    Pr=0;
end

if(f==1)
    Pr=1;
    for r=1:w
        Pr=Pr*P(t,r,i,k);
    end
end