Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement