Помогите с моделированием арбитражного уравнения

[Vladimir_0]

Необходимо написать функцию, которая возвращает значение (−p˙) как функцию от p и (−t).
Уравнение динамики рыночной стоимости акций r*p=d(t)+p ̇
r - рыночная ставка, равна 0,1
p - рыночная стоимость акций
d(t) - дивиденды
В условие задано, что при t>=20 d(t)=0.
Функция d(t) рассчитывается по следующему коду:

function d = dividend(t)
d = 1 * and(t >= 5, t < 10) + 0.02 * t .* (t < 5) + 0.1 * (15 - t) .* and(t >= 10, t < 15);

При этом должно быть удовлетворено условие: lim(|p(t)|)<∞ при t →∞
Какое тогда будет уравнение p(t)?
Пишу следующий код функции:

function ndp = neg_dot_p(neg_t,p)
%NEG_DOT_P
%функции возвращает значение (-pМ‡)
%как функцию от p и (-t)
 
%получаем доступ к глобальным переменным
global r
%используем уравнение r*p=d(t)+pМ‡
if (abs(neg_t)>=20)
    ndp = r*p;
else
    ndp =  r*p - dividend(neg_t); %
end
end

Функция должна проходить текст:

    for i = 1:1000;
      res = res + neg_dot_p(-i/100, i/1000);
    end
    if abs(res - 475) > 1
      fprintf('test not passed\n');
      error('Function "neg_dot_p" is not correct');
    end

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

[Vladimir_0]

Разобрался. Тему можно закрывать