请问牛顿迭代法求解三元非线性方程组，不同迭代初值貌似有很多不同结果，如何求解到满足要求的解，0<z1<z2<18,0.004<t0<0.008,程序如下(尝试了一下割线法也不对)，请教大佬如何做些修改满足如上要求?
FindRoot[{-6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) z1 -
1.3980263449472176`*^-10 (1.438702749038092`*^8 +
1.2840021082598146`*^7 t0) z1 +
0.00011254833333333333` z1^3 - (
0.0005367544235159816` (E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) + (
0.0005367544235159816` (-E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) - t0 == 0,
1.3980263449472176`*^-10 (1.438702749038092`*^8 +
1.2840021082598146`*^7 t0) +
0.000338625` z1^2 - (-6.565159793181527`*^-6 \
(0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 - (
0.00021470176940639266` (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) + (
0.00021470176940639266` (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2))) == 0,
0.00031751999999999996` -
6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 - 9.8`*^-7 z2^2 - (
6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (-E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) - (
6.606735159817352`*^-7 (E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(
0.8` z2)) - (-0.000031536` (-0.65638188` +
0.00011320166666666666` z2^3 - (
0.0005367544235159816` (E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) - (
0.0005334510559360729` (-E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) +
18 (-6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 - 9.8`*^-7 z2^2 - (
6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (-E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (
6.606735159817352`*^-7 (E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2))) +
z2 (-6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 - 9.8`*^-7 z2^2 - (
6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (-E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (
6.606735159817352`*^-7 (E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2))) +
18 (-6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 + 0.000338625` z2^2 - (
6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (
0.00021404109589041095` (-E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
0.00021404109589041095` (E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2))) -
z2 (-6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 + 0.000338625` z2^2 - (
6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (
0.00021404109589041095` (-E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
0.00021404109589041095` (E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2))))) == 0}, {z1, 3, 6}, {z2, 9,
17}, {t0, 0, 0.00837}]
{z1 -> 10.9876, z2 -> 14.529, t0 -> -0.169798}
FindRoot[{-6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) z1 -
1.3980263449472176`*^-10 (1.438702749038092`*^8 +
1.2840021082598146`*^7 t0) z1 +
0.00011254833333333333` z1^3 - (
0.0005367544235159816` (E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) + (
0.0005367544235159816` (-E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) - t0 == 0,
1.3980263449472176`*^-10 (1.438702749038092`*^8 +
1.2840021082598146`*^7 t0) +
0.000338625` z1^2 - (-6.565159793181527`*^-6 \
(0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 - (
0.00021470176940639266` (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) + (
0.00021470176940639266` (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2))) == 0,
0.00031751999999999996` -
6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 - 9.8`*^-7 z2^2 - (
6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (-E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) - (
6.606735159817352`*^-7 (E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(
0.8` z2)) - (-0.000031536` (-0.65638188` +
0.00011320166666666666` z2^3 - (
0.0005367544235159816` (E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 - 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) - (
0.0005334510559360729` (-E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 + 3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2)) +
18 (-6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 - 9.8`*^-7 z2^2 - (
6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (-E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (
6.606735159817352`*^-7 (E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2))) +
z2 (-6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 - 9.8`*^-7 z2^2 - (
6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (-E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (
6.606735159817352`*^-7 (E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2))) +
18 (-6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 + 0.000338625` z2^2 - (
6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (
0.00021404109589041095` (-E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
0.00021404109589041095` (E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2))) -
z2 (-6.565159793181527`*^-6 (0.008850684113090817` -
0.7905181223796485` (0.008850684113090817` - t0)) +
9.8`*^-7 z1^2 + 0.000338625` z2^2 - (
6.606735159817352`*^-7 (-E^(-0.4` z1) + E^(
0.4` z1)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) - (
0.00021404109589041095` (-E^(-0.4` z2) + E^(
0.4` z2)) (-3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 -
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
6.606735159817352`*^-7 (E^(-0.4` z1) + E^(
0.4` z1)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(E^(0.8` z1) - E^(0.8` z2)) + (
0.00021404109589041095` (E^(-0.4` z2) + E^(
0.4` z2)) (3.696` E^(0.4` z1) z1 +
3.696` E^(0.4` z1 + 0.8` z2) z1 +
3.696` E^(0.8` z1 + 0.4` z2) z2 +
3.696` E^(0.4` z2) z2))/(
E^(0.8` z1) - E^(0.8` z2))))) == 0}, {z1, 3, 6}, {z2, 9,
17}, {t0, 0, 0.00837}]
{z1 -> 10.9876, z2 -> 14.529, t0 -> -0.169798}
