% Loesung von Aufgabe 7
m = 350;
c = 1.5e6;
l = 50e-3;
g = 9.81
func = @(y) y.*(1 - 1./sqrt(1 + y.^2));
y = 0:0.01:5;
f = func(y);
plot(y,f)
xlabel('y')
ylabel('f(y)')
F = getframe(gcf);
imwrite(F.cdata, 'bild70.png');
p = m*g/(2*c*l)
y0 = fzero(@(y) func(y)-p, 0)
x0 = l*y0