[MATLAB毕业设计] 一个函数拟合的问题

[ning19890820]

T D K
13.39397 10779.24086 26376.57565
13.18762 10624.98982 28662.55399
12.12091 9529.590057 28863.00213
12.29292 9550.568789 29101.70489
12.77686 9570.945701 30343.15172
12.96479 9503.750018 30721.88902
24.56432 16004.43478 34224.24505
23.86597 16041.78065 34676.95219
22.79274 16801.20476 35004.60145
17.15393 10327.30357 38834.72723
17.19834 10363.02472 39477.20382
23.30241 17202.82461 39725.03034
19.90063 10691.38352 39808.59269
17.66709 17171.67696 40183.9513
14.30726 5943.092718 40368.17276
16.37882 10667.03923 40382.59636
14.26416 11686.96403 40423.08912
17.78141 16943.18455 40532.53862
19.95826 10672.33229 40616.43613
13.65124 10821.18729 40995.32794

上面给出了K T D的数据 要拟合如图这样的函数：
主要是拟合出a c f
还要画出K T D的三维图
求程序……虽然很想自己研究 但是时间不允许了大学四年也没学过 谢谢各位
下面是别人写的 可是我运行不出来 而且也没有画三维图的步骤 能不能稍微改一下呢
data=[13.39397 10779.24086 26376.57565
13.18762 10624.98982 28662.55399
12.12091 9529.590057 28863.00213
12.29292 9550.568789 29101.70489
12.77686 9570.945701 30343.15172
12.96479 9503.750018 30721.88902
24.56432 16004.43478 34224.24505
23.86597 16041.78065 34676.95219
22.79274 16801.20476 35004.60145
17.15393 10327.30357 38834.72723
17.19834 10363.02472 39477.20382
23.30241 17202.82461 39725.03034
19.90063 10691.38352 39808.59269];
T=data(:,1);
D=data(:,2);
x=[T D];
K=data(:,3);
fun=@(a,x)a(1)*x(:,1).^a(2).*x(:,2).^a(3);
a0=[1 0.5 2];
beta = nlinfit(x,K,fun,a0)