Contents

% /home/dimarzio/Documents/working/12303/matlab/nodes.m
% Thu Sep 27 12:08:57 2018
R1=1000;R2=200;R3=400;R4=2000;R5=800;

As I solved it

disp('-----------------------------------');
Is=0.02;Vs=5;
disp(['Is = ',num2str(Is),' A']);

M=[1/R1+1/R4, -1/R1,           0;...
   -1/R1,      1/R2+1/R3+1/R1 -1/R3;...
     0,       -1/R3,          1/R3+1/R5]
y=[Is+Vs/R1;-Vs/R1;-Is]
x=inv(M)*y

V1=x(1),V2=x(2),V3=x(3)

sum1=(V1-Vs-V2)/R1+V1/R4-Is
sum2=V2/R2+(V2-V3)/R3-(V1-Vs-V2)/R1
sum3=(V3-V2)/R3+Is+V3/R5

-----------------------------------
Is = 0.02 A
M =
    0.0015   -0.0010         0
   -0.0010    0.0085   -0.0025
         0   -0.0025    0.0037
y =
    0.0250
   -0.0050
   -0.0200
x =
   16.4865
   -0.2703
   -5.5135
V1 =
   16.4865
V2 =
   -0.2703
V3 =
   -5.5135
sum1 =
   3.4694e-18
sum2 =
  -1.7347e-18
sum3 =
  -1.7347e-18

With 0.2 mA

disp('-----------------------------------');
Is=0.2e-3;Vs=5;
disp(['Is = ',num2str(Is),' A']);

M=[1/R1+1/R4, -1/R1,           0;...
   -1/R1,      1/R2+1/R3+1/R1 -1/R3;...
     0,       -1/R3,          1/R3+1/R5]
y=[Is+Vs/R1;-Vs/R1;-Is]
x=inv(M)*y

V1=x(1),V2=x(2),V3=x(3)

sum1=(V1-Vs-V2)/R1+V1/R4-Is
sum2=V2/R2+(V2-V3)/R3-(V1-Vs-V2)/R1
sum3=(V3-V2)/R3+Is+V3/R5

-----------------------------------
Is = 0.0002 A
M =
    0.0015   -0.0010         0
   -0.0010    0.0085   -0.0025
         0   -0.0025    0.0037
y =
    0.0052
   -0.0050
   -0.0002
x =
    3.2865
   -0.2703
   -0.2335
V1 =
    3.2865
V2 =
   -0.2703
V3 =
   -0.2335
sum1 =
  -1.3553e-19
sum2 =
  -2.1684e-19
sum3 =
   5.4210e-20

Published with MATLAB® R2017a