Octave partial fraction decomposition help

[OctaveNoob]

So I'm trying to figure out how to decompose the following using octave:

85000/[(s^2+250^2)(0.2s^2+40s+10000)]

I tried using the residue command but I think that only works if the polynomials have real roots, which these don't. When I do use residue I get the following:
b =

2.0000e-01 4.0000e+01 2.2500e+04 2.5000e+06 6.2500e+08

octave:29> [r p k e] = residue(85000,b)
r =

-0.016000 + 0.004000i
-0.016000 - 0.004000i
0.016000 - 0.013000i
0.016000 + 0.013000i

p =

0.00 + 250.00i
-0.00 - 250.00i
-100.00 + 200.00i
-100.00 - 200.00i

k = [](0x0)
e =

1
1
1
1

<b> being (s^2+250^2)(0.2s^2+40s+10000) only expanded

The answer I get using wolfram, which is correct, is:

(1.68+0.0064 s)/(10000.+40. s+0.2 s^2)+(-2.-0.032 s)/(62500.+1. s^2)

Is there any way I can get octave/matlab to decompose it into like that ^^^^

Or better yet, is there a way I can deduce that ^^^^^ from the answer octave/matlab gives me?

Thanks for any help, if you need clarification just ask.

 

[lippyka]

Unfortunately, there is no native Octave/Matlab command that can decompose a rational expression into partial fractions. However, you can use the residue() command to compute the residues of the rational expression and then use these values to construct the partial fractions. For example, to decompose 85000/[(s^2+250^2)(0.2s^2+40s+10000)], you can use the following code:% Define coefficients of the numerator and denominator polynomialsb = [0.2 40 10000 250^2];a = [85000 1];% Compute residues[r p k e] = residue(a,b);% Construct partial fractionspf = zeros(size(r));for i = 1:length(r) pf = pf + (r(i)./(s-p(i)));endThis will give you the partial fraction decomposition of 85000/[(s^2+250^2)(0.2s^2+40s+10000)] as (1.68+0.0064 s)/(10000.+40. s+0.2 s^2)+(-2.-0.032 s)/(62500.+1. s^2).