不动点迭代法

Domx1=linspace(-10,10,200)
FX1=((Domx1+1)/2).^(1/3)
plot(Domx1,FX1)
hold on
plot(Domx1,Domx1)
hold on
k=1
x(1,k)=0
for i=1:10
x(1,k+1)=((x(1,k)+1)/2).^(double(1)/double(3))
domx1(1:100)=x(1,k)
y1=linspace(x(1,k),((x(1,k)+1)/2).^(double(1)/double(3)),100)
plot(domx1,y1)
hold on
domx2=linspace(x(1,k),x(1,k+1),100)
y2(1:100)=x(1,k+1)
plot(domx2,y2)
hold on
k=k+1
end
x‘

Domx1 =

  Columns 1 through 15

  -10.0000   -9.8995   -9.7990   -9.6985   -9.5980   -9.4975   -9.3970   -9.2965   -9.1960   -9.0955   -8.9950   -8.8945   -8.7940   -8.6935   -8.5930

  Columns 16 through 30

   -8.4925   -8.3920   -8.2915   -8.1910   -8.0905   -7.9899   -7.8894   -7.7889   -7.6884   -7.5879   -7.4874   -7.3869   -7.2864   -7.1859   -7.0854

  Columns 31 through 45

   -6.9849   -6.8844   -6.7839   -6.6834   -6.5829   -6.4824   -6.3819   -6.2814   -6.1809   -6.0804   -5.9799   -5.8794   -5.7789   -5.6784   -5.5779

  Columns 46 through 60

   -5.4774   -5.3769   -5.2764   -5.1759   -5.0754   -4.9749   -4.8744   -4.7739   -4.6734   -4.5729   -4.4724   -4.3719   -4.2714   -4.1709   -4.0704

  Columns 61 through 75

   -3.9698   -3.8693   -3.7688   -3.6683   -3.5678   -3.4673   -3.3668   -3.2663   -3.1658   -3.0653   -2.9648   -2.8643   -2.7638   -2.6633   -2.5628

  Columns 76 through 90

   -2.4623   -2.3618   -2.2613   -2.1608   -2.0603   -1.9598   -1.8593   -1.7588   -1.6583   -1.5578   -1.4573   -1.3568   -1.2563   -1.1558   -1.0553

  Columns 91 through 105

   -0.9548   -0.8543   -0.7538   -0.6533   -0.5528   -0.4523   -0.3518   -0.2513   -0.1508   -0.0503    0.0503    0.1508    0.2513    0.3518    0.4523

  Columns 106 through 120

    0.5528    0.6533    0.7538    0.8543    0.9548    1.0553    1.1558    1.2563    1.3568    1.4573    1.5578    1.6583    1.7588    1.8593    1.9598

  Columns 121 through 135

    2.0603    2.1608    2.2613    2.3618    2.4623    2.5628    2.6633    2.7638    2.8643    2.9648    3.0653    3.1658    3.2663    3.3668    3.4673

  Columns 136 through 150

    3.5678    3.6683    3.7688    3.8693    3.9698    4.0704    4.1709    4.2714    4.3719    4.4724    4.5729    4.6734    4.7739    4.8744    4.9749

  Columns 151 through 165

    5.0754    5.1759    5.2764    5.3769    5.4774    5.5779    5.6784    5.7789    5.8794    5.9799    6.0804    6.1809    6.2814    6.3819    6.4824

  Columns 166 through 180

    6.5829    6.6834    6.7839    6.8844    6.9849    7.0854    7.1859    7.2864    7.3869    7.4874    7.5879    7.6884    7.7889    7.8894    7.9899

  Columns 181 through 195

    8.0905    8.1910    8.2915    8.3920    8.4925    8.5930    8.6935    8.7940    8.8945    8.9950    9.0955    9.1960    9.2965    9.3970    9.4975

  Columns 196 through 200

    9.5980    9.6985    9.7990    9.8995   10.0000

