clc;clear;close all

T0=0;Q=1672;rho=1000;cp=4.18;

q=[1e-4:1e-4:1];

v0 = T0+Q./q./rho./cp;

figure
plot(q,v0);