作者agody (讲不听)
看板Fortran
标题Re: [问题] 关於微分的写法
时间Thu Jun 4 00:04:45 2009
※ 引述《dana79810 (kicala)》之铭言:
: ※ 引述《agody (讲不听)》之铭言:
: : 随手写的,不知道有没有错
: : program main
: : external y
: : call diff_1st(y,1,f)<---後面(y,l,f) 为什麽和下面附程式名子的不一样?
: 是有其他的意释吗?
你去查一下external的函数使用方法
我可以任意给y,只要格式对,diff_1st中的fun就会去对到y subroutine
: : print *,f
: : end
: : subroutine diff_1st(fun,x,f)
: : external fun
: : dx = 1
: : k=1
: : f_old = (fun(x+dx)-fun(x-dx))/(2*dx)
: : dx = dx/2
: : do while(k) <--K是? 是 do while( K>0 )作回圈的意释吗?
k=1 true 回圈
k=0 f 结束
: : f_new = (fun(x+dx)-fun(x-dx))/(2*dx)
: : dx = dx/2
: : if(f_new.eq.f_old) k = 0
: : f_old = f_new
: : end do
: : f = f_new
: : return
: : end
: : function y(x) <---哪一个步骤是将值带入?
diff_1st中所有fun的函数
: : y = x**3 + x**2 - x + 3
: : return
: : end
: 因为我不太懂阵列的运用 ,所以有很多问题..请各位板大指点!!
: orz
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 122.126.120.172
1F:推 dana79810:谢谢唷!!! 我了解了 06/10 09:45