FX1 =

  Columns 1 through 8

   0.8255 + 1.4298i   0.8224 + 1.4244i   0.8193 + 1.4191i   0.8162 + 1.4136i   0.8130 + 1.4082i   0.8098 + 1.4027i   0.8066 + 1.3971i   0.8034 + 1.3915i

  Columns 9 through 16

   0.8001 + 1.3859i   0.7968 + 1.3802i   0.7935 + 1.3744i   0.7902 + 1.3687i   0.7868 + 1.3628i   0.7834 + 1.3569i   0.7800 + 1.3510i   0.7765 + 1.3450i

  Columns 17 through 24

   0.7731 + 1.3390i   0.7695 + 1.3329i   0.7660 + 1.3267i   0.7624 + 1.3205i   0.7588 + 1.3143i   0.7551 + 1.3079i   0.7514 + 1.3015i   0.7477 + 1.2951i

  Columns 25 through 32

   0.7439 + 1.2886i   0.7401 + 1.2820i   0.7363 + 1.2753i   0.7324 + 1.2686i   0.7285 + 1.2618i   0.7245 + 1.2549i   0.7205 + 1.2480i   0.7165 + 1.2410i

  Columns 33 through 40

   0.7124 + 1.2338i   0.7082 + 1.2267i   0.7040 + 1.2194i   0.6998 + 1.2120i   0.6955 + 1.2046i   0.6911 + 1.1970i   0.6867 + 1.1894i   0.6822 + 1.1816i

  Columns 41 through 48

   0.6777 + 1.1738i   0.6731 + 1.1659i   0.6685 + 1.1578i   0.6637 + 1.1496i   0.6589 + 1.1413i   0.6541 + 1.1329i   0.6492 + 1.1244i   0.6441 + 1.1157i

  Columns 49 through 56

   0.6391 + 1.1069i   0.6339 + 1.0979i   0.6286 + 1.0888i   0.6233 + 1.0796i   0.6179 + 1.0702i   0.6123 + 1.0606i   0.6067 + 1.0508i   0.6009 + 1.0409i

  Columns 57 through 64

   0.5951 + 1.0307i   0.5891 + 1.0204i   0.5830 + 1.0098i   0.5768 + 0.9990i   0.5704 + 0.9880i   0.5639 + 0.9767i   0.5573 + 0.9652i   0.5504 + 0.9534i

  Columns 65 through 72

   0.5434 + 0.9413i   0.5363 + 0.9288i   0.5289 + 0.9160i   0.5213 + 0.9029i   0.5135 + 0.8893i   0.5054 + 0.8754i   0.4971 + 0.8609i   0.4884 + 0.8460i

  Columns 73 through 80

   0.4795 + 0.8305i   0.4702 + 0.8144i   0.4605 + 0.7977i   0.4504 + 0.7802i   0.4399 + 0.7619i   0.4288 + 0.7427i   0.4171 + 0.7224i   0.4047 + 0.7009i

  Columns 81 through 88

   0.3915 + 0.6780i   0.3773 + 0.6535i   0.3620 + 0.6269i   0.3452 + 0.5979i   0.3267 + 0.5658i   0.3057 + 0.5296i   0.2815 + 0.4875i   0.2521 + 0.4366i

  Columns 89 through 96

   0.2135 + 0.3698i   0.1512 + 0.2618i   0.2828 + 0.0000i   0.4177 + 0.0000i   0.4975 + 0.0000i   0.5576 + 0.0000i   0.6070 + 0.0000i   0.6494 + 0.0000i

  Columns 97 through 104

   0.6869 + 0.0000i   0.7207 + 0.0000i   0.7516 + 0.0000i   0.7802 + 0.0000i   0.8068 + 0.0000i   0.8317 + 0.0000i   0.8553 + 0.0000i   0.8776 + 0.0000i

  Columns 105 through 112

   0.8988 + 0.0000i   0.9191 + 0.0000i   0.9385 + 0.0000i   0.9572 + 0.0000i   0.9751 + 0.0000i   0.9924 + 0.0000i   1.0091 + 0.0000i   1.0253 + 0.0000i

  Columns 113 through 120

   1.0410 + 0.0000i   1.0562 + 0.0000i   1.0710 + 0.0000i   1.0855 + 0.0000i   1.0995 + 0.0000i   1.1132 + 0.0000i   1.1265 + 0.0000i   1.1396 + 0.0000i

  Columns 121 through 128

   1.1523 + 0.0000i   1.1648 + 0.0000i   1.1770 + 0.0000i   1.1890 + 0.0000i   1.2007 + 0.0000i   1.2122 + 0.0000i   1.2235 + 0.0000i   1.2346 + 0.0000i

  Columns 129 through 136

   1.2455 + 0.0000i   1.2562 + 0.0000i   1.2667 + 0.0000i   1.2771 + 0.0000i   1.2873 + 0.0000i   1.2973 + 0.0000i   1.3072 + 0.0000i   1.3169 + 0.0000i

  Columns 137 through 144

   1.3265 + 0.0000i   1.3360 + 0.0000i   1.3453 + 0.0000i   1.3545 + 0.0000i   1.3635 + 0.0000i   1.3725 + 0.0000i   1.3813 + 0.0000i   1.3901 + 0.0000i

  Columns 145 through 152

   1.3987 + 0.0000i   1.4072 + 0.0000i   1.4156 + 0.0000i   1.4239 + 0.0000i   1.4321 + 0.0000i   1.4402 + 0.0000i   1.4483 + 0.0000i   1.4562 + 0.0000i

  Columns 153 through 160

   1.4641 + 0.0000i   1.4718 + 0.0000i   1.4795 + 0.0000i   1.4871 + 0.0000i   1.4947 + 0.0000i   1.5021 + 0.0000i   1.5095 + 0.0000i   1.5168 + 0.0000i

  Columns 161 through 168

   1.5241 + 0.0000i   1.5313 + 0.0000i   1.5384 + 0.0000i   1.5454 + 0.0000i   1.5524 + 0.0000i   1.5593 + 0.0000i   1.5662 + 0.0000i   1.5730 + 0.0000i

  Columns 169 through 176

   1.5797 + 0.0000i   1.5864 + 0.0000i   1.5930 + 0.0000i   1.5996 + 0.0000i   1.6061 + 0.0000i   1.6126 + 0.0000i   1.6190 + 0.0000i   1.6254 + 0.0000i

  Columns 177 through 184

   1.6317 + 0.0000i   1.6380 + 0.0000i   1.6442 + 0.0000i   1.6503 + 0.0000i   1.6565 + 0.0000i   1.6626 + 0.0000i   1.6686 + 0.0000i   1.6746 + 0.0000i

  Columns 185 through 192

   1.6805 + 0.0000i   1.6865 + 0.0000i   1.6923 + 0.0000i   1.6982 + 0.0000i   1.7039 + 0.0000i   1.7097 + 0.0000i   1.7154 + 0.0000i   1.7211 + 0.0000i

  Columns 193 through 200

   1.7267 + 0.0000i   1.7323 + 0.0000i   1.7379 + 0.0000i   1.7434 + 0.0000i   1.7489 + 0.0000i   1.7544 + 0.0000i   1.7598 + 0.0000i   1.7652 + 0.0000i

