%
%  rectangular_plate_entry.m  ver 1.2  October 8, 2012
%
function[a,b,h,rho,D,mass]=rectangular_plate_entry(E,rho,mu)
%
disp(' ');
%
length=input(' Enter the length (inch) ');
width=input(' Enter the width (inch) ');
thick=input(' Enter the thickness (inch) ');
area=width*length;
volume=area*thick;
%
a=length;
b=width;
h=thick;
%
smass=(rho*volume)*386.;
disp(' ');
out1=sprintf(' Structural mass = %8.4g lbm',smass);
disp(out1);
disp(' ');
disp(' Add non-structural mass ?');
disp(' 1=yes  2=no ')
imass=input(' ');
%
NSM=0;
if(imass==1)
  disp(' ');
  disp(' Enter added mass(lbm) ');
  NSM=input(' ');
end
%
mass=NSM+smass;
M=mass;
%
disp(' ');
out1=sprintf(' Total mass = %8.4g lbm ',M);
disp(out1);
disp(' ');
out1=sprintf(' Total mass density = %8.4g lbm/in^3 ',M/(a*b*h));
disp(out1);
%
mass=mass/386;
M=M/386;
rho=M/volume;
%
D=E*h^3/(12*(1.-mu^2));
%
disp(' ');
out1=sprintf(' Plate Stiffness Factor D = %8.4g (lbf in) ',D);
disp(out1);
%