matlab 参数 计算矩阵,matlab中exp（含参数矩阵）怎么编程求解

clear

clc

R=1;L=20;h=0.002;y=5;N=1;m=1;n=1;

%p1%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

cc11=6.160e10;cc12=4.744e10;cc13=4.744e10;cc22=15.825e10;cc23=3.151e10;

cc33=15.825e10;cc44=6.337e10;cc55=4.000e10;cc66=4.000e10;pp=7173;

%p0%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

c11=6.110e10;c12=2.079e10;c13=2.079e10;c22=5.858e10;c23=2.502e10;

c33=5.858e10;c44=1.678e10;c55=1.658e10;c66=1.658e10;p=1740;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

syms o

for j=1:5;

C11=c11*(1-((2*j-1)/(2*y))^N)+cc11*((2*j-1)/(2*y))^N;

C12=c12*(1-((2*j-1)/(2*y))^N)+cc12*((2*j-1)/(2*y))^N;

C13=c13*(1-((2*j-1)/(2*y))^N)+cc13*((2*j-1)/