我仍在学习Scilab(5.5.2),因此我正在编写并运行测试代码以熟悉本软件。 要测试的数值微分方程解算器,我开始容易从方程DY/DX = A,其具有作为溶液Y =斧+ C(线方程)。 这是我写的代码: // Function y = A*x+1
function ydot=fn(x, A)
ydot=A
endfunction
A=2;
//Initial condition
我有SCILAB 我如何可以绘制包含函数的一点改进绘图函数,如果和<像 function y = alpha(t)
if (t < 227.8) then
y = 0.75;
elseif (t < 300) then
y = 2.8 - 0.009 .* t;
else
y = 0.1;
end
endfunction
我有这样的循环: Nbcomp=floor(long/(2*dpas)+0.5);
xn=input("Voulez-vous modifier les éléments de compensations si oui tapé 1 sinon 0")
if xn==1 then
for ca=0:Nbcomp
pas_pr1=(long-2*dpas)/(Nbcomp-1