model sistemaTermicoDoisForcado
  parameter Real c = 466;
  parameter Real m1 = 200;
  parameter Real T01 = 400;
  parameter Real m2 = 200;
  parameter Real T02 = 300;
  parameter Real k = 400;
  parameter Real A = 8000e-6;
  parameter Real L = 0.5;
  parameter Real R = 5;
  parameter Real V0 = 220;
  parameter Real t0 = 20000;
  Real T1, T2;
  Real q, P;
  Real i;
  Real U, V;
initial equation
  i = 0;
  T1 = T01;
  T2 = T02;
equation
  m1 * c * der(T1) = -q + P;
  q = k*A/L * (T1 - T2);
  m2 * c * der(T2) = q;
  V = R * i;
  P = V * i;
  U = V;
  U = if time < t0 then 0 else V0;
end sistemaTermicoDoisForcado;