[饿得快]

matlab程序已经写出，但是在与vb混合编程的时候卡住了，搞得焦头烂额，主要是想实现在vb界面上输入参数和显示最后画图，在此请教各位了，如能赐教，不胜感激，matlab程序如下

k=1.0,ds=0.6,tc=0.5,td=0.1,tm=0.24,v=1,syms t,lao=0.3,hao=0.2,
wf=0.09;
wrl=0.23;
ltr=0.32;
lth=0.38;
lsh=0.35;
lan=0.07;
laf=0.2;
lab=0.05;
qb=pi/8;
qf=pi/8;
a=[1.0, 0.0, 0.0, 0.0; 1.0, td, td^2, td^3; 1.0, tm, tm^2, tm^3; 1.0, tc, tc^2, tc^3]
c1=[lan; laf*sin(qb)+lan*cos(qb); hao; lab*sin(qf)+lan*cos(qf)]
b1=pinv(a)*c1
y1=b1(1)+b1(2)*(t-k*tc)+b1(3)*((t-k*tc).^2)+b1(4)*((t-k*tc).^3)


T1 = 0.5:0.01:1;
Result = zeros(1,length(T1));%保存结果的数组(这里只采用二维数组)
%假如你想查看qm(n1,n2)随时间变化的情况

for k = 1:length(T1)
tmp = subs(y1,t,T1(k));%求qm的真实值
Result(k) = tmp;%保存结果
end
plot(T1,Result);

我试过用matrixVB做，但是做到y1的表达式就做不下去了，程序如下
td=0.1
tm=0.24
tc=0.5
lan=0.07
laf=0.2
lab=0.05
qb=3.14/4
hao=0.2
qf=3.14/4
a=creatematrix(1.0, 1.0, 1.0, 1.0, 0, td, tm, tc, 0, td^2, tm^2, tc^2, 0, td^3, tm^3, tc^3)
a=reshape(a,4,4)
c1=creatematrix(lan, laf*sin(qb)+lan*cos(qb), hao, lab*sin(qf)+lan*cos(qf))
c1=reshape(c1,4,1)
c1.show

想试用com编译，但是出很多问题，matlab2009没comtool.matlab6.5又老是出错，我要疯掉了

留个qq11375717,希望懂得的朋友不吝赐教啊~~
+急