{VERSION 5 0 "IBM INTEL LINUX" "5.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 14 0 0 0 0 2 2 2 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "R 3 Font 0" -1 256 1 {CSTYLE "" -1 -1 "Helvetica" 1 14 0 0 0 0 2 1 2 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "R3 Font 2" -1 257 1 {CSTYLE "" -1 -1 "Courier" 1 14 0 0 0 0 2 2 2 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(plots);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 62 "cylinderplot([r,theta,r^4*sin(2*theta)],theta=0..2* Pi,r=0..1);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "g:=(theta)-> sin(5*theta);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "MAX :=20;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "a:=array(0..MAX); \+ b:=array(1..MAX);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "a[0]: =evalf(int(g(theta),theta=0..2*Pi)/(2*Pi));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "for n from 1 by 1 to MAX" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 66 " a[ n]:=evalf(int(cos(n*theta)*g(theta) ,theta=0..2*Pi)/Pi);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " od;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "unassign('n');" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "for n from 1 by 1 to MAX" } }{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 66 " b[n]:=evalf(int(sin(n*theta)*g(theta) ,theta= 0..2*Pi)/Pi);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " od;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "unassign('n');" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 93 "u:=(r,theta)->a[0]+sum(a[n]*r^n*cos(n*theta) ,n=1..MAX)+sum(b[n]*r^n*sin(n*theta) ,n=1..MAX);" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 56 "cylinderplot([r,theta,u(r,theta)],theta=0..2 *Pi,r=0..1);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "g:=(theta)- >sin(5*theta);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 8 "MAX:=20;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "a:=array(0..MAX); b:=array(1..MAX);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "a[0]:=evalf(int(g(theta),theta=0..Pi/5)/(2*Pi)); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "for n from 1 by 1 to MA X" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 66 " a[n]:=evalf(int(cos(n*theta)*g(theta) ,theta= 0..Pi/5)/Pi);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " od;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "unassign('n');" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "for n from 1 by 1 to MAX" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " d o" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 66 " b[n]:=evalf(int(sin(n* theta)*g(theta) ,theta=0..Pi/5)/Pi);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " od;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "unassign('n');" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 93 "u:=(r,theta)->a[0]+sum(a[n]*r^n*cos(n*theta) ,n=1..MAX)+sum(b[n]*r^n*sin(n*theta) ,n=1..MAX);" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 56 "cylinderplot([r,theta,u(r,theta)],theta=0..2 *Pi,r=0..1);" }}}}{MARK "27" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 } {PAGENUMBERS 0 1 2 33 1 1 }