求助，逼近积分

splelvis

u''(x)=f(x),
边界boundary conditions u(a)=0,u(b)=0.
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
( N6 \0 q% @/ g# E6 g7 P) k
maltab code:

, u8 w* x7 `* K) A8 J! U. oclear all
- u+ w" \4 y4 v9 Y% b' e, |1 Ya=0;
% \' s) V. ^0 V9 Eb=1;
) ~* w1 J: J0 }( Gn=10;
2 p. _0 S! I4 P8 y0 c. uh=(b-a)/(n+1);
x_with_boundary=linspace(a,b,n+2)';
- B! d; {. B3 @' O/ d- _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);
- k7 A6 I& @0 f6 S  ^sol=A\rhs;
sol_with_boundary_conds=[0;sol;0];
9 B/ h; o1 F( q: x6 w- y: aplot(x_with_boundary,sol_with_boundary_conds);
& ]6 \% g. U& q) T
( g. U$ z1 p5 b# I& q, q" w
另开一文档，f4,
function y=f4(x)
& H" r2 k6 c) r3 Qy=ones(length(x),1);

: r: o' B4 a( p; G& H* i2 x! u$ y& Z出来后是抛物线。
8 A& I9 P- o) l' l  K, I

2 T% u. b/ l  T. O- |+ ~现在问题是，
u''(x)=sin(2pix);
! ?0 y2 V* v2 _u(-1)=0;
1 N. M' i, g$ D+ ju(1)=0;
7 M2 O+ t; Q" z2 I2 }. @范围[-1;1]
) B+ R+ [0 W+ G: N
* d1 j9 J9 F' D# o4 J: @  L  y请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

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