c++做一个计算两向量夹角的程序

请测试，结果正好是90度

#include<iostream>
#include<cmath>
using namespace std;
#define PI 3.1415926
/*
    向量 A (a,b)  B(c,d) 的夹角为r
cosr= 向量A .  向量B / （向量A的摸 * 向量B的摸）
*/
int main()
{
double a[2]={1,3},b[2]={3,-1};
double ab,a1,b1,cosr;
ab=a[0]*b[0]+a[1]*b[1];
a1=sqrt(a[0]*a[0]+a[1]*a[1]);
b1=sqrt(b[0]*b[0]+b[1]*b[1]);
    cosr=ab/a1/b1;
cout<<acos(cosr)*180/PI<<endl;

}