% check examle case 1 (c1): comparison bo_d_speso and dr_d_wedz

Vs=0.4; TS2=9.4; TSu=2.7;  T0=14.2;   %parameters
c1num=1; c1den=conv([TS2  1],conv([TSu  1],[T0  0]));

% Here case 1 == code 11 or 1
[R2speso__VR_TN_TR0_TV,R2speso__TR1_TR2]=bo_d_speso(11,Vs,TS2,TSu,T0);
VRc1=R2speso__VR_TN_TR0_TV(1,1), 
TNc1=R2speso__VR_TN_TR0_TV(1,2),
TR0c1=R2speso__VR_TN_TR0_TV(1,3),
TVc1=R2speso__VR_TN_TR0_TV(1,4),
TR1c1=R2speso__TR1_TR2(1,1),
TR2c1=R2speso__TR1_TR2(1,2),

% for comparison Naslin Polynomial Method code 13 or 3
R2_dr_xz=dr_d_wedz(13,Vs,TS2,TSu,T0);        %Doppelverhältnisse, allgemein
VRdrX_c1=R2_dr_xz{1,1},  TNdrX_c1=R2_dr_xz{1,2},  TVdrX_c1=R2_dr_xz{1,3}
Z1xc1   =R2_dr_xz{2,1},     Z2xc1=R2_dr_xz{2,2},     T0xc1=R2_dr_xz{2,3}
      
case1_bo_d_sim; sim('case1_bo_d_sim',220),
