如何用蒙特卡洛算简单条件概率？
 

题目：某厂生产的灯泡能用1000小时的概率为0.8, 能用1500小时的概率为0.4 , 求已用1000小时的灯泡能用到1500小时的概率
答案应该为0.5，可是我自己写的程序答案却在0.4附近。请帮我看看是什么地方有问题或者求一个正确的程序，谢谢
function liti3(mm)
randnum1=binornd(1,0.8,1,mm);
randnum2=binornd(1,0.4,1,mm);
randnum=[randnum1;randnum2]
a=0;b=0;
pro=zeros(1,mm);
for i=1:mm;
if randnum(1,i)==1
a=a+1
end
if randnum(1,i)==1&&randnum(2,i)==1
b=b+1
end
pro(i)=b./a;
end
pro=pro
num=1:mm
plot(num,pro)