求助，逼近积分

splelvis

u''(x)=f(x),
  J- W, z' k3 z# n+ x) t% M0 q. v边界boundary conditions u(a)=0,u(b)=0.
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);

3 ?. s- t7 f2 g( ~maltab code:
: q: N$ b, |  ^( @2 |
0 ^! ?2 l% }7 {+ {. s5 vclear all
a=0;
b=1;
n=10;
h=(b-a)/(n+1);
3 I5 p) ^3 G4 ]. ]x_with_boundary=linspace(a,b,n+2)';
x=x_with_boundary(2:n+1);
' H' c! F0 v1 L- _8 g% qA=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
rhs=f4(x);
  s; L9 ~9 Q+ x. y5 m9 K; ysol=A\rhs;
3 P- n" f# e, }: Q2 k  s+ osol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);

0 W( L$ \5 ?$ E/ \% p$ w6 I
另开一文档，f4,
& w# b) T; _+ afunction y=f4(x)
y=ones(length(x),1);

' Y* p! e/ k7 g! k# \; {出来后是抛物线。
0 `* U+ G4 }" t! e) w. Q

6 P' a( n1 G9 }- }5 o" s! \现在问题是，
1 H0 N( Y& g7 _2 e  iu''(x)=sin(2pix);
u(-1)=0;
u(1)=0;
) o4 W6 G# Z2 X) i范围[-1;1]
" k9 O# H# f9 V9 x+ G0 j
请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

没人能帮忙一下吗，
/ q; A/ k3 `+ U7 l8 d还是有什么资料不明确？