# Problem with matrix echelon form

Dear all,

if we define A to be

```
A = matrix([[1,-3,0,2],[-2,4,0,h],[0,-2,1,k],[3,-1,7,1]])
A.echelon_form()
```

gives diag(1,1,1,1). But this is wrong for some values of h and k.In some cases the last 1 could be 0.

How can I fix this? Thanks, N

A possible approach would be defining a function of

`h`

and`k`

.Sorry not clear enough for me … What do you mean?

Like this:

I see. Thanks. But in this case I have to enter the values for h and k. That's not what I want... I'd like Sage to tell me what condition on h and k gives an echelon form like drag(1,1,1,1) and which one gives drag(1,1,1,0).

diag not drag :-)