求解第一类边值问题 源：matlab

matlab爱好者。
求解第一类边值问题的四种方法（完整），包括相互的验证：
任何第一类边值问题都可以用这四种方法快速解决。矩阵法这里只能适用于5个离散点，想得到通解请自己推一遍矩阵法。
我们以一道例题开始：
题目：
解法：①初值解法解边值：通过给dy来求y（end），直到拟合。首先做一个函数，给一个dy,返回一个y（end）。
函数名为findf，
如下:
function s=findf(dy)
t=linspace(0,4,11);
dt=t(2)-t(1);
y(1)=0; v(1)=dy;
for i=1:length(t)-1
y(i+1)=y(i)+v(i)dt;
v(i+1)=v(i)+(2v(i)cos(t(i))-y(i)sin(4t(i))-cos(3t(i)))dt;
end
s=y(end);
程序源代码：
clear
x=linspace(0,2,10);
for i=1:length(x)
ss(i)=findf(x(i));
end
plot(x,ss)
hold on;
o=linspace(0,2,10);
for j=1:length(o)
g(j)=2;
end
plot(o,g,’-’);