matlab练习（11.7）

clear;clc
a=input('请输入a的值：');
b=input('请输入b的值：');
c=input('请输入c的值：');
x=input('请输入x的值：');
if x>=0.5&x<1.5
 y=a*x^2+b*x+c
elseif x<3.5&x>=1.5
 y=a*sin(b)^c+x
elseif x>=3.5&x<5.5
 y=log(abs(b+(c./x)))
else
 y='ERROR!'
end

clc;clear;
A=[1 -1 2 3;5 1 -4 2;3 0 5 2;11 15 0 9]
diag(A)%主对角元素
triu(A)%上三角
tril(A)%下三角
rank(A)%秩
trace(A)%迹
norm(A)%范数
cond(A)%条件

%函数文件
function dy = rigid(t,y)
dy=zeros(3,1);
dy(1)=y(2)*y(3);
dy(2)=-y(1)*y(3);
dy(3)=-0.51*y(1)*y(2);
end
%命令：
clc;clear
[T,Y]=ode45('rigid',[0,20],[0,1,1])
plot(T,Y(:,1),'-',T,Y(:,2),'+',T,Y(:,3),'*');

 (9条消息) Matlab学习——求解微分方程(组)_Nikki_o3o的博客-CSDN博客

%%
%1. 求解微分方程 y ' + 2xy = xe-x2
tic
clear;clc
syms x y;
y=dsolve('Dy+2*x*y=x*exp(-x^2)','x')
toc

%%
%2. 求微分方程 xy ' + y - e x  = 0 在初始条件 y (1) = 2e 下的特解并画出解函数的图形.
syms x y;
y=dsolve('x*Dy+y-exp(x)=0','y(1)=2*exp(1)','x')
ezplot(y)