fortran。谁能帮我看看这个程序哪里错了！！尽快！

显示全部楼层 · 2012-10-20 11:38:56

你的算法我不了解。但是你循环，用 un1 作为精度，要求满足到 eps 的精度。但是 un1=1./(x**2-(n*pi)**2)其中，x 不变，n 不变，pi 也不变。换句话说，永远满足不了精度。你是否忘了 n=n+1 ？？？（只是猜测）...

千问 · 2012-10-20 11:38:56

低级错误太多，建议先回去好好翻翻FORTRAN教程...

千问 · 2012-10-20 11:38:56

注意循环do while(abs(un1)>=eps)
un1=1./(x**2-(n*pi)**2)并没有破坏循环的条件,un1的值实际并没有发生改变，这是一个死循环。...

千问 · 2012-10-20 11:38:56

un1=1./(x**2-(n*pi)**2)un2=flag*un1enddos=(1+2*x**2*sum1)/(1+2*x**2*sum2)end function g都错了...