求助，逼近积分

splelvis

u''(x)=f(x),
1 W6 n* G0 {: f+ R! S" Z边界boundary conditions u(a)=0,u(b)=0.
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);

8 n# {# o' S" L& }, H8 xmaltab code:

clear all
a=0;
% h& U; _# G8 R* y8 ~1 }2 ^b=1;
* O4 r& Y- S9 T" r7 An=10;
. x- G& c% L- Z' C1 Hh=(b-a)/(n+1);
x_with_boundary=linspace(a,b,n+2)';
2 t6 N2 t  ^* _+ m% \/ `+ [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);
sol=A\rhs;
9 r# g2 _) K6 V2 Asol_with_boundary_conds=[0;sol;0];
) y7 e9 K1 \; c8 r0 f, fplot(x_with_boundary,sol_with_boundary_conds);

; s: A7 c' v+ C1 [
3 c: p' L3 L3 J( m8 ~7 h+ ]9 T+ J3 ?另开一文档，f4,
function y=f4(x)
y=ones(length(x),1);
! z4 U( ~( j9 Y! y! M% H  L
  f" ?5 }3 }1 X% J- F) u7 O出来后是抛物线。
7 r/ R8 t& t. x) x' ]( ^

% u* U# M* S+ Y! L* p现在问题是，
u''(x)=sin(2pix);
u(-1)=0;
u(1)=0;
0 V' }% q# a* j, T% g1 }9 n+ M范围[-1;1]
+ p1 h9 {( g( P1 H- m/ z
请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

没人能帮忙一下吗，
" @2 i7 N; U' g5 N0 T还是有什么资料不明确？