CmSc 117.zip_CMSC_beam

function L_b = beamLengths(joints, beams) sizejoints = size(joints,1); sizebeams = size(beams,1); joints xjoints(1) = 0; for i = 1:sizejoints xjoints(i) = joints(i,1); end xjoints yjoints(1) = 0; for j = 1:sizejoints yjoints(j) = joints(j,2); end yjoints beams xbeams(1) = 0; for i = 1:sizebeams xbeams(i) = beams(i,1); end xbeams ybeams(1) = 0; for j = 1:sizebeams ybeams(j) = beams(j,2); end ybeams distance(1) = 0; for i = 1:sizebeams if yjoints(xbeams(i)) == yjoints(ybeams(i)) L_b(i) = xjoints(ybeams(i)) - xjoints(xbeams(i)); elseif xjoints(xbeams(i)) == xjoints(ybeams(i)) L_b(i) = yjoints(ybeams(i)) - yjoints(xbeams(i)); else a = xjoints(ybeams(i)) - xjoints(xbeams(i)); b = yjoints(xbeams(i)) - yjoints(ybeams(i)); L_b(i) = (a.*a + b.*b).^(1/2); end end L_b end