Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: surface area of pdedemo3
Replies: 5   Last Post: Jul 2, 2013 4:16 PM

 Messages: [ Previous | Next ]
 Steven Finch Posts: 9 Registered: 7/17/08
Re: surface area of pdedemo3
Posted: Jul 2, 2013 4:16 PM

Alan_Weiss <aweiss@mathworks.com> wrote in message <kqcjfk\$ffg\$1@newscl01ah.mathworks.com>...
>
> I imagine that it is easier and less error-prone to use Bill Greene's
> suggestion rather than my special-purpose triangle summation function.

Running the following code in Matlab R2013a:

format long
g='circleg';
b='circleb2';
c='1./sqrt(1+ux.^2+uy.^2)';
a=0;
f=0;
rtol=1e-6;
[p,e,t]=initmesh(g);
[p,e,t]=refinemesh(g,p,e,t);
[p,e,t]=refinemesh(g,p,e,t);
[p,e,t]=refinemesh(g,p,e,t);
[p,e,t]=refinemesh(g,p,e,t);
[p,e,t]=refinemesh(g,p,e,t);
[p,e,t]=refinemesh(g,p,e,t);
[p,e,t]=refinemesh(g,p,e,t);
u=pdenonlin(b,p,e,t,c,a,f,'tol',rtol);
tArea = pdetrg(p,t);
surfaceArea = sqrt(1 + ux.^2 + uy.^2)*tArea';
surfaceArea

gives 3.82697... as an estimate for the surface area
(all five digits of which agree with Ken Brakke's
"Surface Evolver" software). I am thankful to
Alan Weiss and Bill Greene for their generous help.

Steve Finch

Date Subject Author
6/24/13 Steven Finch
6/25/13 Alan Weiss
6/25/13 Steven Finch
6/25/13 Bill Greene
6/25/13 Alan Weiss
7/2/13 Steven Finch