求解微分方程ⅆy/ⅆx+ⅇ^(-x)y=x²（y(0)=0）Scilab代码报错求助

错误原因 #

Scilab的ode函数对微分方程右端的函数有明确要求：必须接收两个输入参数（当前自变量x和当前因变量y）。你写的function ydot=f(x)只定义了一个参数x，同时函数内部还引用了未传入的变量y，直接导致参数数量不匹配，触发了Wrong number of input arguments错误。

另外还有两个小问题：

• title函数用法错误：Scilab中title仅用于设置图标题，x/y轴标签需要用xlabel和ylabel单独设置
• 显示解的代码无效：sol='[x;y]'只是定义了一个字符串，并没有实际存储x和y的数值

修正后的完整代码 #

clear
clc
clf

// 修正：函数必须接收x和y两个参数
function ydot=f(x,y)
    ydot = x^2 - exp(-x)*y;
endfunction

x0 = 0; // 初始x值
y0 = 0; // 初始y值
x = 0:.1:5; // x范围
y = ode(y0, x0, x, f);

plot(x, y, '.-')
title('Plot of Y vs X')
xlabel('X axis') // 单独设置x轴标签
ylabel('Y axis') // 单独设置y轴标签

// 正确显示解：直接输出x和y的数值矩阵
disp('The solution is:')
disp([x; y])

关键修正点说明 #

• 函数定义改为function ydot=f(x,y)：满足ode函数对输入参数的要求，确保y能被正确传入计算
• 拆分标题与轴标签：用xlabel和ylabel分别设置轴标签，符合Scilab的语法规范
• 修正解的显示逻辑：用[x; y]生成数值矩阵再输出，能实际看到每个x对应的y值

内容的提问来源于stack exchange，提问作者Rakesh