:Steffensenٵ㷨

ʽ:
[x1,er,l]=SolveImplicitSF(f,x,x0,,LoopSum,Error,,,y1,y2,,yn)
[x1,er,l]=SolveImplicitSF(f,x,x0,,LoopSum,Error,,,y1,y2,)
[x1,er,l]=SolveImplicitSF(f,x,x0,,LoopSum,Error,,,y1,y2)
[x1,er,l]=SolveImplicitSF(f,x,x0,,LoopSum,Error,,,y1)
[x1,er,l]=SolveImplicitSF(f,x,x0,,LoopSum,Error)
[x1,er,l]=SolveImplicitSF(f,x,x0,,LoopSum)
[x1,er,l]=SolveImplicitSF(f,x,x0,)
[x1,er,l]=SolveImplicitSF(f,x,x0)
[x1,er,l]=SolveImplicitSF(f,x)
[x1,er,l]=SolveImplicitSF(f)

f:ű洢x,y1,y2,,ynΪıʽ
x:ҪԱ,ĬΪx
x0:Ӧxĳʼֵ,ֵΪ,ÿԪضӦy1,y2,,ynԪءֵĬϵͳʼ
LoopSum:ѭ,Ĭ3000
Error:Ƶľ,Ĭ1e-12
y1,y2,,yn:ΪʽfԱ,Ϊֵ߾,άһ

x1:ضӦĽ
er:ؽԾ
l:⵽Žѭ

ע:

1ֻһδ֪
2ж,ֻһ
3زʱ,עŵ


ο://ִӦѧֲ.ֵ[M].廪ѧ,,2007:267-269

:
//֪бʽ 6*sin(x+z*y)+5-6*y=0 ,֪xz,yֵ

x = [ 4.00000000000000    5.00000000000000    6.00000000000000
      7.00000000000000    8.00000000000000    9.00000000000000 ]
z = [-5.00000000000000    6.00000000000000    0.00000000000000
      9.00000000000000    3.00000000000000    4.00000000000000 ]

//ϵͳѾxz,ִ
f="6*sin(x+z*y)+5-6*y";
[y,er,l]=SolveImplicitSF(f,"y",,,3000,1e-12,,,x,z)//سõ½(ע,˼)xΪ4zΪ-5ʱӦyΪ1.60443127935459ʱΪ1.7763568394E-15
y =
[ 1.60443127935459    1.35192988769122    0.55391783513440
  1.38101651045136    0.56535672402838    0.91103856432234 ]
er =
[ 1.7763568394E-15    3.5527136788E-15    0.00000000000000
  1.5809575870E-13    5.7731597280E-15    3.5527136788E-15 ]
l =
[ 12.0000000000000    8.00000000000000    2.00000000000000
  13.0000000000000    70.0000000000000    298.000000000000 ]