cosとsinの計算結果が違う。

Question

###実現したいこと
平面上の点を、原点を中心に回転させた後の座標を出力したい。
###試したこと
std::cout<<theta<<"ラジアンは"<<thdo<<"度\n";をしたら60.00..と出てきたのでラジアンの変換はできています。
###問題点
C++による誤差なのかもしれないのですが、回転させた後の座標が違う気がします。
手計算でvoid rotateの中の式x=x*cos(thdo)-y*sin(thdo);を計算してみたところ4になりました。
計算結果とあまりにも違います。
三角関数のプログラムの仕方を学習したばかりなのであまり使い方がわかりません。教えてください。
ちなみに、1.0472radは約60°です。
また、回転させた後の座標を見る際にこちらのサイトを参考にしました。これによると、実行結果は(2,3.4641)になるはずです。
実行結果

平面上の点を1つ入力してください: 4 0
回転させる角度を入力してください(※角度はラジアンを単位として): 1.0472
回転させた後の座標は
(-3.80948,1.16168)

###コード

//点の回転

#include<iostream>
#include<cmath>

//x座標y座標の定義
struct point
{
    double x; //x座標
    double y; //y座標
    
    //原点を中心に反時計回りにθだけ回転させるメンバ関数
    void rotate(double theta){
        double pai=std::acos(-1.0); //円周率π
        double thdo=theta*(180/pai); //ラジアンを度にする
        //std::cout<<theta<<"ラジアンは"<<thdo<<"度\n";
        x=x*cos(thdo)-y*sin(thdo);
        //std::cout<<"x座標"<<x;
        y=x*sin(thdo)+y*cos(thdo); //点(x,y)を原点を中心にθだけ回転させると(xcosθ-ysinθ,xsinθ+ycosθ)となる
    }
};

int main()
{
    point p;
    double kakudo;
    std::cout<<"平面上の点を1つ入力してください: ";
    std::cin>>p.x;
    std::cin>>p.y;
    std::cout<<"回転させる角度を入力してください(※角度はラジアンを単位として): ";
    std::cin>>kakudo;
    p.rotate(kakudo);

    std::cout<<"回転させた後の座標は\n";
    std::cout<<"("<<p.x<<","<<p.y<<")";

    return 0;
}

Answer 1

        x=x*cos(thdo)-y*sin(thdo); // [1]
        y=x*sin(thdo)+y*cos(thdo); // [2]

[1]をやると、xが書き換わってるので[2]の結果が狂います。

Comments

[langhtorn]

xの値が変わらないよう先に別の変数に置いたらできました。ありがとうございます。

[episteme]

これまでの質問を見ると、このテの"不注意"が原因のものが多いようで...

Answer 2

cos,sinの使い方を勘違いしていたようです。
cosθのθの部分の単位はラジアンでした。度ではないみたいです。それを治したら解けました！

//点の回転

#include<iostream>
#include<cmath>

//x座標y座標の定義
struct point
{
    double x; //x座標
    double y; //y座標
    
    //原点を中心に反時計回りにθだけ回転させるメンバ関数
    void rotate(double theta){
        /*double pai=std::acos(-1.0); //円周率π
        double thdo=theta*(180/pai); //ラジアンを度にする
        std::cout<<theta<<"ラジアンは"<<thdo<<"度\n";*/
        double dx=x; //x座標が変わらないように置いておく
        //std::cout<<"cos(thdo)="<<cos(theta)<<"\n";
        x=dx*cos(theta)-y*sin(theta);
        //std::cout<<"x座標"<<x;
        //std::cout<<"dx="<<dx;
        y=dx*sin(theta)+y*cos(theta); //点(x,y)を原点を中心にθだけ回転させると(xcosθ-ysinθ,xsinθ+ycosθ)となる
    }
};

int main()
{
    point p;
    double kakudo;
    std::cout<<"平面上の点を1つ入力してください: ";
    std::cin>>p.x;
    std::cin>>p.y;
    std::cout<<"回転させる角度を入力してください(※角度はラジアンを単位として): ";
    std::cin>>kakudo;
    p.rotate(kakudo);

    std::cout<<"回転させた後の座標は\n";
    std::cout<<"("<<p.x<<","<<p.y<<")";

    return 0;
}

/*実行結果
平面上の点を1つ入力してください: 4 0
回転させる角度を入力してください(※角度はラジアンを単位として): 1.0472
回転させた後の座標は
(1.99999,3.46411)
*/