| 1 | initial version |
This is expected as t=0 makes the denominator zero. You need to not substitute t=0 but find the limit as t tends to 0 here:
P.right_kernel().matrix()[0].apply_map(lambda z: z.limit(t=0))
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.