Question

x=[520 280 156.6 78 39.62 3.1 0] y=[0 -30 -36 -35 -28.44 -9.4 0] for j=1:6 h(j)=x(j+1)-x(j) end for

x=[520 280 156.6 78 39.62 3.1 0]
y=[0 -30 -36 -35 -28.44 -9.4 0]
for j=1:6
h(j)=x(j+1)-x(j)
end
for j=1:6
i=j+1
a(j)=h(i-1)/(h(i-1)+h(1))
b(j)=1-a(j)
c(j)=6*((y(2)-y(1))/h(1)-(y(i)-y(i-1))/h(i-1))/(h(1)+h(i-1))
end
A=2*eye(6)
for i=1:5
A(i,i+1)=b(i)
A(i+1,i)=a(i+1)
end
A(1,6)=a(1)
A(6,1)=b(6)
[L,U]=lu(A)
p=L\c'
N=U\p
M=N'
for i=1:6
t=x(i+1):0.01:x(i)
s=M(i+1)*(t-x(i))^3/(6*h(i))-M(i)*(t-x(i+1))^3/(6*h(i))+(y(i+1)-M(i+1)*h(i)^2/6)*(t-x(i))/h(i)-(y(i)-M(i)*h(i)^2/6)*(t-x(i+1))/h(i)
plot(t,s)
hold on
end
请问我的程序对吗

Answer 1

x=[520 280 156.6 78 39.62 3.1 0] y=[0 -30 -36 -35 -28.44 -9.4 0] for j=1:6 h(j)=x(j+1)-x(j) end for

Answer 2

x=[520 280 156.6 78 39.62 3.1 0] y=[0 -30 -36 -35 -28.44 -9.4 0] for j=1:6 h(j)=x(j+1)-x(j) end for