Warning: Imaginary parts of complex X and/or Y arguments ignored 

k =

     1

x =

     0

x =

         0    0.7937

domx1 =

  Columns 1 through 25

     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0

  Columns 26 through 50

     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0

  Columns 51 through 75

     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0

  Columns 76 through 100

     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0     0

y1 =

  Columns 1 through 15

         0    0.0080    0.0160    0.0241    0.0321    0.0401    0.0481    0.0561    0.0641    0.0722    0.0802    0.0882    0.0962    0.1042    0.1122

  Columns 16 through 30

    0.1203    0.1283    0.1363    0.1443    0.1523    0.1603    0.1684    0.1764    0.1844    0.1924    0.2004    0.2084    0.2165    0.2245    0.2325

  Columns 31 through 45

    0.2405    0.2485    0.2565    0.2646    0.2726    0.2806    0.2886    0.2966    0.3047    0.3127    0.3207    0.3287    0.3367    0.3447    0.3528

  Columns 46 through 60

    0.3608    0.3688    0.3768    0.3848    0.3928    0.4009    0.4089    0.4169    0.4249    0.4329    0.4409    0.4490    0.4570    0.4650    0.4730

  Columns 61 through 75

    0.4810    0.4890    0.4971    0.5051    0.5131    0.5211    0.5291    0.5372    0.5452    0.5532    0.5612    0.5692    0.5772    0.5853    0.5933

  Columns 76 through 90

    0.6013    0.6093    0.6173    0.6253    0.6334    0.6414    0.6494    0.6574    0.6654    0.6734    0.6815    0.6895    0.6975    0.7055    0.7135

  Columns 91 through 100

    0.7215    0.7296    0.7376    0.7456    0.7536    0.7616    0.7696    0.7777    0.7857    0.7937

