求助，逼近积分

splelvis

u''(x)=f(x),
: w' T0 `$ w( Y" t* n边界boundary conditions u(a)=0,u(b)=0.
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);

3 O+ ~: `1 K+ {) \/ l9 _maltab code:

clear all
a=0;
$ g9 Z8 H  P  \- Y) C( v8 Xb=1;
n=10;
h=(b-a)/(n+1);
x_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));
7 k* i+ {3 p0 r9 K* `rhs=f4(x);
1 f& s, a& O3 h, I, nsol=A\rhs;
sol_with_boundary_conds=[0;sol;0];
4 n5 v7 Q' ^" H' v( a0 [/ jplot(x_with_boundary,sol_with_boundary_conds);
% N9 F, g0 c$ J" f0 H- t% G- `

5 n% l! B( Y% U/ A5 c另开一文档，f4,
2 }, ~8 F% Z& x1 c2 v. J. rfunction y=f4(x)
* N/ s8 Q0 d# u- o2 m% yy=ones(length(x),1);

出来后是抛物线。
. ?# I; s0 j6 w4 X+ _1 f. e; Y
8 z8 ?& b1 O6 w- j( _
现在问题是，
u''(x)=sin(2pix);
4 e2 V/ m+ q, u; h" Ju(-1)=0;
! [# H2 ~* o+ |# J) h0 ?+ Ru(1)=0;
范围[-1;1]

; `  v3 C4 s, e$ {: M3 G  R& T$ w请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

没人能帮忙一下吗，
2 E0 @, W7 V, i  S还是有什么资料不明确？