// 5.1 clc; D=8*10^-3; A=D^2; disp(A,'A=') J=8*10^-3; K=16*10^-3; B=4*J*K; disp (B,'B=') disp('since A