為什麼算不出答案哩 ？
井底的小笨蛙^_^
下面的程式如果只算100組數值的話，可以得到答案，可是如果
增加為1000組，出來的答案就只是 y=-1.#IND00e+000（每一組都是）
不知道是為什麼呢？我是認為可能是在 COS 運算時出了問題，因為
可能大於或是太小於 COS 可以算的範圍~~~~~
請幫我想想，或是和我說要如何解決呢？我的頭好大~~~~
#include<stdio.h>
#include<math.h>
#include<iostream.h>
long double p1=0.0,p2=0.0,p3=0.0,y=0.0,
int m=0,n=0,
double d =1.4,
double b =2.1,
float x=0.6;
int s=0;
void main()
{
FILE *fptr;
char name[12];
cout<<"please type the file name of output file\n";
cin >> name;
fptr=fopen(name,"w");
while (m<1000)
{
x+=0.00001;
for (n=-500 ; n<500 ; n++)
{
p1=(-n*(2-d));
p2=pow(b,n)*x;
s=int(p2/pow(2,63));
if (p2>=0)
p3=p2-(s*pow(2,63));
else
p3=p2+2.0;
y=y+pow(b,p1)*(1-cos(p3));
}
m+=1;
cout<<"y="<<y<<"\n";
fprintf(fptr,"y=%e\n",y);
y=0.0;
}
fclose(fptr);
}
--
▏ [1;37m ∩∩∩∩ 我和 [5;32m胖胖 [0m [1;37m 都常 幻想 我們是兩隻 [1;33m 小貓 . [5;34m☆. [5;37m★ [0m
▎ [1;37m ◎◎⊙＜ 我說： [5;31m 胖胖 [0m [1;37m 我們以後要養很多很多 [1;33m 小小貓 [1;34m☆ [0m
▍ [1;37m [5;31m = [0m [1;37m﹀╰┬╯ 她說： 好啊 要養 [5;36m很可愛 [5;35m很可愛 [0m [1;37m的 [1;33m 小小貓 . [5;36m☆ [0m
▌ [1;37m ～△ ▲～ 她靠著我 在滿天的 [5;34;47m星空 [0m [1;37m下 我倆都笑了 [1;37m . ★ [1;34m [0m [1;32m ☆ [1;37m . [0m
▊▁▃▇▇▇▃▂▁ [5;31m ☆ [0m . [5;36m ★ . [0m
█████████████████████████████▆▄▂▁
[m [1;37m※ 來源:‧蛋捲廣場 bbs.tku.edu.tw‧[FROM: 163.13.136.1] [m