求助，逼近积分

splelvis

u''(x)=f(x),
+ B4 |+ Z. j' g1 M5 v边界boundary conditions u(a)=0,u(b)=0.
6 X  _2 i& Z, |  X6 _1 H* P(u(x+h)-2u(x)+u(x-h))/h^2=f(x);

% G5 R( j  ]  x- Y  Hmaltab code:

& B1 S# a+ W/ y7 ]6 l$ b3 }# \clear all
( T0 Y; x* D! M3 a3 A1 Ya=0;
b=1;
# |( k3 }  t. @1 ^# M$ H5 M8 P9 w: [n=10;
h=(b-a)/(n+1);
x_with_boundary=linspace(a,b,n+2)';
5 M4 D& X  y; Y- i3 |7 z! Ix=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));
, g! g( d$ m, J" E4 ^" R( d1 jrhs=f4(x);
sol=A\rhs;
sol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);

' @' _. W2 T! ]
* E6 j; P$ i+ r. }另开一文档，f4,
; C- c, O4 o9 T3 Ofunction y=f4(x)
y=ones(length(x),1);

出来后是抛物线。
, W9 K! v4 R$ D  t( K

) {7 f; D+ b. z0 z现在问题是，
u''(x)=sin(2pix);
u(-1)=0;
/ g( C) s. \+ D% N) ~2 {u(1)=0;
范围[-1;1]

* a+ D! }* l/ T7 s( I* Y请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

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