SOR-Newton method n=2

Which case should be chosen?

F(x(1),x(2)) = grad((x(1)+x(2))**2 + x(2)**2 + exp(x(1)) + exp(x(2))) = 0

Please type the evaluation program of your function here using FORTRAN rules. This function has the form d22usr(x,y,f1,f2) with input x, y and output f1, f2. Your final statement must be

      f1=
      f2=
You may use the constants pi, e(=exp(1)), sqrt2(=1.414...), the integer variables i,j,k, the logicals bool1,bool2,bool3 and the double precision variables sum,h1,h2,h3,h4,u(100),v(100),a(100,100) which are all intialized with zero resp. .false. . never change x or y!. first is a local integer and set 0 before calling the function the first time. You may use this in order to initialize some local data and set it 1 afterwards to avoid multiple such initialization. Your settings of the local variables are preserved during program execution.

you may give here a title for your results
Case

Input of the initial guess x0:
x0,1 = , x0,2 =

Please specify the relaxation parameter omega:
omega= Important: 0 < omega < 2 !

Please specify the required final precision: eps
|F1(xfinal)|, |F2(xfinal)| < = eps

eps = Important: 1.E-8 <= eps <= 1.E-3 !

Please specify the maximum number of steps (k):
maxiter = Important: 50 < = iter <= 1000 !

Please specify the divergence indicator: we terminate if
||F(xk)|| >= fac*||F(x0)|| .

fac = Important: 2<=fac !

Please specify the rotation angles for the 3D coordinate system:
xang= Important: 0 <= xang <= 180,0 <= zang <= 360.
For xang = 0 and zang=0 means x-axis is on screen horizontal and z-axis is vertically to the screen, that is you look vertically down on the (x,y) -plane. zang describes the rotation around the x-axis.
zang=

Click on "evaluate", in order to submit your input.

Back to the theory page

 Back to the top!

25.04.2016