求助，逼近积分

splelvis

u''(x)=f(x),
边界boundary conditions u(a)=0,u(b)=0.
$ Y/ V' ^, b# e  Z1 @' [- p: `(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
3 l2 @% F) S8 |$ f1 g" b
( v9 D$ b. m7 l8 fmaltab code:

clear all
  v+ ?, v. \, a0 r/ ga=0;
3 k& s9 m% e, ub=1;
n=10;
7 l. z7 {# V% @! t+ R1 n) Hh=(b-a)/(n+1);
1 r% b8 z0 @/ s9 K7 C* mx_with_boundary=linspace(a,b,n+2)';
x=x_with_boundary(2:n+1);
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
rhs=f4(x);
" x3 x2 ~+ Z. X2 D( esol=A\rhs;
sol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);
' q8 G  B- Y7 }- ?# `
$ ~5 s2 P" t* b5 j6 X1 [
7 A1 E4 M& L! U$ r# u另开一文档，f4,
; u8 P# K; m. S- w: mfunction y=f4(x)
# W8 l0 e) s8 P" ]y=ones(length(x),1);

; e  a7 O8 L# |$ ^5 |- Z5 |出来后是抛物线。

. w( ]- f. e) c) L8 K: q% @0 K# E+ C0 w
& U0 n3 M+ o- Z  }8 L现在问题是，
u''(x)=sin(2pix);
6 Y% F" x: v" Z1 y* t# U# _u(-1)=0;
9 z0 i, ^" I% w$ ?u(1)=0;
范围[-1;1]
$ u! M3 _5 j* {! ~
请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

没人能帮忙一下吗，
还是有什么资料不明确？