matlab最简单的代码-juliafem:二维Poisson方程的简单FEM代码

matlab最简单的代码该代码是对“大约50行Matlab的注释：简短的有限元实现”的Matlab代码的重写。 用茱莉亚语言。 对于matlab版本，请参见 去做 添加四边形的绘图解决方案。 Matlab代码使用trisurf，可以绘制三角形和四边形。 我在PyPlot中找不到等效的功能。 我通过PyPlot的tricontour函数在Julia代码中使用等高线图。 如果您在“样本”中运行示例，它将通过对它们进行三角剖分来绘制四边形。 如何在“正方形”中运行示例 您需要安装PyPlot软件包以绘制结果。 import Pkg Pkg.add("PyPlot") 首先，生成网格 bash> cd square 通过在Matlab中运行square.m生成网格（julia版本未完成） matlab> square(30) matlab> quit 运行fem代码 bash> julia run.jl 你应该得到这样的解决方案 构造稀疏矩阵 初始化空的nxn稀疏矩阵的最简单方法是 julia> A = spzeros(n,n) 这利用了64个字节的整数和浮点数。 对于较小的问题，使用32个字