多自由度体系中心差分法matlab编程

clc;
clear;
I=[1;1;1];
x=0:0.02:15;
xg=cos（x）;
M=180000*[1 0 0;0 1.5 0;0 0 1.5];
K=98000000*[1 -1 0;-1 3 -2;0 -2 4.5];
C=100000*[0.393 -0.2254 0;-0.2254 0.9275 -0.4508;0 -0.4508 1.2656];
P=[0;0;0];
u（:2）=[0;0;0];
v（:2）=[0;0;0];
w（:2）=[1;1/1.5;1/2];
h=0.02;
a=[1/（h^2）1/（2*h）2/（h^2）（h^2）/2];