累计误差怎么怎 么 这 么 大 ?# Computation - 科学计算
m*n
1 楼
long double u=0, qq=0;
for(j=n;j>=1;j--){
qq=(-n+j-1.)/j*((1.5-j)/(1.5-j+p))*((Mu[0]-1.5+j)/(Mu[0]+j))*(1+qq);
u=(-n+j-1.)*(1.5-j)*(Mu[0]-1.5+j)/(j*(1.5-j+p)*(Mu[0]+j))*(1+u) ;
{printf("%8.15f %8.15f \n",
(double)(qq),(double)
(u));}
}
者 个 程 序 ， 如 果没 有 累 计 误 差 ， 最 后 算 出 的qq和 u 应 该 是 一 样
。， 但 是 我 用 如 下 方 法数 值 算 出 的 结 果 却 差 远 了 。 qq已 经 完
全 偏 离 了 准 确 值 。
long int n=70, p=35;
long double u=0, qq=0, Mu[0]=5.2666495553474018;
我 在 sun公 司 超 级 计 算 机 buster上 采 用 长 双 精 度 计 算。
for(j=n;j>=1;j--){
qq=(-n+j-1.)/j*((1.5-j)/(1.5-j+p))*((Mu[0]-1.5+j)/(Mu[0]+j))*(1+qq);
u=(-n+j-1.)*(1.5-j)*(Mu[0]-1.5+j)/(j*(1.5-j+p)*(Mu[0]+j))*(1+u) ;
{printf("%8.15f %8.15f \n",
(double)(qq),(double)
(u));}
}
者 个 程 序 ， 如 果没 有 累 计 误 差 ， 最 后 算 出 的qq和 u 应 该 是 一 样
。， 但 是 我 用 如 下 方 法数 值 算 出 的 结 果 却 差 远 了 。 qq已 经 完
全 偏 离 了 准 确 值 。
long int n=70, p=35;
long double u=0, qq=0, Mu[0]=5.2666495553474018;
我 在 sun公 司 超 级 计 算 机 buster上 采 用 长 双 精 度 计 算。