subject_55168.htm

来自「vc」· HTM 代码 · 共 24 行

<p>
序号：55168 发表者：*Casper* 发表日期：2003-10-07 13:21:21
<br>主题：这道题的错误是怎么回事?
<br>内容：/*以下是原题： <BR>用牛顿迭代法求解方程x*x*x*x-3*x*x*x+1.5*x*x-4=0的在x=2.0附近的一个根。 <BR>这是我的程序代码，调试的结果没有错误，帮忙修改一下，也可以用你的做法写一个，谢谢*/ <BR><BR>#include &#34;math.h&#34;<BR>#include &#34;stdio.h&#34;<BR><BR>double sum(double m);<BR>int main(void)<BR>{<BR>&nbsp;&nbsp;double x=2,y;<BR>&nbsp;&nbsp; y=sum(x);&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<BR>&nbsp;&nbsp;printf(&#34;%f&#34;,y);<BR>&nbsp;&nbsp; return 0;<BR>&nbsp;&nbsp; }<BR>&nbsp;&nbsp;&nbsp;&nbsp;double sum(double x)<BR>&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp; {<BR>&nbsp;&nbsp;&nbsp;&nbsp;double y,temp;<BR>&nbsp;&nbsp;&nbsp;&nbsp;double a,b;<BR>&nbsp;&nbsp;&nbsp;&nbsp;<BR><BR><BR><BR><BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;a=3*pow(x,4)-6*pow(x,3)+1.5*pow(x,2)+4;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;b=4*pow(x,3)-9*pow(x,2)+3*x;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;y=a/b;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;temp=fabs(x-y);<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if(temp&lt;=0.000001)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; {<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; return (y); <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; }<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; else sum(y);<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<BR>&nbsp;&nbsp;&nbsp;&nbsp;}
<br><a href="javascript:history.go(-1)">返回上页</a><br><a href=http://www.copathway.com/cndevforum/>访问论坛</a></p>
<hr size=1>
<blockquote><p>
回复者：*Casper* 回复日期：2003-10-07 13:23:22
<br>内容：结果是-1.#IND00我晕哦&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 请大虾们帮小弟看一下<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<br>
<a href="javascript:history.go(-1)">返回上页</a><br><a href=http://www.copathway.com/cndevforum/>访问论坛</a></p></blockquote>
<hr size=1>
<blockquote><p>
回复者：saka 回复日期：2003-10-07 13:45:12
<br>内容：double sum(double x)<BR>{<BR>&nbsp;&nbsp;&nbsp;&nbsp;double temp=1，y;<BR>&nbsp;&nbsp;&nbsp;&nbsp;double a,b;<BR>&nbsp;&nbsp;&nbsp;&nbsp;<BR>&nbsp;&nbsp;&nbsp;&nbsp; while(temp&gt;=0.000001)<BR>&nbsp;&nbsp;&nbsp;&nbsp; {<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;a=3*pow(x,4)-6*pow(x,3)+1.5*pow(x,2)+4;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;b=4*pow(x,3)-9*pow(x,2)+3*x;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;y=a/b;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;temp=fabs(x-y);<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;x=y;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;}<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;return y;<BR>}<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<BR> 
<br>
<a href="javascript:history.go(-1)">返回上页</a><br><a href=http://www.copathway.com/cndevforum/>访问论坛</a></p></blockquote>
<hr size=1>
<blockquote><p>
<font color=red>答案被接受</font><br>回复者：水番 回复日期：2003-10-07 13:45:21
<br>内容：#include &#34;math.h&#34;<BR>#include &#34;stdio.h&#34;<BR><BR>double sum(double m);<BR>int main(void)<BR>{<BR>&nbsp;&nbsp;&nbsp;&nbsp;double x=2.0,y;<BR>&nbsp;&nbsp;&nbsp;&nbsp;y=sum(x);&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<BR>&nbsp;&nbsp;&nbsp;&nbsp;printf(&#34;%f\n&#34;,y);<BR>&nbsp;&nbsp;&nbsp;&nbsp;return 0;<BR>}<BR><BR>double sum(double x)<BR>{<BR>&nbsp;&nbsp;&nbsp;&nbsp;double y,temp;<BR>&nbsp;&nbsp;&nbsp;&nbsp;double a,b;<BR>&nbsp;&nbsp;&nbsp;&nbsp;a=3*pow(x,4)-6*pow(x,3)+1.5*pow(x,2)+4;<BR>&nbsp;&nbsp;&nbsp;&nbsp;b=4*pow(x,3)-9*pow(x,2)+3*x;<BR>&nbsp;&nbsp;&nbsp;&nbsp;y=a/b;<BR>&nbsp;&nbsp;&nbsp;&nbsp;<BR>&nbsp;&nbsp;&nbsp;&nbsp;temp=fabs(x-y);<BR>&nbsp;&nbsp;&nbsp;&nbsp;if(temp&lt;=0.000001 &amp;&amp; temp &gt;= -0.000001)<BR>&nbsp;&nbsp;&nbsp;&nbsp;{<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;return y;<BR>&nbsp;&nbsp;&nbsp;&nbsp;}<BR>&nbsp;&nbsp;&nbsp;&nbsp;else<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;return sum(y);<BR>}
<br>
<a href="javascript:history.go(-1)">返回上页</a><br><a href=http://www.copathway.com/cndevforum/>访问论坛</a></p></blockquote>