% 创建一个交互窗口 disp('欢迎使用交互窗口'); disp('------------------');

while true % 显示选项 disp('请选择操作:'); disp('1. 微分方程数值解'); disp('2. 微分方程解析解'); disp('3. 退出'); %输入方程种类

choice = input('输入你的选择: ');
% 根据用户选择执行相应的操作
if choice == 1
    ode1 = input('请输入微分方程:','s');%输入微分方程
    x0 = input('请输入x初始条件:');
    y0 = input('请输入y初始条件:');%初始条件
   
    [x, y] = ode45(ode1, [x0,x0+1], y0);%利用ode45函数求解
    
    disp('计算结果:');
    disp([x, y]);%输出结果
    plot(x, y, 'r-', 'LineWidth', 2); % 数值解
end
    
if choice==2
    ode2=input('请输入微分方程','s');%输入方程
    syms y(x);
    m=dsolve(ode2);%用dsolve函数求解析解
    disp(m);%打印出解析解
    
    


end

end


原文地址: https://www.cveoy.top/t/topic/ptwT 著作权归作者所有。请勿转载和采集!

免费AI点我,无需注册和登录