b1.m

来自「1、 写出此两点边值问题的积分形式； 2、 构造有限元方程； 3、 编程计算」· M 代码 · 共 8 行

%求解有限元方程右端的常数向量的子函数，b1.m
function b=b1(x1,h1)
syms x h kesi
f=h*(-4*(x+h*kesi)*exp(x+h*kesi)-(x+h*kesi)*(exp(1)-exp(-1)))*(1-kesi);
f=int(f,kesi,0,1);
h=h1;
x=x1;
b=eval(f);