求助，逼近积分

splelvis

u''(x)=f(x),
边界boundary conditions u(a)=0,u(b)=0.
- R0 Q- j& f: k) v(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
/ V6 N. T  O6 G: g2 F
: D! E6 s# E- ~& C9 u7 W. H3 {5 omaltab code:
. E2 ?  f7 l8 F) e' ?8 E$ i; \% ?$ [
3 @5 q7 H/ H0 S" Bclear all
( N: o+ s! Z! U# _a=0;
5 b9 @3 v% Y: k: |$ g, qb=1;
n=10;
/ O2 f+ p, _$ t- L, u3 \7 Ah=(b-a)/(n+1);
8 A( n/ q( M, ^6 z/ O$ z  r$ `x_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));
  W8 O6 S2 F* V# K9 v5 h4 xrhs=f4(x);
sol=A\rhs;
( \' o; Y( g# R+ k( Dsol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);
$ C" h. D7 r$ ]+ _* O
6 D1 G/ v2 U3 h% }2 i( H! V/ j
另开一文档，f4,
function y=f4(x)
y=ones(length(x),1);
3 k7 x; M7 j% e8 n* p; {0 w
出来后是抛物线。

( O0 m* ?1 |3 h: D' u- O- h1 W" w
  J8 n; Z0 L4 r- p% B  y8 r* _" l现在问题是，
: U# V0 r+ x0 ^" _8 `u''(x)=sin(2pix);
u(-1)=0;
u(1)=0;
8 k' u, [, ?; u$ O7 g* G范围[-1;1]

/ [/ D! F4 t4 B( d请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

没人能帮忙一下吗，
/ Y. @3 U2 I2 n7 _8 R: ^+ Y, B还是有什么资料不明确？