clear
clc
close all
%x=[-5:.01:5];
L=5; % [m]
xL=-1:0.000001:2; % [-] %xL:=x/L
x=L*xL; % [m]
%https://www.wolframalpha.com/input/?i=FindRoot%5B-x+%2B+Tan%5Bx%5D+%3D%3D+0%2C+%7Bx%2C+4.49%2C+4.50%7D%2C+WorkingPrecision+-%3E+60%5D
%xWolframAlpha=4.49340945790906417530788092728032208221558387229003290473356; %[-]
fun=@(xsym)tan(xsym)-xsym; % function
tanxMx=fzero(fun,4,optimset('TolX',1E-60));
w1=cos(pi*xL/2)-1; % [m]
w2=-sin(pi*xL); % [m]
w3=-1/(2*pi)*(sin(tanxMx*xL)+tanxMx*(1-xL-cos(tanxMx*xL))); % [m]
w4=(cos(2*pi*xL)-1)/2; % [m]
w5=(cos(pi*xL)-1)/2; % [m]
w0=0*xL; % [m]
kl3=pi/tanxMx; % ~0.699
k3=.5/kl3; %Steigung
w3Linie=k3*xL-k3; %Wirkungslinie der Kraft
plot(x,w1,x,w2,x,w3,x,w4,x,w5,x,w0,x,w3Linie)
grid on
daspect([1 1 1])
print('-dsvg','Eulerfalle')