# -*- coding: utf-8 -*-
"""
Created on 20.01.2022

@author: stovneng
"""

# import matplotlib.pyplot as plt
import numpy as np

rho = 1.2 #massetetthet luft kg/m^3
g = 9.81  #tyngdens akselerasjon m/s^2
Cd = 0.5  #drag-koeffisient kuler

#Data for: bordtennis, golf, tennis, fotball, typisk regndråpe
r = np.asarray([20,22,34,110,1])*0.001 # radius i m
m = np.asarray([2.7,46,58,430,0.0042])*0.001 # masse i kg

A = np.pi*r*r # areal på tvers av fartsretningen m^2

#Newtons 1. lov gir f = mg, som med f = 0.5*rho*A*Cd*vt*vt gir:
vt = np.sqrt(2*m*g/(rho*A*Cd)) # terminalfart m/s

np.set_printoptions(precision=1)
print('Terminalfart for følgende baller:')
print('Bordtennis, golf, tennis, fotball, typisk regndråpe')
print(vt,'m/s')