domx2 =

  Columns 1 through 15

         0    0.0080    0.0160    0.0241    0.0321    0.0401    0.0481    0.0561    0.0641    0.0722    0.0802    0.0882    0.0962    0.1042    0.1122

  Columns 16 through 30

    0.1203    0.1283    0.1363    0.1443    0.1523    0.1603    0.1684    0.1764    0.1844    0.1924    0.2004    0.2084    0.2165    0.2245    0.2325

  Columns 31 through 45

    0.2405    0.2485    0.2565    0.2646    0.2726    0.2806    0.2886    0.2966    0.3047    0.3127    0.3207    0.3287    0.3367    0.3447    0.3528

  Columns 46 through 60

    0.3608    0.3688    0.3768    0.3848    0.3928    0.4009    0.4089    0.4169    0.4249    0.4329    0.4409    0.4490    0.4570    0.4650    0.4730

  Columns 61 through 75

    0.4810    0.4890    0.4971    0.5051    0.5131    0.5211    0.5291    0.5372    0.5452    0.5532    0.5612    0.5692    0.5772    0.5853    0.5933

  Columns 76 through 90

    0.6013    0.6093    0.6173    0.6253    0.6334    0.6414    0.6494    0.6574    0.6654    0.6734    0.6815    0.6895    0.6975    0.7055    0.7135

  Columns 91 through 100

    0.7215    0.7296    0.7376    0.7456    0.7536    0.7616    0.7696    0.7777    0.7857    0.7937

y2 =

  Columns 1 through 15

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 16 through 30

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 31 through 45

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 46 through 60

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 61 through 75

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 76 through 90

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 91 through 100

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

k =

     2

x =

         0    0.7937    0.9644

domx1 =

  Columns 1 through 15

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 16 through 30

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 31 through 45

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 46 through 60

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 61 through 75

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 76 through 90

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

  Columns 91 through 100

    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937    0.7937

y1 =

  Columns 1 through 15

    0.7937    0.7954    0.7971    0.7989    0.8006    0.8023    0.8040    0.8058    0.8075    0.8092    0.8109    0.8127    0.8144    0.8161    0.8178

  Columns 16 through 30

    0.8196    0.8213    0.8230    0.8247    0.8265    0.8282    0.8299    0.8316    0.8333    0.8351    0.8368    0.8385    0.8402    0.8420    0.8437

  Columns 31 through 45

    0.8454    0.8471    0.8489    0.8506    0.8523    0.8540    0.8558    0.8575    0.8592    0.8609    0.8627    0.8644    0.8661    0.8678    0.8695

  Columns 46 through 60

    0.8713    0.8730    0.8747    0.8764    0.8782    0.8799    0.8816    0.8833    0.8851    0.8868    0.8885    0.8902    0.8920    0.8937    0.8954

  Columns 61 through 75

    0.8971    0.8989    0.9006    0.9023    0.9040    0.9058    0.9075    0.9092    0.9109    0.9126    0.9144    0.9161    0.9178    0.9195    0.9213

  Columns 76 through 90

    0.9230    0.9247    0.9264    0.9282    0.9299    0.9316    0.9333    0.9351    0.9368    0.9385    0.9402    0.9420    0.9437    0.9454    0.9471

  Columns 91 through 100

    0.9488    0.9506    0.9523    0.9540    0.9557    0.9575    0.9592    0.9609    0.9626    0.9644

