- 积分
- 7
- 注册时间
- 2002-9-10
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2010-6-17 22:48:15
|
显示全部楼层
来自 山东淄博
将第二个式子带入第一个可简化求解,Forcal代码:
- c: z2(z1) = (1+i*(3*i)^(1/3))*(1/2)-z1;
- c: cz1(z1)= z1^z2(z1) - z1*(z2(z1)+i*z1);
- c: absz1(z1)= abs[cz1(z1)];
- f(x1,x2,y1,y2)= //函数定义
- {
- cz1(x1,x2),
- GetCalByte(&y1,&y2)
- };
- mvar:
- a=1, b=1, //变换初值,得到不同的解
- array=new[rtoi(real_s),rtoi(4),rtoi(EndType),a,b].free(),
- i=XSLSF::netn[HFor("f"),1e-12,0.1,0.1,array,100], //也可以变换0.1,0.1两个值尝试
- array.sys::GA[0,&x,&y],
- printff{"\r\n初值a={1,r}, 初值b={2,r}, 实际迭代次数={3,i},误差={4,r}\r\nz1x={5,r},z1y={6,r},z2x={7,r},z2y={8,r}", a, b, i, absz1(x,y), x, y, z2(x,y), [GetCalByte(&x,&y),y]};
复制代码
几组解:
初值a=1., 初值b=1., 实际迭代次数=10,误差=2.7755575615628914e-016
z1x=-0.62095357379878169,z1y=0.35627563689818542,z2x=0.76039118122192972,z2y=0.2682367463435178
初值a=0., 初值b=0., 实际迭代次数=0,误差=0.
z1x=0.,z1y=0.,z2x=0.13943760742314798,z2y=0.62451238324170322
初值a=1., 初值b=-1., 实际迭代次数=5,误差=7.1993148024852709e-014
z1x=0.84598756933692421,z1y=-0.80817613912306552,z2x=-0.70654996191377628,z2y=1.4326885223647687
初值a=1., 初值b=3., 实际迭代次数=9,误差=4.5498153246749828e-014
z1x=0.45247348603451981,z1y=2.0021798049449022,z2x=-0.31303587861137183,z2y=-1.377667421703199
初值a=1., 初值b=5., 实际迭代次数=6,误差=7.7475686966118219e-013
z1x=1.2689825765920968,z1y=4.5996065256791203,z2x=-1.1295449691689488,z2y=-3.9750941424374169
初值a=1., 初值b=7., 实际迭代次数=8,误差=4.9769702052493514e-014
z1x=1.9680856685532386,z1y=6.6375164068950596,z2x=-1.8286480611300906,z2y=-6.0130040236533562
初值a=1., 初值b=9., 实际迭代次数=11,误差=1.1457157353758233e-013
z1x=2.5858464044055864,z1y=8.4819624180205597,z2x=-2.4464087969824386,z2y=-7.8574500347788563
初值a=-3., 初值b=1., 实际迭代次数=8,误差=1.8820163294834461e-014
z1x=-2.6407595502911008,z1y=0.44564314682191952,z2x=2.7801971577142486,z2y=0.17886923641978369 |
|