|
Revision 706, 0.5 kB
(checked in by smidl, 16 years ago)
|
|
eol-native
|
-
Property svn:eol-style set to
native
|
| Line | |
|---|
| 1 | |
|---|
| 2 | itload('../merger_2d_test.it'); |
|---|
| 3 | XG = reshape(Grid(1,:),Npoints,Npoints); |
|---|
| 4 | YG = reshape(Grid(2,:),Npoints,Npoints); |
|---|
| 5 | |
|---|
| 6 | figure(2); |
|---|
| 7 | M1 = reshape(exp(Res1),Npoints,Npoints); |
|---|
| 8 | contour(XG,YG,M1,7); |
|---|
| 9 | |
|---|
| 10 | figure(3); |
|---|
| 11 | hold off |
|---|
| 12 | mm = max(max(Res2)); |
|---|
| 13 | for i=1:size(Res2,1) |
|---|
| 14 | M2 = reshape(Res2(i,:),Npoints,Npoints); |
|---|
| 15 | contour(XG,YG,M2,[0:mm/7:mm]); |
|---|
| 16 | hold on |
|---|
| 17 | end |
|---|
| 18 | S1g=reshape(exp(S1),Npoints,Npoints); |
|---|
| 19 | S2g=reshape(exp(S2),Npoints,Npoints); |
|---|
| 20 | |
|---|
| 21 | contour(XG,YG,M1,[0:mm/7:mm],'m:'); |
|---|
| 22 | contour(XG,YG,S1g,[0:mm/7:mm],'g:'); |
|---|
| 23 | contour(XG,YG,S2g,[0:mm/7:mm],'r:'); |
|---|
