求助，逼近积分

splelvis

u''(x)=f(x),
边界boundary conditions u(a)=0,u(b)=0.
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
7 q* V8 E" _1 J! c8 g
maltab code:
$ x: q" c: {1 H: `! V; I' E, u
clear all
2 z, H0 u6 G* h( k4 _/ Z( D' m/ |a=0;
( Q. l+ w( i6 P2 F( p5 `b=1;
n=10;
h=(b-a)/(n+1);
x_with_boundary=linspace(a,b,n+2)';
( X0 M3 G3 Z" b) S, Ox=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;
sol_with_boundary_conds=[0;sol;0];
% a9 n* L; C% s, G3 Qplot(x_with_boundary,sol_with_boundary_conds);

) ^+ d9 o% I/ B) S1 Q+ J4 N- Q
5 m" U& c" D2 d( n/ r$ P另开一文档，f4,
function y=f4(x)
y=ones(length(x),1);
9 j, i- F4 P# A4 ]$ P7 ]
0 i' p; J% q2 j出来后是抛物线。


7 l0 K3 d% E' {; M- X/ i- j现在问题是，
' S, j# _2 V6 l3 Q) m( Hu''(x)=sin(2pix);
3 t+ v2 x/ f" }9 Pu(-1)=0;
5 E% B6 c* A6 N8 }9 h8 `' Su(1)=0;
- t" B$ L* `8 l1 I  G) w) e范围[-1;1]

请问各位要怎么把问题代入原题目中，哪些数值需要修改。

splelvis

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