Решить систему дифференциальных уравнений

@coldunox · Регистрация: 03.11.2016

Author24 — интернет-сервис помощи студентам

Система состоит из 14 уравнений (скрин 1)
Есть решение в MATLAB
Как его переписать в MathCAD?
Мои попытки (скрины 2 - 3)

Matlab M

function dx=calculate(t,x)
dx=zeros(12,1);
 
res = 1+0.4*sin(t);
price = 0.25*t*t-1.1*t+2+0.5*t*sin(t);
foreign = 0.25*t*t-0.2*t+1;
obs = 1+0.05*t;
exper = 1+0.7*cos(t);
port = 1+0.3*sin(t);
gov = 1+0.4*sin(t);
 
%dx(1) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
%dx(2) = exper + 1.1 * x(1) + 0.01 - ((2/(x(4) - 0.001) + 0.0001)) * ((2/(x(6) - 0.001) + 0.0001)) * ((2/(x(5) - 0.001) + 0.0001));  
dx(1) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(2) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(3) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(4) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(5) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(6) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(7) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(8) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(9) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(10) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(11) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(12) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(13) = exper - ((2/(x(4) - 0.001) + 1.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
dx(14) = exper - ((2/(x(4) - 0.001) + 0.0001) * (2/(x(5) - 0.001) + 0.0001) * (2/(x(6) - 0.001) + 0.0001));
end

@mathmichel · 21.09.2017, 08:19

Кроме самих дифференциальных уравнений ещё надо задать начальные условия для всех 14 интегрируемых функций. Судя по скринам, Вам необходимо все уравнения и выражения вставить между Given и Find. Уравнения задаются не через знак присвоения ":=", а через жирный знак равенства <ctrl>+<=>. Ваши попытки в Mathcad желательно заархивировать и загрузить сюда.

@mathmichel $Эксперт по математике/физике$ 10452 / 6933 / 3772 Регистрация: 14.01.2014 Сообщений: 15,920
	21.09.2017, 08:19	2
	Кроме самих дифференциальных уравнений ещё надо задать начальные условия для всех 14 интегрируемых функций. Судя по скринам, Вам необходимо все уравнения и выражения вставить между Given и Find. Уравнения задаются не через знак присвоения ":=", а через жирный знак равенства <ctrl>+<=>. Ваши попытки в Mathcad желательно заархивировать и загрузить сюда. 1