clear all
clc
% Detta program beraknar luftmostandet hos ett objekt

Fd=20000;
rho=1e-6;
v=100;
A=1;

Cd=2*Fd/(rho*v^2*A);


%% Uppgift 2.
% Vektor med alla hastigheter
v_vec=[10:10:50 100]';

% Anvand .-notationen

Fd_vec=Cd*rho*v_vec.^2*A/2;

% Detta kommando skriver ut en tabell pa skarmen. 
% Inget ny behover kunna for denna kursen men kan vara
% bra att veta hur man gor.

T=table(v_vec,Fd_vec,'VariableNames',{'v', 'Fd'})

% Plotta luftmotstandet som funktion av v
figure(1)
plot(v_vec,Fd_vec,'-ro','LineWidth',2)
title('Luftmotstand')
xlabel('v m/s')
ylabel('Fd kgm/s2')

