Adam Majewski <adammaj1 <at> o2.pl> writes:
> Curve3b:GiveListAngle(3,1,center3[2],50)
> "mnewton: the process doesn't converge or it converges too slowly."
Well, trace(mnewton) shows that it processes several problems
successfully but then it fails on this one:
mnewton ([z^4+2*c*z^2-z+c^2+c, 4*z^3+4*c*z - %i*sin(2.0*%pi) - cos(2.0*%pi)],
[c, z], [-1, 0]);
=> mnewton: the process doesn't converge or it converges too slowly.
But I find that if you just change the initial point, it succeeds:
mnewton ([z^4+2*c*z^2-z+c^2+c, 4*z^3+4*c*z - %i*sin(2.0*%pi) - cos(2.0*%pi)],
[c, z], [0, 1]);
=> [[c = 7.37444285991106E-34*%i+0.25,z = 0.5-6.122503540205582E-17*%i]]
Maybe the solution of the preceding problem is a good starting
point for the next one? I don't know.
As a separate issue, maybe you want to make the increment of each
step a rational number instead of float, so that you would get
sin(2*%pi) and cos(2*%pi) which will simplify. I don't think it
will make much practical difference.
Hope this helps,
Robert Dodier