请问关于解微分方程组的问题

fmx2000

用MATLAB求解该方程组时，我编写了m文件和工作窗口调用命令：
%创建状态方程
function dy=generator(t,y);
%定义常量
Ws=377;
Tj=6.4;
D=2;
Pm=0.7516;
dy==zeros(2,1); %初始化
dy(1)=[Pm-Pe-D*y(1)]/Tj;%pe是变量
dy(2)=Ws*y(1);

%窗口调用命令
%计算输出
y0 = [0; 19.25];           % 取初值
options = optimset('Display','off');
t=0:0.001:2;             % 变量取值范围[0 2]
for i=1:1:723            %仿真输出的pe是大小为723的Array存储在workplace里
pe=pe(:,2);                        
x = fsolve(@generator, y0, options,pe);%求解非线性方程组
y1(i)=y(1);
y2(i)=y(2);
end
plot(k,y1,'-b',k,y2,'-r');
xlabel('t')
legend('y1','y2')
simulink仿真时采取的是变步长ode23tb法，仿真得到的pe是大小为723的Array。请问各位朋友：在用fsolve求解程序中应该如何调用Pe的值，以及如何控制t的步长，能否使其做到与仿真的变步长一样来求解？用上面程序出现维数不匹配问题。谢谢指导。