大家好 打給賀!
我想寫一條通式,讓fortran可以自動對變數微分
變數[x]=x1,x2,...,xn [y]=y1,y2,...,yn
我大概的想法是這樣:
__
do i=1,n
do j=1,n
do k=1,n
if (i/=j .and. i/=k .and. j<k) then
[x(j)-x(i)][x(j)-x(k)]+[y(j)-y(i)][y(j)-y(k)]
ceta(i,j,k)=acos(
我有想過 DERIV(F,KORDER,X,BGSTEP,TOL) 這語法,但我的電腦資料庫好像沒有這個函式@@