matlab求解拉普拉斯方程完整程序,急求用matlab编写解拉普拉斯方程的程序

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

这么处理狄拉克边界条件，我的精确解和数值解的误差很大部只哪错了

附程序

% examp1a2.m

clear all

clc

N=5;M=5;n=N*M;

L=2;h=2;

tic

c=sqrt((L/(N-1))^2+(h/(M-1))^2);

dm=2.5*c;

mu=0.25;

E=1e5;

beta=1*10^5;

for j=1:M

for i=1:N

x(i+(j-1)*N)=L/(N-1)*(i-1);

y(i+(j-1)*N)=h/(M-1)*(j-(M+1)/2);

end

end

t=[-0.8611363;-0.339880;0.339880;0.8611363];

G=[0.3478548;0.6521452;0.6521452;0.3478548];

for j=1:(M-1)

for i=1:(N-1)

for ky=1:4

for kx=1:4

xg(16*(i-1)+16*(j-1)*(N-1)+4*(ky-1)+kx)=(x(i+1+(j-1)*N)-x(i+(j-1)*N))/2*t(kx)+(x(i+1+(j-1)*N)+x(i+(j-1)*N))/2;

yg(16*(i-1)+16*(j-1)*(N-1)+4*(ky-1)&