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

[IRENN_Ka]

Здравствуйте! Помогите, пожалуйста подкорректировать программку для решения системы ДУ с граничными условиями на отрезке [7.9995, 8] Н2(7.9995)=0,Н3(7.9995)=0, H2(8)=0, H3(8)=0, функция Н1 изменяется произвольно. Код matlab:

function main
% инициализация:
solinit = bvpinit(linspace(7.9995,0.00002,8),[0 0]);
% linspace(7.9995,0.00002,8) - r изменяется в диапазоне от 7,9995 до 8
% [0 0] начальное приближение

% решение:
sol =bvp4c(@diffsys, @grn, solinit);

% графики:
plot(sol.R,sol.H);
legend('y(x)','dy/dx')
grid on
end

function dhdr = diffsys(r, h)
% обозначения:
% h(1) -> H1
% h(2) -> H2
% h(3) -> H3
% h(4) -> dH1/dr
% h(5) -> dH2/dr
% h(6) -> dH3/dr

% dhdr(1) -> dH1/dr
% dhdr(2) -> dH2/dr
% dhdr(3) -> dH3/dr
% dhdr(4) -> d2H1/dr2
% dhdr(5) -> d2H2/dr2
% dhdr(6) -> d2H3/dr2

w = 500000000;
gamma=1*10^7;
mu=(1.4/64000)*exp(-i*30);
k=1i*gamma*w*mu;
l=gamma*i*w;
a1 = (1.4/(64000^2))*exp(-i*30);
a2 = (1.4/(64000^3))*exp(-i*30);



dhdr = zeros(6,1);
dhdr(1) = h(4);
dhdr(2) = h(5);
dhdr(3) = h(6);
dhdr(4) = k*h(1) - h(4)/r;
dhdr(5) = 2*k*h(2) + 2*l*a1*h(1)^2 - h(5)/r;
dhdr(6) = 3*k*h(3) + 6*l*a1*h(1)*h(2) + 3*l*a2*h(1)^3 - h(6)/r;

end

function res = grn(ha , hb)
% граничные условия
grn = [ ha(2) % H2(7.9995)=0
hb(2) % H2(8)=0
ha(3) % H3(7.9995)=0
hb(3)]; %H3(8)=0

За ранее спасибо