求助，逼近积分

splelvis

u''(x)=f(x),
$ I! ?8 v( J* S( S! U边界boundary conditions u(a)=0,u(b)=0.
9 r1 p" g. ~4 e(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
( B6 I5 ]8 D+ ^
maltab code:
! c5 F( r& ~, V7 @9 X
- o' u3 G$ O9 e5 ~% W. rclear all
# M3 l5 m" z* r: H1 D. }a=0;
- k! T) b+ d1 q4 h7 jb=1;
n=10;
h=(b-a)/(n+1);
x_with_boundary=linspace(a,b,n+2)';
/ j9 {" G- }% i0 b1 S3 A) Vx=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));
3 Z: k$ a6 y& a' I- R) X0 Prhs=f4(x);
sol=A\rhs;
sol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);

0 D/ g- e4 K6 C" f. L' S& f
% R  V' S- y5 U9 i6 L另开一文档，f4,
1 c$ B$ {2 ~- W( ~3 |function y=f4(x)
y=ones(length(x),1);

! o' O2 s( m7 [4 O% I- M出来后是抛物线。


现在问题是，
: k% P7 E. ]2 E% r( N$ }$ M% Gu''(x)=sin(2pix);
u(-1)=0;
u(1)=0;
范围[-1;1]

$ \+ |4 z) Y: T! b请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

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