N?=?4;
dag?=?zeros（?N?N?）;
C?=?1;?
S?=?2;
R?=?3;
W?=?4;
dag（?C?[R?S]?）?=?1;
dag（?R?W?）?=?1;
dag（?S?W?）?=?1;
?
node_sizes?=?2*ones（1N）;
%node_sizes?=?[?4?2?3?5];
?
%???????
onodes?=?[];
bNet?=?mk_bnet（?dag?node_sizes?）;
?
%???????
bNet.CPD{C}?=?tabular_CPD（?bNet?C?‘CPT‘?[0.5?0.5]?）;
bNet.CPD{R}?=?tabular_CPD（?bNet?R??‘CPT‘[0.8?0.2?0.2?0.8?]?）;
bNet.CPD{S}?=?tabular_CPD（?bNet?S??‘CPT‘?[0.5?0.9?0.5?0.1]?）;
bNet.CPD{W}?=?tabular_CPD（?bNet?W?‘CPT‘?[1?0.1?0.1?0.01?0?0.9?0.9?0.99?]?）;