[課業] 請問各位一個國考程式題目

[fjf1980 (Yes We Can!)]

※ [本文轉錄自 Examination 看板]

作者: fjf1980 (Yes We Can!) 看板: Examination
標題: [課業] 96地特/程式語言/泰勒展開
時間: Mon Oct 26 19:33:45

4.想問的內容：
題目：f(x)=1- x/2! + x/4! -x/6!...
      每項算法為si = (-1)^i * x^2i / 2i!
題目要求設計把引數x先模數到 0<= x <2*pi, 然後用級數計算f(x)相加
準確到項數|Sn|<10^-8


5.想法：下面是我在DEV-C++實際跑程式，可是不管我x換多少，
我的答案都是1.000000   = = 買尬！ 為什麼？
#include <stdio.h>
#include <stdlib.h>
#define pi 3.14159
double sn(int x)
{
    double result=0.0;
    double temp=0.0;
    int i=0;
    while(temp<0.00000001 && temp >-0.00000001)
    {
           temp = pow(-1,i)*pow(x*pi,2*i)/fact(2*i);
           i++;
           result += temp;
    }
    return result;
}
int fact(int x)
{
     if(x==0 || x==1)
        return 1;
     else
        return x*fact(x-1);
}
int main(int argc, char *argv[])
{
  double x=1.1;
  double input_x=0.0;
  x = x-((int)(x/2.0)*2.0);
  printf("%f\n",sn(x));

  system("PAUSE");
  return 0;
}

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.230.197.214

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.230.197.214
→ fjf1980:請教高手們  我這樣寫有問題嗎？
※ 編輯: fjf1980         來自: 61.230.197.214       (10/26 20:03)
推 POSIX:sn的argument 是int? main裡面給double?
→ fjf1980:我改double sn(double x) 結果還是一樣耶
推 VictorTom:實作演算法有沒有錯呢??i從0開始代的話, -1^0是1, x*pi
→ VictorTom:只要不是0, 代入2*i也相當於是^0, fact傳入0回傳1,
→ VictorTom:(上面的^暫當作次方描述, 不是C/C++的xor op^^||)
→ VictorTom:temp在第一輪迴圈會算出1, 超過你的while執行範圍, 就算
→ VictorTom:出1跳回去了.
→ VictorTom:既然是誤差小於0.00000001, 所以應該是它和某個值相減以
→ VictorTom:後的絕對值小於它, 而您temp的值看起來只是一項(第一項)
→ VictorTom:所以您可能要重新檢查您的實作與演算法:)
→ VictorTom:試著step by step/debug幾個測試資料就會發現這問題了:)
→ VictorTom:另外, 像2F大大說的型別問題最好也注意一下, 還有printf
→ VictorTom:要印double時要用%lf印:)
推 PsMonkey:這樣推文成章甘好？
推 VictorTom:小弟的推文常常是想到一點就推一點, 上面其實是分開想到
→ VictorTom:的三件事Orz 要是一開始就想完了, 就會回文騙P幣了XD
→ fjf1980:甘溫唷！
推 cismjmgoshr:fact(int x)也會有問題..x超過13的話會造成整數溢位
→ cismjmgoshr:不過我覺得直接用math.h裡面的cos(x)比較快 :p
推 VictorTom:(大驚)原來這是算cos用的展開示啊....Orz
→ fjf1980:越來越複雜了  數學果然煩人