Chinaunix

标题: 继续问很弱的问题 [打印本页]

作者: susu1732 时间: 2008-09-03 15:50
标题: 继续问很弱的问题
double WS(double radius, double z)
{

  //cout <<"R= " << R << endl;

  double temp;
  temp=radius*radius+z*z;
  temp=sqrt(temp);
  temp-=Rnuc;
  temp/=anuc;
  return 1/(1+exp(temp));
}

//overlap function T_A
double TA(double radius)
{
  double temp=0;
  for (int i=0;i<1200;i++)

temp+=WS(radius,(i+0.5)*Rnuc/400.)*Rnuc/400.;

//return result normalized to a gold nucleus
return 2*temp*197./1175.22;
}
表情中间那段输出的是一个什么样的结果？
Rnhc，和anuc是常数。。。

作者: 5毛党党员 时间: 2008-09-03 16:01
表情那段没有输出。。。还有lz最好说明一下这是什么函数

作者: nicolas.shen 时间: 2008-09-03 16:27
椭圆的面积？

作者: susu1732 时间: 2008-09-03 16:35
有点迷糊，应该是表情那段的结果会是个什么样子，呵呵
这是关于一个流体的计算

作者: 5毛党党员 时间: 2008-09-03 16:37
400. 和400.0的效果是一样的，都是为了防止编译器把double退化成int

欢迎光临 Chinaunix (http://bbs.chinaunix.net/)