function [ J ] = Diode3( z )
 
   V1 = z(1);
   V2 = z(2);
   V3 = z(3);
 
   Idss = 1e-8;
   nVt = 0.052;
   
   Id1 = Idss* exp( (V1 - V2)/nVt - 1 );
   Id2 = Idss* exp( (V2 - V3)/nVt - 1 );
   Id3 = Idss* exp( (V3 - 0)/nVt - 1 );
   
   e1 = (V1 - 10)/100 + Id1 + (V1-V2)/200 + (V1-V3)/300 + (V1/400);
   e2 = (V2-V1)/200 - Id1 + Id2;
   e3 = (V3-V1)/300 - Id2 + Id3;
   
   J = (e1)^2 + (e2)^2 + (e3)^2;

   disp([V1, V2, V3, log10(J)]);
   pause(0.1);

   end