c++-用C++計算erf的導函數，算不出結果。

用C++計算erf的導函數，算不出結果。

做畢設的時候要計算一個函數,公式如下：

想用c++算出一個erf導函數的結果的近似值（k=100）。

 #include <stdio.h>
#include <math.h>
#define PI 3.1415926535898


long double c_K(int k){
    if (k == 0){
        return (float)1;
    }
    else  {
        long double num=0.0;
        for(int m=0;m<k;m++){
            num += ((c_K(m)*c_K(k-1-m))/((m+1)*(2*m+1)));}
        return num;
    }
}

long double e_F(int n){
    long double num=0.0;
    for(int k=0;k<n;k++){
        num += (c_K(k)*(pow((PI*n/2),(2*k+1))/(2*k+1)));
    } 
    return num;
}


int main(){
    long double erf_D=e_F(100);
    printf("erf_D=%f\n",erf_D);
    return 0;
} 

但是怎麼都算不出來啊，參數取10以內要算幾分鐘，再大點就算不出來了，怎麼辦。。。
我這代碼寫得有沒有問題呀，ck我測試了一下，對照現有的參數表都正確。總覺得有更好的方法做這種計算，但是導師什麼也沒說。。。

最佳回答：

樓主，程序性能主要耗在算cK(k)函數上了:計算cK(9)，要把cK(8)----cK(0)重算很多遍
是否可以做一個數組記錄計算出來的cK(k)的值，這樣每次每次先查表，沒有的話再計算