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

ctop=-95.6;
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= 5:-0.25:3;
x_top = [-95.65217391,-75,-56.52173913,-44.56521739,-34.7826087,-23.91304348,-16.30434783,-8.695652174,-6.52173913];

y_bot= -5:0.25:-3;
x_bot = [-95.65217391,-75,-60.86956522, -46.73913043,-36.95652174,-27.17391304,-18.47826087,-9.782608696,-4.347826087];



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