子程序ex713f 函数程序应另存成一个文件ex713f.m function zprime=ex713f(t,z) global vt vm zprime-=[0;0];%给出t0之前zprime初值 zprime(1)=-vt-vm*z(1)/sqrt(z(1)2+z(2)^2); zprime(2)=-vm*z(2)/sqrt(z(1)^2+z(2)^2): %上面两句可换成一个矩阵语句： zprime=-vt*[1;0]-vm*z/sqrt(Z(1)^2+z(2)^2) 子程序ex713f 函数程序应另存成一个文件ex713f.m function zprime=ex713f(t,z) global vt vm zprime=[0;0]; % 给出t0之前zprime初值 zprime(1)=-vt-vm*z(1)/sqrt(z(1)^2+z(2)^2); zprime(2)=-vm*z(2)/sqrt(z(1)^2+z(2)^2); %上面两句可换成一个矩阵语句： zprime=-vt*[1;0]-vm*z/sqrt(z(1)^2+z(2)^2);