domx2 =

  Columns 1 through 15

    0.7937    0.7954    0.7971    0.7989    0.8006    0.8023    0.8040    0.8058    0.8075    0.8092    0.8109    0.8127    0.8144    0.8161    0.8178

  Columns 16 through 30

    0.8196    0.8213    0.8230    0.8247    0.8265    0.8282    0.8299    0.8316    0.8333    0.8351    0.8368    0.8385    0.8402    0.8420    0.8437

  Columns 31 through 45

    0.8454    0.8471    0.8489    0.8506    0.8523    0.8540    0.8558    0.8575    0.8592    0.8609    0.8627    0.8644    0.8661    0.8678    0.8695

  Columns 46 through 60

    0.8713    0.8730    0.8747    0.8764    0.8782    0.8799    0.8816    0.8833    0.8851    0.8868    0.8885    0.8902    0.8920    0.8937    0.8954

  Columns 61 through 75

    0.8971    0.8989    0.9006    0.9023    0.9040    0.9058    0.9075    0.9092    0.9109    0.9126    0.9144    0.9161    0.9178    0.9195    0.9213

  Columns 76 through 90

    0.9230    0.9247    0.9264    0.9282    0.9299    0.9316    0.9333    0.9351    0.9368    0.9385    0.9402    0.9420    0.9437    0.9454    0.9471

  Columns 91 through 100

    0.9488    0.9506    0.9523    0.9540    0.9557    0.9575    0.9592    0.9609    0.9626    0.9644

y2 =

  Columns 1 through 15

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 16 through 30

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 31 through 45

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 46 through 60

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 61 through 75

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 76 through 90

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 91 through 100

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

k =

     3

x =

         0    0.7937    0.9644    0.9940

domx1 =

  Columns 1 through 15

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 16 through 30

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 31 through 45

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 46 through 60

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 61 through 75

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 76 through 90

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

  Columns 91 through 100

    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644    0.9644

y1 =

  Columns 1 through 15

    0.9644    0.9647    0.9650    0.9653    0.9656    0.9659    0.9662    0.9665    0.9668    0.9671    0.9674    0.9677    0.9680    0.9683    0.9686

  Columns 16 through 30

    0.9689    0.9692    0.9695    0.9698    0.9701    0.9704    0.9707    0.9710    0.9713    0.9716    0.9719    0.9722    0.9725    0.9728    0.9731

  Columns 31 through 45

    0.9734    0.9737    0.9739    0.9742    0.9745    0.9748    0.9751    0.9754    0.9757    0.9760    0.9763    0.9766    0.9769    0.9772    0.9775

  Columns 46 through 60

    0.9778    0.9781    0.9784    0.9787    0.9790    0.9793    0.9796    0.9799    0.9802    0.9805    0.9808    0.9811    0.9814    0.9817    0.9820

  Columns 61 through 75

    0.9823    0.9826    0.9829    0.9832    0.9835    0.9838    0.9841    0.9844    0.9847    0.9850    0.9853    0.9856    0.9859    0.9862    0.9865

  Columns 76 through 90

    0.9868    0.9871    0.9874    0.9877    0.9880    0.9883    0.9886    0.9889    0.9892    0.9895    0.9898    0.9901    0.9904    0.9907    0.9910

  Columns 91 through 100

    0.9913    0.9916    0.9919    0.9922    0.9925    0.9928    0.9931    0.9934    0.9937    0.9940

domx2 =

  Columns 1 through 15

    0.9644    0.9647    0.9650    0.9653    0.9656    0.9659    0.9662    0.9665    0.9668    0.9671    0.9674    0.9677    0.9680    0.9683    0.9686

  Columns 16 through 30

    0.9689    0.9692    0.9695    0.9698    0.9701    0.9704    0.9707    0.9710    0.9713    0.9716    0.9719    0.9722    0.9725    0.9728    0.9731

  Columns 31 through 45

    0.9734    0.9737    0.9739    0.9742    0.9745    0.9748    0.9751    0.9754    0.9757    0.9760    0.9763    0.9766    0.9769    0.9772    0.9775

  Columns 46 through 60

    0.9778    0.9781    0.9784    0.9787    0.9790    0.9793    0.9796    0.9799    0.9802    0.9805    0.9808    0.9811    0.9814    0.9817    0.9820

  Columns 61 through 75

    0.9823    0.9826    0.9829    0.9832    0.9835    0.9838    0.9841    0.9844    0.9847    0.9850    0.9853    0.9856    0.9859    0.9862    0.9865

  Columns 76 through 90

    0.9868    0.9871    0.9874    0.9877    0.9880    0.9883    0.9886    0.9889    0.9892    0.9895    0.9898    0.9901    0.9904    0.9907    0.9910

  Columns 91 through 100

    0.9913    0.9916    0.9919    0.9922    0.9925    0.9928    0.9931    0.9934    0.9937    0.9940

