Спасибо,
ewert. К сожалению, пока техническая проблема не решается.
Опишу ее подробнее. Быть может совместными усилиями удастся придумать решение.
Есть код:
Код:
function [u,v,p] = Navier_Stokes_thoer( x, y, t, Rey )
ut = inline('-exp(-t*2/Rey)*cos(x)*sin(y)','x','y','t','Rey');
vt = inline('0','x','y','t','Rey');
pt = inline('1','x','y','t','Rey');
m = length(x);
n = length(y);
u = zeros(m,n);
v = zeros(m,n);
p = zeros(m,n);
for k = 1:n
u(:,k) = ut(x,y(k),t,Rey);
v(:,k) = vt(x,y(k),t,Rey);
p(:,k) = pt(x,y(k),t,Rey);
end
Здесь записана функция, которая получает
x,y, t, Rey и выдает
u,v,p. Вместо
ut я хочу вставить разрывную функцию типа той, что описана выше. В инлайне такая запись не работает. Править и менять эту часть скрипта, к сожалению, не вариант, потому что векторный синтаксис этой маленькой функции используется все время дальше в большом скрипте...