% 4 mm annealed glass (Achintha & Balan)
clear all
h = 4;
h1 = h/2;
%hzero= 0.288675*h;

ctop=-100;
step=0.2;
% c2=0;

for n=1:(((h/2)/step)+1)
    y(n) = (n-1)*step;
    y_(n) = -y(n);
    sigma(n) = 6*ctop/(h^2)*y(n)^2 - ctop/2;
    
end

y_top=2:-0.25:0;
x_top = [-100, -52.5, -33.75, -10, 12.5,32.5, 43.75, 50, 52.5];


y_bot = -2:0.25:0;
x_bot = [-100,-61.25,-33.75,-7.5,18.125,31.25, 43.75,50,52.5];


plot(sigma,y, sigma, y_, x_top, y_top, x_bot, y_bot);