clear;?
U=input（‘please?input?U:?‘）;?
X=input（‘please?input?X:?‘）;
Y=input（‘please?input?Y:?‘）;?
Z=input（‘please?input?Z:?‘）;
W=input（‘please?input?W:?‘）;?
for?i=7:length（W）?
M{i-3}（14）=0;M{i-3}（34）=0;M{i-3}（24）=0;?
M{i-3}（11）=（U（i+1）-U（i））.^2/（U（i+1）-U（i-1））./（U（i+2）-U（i-1））;?
M{i-3}（13）=（U（i）-U（i-1））.^2/（U（i+1）-U（i-1））./（U（i+2）-U（i-1））;?
M{i-3}（23）=3.*（U（i）-U（i-1））.*（U（i+1）-U（i-1））./（U（i+1）-U（i-1））./（U（i+2）-U（i-1））;?
M{i-3}（33）=3.*（U（i+1）-U（i））.^2/（U（i+1）-U（i-1））./（U（i+2）-U（i-1））;?
M{i-3}（44）=（U（i+1）-U（i））.^2/（U（i+2）-U（i））./（U（i+3）-U（i））;?

M{i-3}（12）=1-M{i-3}（11）-M{i-3}（13）;??
M{i-3}（21）=-3.*M{i-3}（11）;?
M{i-3}（22）=3.*M{i-3}（11）-M{i-3}（23）;?
M{i-3}（31）=3.*M{i-3}（11）;?
M{i-3}（32）=-3.*（M{i-3}（21）+M