y2 =

  Columns 1 through 15

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 16 through 30

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 31 through 45

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 46 through 60

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 61 through 75

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 76 through 90

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 91 through 100

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

k =

     4

x =

         0    0.7937    0.9644    0.9940    0.9990

domx1 =

  Columns 1 through 15

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 16 through 30

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 31 through 45

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 46 through 60

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 61 through 75

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 76 through 90

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

  Columns 91 through 100

    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940    0.9940

y1 =

  Columns 1 through 15

    0.9940    0.9941    0.9941    0.9942    0.9942    0.9943    0.9943    0.9944    0.9944    0.9945    0.9945    0.9946    0.9946    0.9947    0.9947

  Columns 16 through 30

    0.9948    0.9948    0.9949    0.9949    0.9950    0.9950    0.9951    0.9951    0.9952    0.9952    0.9953    0.9953    0.9954    0.9954    0.9955

  Columns 31 through 45

    0.9955    0.9956    0.9956    0.9957    0.9957    0.9958    0.9958    0.9959    0.9959    0.9960    0.9960    0.9961    0.9961    0.9962    0.9962

  Columns 46 through 60

    0.9963    0.9963    0.9964    0.9964    0.9965    0.9965    0.9966    0.9966    0.9967    0.9967    0.9968    0.9968    0.9969    0.9969    0.9970

  Columns 61 through 75

    0.9970    0.9971    0.9971    0.9972    0.9972    0.9973    0.9973    0.9974    0.9974    0.9975    0.9975    0.9976    0.9976    0.9977    0.9977

  Columns 76 through 90

    0.9978    0.9978    0.9979    0.9979    0.9980    0.9980    0.9981    0.9981    0.9982    0.9982    0.9983    0.9983    0.9984    0.9984    0.9985

  Columns 91 through 100

    0.9986    0.9986    0.9987    0.9987    0.9988    0.9988    0.9989    0.9989    0.9990    0.9990

domx2 =

  Columns 1 through 15

    0.9940    0.9941    0.9941    0.9942    0.9942    0.9943    0.9943    0.9944    0.9944    0.9945    0.9945    0.9946    0.9946    0.9947    0.9947

  Columns 16 through 30

    0.9948    0.9948    0.9949    0.9949    0.9950    0.9950    0.9951    0.9951    0.9952    0.9952    0.9953    0.9953    0.9954    0.9954    0.9955

  Columns 31 through 45

    0.9955    0.9956    0.9956    0.9957    0.9957    0.9958    0.9958    0.9959    0.9959    0.9960    0.9960    0.9961    0.9961    0.9962    0.9962

  Columns 46 through 60

    0.9963    0.9963    0.9964    0.9964    0.9965    0.9965    0.9966    0.9966    0.9967    0.9967    0.9968    0.9968    0.9969    0.9969    0.9970

  Columns 61 through 75

    0.9970    0.9971    0.9971    0.9972    0.9972    0.9973    0.9973    0.9974    0.9974    0.9975    0.9975    0.9976    0.9976    0.9977    0.9977

  Columns 76 through 90

    0.9978    0.9978    0.9979    0.9979    0.9980    0.9980    0.9981    0.9981    0.9982    0.9982    0.9983    0.9983    0.9984    0.9984    0.9985

  Columns 91 through 100

    0.9986    0.9986    0.9987    0.9987    0.9988    0.9988    0.9989    0.9989    0.9990    0.9990

y2 =

  Columns 1 through 15

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 16 through 30

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 31 through 45

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 46 through 60

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 61 through 75

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 76 through 90

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 91 through 100

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

k =

     5

x =

         0    0.7937    0.9644    0.9940    0.9990    0.9998

domx1 =

  Columns 1 through 15

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 16 through 30

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 31 through 45

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 46 through 60

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 61 through 75

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 76 through 90

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

  Columns 91 through 100

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990

