[½»Á÷] ¡¾ÇóÖú¡¿macµÄ±àÒë»·¾³xcodeµÄʹÓÃÎÊÌâ

#include
#include
#include
#include

using namespace std;
class ffts
{
private:
        double a,b,eps,integ;
public:
        ffts(double aa,double bb,double es)
        {a=aa;b=bb;eps=es;}
        void integration();
        void output();
        double func (double);
};

void ffts::integration()
{
        int n,k;
        double fa,fb,h,t1,p,s,x,t;
        fa=func(a);fb=func(b);
        n=1;h=b-a;
        t1=h*(fa+fb)/2.0;
        p=eps+1.0;
        while(p>=eps)
        {
                s=0.0;
                for(k=0;k                 {
                        x=a+(k+0.5)*h;
                        s=s+func(x);
                }
                t=(t1+h*s)/2.0;
                p=fabs(t1-t);
                t1=t;n=n+n;h=h/2.0;
        }
        integ=t;
}

void ffts::output()
{
        char str2[20];
        cout<<"\nÊä³öÎļþÃû:"£»
        cin>>str2;
        ofstream fout(str2);
        if(!fout)
        {cout<<"\n²»ÄÜ´ò¿ªÕâ¸öÎļþ"<         fout<         fout.close();
}

double ffts::func(double x)
{ return exp(-x*x);£ý
       
void main()
{
        ffts solution(0.0,1.0,0.000001);
        solution.integration();
        solution.output();
}

Originally posted by zyj8119 at 2011-04-05 20:29:11:
[code]#include<iostream>
#include<fstream>
#include<cmath>

using namespace std;
class ffts
{
private:
        double a,b,eps,integ;
public:
        ffts(double aa ...

Originally posted by sudo at 2011-04-05 16:26:47:
Ò»ÑÛ¿´¹ýÈ¥¾Í·¢ÏÖÓÐÖÐÎÄ×Ö·û´íÎóºÍȱ´óÀ¨ºÅ´íÎó£¬ÇëÂ¥Ö÷ÏÈÈÏÕæ¼ì²éÔÙÌù¹ýÀ´...

Áí£º
1. void main()²»ºÏ±ê×¼£¬lzÊÜVC6»òÒÔÇ°°æ±¾¶¾º¦ÁË

2. ¹ØÓÚÍ·Îļþ£º
#include <cstdio>
#include <iostream ...

#include
#include
#include

using namespace std;
class ffts
{
private:
        double a,b,eps,integ;
public:
        ffts(double aa,double bb,double es)
        {a=aa;b=bb;eps=es;}
        void integration();
        void output();
        double func (double);
};

void ffts::integration()
{
        int n,k;
        double fa,fb,h,t1,p,s,x,t;
        fa=func(a);fb=func(b);
        n=1;h=b-a;
        t1=h*(fa+fb)/2.0;
        p=eps+1.0;
        while(p>=eps)
        {
                s=0.0;
                for(k=0;k                 {
                        x=a+(k+0.5)*h;
                        s=s+func(x);
                }
                t=(t1+h*s)/2.0;
                p=fabs(t1-t);
                t1=t;n=n+n;h=h/2.0;
        }
        integ=t;
}

void ffts::output()
{
        char str2[20];
        cout<<"\nÊä³öÎļþÃû:";
        cin>>str2;
        ofstream fout(str2);
        if(!fout)
        {cout<<"\n²»ÄÜ´ò¿ªÕâ¸öÎļþ"<         fout<         fout.close();
}

double ffts::func(double x)
{ return exp(-x*x);}
        
void main()
{
        ffts solution(0.0,1.0,0.000001);
        solution.integration();
        solution.output();
}