总结代码排版给我

好的，以下是使用ode45求解该系统的ODE方程的完整代码，已经按照MATLAB的代码排版规范排版：

clear
clc

% 定义参数和变量
syms t u(t) v(t)
z = 0;
R1 = 1.2e-2;
R2 = 9.2e-3;
Cin = 1.1e6/60;
Cwall = 1.86e8/60;
PN = 8000;
qin = 20;
qout = -15;
u1 = diff(u);
v1 = diff(v);

% 定义ODE方程
eq0 = Cin*u1 == - (u - v)/R1;
eq1 = Cin*u1 == PN - (u - v)/R1;
eq2 = Cwall*v1 == (u - v)/R1 - (v - qout)/R2;

% 定义初始条件
eq3 = u(0) == qin;
eq4 = v(0) == z;

% 使用dsolve求解ODE方程
[uSol1(t), vSol1(t)] = dsolve(eq1, eq2, eq3, eq4);
[uSol2(t), vSol2(t)] = dsolve(eq0, eq2, eq3, eq4);

% 定义初始条件和时间范围
tspan = [0, 60];
y0 = [qin; z; 0; 0];

% 定义ODE方程
function dydt = myODEfun(t, y)
    % 定义参数和变量
    R1 = 1.2e-2;
    R2 = 9.2e-3;
    Cin = 1.1e6/60;
    Cwall = 1.86e8/60;
    PN = 8000;
    qin = 20;
    qout = -15;
    u = y(1);
    v = y(2);

    % 定义ODE方程
    du_dt = y(3);
    dv_dt = y(4);
    dudt = Cin*du_dt + (u - v)/R1;
    dvdt = Cwall*dv_dt - (u - v)/R1 + (v - qout)/R2 - PN;

    % 返回ODE方程的结果
    dydt = [du_dt; dv_dt; dudt; dvdt];
end

% 使用ode45求解ODE方程
[t, y] = ode45(@myODEfun, tspan, y0);

% 绘制结果
figure
subplot(2,2,1), plot(t, y(:,1)), xlabel('时间/min'),ylabel('开启时室内温度')
subplot(2,2,2), plot(t, y(:,2)), xlabel('时间/min'),ylabel('开启时室内温度')
subplot(2,2,3), plot(t, y(:,3)), xlabel('时间/min'),ylabel('开启时室内温度')
subplot(2,2,4), plot(t, y(:,4)), xlabel('时间/min'),ylabel('开启时室内温度')

希望这能帮到你！