求助，逼近积分

splelvis

u''(x)=f(x),
1 x9 Y" ~5 k8 I  X边界boundary conditions u(a)=0,u(b)=0.
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
0 X8 v& ]: B) f, [; w0 c8 p4 F
$ }; w7 N8 J7 e3 e; ^& t  Amaltab code:

+ |, f3 L* G1 nclear all
! c2 u: J" E8 E/ ja=0;
b=1;
n=10;
9 C- J% ]% j# t8 `: Mh=(b-a)/(n+1);
2 S4 X. e4 H: ~- B; N" s, qx_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));
rhs=f4(x);
1 y0 n5 V0 g& @) Xsol=A\rhs;
sol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);


另开一文档，f4,
function y=f4(x)
y=ones(length(x),1);
6 _5 C9 @# r( \3 U2 j7 Q! [" g1 ]: _
* w+ k! e6 [! y3 d5 ~% c出来后是抛物线。


" g- o& F- `) a9 r4 C( l现在问题是，
u''(x)=sin(2pix);
: O, R; Q- u, w. t* L8 Xu(-1)=0;
3 J: c2 z7 a8 g7 ^u(1)=0;
范围[-1;1]
4 W9 s- u$ `9 h& O/ h4 ]
请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

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