求助，逼近积分

splelvis

u''(x)=f(x),
( g; |5 ?' h% ]5 _# }边界boundary conditions u(a)=0,u(b)=0.
5 K/ n) q4 s# L/ W/ z& j  _3 ^(u(x+h)-2u(x)+u(x-h))/h^2=f(x);

maltab code:

clear all
a=0;
b=1;
1 l2 ?4 I( w+ mn=10;
! j. d$ t4 u2 G% ?" c! ih=(b-a)/(n+1);
, L8 g. H7 N; w  t; ]: c  Y3 T& C/ N8 yx_with_boundary=linspace(a,b,n+2)';
0 f6 W+ v. v! O1 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;
sol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);
8 z( e% K9 J  ^% u3 T- t* P
( v" v8 A9 T# l8 K: Z
另开一文档，f4,
function y=f4(x)
# Z! M$ |8 ^$ J2 c5 N; x+ Vy=ones(length(x),1);

$ F8 O1 d* |" ^5 c出来后是抛物线。


现在问题是，
u''(x)=sin(2pix);
u(-1)=0;
u(1)=0;
范围[-1;1]

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

splelvis

没人能帮忙一下吗，
( t) {3 ]8 ?8 J. J: {还是有什么资料不明确？