clc
clear
t11=input（‘t11?‘）;
t12=input（‘t12?‘）;
t21=input（‘t21?‘）;
t22=input（‘t22?‘）;
P1=input（‘p1?‘）;
P2=input（‘p2?‘）;
M1=input（‘M1?‘）;


tm1=（t11+t12）/2;
cp1=input（‘cp1?單位kJ/kg.K?‘）;
cp1=cp1*1000;
p1=input（‘p1密度?‘）;%密度
u1=input（‘u1?‘）;
lm1=input（‘lm1?‘）;
Pr1=（u1*cp1）/lm1;
tm2=（t21+t22）/2;
cp2=input（‘cp2?單位kJ/kg.K‘）;
cp2=cp2*1000;
p2=input（‘p2?密度?‘）;%密度



lm2=input（‘lm2?‘）;
u2=input（‘u2?‘）;
Pr2=（u2*cp2）/lm2;

nL=0.98;??????%取用
Q=M1*cp1*（t11-t12）*nL;
M2=Q/（cp2*（t22-t21））;
deta_t1mc=abs（（t11-t12）-（t22-t21））/（log（abs（（t11-t12）/（t22-t21））））;
P=（t22-t21）/（t11-t21）;
R=（t11-t12）/（t22-t21）;


%（1-2）型，具體見p18
%w1=sqrt（R^2+1）/（R-1）/2;
%w2=log（（1-R）/（1-P*R））;
%w3=2*（1+sqrt（（1-P）*（1-P*R）））-P*（1+R-sqrt（R^2+1））;
%w4=2*（1+sqrt（（1-P）*（1-P*R）））-P*（1+R+sqrt（R^2+1））;