#!/usr/bin/octave -qf

#zaokrouhleni
output_precision = 15;
#Úloha 3
disp ("\n \n |---------------------------------|\n \n")
disp ("Script pro výpočet úloly 3 do")
disp ("VB005 - Panorama fyziky")
disp ("\n \n |---------------------------------|\n \n")
format short g

disp ("Zadejte prosím několik vstupních dat pro:\n")

m_1 =  input ("\n Vlož hmotnost první koule v kilogramech: ");
m_2 =  input ("\n Vlož hmotnost druhé koule v kilogramech: ");
l = input ("\n Vlož vzdálenost mezi koulemi v metrech: ");
delta_l = input ("\n Vlož změnu vzdálenosti v metrech: ");

#gravitacni konstanta
G = 6.672e-11 #N.m2.kg-2;

#Sila mezi koulemi
F = G * (m_1 * m_2) / (l**2);

#cas
a = F / (m_1 + m_2);
t_sec = sqrt((2*delta_l) / (a));
t_min = t_sec/60;

disp ("\n \n |---------------------------------|\n \n")
disp ("Výsledky:")
printf ("\n Na koule působí síla %g [N]\n", F)
printf ("\n Tato síla vyvolá zrychlení %g [ms^(-2)]\n", a)
printf ("\n Na pozadovanou vzdálenost se koule dostanou za %g [sec], resp. %g [min]\n", t_sec, t_min)
disp ("\n \n |---------------------------------|\n \n")