y1 =

  Columns 1 through 15

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9991    0.9991    0.9991    0.9991    0.9991    0.9991    0.9991    0.9991    0.9991

  Columns 16 through 30

    0.9991    0.9991    0.9991    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992

  Columns 31 through 45

    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9994    0.9994    0.9994

  Columns 46 through 60

    0.9994    0.9994    0.9994    0.9994    0.9994    0.9994    0.9994    0.9994    0.9994    0.9995    0.9995    0.9995    0.9995    0.9995    0.9995

  Columns 61 through 75

    0.9995    0.9995    0.9995    0.9995    0.9995    0.9995    0.9996    0.9996    0.9996    0.9996    0.9996    0.9996    0.9996    0.9996    0.9996

  Columns 76 through 90

    0.9996    0.9996    0.9996    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997

  Columns 91 through 100

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

domx2 =

  Columns 1 through 15

    0.9990    0.9990    0.9990    0.9990    0.9990    0.9990    0.9991    0.9991    0.9991    0.9991    0.9991    0.9991    0.9991    0.9991    0.9991

  Columns 16 through 30

    0.9991    0.9991    0.9991    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992    0.9992

  Columns 31 through 45

    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9993    0.9994    0.9994    0.9994

  Columns 46 through 60

    0.9994    0.9994    0.9994    0.9994    0.9994    0.9994    0.9994    0.9994    0.9994    0.9995    0.9995    0.9995    0.9995    0.9995    0.9995

  Columns 61 through 75

    0.9995    0.9995    0.9995    0.9995    0.9995    0.9995    0.9996    0.9996    0.9996    0.9996    0.9996    0.9996    0.9996    0.9996    0.9996

  Columns 76 through 90

    0.9996    0.9996    0.9996    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997    0.9997

  Columns 91 through 100

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

y2 =

  Columns 1 through 15

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 16 through 30

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 31 through 45

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 46 through 60

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 61 through 75

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 76 through 90

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 91 through 100

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

k =

     6

x =

         0    0.7937    0.9644    0.9940    0.9990    0.9998    1.0000

domx1 =

  Columns 1 through 15

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 16 through 30

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 31 through 45

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 46 through 60

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 61 through 75

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 76 through 90

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

  Columns 91 through 100

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998

y1 =

  Columns 1 through 15

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9999    0.9999    0.9999

  Columns 16 through 30

    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999

  Columns 31 through 45

    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999

  Columns 46 through 60

    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999

  Columns 61 through 75

    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999

  Columns 76 through 90

    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

domx2 =

  Columns 1 through 15

    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9998    0.9999    0.9999    0.9999

  Columns 16 through 30

    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999

  Columns 31 through 45

    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999

  Columns 46 through 60

    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999

  Columns 61 through 75

    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999

  Columns 76 through 90

    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    0.9999    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

y2 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

k =

     7

x =

         0    0.7937    0.9644    0.9940    0.9990    0.9998    1.0000    1.0000

domx1 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

y1 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

domx2 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

y2 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

k =

     8

x =

         0    0.7937    0.9644    0.9940    0.9990    0.9998    1.0000    1.0000    1.0000

domx1 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

y1 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

domx2 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

y2 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

k =

     9

x =

         0    0.7937    0.9644    0.9940    0.9990    0.9998    1.0000    1.0000    1.0000    1.0000

domx1 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

y1 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

domx2 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

y2 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

k =

    10

x =

         0    0.7937    0.9644    0.9940    0.9990    0.9998    1.0000    1.0000    1.0000    1.0000    1.0000

domx1 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

y1 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

domx2 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

y2 =

  Columns 1 through 15

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 16 through 30

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 31 through 45

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 46 through 60

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 61 through 75

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 76 through 90

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

  Columns 91 through 100

    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000    1.0000

k =

    11

ans =

         0
    0.7937
    0.9644
    0.9940
    0.9990
    0.9998
    1.0000
    1.0000
    1.0000
    1.0000
    1.0000

上图为迭代图像，可见最终将会收敛于一点。

时间： 2024-10-19 07:33:05

不动点迭代法的相关文章

MATLAB用二分法、不动点迭代法及Newton迭代（切线）法求非线性方程的根

一.实验原理 二.实验步骤 三.实验过程 1.(程序) (1)二分法:求   在区间(1,2)之间的根,取 (a)bipart.m: function [x,m]=bipart(fun,a0,b0,tol) a=a0;b=b0; m=1+round(round(log((b-a)/tol))/log(2)); for k=1:m p=(a+b)/2; if fun(p)*fun(b)<0 a=p; else b=p; end x=p; end end (b)fun1.m: function f=

