求助，逼近积分

splelvis

u''(x)=f(x),
边界boundary conditions u(a)=0,u(b)=0.
/ e/ G# k' i3 v( D% Q8 j(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
" p! R  \! Y9 T, Q, C+ S+ j
  c) ?& [! N* Q( w/ u: O; ]maltab code:

# G: H  g3 ?8 G$ jclear all
a=0;
9 W. i! u- s( a8 Bb=1;
n=10;
h=(b-a)/(n+1);
x_with_boundary=linspace(a,b,n+2)';
: k- B- _/ M& g9 O# R) I( Sx=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);
6 D4 [1 i/ s6 p5 b. Osol=A\rhs;
sol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);
$ Y. Z7 }" L* N$ A% l

另开一文档，f4,
function y=f4(x)
y=ones(length(x),1);

" o) {% z4 z' {1 z7 M. m* Q出来后是抛物线。


- g1 w0 d2 f( @9 w7 q) W8 Q现在问题是，
8 I5 u: C! ]- g/ c$ Nu''(x)=sin(2pix);
& u, z! |9 D+ C" vu(-1)=0;
( G4 L6 ?9 a1 x9 lu(1)=0;
范围[-1;1]
$ b( S- k3 ^( K' }+ d
. `2 j: s) g7 J0 L6 L9 K! L1 y请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

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