clc;
clear;
P_box1=0.25; //P(box1)
P_box2=0.25; //P(box2)
P_box3=0.25; //P(box3)
P_box4=0.25; //P(box4)

Pdef_1=0.05;  //P(defective/box1)
Pdef_2=0.4;  //P(defective/box2)
Pdef_3=0.10;  //P(defective/box3)
Pdef_4=0.10;  //P(defective/box4)

//a)

Pcomp_def=(P_box1*Pdef_1)+(P_box2*Pdef_2)+(P_box3*Pdef_3)+(P_box4*Pdef_4); //Theoram of total probability

//b)
Pbox2_def=(P_box2*Pdef_2)/((P_box1*Pdef_1)+(P_box2*Pdef_2)+(P_box3*Pdef_3)+(P_box4*Pdef_4)); //Bayes theoram

disp(Pcomp_def," a) P(component is defective)=");

disp(Pbox2_def," b) P(box2|defective)=");