层合板刚度矩阵求解matlab

function [ output_args ] = Untitled3（ input_args ）
%TEST.m
%MECHANICS OF COMPOSITE MATERIALS EX1
%清屏，清除原有变量，程序开始
echo off;
clc;
clear;
format short;
%…………………………………………………………………………………………………………
%1、预处理，定义已知材料性能指标，输入相关参数

%N = input（‘Input N:‘）;
N  = 4;
%E1 = input（‘Input E1:‘）;
E1 = 140;
%E2 = input（‘Input E2:‘）;
E2 = 10;
%G12 = input（‘Input G12:‘）;
G12 = 5;
%u21 = input（‘Input u12:‘）;
u21 = 0.3;
%第一层到最后一层与主方向夹角，单位°
    Agl = zeros（1N）;
    for n = 1:N
    Agl（1n） = input（[‘Input No‘num2str（n）‘‘s angle:‘]）;
    end;


%对应每层单层板的厚度，单位mm
str = input（‘Is tk changeless? please input y or n :‘‘s‘）;
for ii = 1:length（str）
if uint8（str（ii））==89 || uint8（str（ii））==121
    tk_chs = input（‘Input tk:‘）;
    tk = zeros（1N）;
    for n = 1:N
    tk（1n） = tk_chs;