Матлаб - построить поверхность (пользов. функция)

[vedro-compota]

Здравствуйте уважаемые участники научного форума)

Пожалуйста - подскажите как решить такую проблему:
Мне необходимо построить поверхность (3-ёхмерный график) на котором была бы показана зависимость значения функции (пользовательской) от двух значений.

Так как наборы моих значений не равны по длине (1-ое изменяется , например так: [2,4,6....32], а второе от нуля до единицы изменяясь каждый раз на 0.0001), то - изрядно поискав в интернете я решил использовать функцию
в таком ключе:

clear();
  figure
  [S, n] = meshgrid(1::5, 1:10);
  Z = ppscq(S,n);
  surf(S, n, Z);

но тут с построением поверхности возникает проблема -
пользовательская функция ppscq() =

function R = ppscq(S, n)
   
    U = [];
    for k = 1:1:10
           disp 'ля-ля-ля';
    disp (S);
        S = (1/pi)*acos(cos(10^n)*S); U(k)=S; 
    end;
 
    R2 = S + n -S - n + preperiod(U);
    R = R2;
end

работает со значениями ,которые получаются на основе предыдущего -
а так как получается, что передаётся в качестве параметра матрица (судя по всему из-за этого), то происходит ошибка в строке:

 S = (1/pi)*acos(cos(10^n)*S); U(k)=S; 

Пожалуйста, подскажите, как можно исправить эту ошибку?

Заранее благодарю.