timour_matlab/mscripts/forcing.m

Heat=1.0;
Nx=10;
Ny=20;
for i=1:Nx
for j=1:Ny
heating(i,j)=-Heat*(sin(pi*(2*j/Ny)));;
end
end
Stress=1.0;
Nx=10;
Ny=20;
for i=1:Nx
for j=1:Ny
stress(i,j)=-Stress*(cos(pi*(2*j/Ny)));
end
end
for i=1:Nx
for j=1:Ny
stressy(i,j)=-Stress*(cos(pi*(2*j/Ny)));
end
end
for i=1:Nx
for j=1:Ny
stressx(i,j)=-Stress*(cos(pi*(2*i/Nx)));;
end
end
stress=stressx+stressy;
1	Heat=1.0;
2	Nx=10;
3	Ny=20;
4	for i=1:Nx
5	for j=1:Ny
6	heating(i,j)=-Heat(sin(pi(2*j/Ny)));;
7	end
8	end
9	Stress=1.0;
10	Nx=10;
11	Ny=20;
12	for i=1:Nx
13	for j=1:Ny
14	stress(i,j)=-Stress(cos(pi(2*j/Ny)));
15	end
16	end
17	for i=1:Nx
18	for j=1:Ny
19	stressy(i,j)=-Stress(cos(pi(2*j/Ny)));
20	end
21	end
22	for i=1:Nx
23	for j=1:Ny
24	stressx(i,j)=-Stress(cos(pi(2*i/Nx)));;
25	end
26	end
27	stress=stressx+stressy;