clear;
echo off;

sym Rs;
sym Vb;
sym N;
sym Vd;
sym Rl;
sym x;
syms S;

Cj = 1e-13;
Omega = 2*pi*9e8;

  S(1) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.5','Rl=3e4','tan(x)-x=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','x');
  S(2) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.5','Rl=6e4','tan(x)-x=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','x');
  S(3) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.5','Rl=9e4','tan(x)-x=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','x');
  S(4) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.5','Rl=12e4','tan(x)-x=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','x');
  S(5) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.5','Rl=15e4','tan(x)-x=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','x');
  S(6) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.5','Rl=18e4','tan(x)-x=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','x');
  S(7) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.5','Rl=21e4','tan(x)-x=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','x');
  S(8) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.5','Rl=24e4','tan(x)-x=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','x');
  S(9) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.5','Rl=27e4','tan(x)-x=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','x');
  S(10) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.5','Rl=30e4','tan(x)-x=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','x');

for i=1:1:10;

%  S(i) = solve ('Rs=10','Vb=0.7','N=2','Vd=1.5','Rl=5e4','tan(x)-x=pi*Rs/(Rl*(0.5/N+Vb/Vd))','Rs','Vb','N','Vd','Rl','x');
 

  ThetaOn(i) = S(i).x;
  Rs(i) = S(i).Rs;
  Vb(i) = S(i).Vb;
  N(i) = S(i).N;
  Vd(i) = S(i).Vd;
  Rl(i) = S(i).Rl;
  
  Eta1(i) = (N(i)*Rl(i)/pi/Rs(i))*((1/2/N(i)+Vb(i)/Vd(i))^2)*(ThetaOn(i)*(2+1/(cos(ThetaOn(i))^2))-3*tan(ThetaOn(i)));
  Eta2(i) = (N(i)*Rs(i)*Rl(i)*(Cj^2)*(Omega^2)/pi)*((1/2/N(i)+Vb(i)/Vd(i))^2)*((pi-ThetaOn(i))/((cos(ThetaOn(i)))^2)+tan(ThetaOn(i)));
  Eta3(i) = (2*N(i)*Rl(i)/pi/Rs(i))*(1/2/N(i)+Vb(i)/Vd(i))*(Vb(i)/Vd(i))*(tan(ThetaOn(i))-ThetaOn(i));
  Eta(i)  = 1/(1+Eta1(i)+Eta2(i)+Eta3(i));

end

plot (Rl, Eta);
   
echo on;