function sum = func_integrate_viscdisc(n, mu, h, omega, rmin, rmax)
 sum = 0;
 rlen = rmax - rmin;
 dr = rlen/n;
 for idx = 1:n;
     r = rmin + dr*(idx - 0.5);
     ig = ((2*pi*mu*omega)/h)*(r^3);
     sum = sum + dr*ig;
 end


end