#include?
#include?
void?main（）
{?int?ijkln;
?float?ha[100][100]b[100]x[100]max;
?printf（“輸入區間等分數n（步長h=1.0000/n）\n“）;
?scanf（“%d“&n）;
?h=1.00000/n;
b[0]=0;b[1]=2*h;
for（i=2;i{b[i]=（i-1）*h*h*h;
}
for（i=0;i{for（j=0;ja[i][j]=0;
}
a[0][0]=-3;a[0][1]=4;a[0][2]=-1;
a[1][n-2]=1;a[1][n-1]=-4;a[1][n]=3;
for（i=2;i{a[i][i-2]=1;
a[i][i-1]=-2-h*h;
a[i][i]=1;
}
printf（“輸出矩陣a[n+1][n+1]的值\n“）;
for（i=0;i<=n;i++）
{for（j=0;j<=n;j++）
printf（“%a[%d][%d]=%f\n“ija[i][j]）;
}
printf（“輸出數組b[n+1]的值\n“）;
for（i=0;i<=n;i++）
printf（“%b[%d]=%