matlab解决矩量法的问题

问题是求解一个圆柱导体的电荷分布，使用矩量法，讨论单元数对电势的影响，请问其中的矩阵计算怎么编程？没有头绪啊 。。跪求大神指导