Как построить график функции Бесселя

[IRENN]

После символьного решения дифференциального уравнения получается сумма модифицированных уравнений Бесселя. Но график почему-то не строится =(

H1=dsolve('D2H1+(1/r)*DH1=H1','r')
ezplot('H1', [-100 100])
hold on
grid on

ezplot('imag(H1)', [-100 100]),
ezplot('abs(H1)', [-100 100])
xlabel('r')
ylabel('H1')
legend('real','image','abs')