R语言求根

求根是数值计算的一个基本问题,一般采用的都是迭代算法求解,主要有不动点迭代法.牛顿-拉富生算法.割线法和二分法. 不动点迭代法 所谓的不动点是指x=f(x)的那些点,而所谓的不懂点迭代法是指将原方程化为x=f(x)形式之后,下一步所用的x值为这一步的f(x),这样的话就可以一直逼近我们需         要的x,即方程的根,但是这种方法可能不会收敛到方程的根,随着初始值选定的大小,可能会有发散的情况,因此需要谨慎使用. ###不动点迭代法 func1 <- function(x){return(

牛顿法及其收敛性

目录 0. 引论 2. 牛顿法 2.1 单变量牛顿法 2.2 多变量牛顿法 3. 简化牛顿法 3.1 平行弦法 3.2 牛顿下山法 0. 引论 在工程中,很多的问题都可以归结为求解一组偏微分方程组.描述如下:在某一个区域\(\Omega\)中,物理量\(u={u_1,u_2,\cdots,u_n}\)满足一组控制方程,并且在边界上满足若干的边界条件: 控制方程利用两个微分算子来表示\(A(u),B(u)\). \[ A(u)=\begin{cases} A_1(u)& =0 \A_2(u)&

线性方程组之迭代法篇

不管哪一种数值算法,其设计原理都是将复杂转化为简单的重复,或者说,通过简单的重复生成复杂,在算法设计和算法实现过程中,重复就是力量[1].                                 ----题记 一般地,线性方程组可以表达为                                                                                   Ax = b其中,A称为系数矩阵,b称为右端项,x为待求的未知数向量.       迭代

牛顿迭代法求解平方根

牛顿迭代法求解平方根 2015-05-16 10:30 2492人阅读 评论(1) 收藏 举报 版权声明:本文为博主原创文章,未经博主允许不得转载. 目录(?)[+] 一个实例 迭代简介 牛顿迭代法 牛顿迭代法简介 简单推导 泰勒公式推导 延伸与应用 一个实例 //java实现的sqrt类和方法 public class sqrt { public static double sqrt(double n) { if (n<0) return Double.NaN; double err = 1e

C#跳转语句 迭代法 穷举法

一.跳转语句 break & continue break:跳出循环,终止此循环,不管下面还有多少次,全部跳过. string a=" ", for (int i=1;i<=10;I++) { if(i==5) { break; } a += i +",": } Console.WriteLine(a); 输出结果为 1,2,3,4,5 continue:终止此次循环,直接开始下次循环. string a=" ", for (int

经典算法：牛顿迭代法求平方根

//牛顿迭代法求平方根 1 double mysqrt(double num) 2 { 3 double x = num/2; 4 double y = 0; 5 do{ 6 x = x/2+num/(2*x); 7 y = x*x-num; 8 if(y<0) y = -y; 9 } while(y>0.0001); 10 return x; 11 } 12 int main(int argc, char* argv[]) 13 { 14 printf("%.3f",my

解方程——不动点迭代

定义:如果g(r)=r,那么实数r是函数g的一个不动点. 如果我们有方程f(x)=0,表示为不动点问题时,有:g(x)=x   (注:f(x)=g(x)-x) 不动点r是方程f(x)=0的一个根,几何表示为y=g(x)与y=x的交点就是不动点r 算法分析:x0 = 初始设定值 x1 = g(x0) x2 = g(x1) x3 = g(x2) ...    ... x(k+1) = g(x(k)) 直到收敛至g(r) = g(lim x(i)) = lim x(i+1) = r matlab代码实

51nod 1166 大数开平方(高精度+牛顿迭代法)

分析:直接用二分还是会T,用更快的牛顿迭代法.把问题转化为求x^2-n=0的根,假设解为x0,当前解为x且x^2-n>0,在(x,x^2-n)处作切线,与x轴交点横坐标为新的x,然后迭代即可,比二分法快,但是貌似只能用在凹函数或凸函数上.. java水高精度真是666... 1 import java.io.*; 2 import java.util.*; 3 import java.math.BigInteger; 4 public class Main { 5 public static v