%clear all;clc;
open('SchieferWurf02.slx');

% Angaben für den schiefen Wurf
% Angabe mittels Vektoren

rho = 7874; % [kg/m^3] Eisen
rhostyro = 1090; % [kg/m^3] Styropor

%rhoLuft = 1.3; % [kg/m^3] wenn Luftwiderstand
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

% Berechnung für Styroporkugel
%m=mstyro; 

alpha = 90;% Abwurfwinkel
alpharad = (alpha*pi)/180;% Umwandlung in Radiant
v0=5;%[m/s] Anfangsgeschw.
vektor=[cos(alpharad) sin(alpharad)];% Richtungsvektor
h0=1.5;%[m] Anfangshöhe
Startpunkt = [0 h0];
g=9.81;%[m/s^2] Erdbeschleunigung
G = [0 -1];