2625,2627d2624
<   if (numel () == 1)
<     return ComplexMatrix (1, 1, exp (m(0)));
< 
2642c2639,2643
<     trshift = trshift.imag ();
---
>     {
>       trshift = trshift.imag ();
>       if (trshift.real () > 709.0)
> 	trshift = 709.0;
>     }
2727,2728c2728,2734
<       npp = m * npp + m * padec[j];
<       dpp = m * dpp + m * (minus_one_j * padec[j]);
---
>       for (octave_idx_type i = 0; i < nc; i++)
> 	{
> 	  npp.elem (i, i) += padec[j];
> 	  dpp.elem (i, i) += minus_one_j * pade0c[j];
> 	}      
>       npp = m * npp;
>       dpp = m * dpp;