difcascad

Ica = 15*10^(-3);
Iba=0.3*10^(-3);
Iea=15.3*10^(-3);
Uca=8;
Uba=0.65;
Uea=Uba+Uca;
Ica1=Ica;
Iea1=Iea;
Iba1=Iba;
Ucea1=8.65;
Ucba1=8;
Ubea1=0.65;
Iea2=Ica1;
Ica2=(158/159)*Iea2;
Iba2=Iea2-Ica2;
Ucea1=Uca;
Ucea2=Ucea1;
Ubea1=Uba;
Ubea2=Ubea1;
Ec=30;
tau=5.31*10^(-4);
Rg=50;
Rn=1000;
rb=16.7;
B=158;
re=3.48;
R0=10;
Ica1=Ica;
Ica2=Ica1;
Iba1=Iba;
Iba2=Iba1;
Iea1=Iea;
Iea2=Iea1;
Ucea1=Uca;
Ucea2=Ucea1;
Ubea1=Uba;
Ubea2=Ubea1;
Ucba1=Uca;
Ucba2=Ucba1;
Rin=rb+(1+B)*(re+R0);
disp("Rin");
disp(Rin);
disp("Rb1");
Rb1=3*Rin;
disp(Rb1);
Ue1=0.2*Ec;
Re=(((Ue1/Iea1) - R0)/2);
disp("Re");
disp(Re);
Rc1=((Ec-Ucea1-Ue1)/Ica1);
disp("Rc1");
disp(Rc1);
R1=((Ec*Rb1)/(Ubea1+Ue1+Iba1*Rb1));
disp("R1");
disp(R1);
R2=((Rb1*R1)/(R1-Rb1));
disp("R2");
disp(R2);
Rin2=Rin;
Rb2=2*Rin2;
disp("Rb2");
disp(Rb2);
Rc2=((Ec-Ucea2-Uea)/(Ica2));
disp("Rc2");
disp(Rc2);
R3=((Ec*Rb2)/(Ubea2+Uea+ Iba1*Rb2));
disp("R3");
disp(R3);
R4=((Rb2*R3)/(R3-Rb2));
disp("R4");
disp(R4);
Rind=2*Rin;
disp("Rind");
disp(Rind);
Rinc=rb+(1+B)*(re+Re+R0);
disp("Rinc");
disp(Rinc);
Routd=2*Rc1;
disp("Routd");
disp(Routd);
Routc=Rc1;
disp("Routc");
disp(Routc);
Kud=(-(-58)/(-58*0.933 + 1))*(((Rc1*Rn/2)/(Rc1+Rn/2))/(Rg+rb+re+R0));
disp("Kud");
disp(Kud);
Kuj=(-Rc1)/R0;
disp("Kuj");
disp(Kuj);
Kus=(-Rc1)/(2*Re);
disp("Kus");
disp(Kus);
Kp=(Kud)/(Kus);
disp("Kp");
disp(Kp);
Iinavg=Iba1;
disp("Iinavg");
disp(Iinavg);