作者agody (讲不听)
看板Fortran
标题Re: [问题] 关於微分的写法
时间Thu Mar 19 00:52:04 2009
※ 引述《latinboy (昵称)》之铭言：
: ※ 引述《comicgirl (一定要没事!!!拜托!!!)》之铭言：
: : 请问各位强者
: : 怎麽将y',y''写成fortran的语法
: : 直接写y',y''似乎不被接受
: : 之前有看过写成y1.y2,不过我也run不出来Q____________Q
: : 我是fortran的初学者,学的好辛苦
: : 希望大家不吝指教,谢谢^^
: 你这样写当然不行......
: 有些语言能够直接作微分 可惜 fortran没有....
: 自己写也很简单 就是自己微阿
: 随便举个例子 假设 y = x^3 + x^2 - x + 3
: fortran写法 y = x**3 + x**2 - x + 3
: dy = 3*(x**2) + 2*x - 1
: d2y = 6*x + 2
随手写的,不知道有没有错
program main
external y
call diff_1st(y,1,f)
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)
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)
y = x**3 + x**2 - x + 3
return
end
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 114.32.182.191
1F：→ agody:没人理我QQ 03/20 21:18