魔方阵matlab

clear all
clc
M=magic(4)
M([1,5,9])=[NaN NaN NaN]
y=nanstd(M)

先生成魔方阵，然后替换1，5，9位置，然后忽略，之后再求忽略nan之后的标准差。

M =

    16     2     3    13
     5    11    10     8
     9     7     6    12
     4    14    15     1


M =

   NaN   NaN   NaN    13
     5    11    10     8
     9     7     6    12
     4    14    15     1


y =

    2.6458    3.5119    4.5092    5.4467