# # A total of EIGHT lines below include the "fopen(...)" command. # These lines must be changed by the user # # to reflect the address of the file being opened. In general, the # string between the quotation marks ("...string...") # # must be changed to reflect the file location... > restart: > with(linalg): Warning, the protected names norm and trace have been redefined and unprotected > T := fopen( "C:/Text/FortranPrograms/Symm/Type.txt",READ): > M := fopen( "C:/Text/FortranPrograms/Symm/M.txt",READ): > DIM := fopen( "C:/Text/FortranPrograms/Symm/Dim.txt",READ): > TT := readdata(T,integer): > TTT := TT[1]: > MM3 := readdata(M,integer): > M3 := MM3[1]: > N3 := M3: > Dimen := readdata(DIM,float): > a := Dimen[1]: > b := Dimen[2]: > coe := array(1..M3,1..N3): > c := array(1..M3): > if (TTT < 2) then > Coef := fopen("C:/Text/FortranPrograms/Symm/Coeffs.txt",READ): > c:= readdata(Coef,float): > i := 1: > for m from 1 to M3 do > for n from 1 to N3 do > coe[m,n] := c[i]: > i := i + 1 > end do > end do: > w := add((add(coe[m,n]*sin(m*Pi*x/a)*sin(n*Pi*y/b),n=1..N3)),m = > 1..M3): > plot3d(w,x=0...a,y=0...b,scaling = constrained, axes = normal,style = > contour); > else > Coef := fopen("C:/Text/FortranPrograms/Symm/Eivec1.txt",READ): > c:= readdata(Coef,float): > i := 1: > for m from 1 to M3 do > for n from 1 to N3 do > coe[m,n] := c[i]: > i := i + 1 > end do > end do: > w := add((add(coe[m,n]*sin(m*Pi*x/a)*sin(n*Pi*y/b),n=1..N3)),m > =1..M3): > plot3d(w,x=0...a,y=0...b,scaling = constrained, axes = normal,style = > contour, title = "Mode 1"); > Coef := fopen("C:/Text/FortranPrograms/Symm/Eivec2.txt",READ): > c:= readdata(Coef,float): > i := 1: > for m from 1 to M3 do > for n from 1 to N3 do > coe[m,n] := c[i]: > i := i + 1 > end do > end do: > w := add((add(coe[m,n]*sin(m*Pi*x/a)*sin(n*Pi*y/b),n=1..N3)),m = > 1..M3): > plot3d(w,x=0...a,y=0...b,scaling = constrained, axes = normal,style = > contour, title = "Mode 2"); > Coef := fopen("C:/Text/FortranPrograms/Symm/Eivec3.txt",READ): > c:= readdata(Coef,float): > i := 1: > for m from 1 to M3 do > for n from 1 to N3 do > coe[m,n] := c[i]: > i := i + 1 > end do > end do: > w := add((add(coe[m,n]*sin(m*Pi*x/a)*sin(n*Pi*y/b),n=1..N3)),m = > 1..M3): > plot3d(w,x=0...a,y=0...b,scaling = constrained, axes = normal,style = > contour, title = "Mode 3"); > Coef := fopen("C:/Text/FortranPrograms/Symm/Eivec4.txt",READ): > c:= readdata(Coef,float): > i := 1: > for m from 1 to M3 do > for n from 1 to N3 do > coe[m,n] := c[i]: > i := i + 1 > end do > end do; > w := add((add(coe[m,n]*sin(m*Pi*x/a)*sin(n*Pi*y/b),n=1..N3)),m = > 1..M3): > plot3d(w,x=0...a,y=0...b,scaling = constrained, axes = normal,style = > contour, title = "Mode 4"); > end if; Coef := 3 -7 -7 -8 c := [-0.35861799 10 , 0.12840639 10 , -0.51064912 10 , -8 -8 -9 0.25172233 10 , -0.13549625 10 , 0.83194263 10 , -9 -9 -9 -0.52528415 10 , 0.36609543 10 , -0.25431671 10 , -9 -7 -8 0.19572426 10 , -0.10816962 10 , 0.67420425 10 , -8 -8 -8 -0.41445034 10 , 0.21963991 10 , -0.13715762 10 , -9 -9 -9 0.81596785 10 , -0.56299698 10 , 0.36796971 10 , -9 -9 -8 -0.27522079 10 , 0.18802301 10 , -0.27585452 10 , -8 -8 -8 0.25475522 10 , -0.18003550 10 , 0.10475060 10 , -9 -9 -9 -0.67637734 10 , 0.39755804 10 , -0.27429162 10 , -9 -9 -10 0.16632433 10 , -0.12458244 10 , 0.60680946 10 , -9 -9 -9 -0.59103683 10 , 0.28673375 10 , -0.15879852 10 , -10 -10 -9 -0.32571244 10 , 0.62773439 10 , -0.11171444 10 , -10 -9 -10 0.92920477 10 , -0.10291816 10 , 0.80060659 10 , -10 -10 -9 -0.98184814 10 , 0.51158376 10 , -0.27433578 10 , -9 -9 -9 0.39650519 10 , -0.43791221 10 , 0.39763570 10 , -9 -9 -9 -0.34691736 10 , 0.28365746 10 , -0.23520799 10 , -9 -9 -9 0.18992130 10 , -0.16222594 10 , 0.11492352 10 , -9 -9 -9 -0.34865077 10 , 0.42745671 10 , -0.44549528 10 , -9 -9 -9 0.40232395 10 , -0.34275557 10 , 0.28346342 10 , -9 -9 -9 -0.22424715 10 , 0.18355702 10 , -0.13027586 10 , -10 -9 -9 0.77005013 10 , -0.18737169 10 , 0.24231203 10 , -9 -9 -9 -0.25026875 10 , 0.22546237 10 , -0.18465475 10 , -9 -9 -10 0.14896485 10 , -0.10648506 10 , 0.83720614 10 , -10 -11 -10 -0.33535327 10 , 0.66899923 10 , -0.68591001 10 , -10 -10 -10 0.74391264 10 , -0.72208711 10 , 0.49933557 10 , -10 -10 -11 -0.30287821 10 , 0.11091294 10 , 0.67327003 10 , -10 -10 -10 -0.15852255 10 , 0.42785771 10 , -0.19213384 10 , -10 -10 -10 0.13895564 10 , -0.24045151 10 , 0.34417084 10 , -10 -10 -10 -0.53951371 10 , 0.63596774 10 , -0.74633029 10 , -10 -10 -10 0.76202475 10 , -0.78716325 10 , 0.79759477 10 , -10 -10 -10 -0.32395166 10 , 0.31015673 10 , -0.51626911 10 , -10 -10 -10 0.62722404 10 , -0.78663222 10 , 0.79857995 10 , -10 -10 -10 -0.83501754 10 , 0.73902912 10 , -0.71076492 10 , -10 0.51536955 10 ] i := 1 -7 -7 w := -0.35861799 10 %20 %11 + 0.12840639 10 %20 %10 -8 -8 - 0.51064912 10 %20 %9 + 0.25172233 10 %20 %8 -8 -9 - 0.13549625 10 %20 %7 + 0.83194263 10 %20 %6 -9 -9 - 0.52528415 10 %20 %5 + 0.36609543 10 %20 %4 -9 -9 - 0.25431671 10 %20 %3 + 0.19572426 10 %20 %1 -7 -8 - 0.10816962 10 %19 %11 + 0.67420425 10 %19 %10 -8 -8 - 0.41445034 10 %19 %9 + 0.21963991 10 %19 %8 -8 -9 - 0.13715762 10 %19 %7 + 0.81596785 10 %19 %6 -9 -9 - 0.56299698 10 %19 %5 + 0.36796971 10 %19 %4 -9 -9 - 0.27522079 10 %19 %3 + 0.18802301 10 %19 %1 -8 -8 - 0.27585452 10 %18 %11 + 0.25475522 10 %18 %10 -8 -8 - 0.18003550 10 %18 %9 + 0.10475060 10 %18 %8 -9 -9 - 0.67637734 10 %18 %7 + 0.39755804 10 %18 %6 -9 -9 - 0.27429162 10 %18 %5 + 0.16632433 10 %18 %4 -9 -10 - 0.12458244 10 %18 %3 + 0.60680946 10 %18 %1 -9 -9 - 0.59103683 10 %17 %11 + 0.28673375 10 %17 %10 -9 -10 - 0.15879852 10 %17 %9 - 0.32571244 10 %17 %8 -10 -9 + 0.62773439 10 %17 %7 - 0.11171444 10 %17 %6 -10 -9 + 0.92920477 10 %17 %5 - 0.10291816 10 %17 %4 -10 -10 + 0.80060659 10 %17 %3 - 0.98184814 10 %17 %1 -10 -9 + 0.51158376 10 %16 %11 - 0.27433578 10 %16 %10 -9 -9 + 0.39650519 10 %16 %9 - 0.43791221 10 %16 %8 -9 -9 + 0.39763570 10 %16 %7 - 0.34691736 10 %16 %6 -9 -9 + 0.28365746 10 %16 %5 - 0.23520799 10 %16 %4 -9 -9 + 0.18992130 10 %16 %3 - 0.16222594 10 %16 %1 -9 -9 + 0.11492352 10 %15 %11 - 0.34865077 10 %15 %10 -9 -9 + 0.42745671 10 %15 %9 - 0.44549528 10 %15 %8 -9 -9 + 0.40232395 10 %15 %7 - 0.34275557 10 %15 %6 -9 -9 + 0.28346342 10 %15 %5 - 0.22424715 10 %15 %4 -9 -9 + 0.18355702 10 %15 %3 - 0.13027586 10 %15 %1 -10 -9 + 0.77005013 10 %14 %11 - 0.18737169 10 %14 %10 -9 -9 + 0.24231203 10 %14 %9 - 0.25026875 10 %14 %8 -9 -9 + 0.22546237 10 %14 %7 - 0.18465475 10 %14 %6 -9 -9 + 0.14896485 10 %14 %5 - 0.10648506 10 %14 %4 -10 -10 + 0.83720614 10 %14 %3 - 0.33535327 10 %14 %1 -11 -10 + 0.66899923 10 %13 %11 - 0.68591001 10 %13 %10 -10 -10 + 0.74391264 10 %13 %9 - 0.72208711 10 %13 %8 -10 -10 + 0.49933557 10 %13 %7 - 0.30287821 10 %13 %6 -10 -11 + 0.11091294 10 %13 %5 + 0.67327003 10 %13 %4 -10 -10 - 0.15852255 10 %13 %3 + 0.42785771 10 %13 %1 -10 -10 - 0.19213384 10 %12 %11 + 0.13895564 10 %12 %10 -10 -10 - 0.24045151 10 %12 %9 + 0.34417084 10 %12 %8 -10 -10 - 0.53951371 10 %12 %7 + 0.63596774 10 %12 %6 -10 -10 - 0.74633029 10 %12 %5 + 0.76202475 10 %12 %4 -10 -10 - 0.78716325 10 %12 %3 + 0.79759477 10 %12 %1 -10 -10 - 0.32395166 10 %2 %11 + 0.31015673 10 %2 %10 -10 -10 - 0.51626911 10 %2 %9 + 0.62722404 10 %2 %8 -10 -10 - 0.78663222 10 %2 %7 + 0.79857995 10 %2 %6 -10 -10 - 0.83501754 10 %2 %5 + 0.73902912 10 %2 %4 -10 -10 - 0.71076492 10 %2 %3 + 0.51536955 10 %2 %1 %1 := sin(1.000000000 Pi y) %2 := sin(0.5000000000 Pi x) %3 := sin(0.9000000000 Pi y) %4 := sin(0.8000000000 Pi y) %5 := sin(0.7000000000 Pi y) %6 := sin(0.6000000000 Pi y) %7 := sin(0.5000000000 Pi y) %8 := sin(0.4000000000 Pi y) %9 := sin(0.3000000000 Pi y) %10 := sin(0.2000000000 Pi y) %11 := sin(0.1000000000 Pi y) %12 := sin(0.4500000000 Pi x) %13 := sin(0.4000000000 Pi x) %14 := sin(0.3500000000 Pi x) %15 := sin(0.3000000000 Pi x) %16 := sin(0.2500000000 Pi x) %17 := sin(0.2000000000 Pi x) %18 := sin(0.1500000000 Pi x) %19 := sin(0.1000000000 Pi x) %20 := sin(0.05000000000 Pi x) > > > > >