clear all
h = 6;
h1 = h/2;
%hzero= 0.288675*h;

ctop=5;
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
plot(sigma,y, sigma, y_);