c++问题 大虾们帮帮忙…………不知道哪儿出问题了

#include<math.h>
#include<iostream>
using namespace std;
double main()
{double l,n2,u,d,p,T1,n1,a,d1,k,fai,Ze,Bh,N;
double pi;
float Zh=2.5;
int i=7;
cin>>u>>d>>fai>>k>>Ze>>Bh;
p=625*40/60;
l=3;
pi=3.141592;
N=(1-u+sqrt((u-1)*(u-1)-4*3500/3/pi))/2; //求层数
n2=40/pi/(d*u+2*(N-1)*d); //大齿轮的转速
n1=7*n2; //小齿轮的转速
T1=9.55*pow(10,6)*p/2/n1;
a=2*k*T1/fai*(1+1/i)*pow(Ze*Zh/Bh,2);
d1=pow(a,0.33);
cout <<"u="<<u<<"\n"<<"d="<<d<<"\n"<<"fai="<<fai
<<"\n"<<"k="<<k<<"\n"<<"Ze="<<Ze<<"\n"<<"Bh="<<Bh<<"\n"
<<"n2="<<n2<<"\n"<<"n1="<<n1<<"\n"<<"p="<<p<<"\n"
<<"T1="<<T1<<"\n"<<"d1="<<d1<<"N="<<N;
}

检查下公式，公式错误了吧，我记得齿轮设计公式a=2*k*T1/fai*(1+1/i)*pow(Ze*Zh/Bh,2);中应该是a=2*k*T1/fai*(1+1/u)*pow(Ze*Zh/Bh,2);而d1=pow(a,0.33);应该是d1=pow(a,2/3);其他地方自己对照下公式，应该是有些地方化简有问题，这种公式化想用程序表达就要有耐心和细心。学习也是这样，加油！