求助，逼近积分

splelvis

u''(x)=f(x),
( i, {2 c& P! r) }2 Y- T$ k边界boundary conditions u(a)=0,u(b)=0.
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
; Y; D6 a) q+ ?1 b
maltab code:
4 b9 H) S4 s( {) ?# H
/ Z9 y( q) H5 ^, k; q& I* Yclear all
* n! v, v9 u  Y9 Da=0;
) M3 G  _$ @2 ?b=1;
, i0 J; O. V& D1 A2 j3 In=10;
; h$ F: Q! f% S* O; q7 v5 C. }. R) Qh=(b-a)/(n+1);
6 z( j2 e  Z' G$ ?2 k+ gx_with_boundary=linspace(a,b,n+2)';
x=x_with_boundary(2:n+1);
5 H% r7 r' E  U8 k1 o# }: D* wA=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
. a$ ]3 d/ f+ q6 X% |: [' u- hrhs=f4(x);
sol=A\rhs;
! k2 A% r) h; A( J; k7 T3 M8 Gsol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);
1 d' Z! i6 s4 c* a$ J# d$ g/ s

另开一文档，f4,
/ @6 V1 X) o) H9 \function y=f4(x)
y=ones(length(x),1);
4 U1 b# T" {5 z0 _
出来后是抛物线。
- o2 q8 o6 V/ i1 X

- O4 S$ F0 L" o/ U" _# E! L; k现在问题是，
) @( ~6 ~5 P# x5 G' U) fu''(x)=sin(2pix);
u(-1)=0;
u(1)=0;
范围[-1;1]

' W/ R) g4 l2 I0 x9 i请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

没人能帮忙一下吗，
8 z# A6 ]: o2 t* U! G还是有什么资料不明确？