mallanodal1.m

clc

u=zeros(6,6);

for i=2:5
    u(1,i)=10;
    u(i,1)=75;
    u(6,i)=100;
    u(i,6)=50;
    u(1,1)=(10+75)/2;
    u(1,6)=(10+100)/2;
    u(6,6)=(50+100)/2;
    u(6,1)=(75+50)/2;
end
 
r=[-u(1,2)-(B*u(1,2));-B*u(3,1);-B*u(4,1);-u(6,2)-B*u(5,1);-u(1,3)];