@hooshmand 

Certainly the last command was intended to be implicitplot, not plots. But for the plot to work, ellips would need to be a function of x and y. However, ellips is defined on the second-to-last line and it doesn't depend on any variable.

Perhaps we could help you if you could describe in words what you're trying to do.

@testht06 

There are several ways in Maple to represent the field. The way that I used uses RootOf instead of GF.

alias(x= RootOf(y^8+y^4+y^3+y+1)):
M:= Matrix([[0,1,0,0],[0,0,1,0],[0,0,0,1],[1,x,1,x^2]]):
eig:= Eigenvalues(M) mod 2;

Eigenvector(M, eig[1,1]) mod 2;

Eigenvector(M, eig[2,1]) mod 2;

Where does this code come from? It requires significant knowledge of Maple to write code this bad.

@tomleslie 

The fact that polynomials of degree greater than four over fields of characteristic zero are in general insoluble is irrelevant. All roots of any polynomial over any finite field can be computed, and there are efficient algorithms to do it. Indeed, the complete factorization into irreducibles of any polynomial over any finite field can be computed, and this is the basis for algorithms that factor polynomials over the rationals.

Newton's method works in modular arithmetic, and it converges to the exact answer.

@reinhardsiegfried 

The {x,y} is the second argument to coeffs. If map is called with three arguments, as in map(f, a, b), then b is used as the second argument to f. See ?coeffs and ?map.

@zia9206314 The values of the roots computed by NeztZero in Maple 17 and Maple 18 are identical down to the last digit. It is only the computation of the residuals that is different. This is probably because of a bug in Maple 17 that was fixed in Maple 18.

No, I can't send you Maple 18.

evaluation.mw

@ecterrab I didn't say that there were several bugs; I said that I'd seen several bug reports. All the reports may have been about the same bug.

It must've been Maple 18 where that bug was fixed. Executing the worksheet evaluation.mw (from the head of this subthread) in Maple 18 and Maple 17 shows huge differences in evaluation of a HeunB function such that the residuals from both root finders (NextZero and fsolve) make no sense in Maple 17. Since the OP is using Maple 17, I think that the Heun functions can be blamed.

@zia9206314 

I could not duplicate your results. I executed your worksheet in Maple 18, and all 14 roots found by NextZero have residuals less than 5*10^(-Digits). Try executing the worksheet again. And what version of Maple are your using?

That being said, I have seen several bug reports about numerical evaluation of Heun functions. So, if there is actually a discrepancy with the residuals, I'd blame that instead of blaming NextZero or fsolve.

By the way, guardDigits= 22 is excessive. At Digits=15, guardDigits= 4 is enough to get the first 14 roots.

@Markiyan Hirnyk 

It takes guardDigits= 9 to get the first 20 roots.

I think that the OP was just blindly copying my code from an earlier problem and didn't really want 20 roots.

@ecterrab 

I don't see how evalb makes any difference. It's the default boolean evaluator anyway, and of course it's idempotent. But perhaps it's fine if it errors out from a FAIL---it's a sign that the parameters to NextZero need to be adjusted.

If your Matrix has 8 bytes per element, that would be 8*63001^2 ~ 32 GB. Does your computer have that much free memory? Since you're only using Digits = 5, you can use 4-byte floats in the matrix by declaring

m_Gi:= Matrix(num_Strat,  num_Strat, datatye= float[4]): 

Since the columns in your matrix can be generated on the fly, you may be able to use black-box methods, i.e., the columns of the matrix are generated as needed and discarded. How do you plan on using this matrix?

@shzan I already told you how to make it constant! Use convert(..., D).

Alejandro's Answer above explains the source of the problem, not how to correct it. If you are a new Maple user, I wouldn't try to understand it.

@Mac Dude I believe that the OP's problem is that the eval(diff(...)) form is not treated as a constant with respect to integration.

Compare

int(diff(f(x), x), x= -1), t)

with

int(D(f)(-1), t)

@shzan If Ans is the answer returned by the program, do convert(Ans, D).

@Markiyan Hirnyk 

That's very odd because I wrote the code in Maple 2015 Standard, mode= worksheet. What happens for you when you run it? Are you using 1D or 2D input?

What about you, Rouben? Can you run it?