%FemReadDat

N=22;
str='dat.txt';
fid=fopen(str,'r');

M=fscanf(fid,'%f',[N N]);
fclose(fid);

nx=N-1;ny=N-1;
dx=1/nx;
dy=1/ny;
x=0:(1/nx):1;
y=0:(1/ny):1;
[X,Y]=meshgrid(x,y);
subplot(2,2,1)
Mesh(X,Y,M)
xlabel('X-axis'),ylabel('Y-axis'),zlabel('溫度T-axis')
title('有限元法解二維熱傳導問題計算結果數據可視化 圖 1')

subplot(2,2,2)
Mesh(X,Y,M)
xlabel('X-axis'),ylabel('Y-axis'),zlabel('溫度T-axis')
title('圖 2')
view(30,20)
subplot(2,2,3)
Mesh(X,Y,M)
xlabel('X-axis'),ylabel('Y-axis'),zlabel('溫度T-axis')
title('圖 3')
view(-30,45)
subplot(2,2,4)
surfl(X,Y,M)
view(-60,45)
shading interp
xlabel('X-axis'),ylabel('Y-axis'),zlabel('溫度T-axis')

title('圖 4')