clear all;
clc;
open('SchieferWurf04.slx');

% Anfangswerte für den schiefen Wurf
rho = 7874; % [kg/m^3] Eisen
rhostyro = 1090; % [kg/m^3] Styropor
rhoLuft = 1.3; % [kg/m^3] Luftdichte
%rhoLuft = 0; % kein Luftwiderstand
cw = 0.45; % für die Kugel bei Re<1.7*10^5
m = 5; % für die Eisenkugel
V = m/rho; %[m^3]
r = ((3*V)/(4*pi))^(1/3);% Radius der Kugel
A = r^2*pi; % Fläche der Kugel
mstyro=V*rhostyro;% Styroporkugel mit gleichem Volumen und Fläche

%m=mstyro;

alpha = 30;% Abwurfwinkel
alpharad = (alpha*pi)/180;% Umwandlung in Radiant
v0=5;%[m/s] Anfangsgeschw.
h0=1.5;%[m] Anfangshöhe
g=9.81;%[m/s^2] Erdbeschleunigung