b=2
h=38
L=574
E=210*(10^3)
p=2700/(10^9)
S=b*h
m=p*S*h
I=(b*h^3)/12
w=[]
F=[]
for a=0:0.25*L:L
    w_n=sqrt((3*E*I)/(m*a^3))
    w=[w w_n]
end
x=0:L
y=[]
for k=1:length(x)
    y_n=2*(x(k)^3-3*L*(x(k))^(2))/(6*E*I)
    F_n=(y_n*3*E*I)/(x(k)^(3))
    y=[y y_n]
    
    F=[F F_n]
end
plot(x,y)
F_n=(0.18*3*E*I)/(a+b)^(3)