→ kikiqqp:歐歐歐 這用Fortran正好寫 10/02 17:37
→ james732:2的n0、3的compute p_i、4的TOL不懂指的是什麼 orz 10/02 19:21
噓 MOONRAKER:找什麼藉口，這跟寫完也沒太大差別 10/02 19:47
→ MOONRAKER:自己寫一寫，沒有什麼難的 10/02 19:47
恩.我把algorithm拉出來給大家想想
f(x)=0 .x介於[p0.p1]
step1. i=2;q0=f(p0),q1=f(p1)
step2.while i小於等於n0 do step 3-7
step3.set p=p1-q1(p1-p0)/(q1-q0) . [compute p_i]
step4. if abs[p-p_i] < TOL then OUTPUT(P); STOP
step5.set i=i+1 ; q=f(p)
step6 if q*q1<0 then set p0 =p1;q0=q1.
step7 set p1=p;q1=q.
step8 output;stop
各位大大要的是這個嗎 @@?
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 113.61.179.15
※ 編輯: coolbetter33 來自: 113.61.179.15 (10/02 17:37)