Переделать код под С++

t=0:0.1:1; 
f=zeros(2,11); 
xt=exp(t)+exp(-3*t)+t.^2-t+9;
f(1,1)=11; f(2,1)=-3;
h=0.1;
k1=@(t,x,y)h*([y;-2*y+3*x-3*(t^2)+7*t-27]);
for i=1:1:11
a1=k1(t(i),f(1,i),f(2,i));
a2=k1(t(i)+h/2,f(1,i)+a1(1)/2,f(2,i)+a1(2)/2);
a3=k1(t(i)+h/2,f(1,i)+a2(1)/2,f(2,i)+a2(2)/2);
a4=k1(t(i)+h,f(1,i)+a3(1),f(2,i)+a3(2));
f(:,i+1)=f(:,i)+(a1+2*a2+2*a3+a4)/6;
end
hold on; grid on;
plot(t(1:11),f(1,1:11),'o',t,xt,'-')
legend('Метод Рунге-Кутта','Точное решение')

t=0:0.1:1; // параметр, с шагом 0.1
f=zeros(2,11);
xt=exp(t)+exp(-3*t)+t.^2-t+9; // функция (решение дифференциального уравнения)
f(1,1)=11; f(2,1)=-3; // (начальные условия (f(1,1) — значение функции в нуле, f(2,1) — значение производной функции в нуле))
h=0.1; // шаг 
k1=@(t,x,y)h*([y;-2*y+3*x-3*(t^2)+7*t-27]); //функция для расчета приближенного значения задачи Коши
for i=1:1:11
a1=k1(t(i),f(1,i),f(2,i)); // 1 (ниже 2, 3 и 4 коэффициенты (Рунге-Кутта 4 порядка применяется))
a2=k1(t(i)+h/2,f(1,i)+a1(1)/2,f(2,i)+a1(2)/2);
a3=k1(t(i)+h/2,f(1,i)+a2(1)/2,f(2,i)+a2(2)/2);
a4=k1(t(i)+h,f(1,i)+a3(1),f(2,i)+a3(2));
f(:,i+1)=f(:,i)+(a1+2*a2+2*a3+a4)/6; // нахождение точек функции при конкретных значениях параметра и х
end
hold on; grid on;
plot(t(1:11),f(1,1:11),'o',t,xt,'-') // построение графиков
legend('Метод Рунге-Кутта','Точное решение') // легенда