disp('  ');
disp('  fixed_fixed_fixed_fixed_plate.m   ver 1.1  July 26, 2012 ');  
disp('  by Tom Irvine  Email: tomirvine@aol.com ');
disp('  ');
disp('  This program calculates the fundamental bending frequency ');
disp('  of a fixed-fixed-fixed_fixed plate via the Rayleigh method. ');
disp(' ');
%
clear all;
%
clear E;
clear rho;
clear mu;
clear D;
clear rho;
clear h;
clear a;
clear b;
clear mass_per_area;
% 
disp(' ');  
%
fig_num=1;
%
iu=1;
% 
[E,rho,mu]=materials(iu);
%
[a,b,h,rho,D,mass]=rectangular_plate_entry(E,rho,mu);
% 
mass_per_area=rho*h;
%
[fn,zzr,fig_num,P,dPdx,d2Pdx2,W,dWdy,d2Wdy2,norm,PF]=...
fixed_fixed_fixed_fixed_plate_Rayleigh(D,rho,h,a,b,mu,mass_per_area,fig_num);
%
disp(' ');
disp(' Apply uniform pressure?  1=yes 2=no ');
iup=input(' ');
%
if(iup==1)
%    
    [disp_transfer,von_Mises_stress_transfer]=...
        fixed_fixed_fixed_fixed_plate_pressure(E,rho,mu,a,b,h,D,...
        fn,zzr,fig_num,P,dPdx,d2Pdx2,W,dWdy,d2Wdy2,norm,PF);
%    
end