$\frac{d{x}_1}{dt}=3x^2+2x_2+u_1=f(x_1,x_2,n_i)$

$\frac{d{x}_3}{dt}=-1\mathrm{.}5x_1+5x_2^2+H_t=f_t(x_1,x_2,H_2)$

$k{x}_1=$I,$f(x_1,x_2,H_1)$

$k{x}_2=T{f}_1(x_1,x_2,n_2)$

$k_2{x}_1=$I,$f_f(x_1+\frac{k{x}_1}{2},x_2+\frac{k{x}_2}{2},x_1)$

$k_1{x}_2=I_2{f}_2(x_1+\frac{k{x}_1}{2},x_2+\frac{k{x}_2}{2},{}_2)$

$k{x}_1=I{f}_f(x_1+\frac{k{x}_1}{2},x_2+\frac{k_{x_2}}{2},u_1)$

$k_2{x}_2=I_1{f}_2(x_1+\frac{k_1{x}_1}{2},x_2+\frac{k_2{x}_2}{2},u_2)$

$k_1{x}_1=I_{f_1}(x_1+k_1{x}_1,x_3+k_1{x}_2,u_1)$

$k_1{x}_1=I{f}_2(x_1+k_3{x}_1,x_2+k{x}_2,u_2)$

$x_1[n+1]=x_1[n]+\frac{1}{6}(k{x}_i+2k_2{x}_i+2k_3{x}_i+k_3{x}_i)$

$x_3[n+1]=x_1[n]+\frac{1}{6}(k_1+2k_2{x}_2+2k_3{x}_2+k_2{x}_2)$

$\frac{d{x}_1}{dt}=3x_1^2+2x_2+u_4=f(x_1,x_2,u_i)$

$\frac{d{u}_2}{dt}=-1\mathrm{.}5x_1+5x_2^2+u_2=f_2(x_1,x_2,u_2)$

$x_1(1)=0\mathrm{.}5$

$x_2(1)=0\mathrm{.}5$

$u_1(1$:$1000)=1\mathrm{.}2$;

$u_2(1$:$1000)=-1\mathrm{.}2$;

for $n=1$:$1000$

$k_1{x}_1(n)=T_1{f}_1(x_1(n),x_2(n{)}_n{u}_1(n))$

$k_1{x}_2(n)=I_1{f}_2(x_1(n),x_2(n),u_2(n))$

$k_2{x}_1(n)=I{f}_1(x_1(n)+\frac{k_1{x}_1(n)}{2},x_2(n)+\frac{k_1{x}_2(n)}{2},n_1(n))$

$k_2{x}_2(n)=I{f}_2(x_1(n)+\frac{k_1{x}_1(n)}{2},x_2(n)+\frac{k_1{x}_2(n)}{2},u_2(n))$

$21$

$k_1{x}_1(n)=T{f}_1(x_1(n)+\frac{k_1{x}_1(n)}{2},x_2(n)+\frac{k_2{x}_2(n)}{2},u_1(n))$

$k_1{x}_2(n)=T_n{f}_2(x_1(n)+\frac{k_2{x}_1(n)}{2},x_2(n)+\frac{k_1{x}_2(n)}{2},n_2(n))$

$k_1{x}_1(n)=T{f}_1(x_1(n)+k_1{x}_1(n),x_2(n)+k_1{x}_2(n),u_1(n))$

$k_4{x}_3(n)=T{f}_2(x_1(n)+k_1{x}_1(n),x_2(n)+k_1{x}_2(n),u_2(n))$

$x_i(n+1)=x_i(n)+\frac{1}{6}(k_1{x}_i(n)+2k_x{x}_i(n)+2k_i{x}_i(n)+k_i{x}_i(n))$

$x_3(n+1)=x_3(n)+\frac{1}{6}(k_1{x}_2(n)+2k_2{x}_2(n)+2k_3{x}_2(n)+k_3{x}_2(n))$

end

function $[dx1]={}_{-}$finction $(x1,x2,u1)$

$[$:$1=3^{**}{1}^{**}1+2^{**}2+u$