请问如何用龙格库塔法求解下面的微分方程组：
dx1/dt=3*x1*x2+x1
dx2/dt=8*x2*x1+x2
t的范围从0-6
初值：x2(6)=0.6*x1(6)
x1(0)=0.8*x2(0)
可以通过龙哥库塔方法求数值解吗，可以的话怎么用matlab编程序 我不会编
谢谢