%% 脚本和实时脚本
% 脚本就是只能写代码的文件，后缀.m；实时脚本就是像这样，既可以写代码也可以写文字、放图片等，后缀.mlx，下图显示了如何创建脚本和实时脚本。接下来我们用实时脚本来讲解。
%
%
% 分节：为了区分每一部分内容的不同，在运行左侧有分节符的按钮，可以分节。
%
%% 简单操作
% ctrl+E：实时脚本当中切换为代码块模式

3

% 它可以是命令窗口那种，直接输入式子就能返回结果

a=3+3

% 当然也可以赋值，这里要引入分号的概念。分号可以让它不输出。

a=3+3;

%
%% 现在开始写代码：实现梯形法积分
% clear：可以清除工作区的所有变量

clear;

% clc：可以清除命令行窗口的所有文本，让屏幕更干净

clc;

% disp：在命令行输出一些东西

disp("Hello World!");

% 用梯形法计算积分区域面积（下面用插入公式获得）
% \int_{3}^{4} x^2 dx = \frac{4^3-3^3}{3} = 37/3

a = 3;
b = 4;
n = 20; % 梯形分割数量
h = (b-a)/n;
s = 0;
for k=1:1:(n-1)
   x= a + h*k;
   s = s + x^2;
end
s=(h*(a^2+b^2)/2+h*s)/(b-a)

%
%% 序列（向量）
% 初始值：步长：终值，本质是向量，这也可以用在 for 循环里

x = 1:1:10
x1 = linspace(0,10,21)

% 为了逐元素相乘，我们需要在符号前加 .

y = x.^2

%
%% 画图
% 计算机画图其实也是要连点成线的，于是刚刚讲的序列就有用了，y=x.^2 就是针对 y=x^2 用 x 序列采样后得到的结果。



%
%% 画子图

subplot(2,1,1);
x = linspace(0,10);
y1 = sin(x);
plot(x,y1)



%
% 假设不是 y=x^2 这个函数呢？于是我们可以新建一个.m文件来存放这个函数，见 fun.m
% 于是我们就可以整合得到一个函数了