Any ideal how combine step() and routh with array

So here my code

Kpid = [50,1,1,70,10];

Ti = [2,50,50,150,100];

Td = [0.5,5,2,2,2];

a = [0.05,0.05,0.05,0.05,0.05];

Kdt = [5,4,20,5,5];

T1 = [1,200,100,10,80];

T2 = [0.2,40,50,2,40];

Wpid = arrayfun(@(Kpid,Ti,Td,a,Kdt,T1,T2) tf([a*Kpid*Ti*Td Kpid*Ti+a*Td*Kpid+Ti*Td*Kpid Kpid],conv([Ti 0],[a*Td 1])),Kpid,Ti,Td,a,Kdt,T1,T2);

Wdt = arrayfun(@(Kdt,T1,T2) tf([Kdt],conv([T1 1],[T2 1])),Kdt,T1,T2);

Wpd = arrayfun(@(Wpid,Wdt) series(Wpid,Wdt),Wpid,Wdt);

W1 = [1,1,1,1,1];

W = arrayfun(@(Wpd,W1) feedback(Wpd,W1),Wpd,W1)

and then i want to use step(W) and stepinfo(W) but isnt working? I still dont know why

second code:

K1 = [25,25,25,25,25];

K2 = [8, 7, 9, 5, 6.5];

T1 = [1, 2, 0.5, 4, 5];

T2 = [0.1, 0.2, 0.4, 0.8, 0.5];

T3 = [0.4, 0.8, 0.05, 0.6, 0.2];

syms K ;

a = arrayfun(@(T1,T2,T3) (T1*T2+T1*T3+T2*T3)/(T1*T2*T3),T1,T2,T3);

b = arrayfun(@(T1,T2,T3) (T1+T2+T3)/(T1*T2*T3) ,T1,T2,T3);

c = arrayfun(@(T1,T2,T3,K1) K1/(T1*T2*T3) ,K1,T1,T2,T3);

d = arrayfun(@(T1,T2,T3) 1/(T1*T2*T3) ,T1,T2,T3);

how can i use routh in that situation with constant K

thanks everyone