% 8 mm Temp glass (Ounapuu) clear all h = 8; h1 = h/2; %hzero= 0.288675*h; ctop=-90; 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 x_top = [-90,-74.14285714,-54.85714286,-38.57142857,-23.57142857,-17.14285714,-4.285714286,2.142857143,15,22.28571429,28.71428571,34.28571429,37.71428571,38.57142857,42.85714286,43.71428571,44.14285714]; y_top = 4:-0.25:0; y_bot = -4:0.25:-0.25; x_bot = [-89.8,-74.7,-54.236,-39.3142857,-23.87142857,-16.94285714,-4.45714286,2.452857143,15.5,22.4571429,28.91428571,34.48571429,37.91428571,38.67142857,42.95714286,43.81428571]; plot(sigma,y, sigma, y_, x_top, y_top, x_bot, y_bot);