function [ J ] = Diode1( z )
 
   Vd = z(1);
 
   Idss = 1e-8;
   nVt = 0.052;
   
   Id = Idss* exp( Vd/nVt - 1 );
   
   e1 = Id + (Vd - 10)/100;
   
   J = (e1)^2;

   disp([Vd, log10(J)]);
   pause(0.1);
   
   end