• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            隨筆-48  評論-259  文章-1  trackbacks-0

            #include"iostream.h"
            #include"fstream.h"

            #define N 100
            float X[N],Y[N],M[N],a,b;
            double f(int i,int j)
            {
              double sum1=0;
              if(j-i==1) return (Y[j]-Y[i])/(X[j]-X[i]);
              else return (f(i+1,j)-f(i,j-1))/(X[j]-X[i]);
            }


            void main()
            {

              ifstream istrm("E:\\data.txt");//這里是文件的輸入,請相應改成普通鍵盤輸入
              int n=3;
              for(int i=0;i<=n;i++)
               istrm>>X[i]>>Y[i];
              istrm>>a>>b;
              float h[N],u[N],q[N],d[N];
              for(i=0;i<n;i++)
               h[i]=X[i+1]-X[i];
              q[0]=1;q[n]=0;
              for(i=1;i<n;i++)
                  q[i]=h[i]/(h[i-1]+h[i]);
              for(i=0;i<=n;i++)
               u[i]=1-q[i];
              d[0]=(6*(f(0,1)-a))/h[0];
              d[n]=(6*(b-f(n-1,n)))/h[n-1];
              for(i=1;i<n;i++)
                 d[i]=6*f(i-1,i+1);
              float m[N];
              for(i=0;i<=n;i++) m[i]=2;
              for(i=1;i<=n;i++)
              {
                m[i]=m[i]-(u[i]/m[i-1])*q[i-1];
             d[i]=d[i]-(u[i]/m[i-1])*d[i-1];
              }
              M[n]=d[n]/m[n];
              cout<<M[n]<<endl;
              for(i=n-1;i>=0;i--)
              {
                M[i]=(d[i]-q[i]*M[i+1])/m[i];
             cout<<M[i]<<endl;
              }
            }

            posted on 2007-06-09 12:16 星夢情緣 閱讀(6124) 評論(4)  編輯 收藏 引用 所屬分類: 數據結構的所有實現程序

            評論:
            # re: 數值分析--三次樣條插值公式 2007-06-09 19:54 | 隨便
            寫這個差點讓我崩潰了!我快放棄寫這了!  回復  更多評論
              
            # re: 數值分析--三次樣條插值公式 2007-06-09 22:17 | 星夢情緣
            嘿嘿,你是誰啊?我認識嗎  回復  更多評論
              
            # re: 數值分析--三次樣條插值公式 2008-03-23 22:44 | 路過
            恩,怎么沒有一點注釋呢。
            代碼好亂的。  回復  更多評論
              
            # re: 數值分析--三次樣條插值公式 2009-02-05 17:39 | sdfabc027
            不介意做個廣告的吧
            數值分析--三次樣條插值函數http://hi.baidu.com/sdfabc027/blog/item/769d09836d4cbfa50df4d2ca.html  回復  更多評論
              
            日本人妻丰满熟妇久久久久久| 26uuu久久五月天| 亚洲国产成人精品无码久久久久久综合| 要久久爱在线免费观看| 久久久无码精品亚洲日韩蜜臀浪潮| 久久精品国产日本波多野结衣| 久久夜色精品国产噜噜噜亚洲AV| 品成人欧美大片久久国产欧美...| 色天使久久综合网天天| 国产精品亚洲美女久久久| 中文精品99久久国产 | 天天久久狠狠色综合| 久久精品国产亚洲AV香蕉| 久久黄视频| 国产精品综合久久第一页| 性色欲网站人妻丰满中文久久不卡| 久久久久亚洲AV无码专区网站 | 久久亚洲中文字幕精品一区| 国产精品一久久香蕉产线看| 久久久久久极精品久久久| 久久婷婷五月综合97色| 久久久久久曰本AV免费免费| 国产亚洲精午夜久久久久久| 国产亚州精品女人久久久久久| 国产高潮国产高潮久久久| 久久精品国产清高在天天线| 久久久久亚洲AV片无码下载蜜桃| 久久一区二区免费播放| 四虎国产精品免费久久5151| 国产一级持黄大片99久久 | 93精91精品国产综合久久香蕉| 囯产极品美女高潮无套久久久| 影音先锋女人AV鲁色资源网久久 | avtt天堂网久久精品| 欧洲人妻丰满av无码久久不卡 | 99久久精品国产综合一区| 青青草国产精品久久| 99久久成人18免费网站| 91麻豆精品国产91久久久久久 | 女人香蕉久久**毛片精品| 成人精品一区二区久久|