数值分析实验 C++_算法_c++_其他

<html xmlns="http://www.w3.org/1999/xhtml"> <head> <meta charset="utf-8"> <meta name="generator" content="pdf2htmlEX"> <meta http-equiv="X-UA-Compatible" content="IE=edge,chrome=1"> <link rel="stylesheet" href="https://static.pudn.com/base/css/base.min.css"> <link rel="stylesheet" href="https://static.pudn.com/base/css/fancy.min.css"> <link rel="stylesheet" href="https://static.pudn.com/prod/directory_preview_static/62f2f84f5425817ffc27b1dd/raw.css"> <script src="https://static.pudn.com/base/js/compatibility.min.js"></script> <script src="https://static.pudn.com/base/js/pdf2htmlEX.min.js"></script> <script> try{ pdf2htmlEX.defaultViewer = new pdf2htmlEX.Viewer({}); }catch(e){} </script> <title></title> </head> <body> <div id="sidebar" style="display: none"> <div id="outline"> </div> </div> <div id="pf1" class="pf w0 h0" data-page-no="1"><div class="pc pc1 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/62f2f84f5425817ffc27b1dd/bg1.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x1 h3 y2 ff1 fs0 fc0 sc0 ls0 ws0">实验<span class="_ _0"> </span><span class="ff2">4</span>：</div><div class="t m0 x1 h3 y3 ff1 fs0 fc0 sc0 ls0 ws0">程序代码：</div><div class="t m0 x1 h4 y4 ff2 fs1 fc1 sc0 ls0 ws0">#include<span class="fc0"> <span class="fc2"><iostream></span></span></div><div class="t m0 x1 h4 y5 ff2 fs1 fc1 sc0 ls0 ws0">using<span class="fc0"> </span>namespace<span class="fc0"> std;</span></div><div class="t m0 x1 h4 y6 ff2 fs1 fc1 sc0 ls0 ws0">struct<span class="fc0"> Point{</span></div><div class="t m0 x2 h5 y7 ff2 fs1 fc1 sc0 ls0 ws0">int<span class="fc0"> n;<span class="_ _1"> </span><span class="fc3">//<span class="ff1">存储数据数</span></span></span></div><div class="t m0 x2 h5 y8 ff2 fs1 fc1 sc0 ls0 ws0">double<span class="fc0"> * p;<span class="_ _2"> </span><span class="fc3">//<span class="ff1">指向申请空间位置</span></span></span></div><div class="t m0 x2 h5 y9 ff2 fs1 fc1 sc0 ls0 ws0">struct<span class="fc0"> Point * next; <span class="fc3">//<span class="ff1">指向下一个节点</span></span></span></div><div class="t m0 x2 h5 ya ff2 fs1 fc1 sc0 ls0 ws0">struct<span class="fc0"> Point * prev; <span class="fc3">//<span class="ff1">指向上一个节点</span></span></span></div><div class="t m0 x1 h4 yb ff2 fs1 fc0 sc0 ls0 ws0">};</div><div class="t m0 x1 h4 yc ff2 fs1 fc1 sc0 ls0 ws0">int<span class="fc0"> main(){</span></div><div class="t m0 x2 h4 yd ff2 fs1 fc0 sc0 ls0 ws0">Point * Romberg();</div><div class="t m0 x2 h4 ye ff2 fs1 fc0 sc0 ls0 ws0">Point * head;</div><div class="t m0 x2 h4 yf ff2 fs1 fc0 sc0 ls0 ws0">head=Romberg();</div><div class="t m0 x2 h4 y10 ff2 fs1 fc1 sc0 ls0 ws0">while<span class="fc0">(head){</span></div><div class="t m0 x3 h4 y11 ff2 fs1 fc1 sc0 ls0 ws0">if<span class="fc0">(head->p){</span></div><div class="t m0 x4 h4 y12 ff2 fs1 fc1 sc0 ls0 ws0">for<span class="fc0">(</span>int<span class="fc0"> i=0;i<head->n;i++){</span></div><div class="t m0 x5 h4 y13 ff2 fs1 fc0 sc0 ls0 ws0">cout<<head->p[i]<<<span class="fc2">" "</span>;</div><div class="t m0 x4 h4 y14 ff2 fs1 fc0 sc0 ls0 ws0">}</div><div class="t m0 x4 h4 y15 ff2 fs1 fc0 sc0 ls0 ws0">cout<<endl;</div><div class="t m0 x3 h4 y16 ff2 fs1 fc0 sc0 ls0 ws0">}</div><div class="t m0 x3 h4 y17 ff2 fs1 fc0 sc0 ls0 ws0">head=head->next;</div><div class="t m0 x2 h4 y18 ff2 fs1 fc0 sc0 ls0 ws0">}</div><div class="t m0 x2 h4 y19 ff2 fs1 fc1 sc0 ls0 ws0">return<span class="fc0"> 0;</span></div><div class="t m0 x1 h4 y1a ff2 fs1 fc0 sc0 ls0 ws0">}</div><div class="t m0 x1 h5 y1b ff2 fs1 fc3 sc0 ls0 ws0">//Romberg <span class="ff1">函数</span></div><div class="t m0 x1 h4 y1c ff2 fs1 fc1 sc0 ls0 ws0">#include<span class="fc0"> <span class="fc2"><cmath></span></span></div><div class="t m0 x1 h4 y1d ff2 fs1 fc1 sc0 ls0 ws0">double<span class="fc0"> func(</span>double<span class="fc0"> x);</span></div><div class="t m0 x1 h4 y1e ff2 fs1 fc0 sc0 ls0 ws0">Point * Romberg(<span class="fc1">void</span>){</div><div class="t m0 x2 h5 y1f ff2 fs1 fc3 sc0 ls0 ws0">//<span class="ff1">数据输入</span></div><div class="t m0 x2 h5 y20 ff2 fs1 fc1 sc0 ls0 ws0">double<span class="fc0"> a,b;<span class="_ _3"> </span><span class="fc3">//<span class="ff1">积分上下限</span></span></span></div><div class="t m0 x2 h4 y21 ff2 fs1 fc1 sc0 ls0 ws0">double<span class="fc0"> h;<span class="_ _4"> </span><span class="fc3">//h=b-a</span></span></div><div class="t m0 x2 h4 y22 ff2 fs1 fc1 sc0 ls0 ws0">int<span class="fc0"> k;<span class="_ _5"> </span><span class="fc3">//T(k,_ )</span></span></div><div class="t m0 x2 h5 y23 ff2 fs1 fc1 sc0 ls0 ws0">double<span class="fc0"> err<span class="_ _6"></span>,eps;<span class="_ _7"> </span><span class="fc3">//<span class="ff1">精度控制，</span></span></span></div><div class="t m0 x2 h4 y24 ff2 fs1 fc1 sc0 ls0 ws0">double<span class="fc0"> x=0.0;</span></div><div class="t m0 x2 h5 y25 ff2 fs1 fc3 sc0 ls0 ws0">//<span class="ff1">数据初始化</span></div><div class="t m0 x2 h5 y26 ff2 fs1 fc0 sc0 ls0 ws0">cout<<<span class="fc2">"<span class="ff1">输入积分上下限</span>a<span class="ff1">和</span>b<span class="ff1">：</span>"</span><<endl;</div><div class="t m0 x2 h4 y27 ff2 fs1 fc0 sc0 ls0 ws0">cin>>a>>b;</div><div class="t m0 x2 h5 y28 ff2 fs1 fc0 sc0 ls0 ws0">cout<<<span class="fc2">"<span class="ff1">输入精度控制数</span>eps"</span><<endl;</div><div class="t m0 x2 h4 y29 ff2 fs1 fc0 sc0 ls0 ws0">cin>>eps;</div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div> </body> </html>