math2930_2011fa_week14&15_PS13Sol_Zehnder

math2930_2011fa_week14&15_PS13Sol_Zehnder - 'r' )...

Info iconThis preview shows pages 1–7. Sign up to view the full content.

View Full Document Right Arrow Icon
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 2
%HW 10.8 #6 a=2; nmax=5; [R,T]=meshgrid(0:a/100:a,0:pi/100:pi); X=R.*cos(T); Y=R.*sin(T); U=0*X; for n=1:nmax; kn=4*(1-cos(n*pi))/(a^n*pi*n^3); U=U+kn*R.^n.*sin(n*T); end figure(1) h=gca; set(h, 'FontSize' ,14, 'FontWeight' , 'bold' ); surf(X,Y,U) xlabel( 'x' ) ylabel( 'y' ) zlabel( 'u(r,\theta)' ) shading interp figure(2) h=gca; set(h, 'FontSize' ,14, 'FontWeight' , 'bold' ); plot(R(1,:),U(10,:), R(1,:), U(20,:), R(1,:), U(30,:), R(1,:), U(40,:), R(1,:), U(50,:)) xlabel(
Background image of page 3

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 4
Background image of page 5

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 6
Background image of page 7
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: 'r' ) ylabel( 'u(r) for \theta=[.1,.2,.3,.4,.5]*\pi' ) figure(3) h=gca; set(h, 'FontSize' ,14, 'FontWeight' , 'bold' ); plot(T(:,1),U(:,20),T(:,1),U(:,40),T(:,1),U (:,60),T(:,1),U(:,80),T(:,1),U(:,101) ) xlabel( '\theta' ) ylabel( 'u(\theta) for r=[.2,.4,.6,.8,1.]*a' ) 1 2 3 4 0.5 1 1.5 2 2.5 u( ) for r=[.2,.4,.6,.8,1.]*a 0.5 1 1.5 2 0.5 1 1.5 2 2.5 r u(r) for =[.1,.2,.3,.4,.5]*...
View Full Document

Page1 / 7

math2930_2011fa_week14&15_PS13Sol_Zehnder - 'r' )...

This preview shows document pages 1 - 7. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online