Вы можете написать функцию-обертку, которая будет принимать вектор в качестве входного аргумента и «распаковывать» этот вектор перед вызовом вашей функции.
f(y1,y2,y3,y4,z1) = (-(z1^3)/ ...; % your function
fv(v) = f(v(1),v(2),v(3),v(4),v(5));
Или, альтернативно, заставьте f ссылаться на элементы входного вектора, используя соответствующие индексы:
f(v) = (-(v(5)^3)/(v(3)^2)); ... % etc.