진자 운동 코드

#include <cmath>
//M_PI 상수 포함
using namespace std;
void printpos(const double maxangle, const int length){
	//maxangle은 진자운동의 진폭, length는 줄의 길이
	//줄의 시작점은 (0, 0)으로 가정
	//오른쪽 정점부터 시작 가정
	double Angle=-maxangle;
	double VAngle=0;
	for(int i=0;i<100;i++){
		cout<<i+1<<". 현재 각도(라디안) : "<<2*M_PI*Angle<<", 구의 좌표(x, y) : ("<<
		-sin(2*M_PI*Angle)*length<<", "<<cos(2*M_PI*Angle*length)<<")\n";
		if(abs(VAngle*2*M_PI) <= 0.01){
			cout<<"진자 운동의 정점입니다.\n";
			VAngle = 0;
		}
		cout<<"현재 각속도 : "<<VAngle*2*M_PI<<'\n';
		VAngle-=1.5/length*sin(2*M_PI*Angle);
		Angle+=VAngle;
	}
}
int main() {
	// your code goes here
	double ma;
	int len;
	cin>>ma>>len;
	printpos(ma, len);
	return 0;
}