Finite Element Toolbox 2.1: meshtest.m - File Exchange

No BSD License

Highlights from
Finite Element Toolbox 2.1

[A,B,b]=genmat1(xn,a,b,s) GENMAT1 Generates FEM1 matrices.
[A,B,b]=genmat2(xy,el,a,b,s) GENMAT2 Generates FEM2 matrices.
[A,B,b]=genmat3(xyz,el,a,b,s) GENMAT3 Generates FEM3 matrices.
[nod2xy,el2nod,bound,W,G]=imp...
[u,x,y,c]=fem2(nod2xy,el2nod,... FEM2 Finite element solver in 2D.
[xnh,ah,bh,sh]=adapt1(xn0,tol... ADAPT1 Generate adaptively refined grid.
el2nod=genmesh2(nod2xy,boundp... GENMESH2 Generate 2D mesh from grid points.
ex1() This is the machine-generated representation of a Handle Graphics object
ex2() This is the machine-generated representation of a Handle Graphics object
plotgrid2(nod2xy,el2nod,fs) PLOTGRID Plots FE nodes and mesh.
u=fem1(xn, alpha,beta,s, bd) FEM1 Finite element solver in 1D.
x=quadspace(a,b,n) QUADSPACE Quadratically spaced vector.
xn=refine1(pts,n) REFINE1 Refine grid uniformly around critical points.
y=vmean(x) VMEAN Mean value of vector elements.
meshtest.m
test1.m
test2.m
View all files

from Finite Element Toolbox 2.1 by Rasmus Anthin
Piecewise linear Galerkin method with piecewise constant coeffs for static PDE probelms.

meshtest.m

x0=linspace(-1,1,10);
alpha=inline('x | 1');
beta=inline('x | 1');
source=inline('x & 0');
bd=[-1 nan nan;1 nan nan];

u=inline('exp(1)/(1-exp(2))*(exp(-x)-exp(x))');
tol=input('tol: ');
[x1,a1,b1,s1]=adapt1(x0,tol,alpha,beta,source,bd);
U1=fem1(x1,a1,b1,s1,bd).';
figure,plot(x1,U1,'.-',x1,u(x1))
figure,plot(x1,abs(u(x1)-U1));
figure,plot(diff(x1),'.-')

Highlights from Finite Element Toolbox 2.1

Highlights from
Finite Element Toolbox 2.1