小编urs*_*oth的帖子

我想翻转我的矩阵。使得 T(1,1)=C(2,1)

我已经制作了这个程序，我发现了一个应该在线完成的代码，C=T(2:1:-1, :)但是当尝试获得应该是 3 的值 C(1,1) 时，我得到 1.3533635457363350E-306。你如何翻转矩阵使得向上变成向下？

program main


implicit none
  integer iMax, jMax
  double precision, dimension(:,:), allocatable :: T,C

double precision x, dx,f,L2old,L2norm

integer i, j,n


 allocate(T(0:2, 0:2))
 allocate(C(0:2, 0:2))


T(1,1)=1
T(1,2)=2
T(2,1)=3
T(2,2)=4

write(*,*) T(2,2)

C=T(2:1:-1, :)

Write(*,*) C(1,2)


end program main

我在 matlab 中做了一个等高线图（参见代码）。我想找到值等于 1 的等高线。现在我刚刚找到它大约在等高线图之间： 这可以做到吗？例如，如果我想从值 0 到 1 绘制 5 条等高线

更新我设法绘制等于 1 的轮廓线，但我希望轮廓线在内部，而不是在轮廓线 =1 外部，正如我使用此代码得到的那样。

[x,y] = meshgrid(-3 : 0.01: 3, -3 : 0.01: 3);
s = x + i*y;
z=abs(1+s+((s.^2)/2)+((s.^3)/6));
figure;


[C,h] = contour(x,y,z,[1 1]);
clabel(C,h)
hold on;
[R,k] = contour(x,y,z,25);
clabel(R,k)