-
Notifications
You must be signed in to change notification settings - Fork 0
/
surfacegraph.m
70 lines (66 loc) · 2.42 KB
/
surfacegraph.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
% creates a titled, single surface plot of each state variable
% then saves the result according to t = T
hf = figure;
set(hf,'position',[10 190 640 500]);
filenameS = strcat('Sr', num2str(paramcond), '_', strrep(num2str(t), '.', '_'), '.eps');
surf(XX,YY,S','FaceColor','interp',...
'EdgeColor','k',...
'FaceLighting','phong');
title(strcat('Surface Plot of S at T = ', num2str(t)));%,'Interpreter','LaTex');
xlabel('x');%,'Interpreter','LaTex');
ylabel('y');%,'Interpreter','LaTex');
zlabel(strcat('S(x, y, ', num2str(t), ')'));%,'Interpreter','LaTex');
%daspect([5 5 0.1]);
axis([0, 1, 0, 1, 0, tmaxS])
view(-50,30);
camlight left;
saveas(gcf,filenameS, 'psc2')
closereq
hf = figure;
set(hf,'position',[10 190 640 500]);
filenameX = strcat('Xr', num2str(paramcond), '_', strrep(num2str(t), '.', '_'), '.eps');
surf(XX,YY,X','FaceColor','interp',...
'EdgeColor','k',...
'FaceLighting','phong');
title(strcat('Surface Plot of X at T = ', num2str(t)));%,'Interpreter','LaTex');
xlabel('x');%,'Interpreter','LaTex');
ylabel('y');%,'Interpreter','LaTex');
zlabel(strcat('X(x, y, ', num2str(t), ')'));%,'Interpreter','LaTex');
%daspect([5 5 0.1]);
axis([0,1, 0, 1, 0, tmaxX])
view(-50,30);
camlight left;
saveas(gcf,filenameX, 'psc2')
closereq
hf = figure;
set(hf,'position',[10 190 640 500]);
filenameI = strcat('Ir', num2str(paramcond), '_', strrep(num2str(t), '.', '_'), '.eps');
surf(XX,YY,I','FaceColor','interp',...
'EdgeColor','k',...
'FaceLighting','phong');
title(strcat('Surface Plot of I at T = ', num2str(t)));%,'Interpreter','LaTex');
xlabel('x');%,'Interpreter','LaTex');
ylabel('y');%,'Interpreter','LaTex');
zlabel(strcat('I(x, y, ', num2str(t), ')'));%,'Interpreter','LaTex');
%daspect([5 5 0.1]);
axis([0,1, 0, 1, 0, tmaxI])
view(-50,30);
camlight left;
saveas(gcf,filenameI, 'psc2')
closereq
hf = figure;
set(hf,'position',[10 190 640 500]);
filenameE = strcat('Er', num2str(paramcond), '_', strrep(num2str(t), '.', '_'), '.eps');
surf(XX,YY,E','FaceColor','interp',...
'EdgeColor','k',...
'FaceLighting','phong');
title(strcat('Surface Plot of E at T = ', num2str(t)));%,'Interpreter','LaTex');
xlabel('x');%,'Interpreter','LaTex');
ylabel('y');%,'Interpreter','LaTex');
zlabel(strcat('E(x, y, ', num2str(t), ')'));%,'Interpreter','LaTex');
%daspect([5 5 0.1]);
axis([0,1, 0, 1, 0, tmaxE])
view(-50,30);
camlight left;
saveas(gcf,filenameE, 'psc2')
closereq