Колебания маятника в С++

Известно, что математический маятник совершает колебания по закону $$x\left(t \right)=Asin\left(\omega t+\varphi _{0} \right)$$ будем считать, что частота колебаний $\omega$ известна. Предполагается, что известно также, что в начальный момент координата маятника положительна и в k раз меньше амплитуды А, а в момент времени $t_{1}$ значение координаты маятника равно $A_{1}$. Ниже приведена программа, в которой определяется амплитуда колебаний А.

#include
#include 
using namespace std;
int main(){
//Частота колебаний:
double omega=0.2;
//Параметры задачи:
double A1,t1,k;
//Амплитуда и начальная фаза:
double A,phi0;
//Ввод параметров:
cout << "Enter t1 = ";
cin >> t1;
cout << "Enter A1 = ";
cin >> A1;
cout << "Enter k = ";
cin >> k;
phi0=asin(1/k);
A=A1/sin(omega*t1+phi0);
cout << "Amplitude A = " << A << "\n";
return 0;
}