用C++计算erf的导函数，算不出结果。

Question

做毕设的时候要计算一个函数,公式如下：

想用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我测试了一下，对照现有的参数表都正确。总觉得有更好的方法做这种计算，但是导师什么也没说。。。

Answer 1

楼主，程序性能主要耗在算cK(k)函数上了:计算cK(9)，要把cK(8)----cK(0)重算很多遍
是否可以做一个数组记录计算出来的cK(k)的值，这样每次每次先查表，没有的话再计算