这是一个关于换热器计算的程序，我觉得不对，可我水平有限，谁能帮我改一下？


clear all; clc
global T1 T2 G t1 JA beta K theta JW CW Cc Q
T1 = 135; T2 = 40; G = 4e4; t1 = 30;
JA = 400; beta = 0. 15; K = 840;
Theta = 7900; Jw = 0. 1; Cw = 4. 184;
Cc = 2. 092;Q = G3 (T1 - T2) ;
T0 = 50;
T2 = fiminsearch (@Totalfee, t0) ;
Fp rintf (‘优化结果: /n /n’)
Fp rintf (‘换热器最优出口温度: %. 2% s/n’, t2, ’℃’)
Allfee = totalfee ( t2 ) ; fp rintf (‘最小年费用为: %. 3f 元/n’, allfee)
[AW ] =Area_water ( t2) ; frintf (‘换热器传热面积为: %.3fm^2 /n’,A)
fp rintf (‘每小时用水量为: %. 1fkg/h /n’,W)
fee1 = JA3 A3 beta; fee2 = JW3 theta3 W /1000; %
function J = totalfee ( t2)
global T1 T2 G t1 JA beta K theta JW CW Cc Q
[AW ] =Area_water( t2) ; J = JA3 A3 beta +JW3 theta3 W /1000%
function [A,W ] =Area_water ( t2)
global T1 T2 G t1 JA beta K theta JW CW Cc Q
var1 = T1 - t2; var2 = T2 - t1; dtm = ( var1 - var2 ) / log( var1 /var2) ;
A =Q / (K3 dtm) ;W =Q /Cw/ ( t2 - t1) ;