# \donttest{
data(triMesh)
data(fea_result)
meshP = triMesh$MeshPts$p
meshT = triMesh$MeshPts$T
PlotVal = abs(fea_result$Stress[,1])
Oc = "slateblue"; ac = "steelblue2"; bc = "cyan2"; cc = "palegreen2";
dc = "darkolivegreen1"; ec = "lemonchiffon"; fc = "lightgoldenrod1"; gc = "gold";
hc= "lightsalmon"; ic= "tomato"; jc= "firebrick3"
a = 1e5; b = 5e5; c = 1e6; d = 5e6; e = 1e7; f = 5e7; g = 1e8; h = 5e8; i = 1e9; j =5e9
PlotSystem.2d(meshP, meshT, PlotVal, a, b, c, d, e, f, g, h, i, j,
Oc, ac, bc, cc, dc, ec, fc, gc, hc, ic, jc)
# }
Run the code above in your browser using DataLab