@@ -0,0 +1,2265 @@
+import pyinduct as pi%0A%0Av = 10%0Ac1, c2, c3 = %5B1, 1, 1%5D%0Al = 5%0AT = 5%0Aspat_domain = pi.Domain(bounds=(0, l), num=51)%0Atemp_domain = pi.Domain(bounds=(0, T), num=100)%0A%0A_, init_funcs1 = pi.cure_interval(pi.LagrangeSecondOrder,%0A                                  spat_domain.bounds,%0A                                  node_count=51)%0A_, init_funcs2 = pi.cure_interval(pi.LagrangeFirstOrder,%0A                                  spat_domain.bounds,%0A                                  node_count=30)%0Api.register_base(%22x1_funcs%22, init_funcs1)%0Api.register_base(%22x2_funcs%22, init_funcs2)%0A%0Au = pi.SimulationInputSum(%5B%0A    pi.SignalGenerator('square', temp_domain, frequency=.03,%0A                       scale=2, offset=4, phase_shift=1),%0A%5D)%0A%0Ax1 = pi.FieldVariable(%22x1_funcs%22)%0Apsi1 = pi.TestFunction(%22x1_funcs%22)%0Ax2 = pi.FieldVariable(%22x2_funcs%22)%0Apsi2 = pi.TestFunction(%22x2_funcs%22)%0A%0Aweak_form1 = pi.WeakFormulation(%0A    %5B%0A        pi.IntegralTerm(pi.Product(x1.derive(temp_order=1), psi1), limits=spat_domain.bounds),%0A        pi.IntegralTerm(pi.Product(x1, psi1.derive(1)), limits=spat_domain.bounds, scale=-v),%0A        pi.ScalarTerm(pi.Product(x1(l), psi1(l)), scale=v),%0A        pi.ScalarTerm(pi.Product(pi.Input(u), psi1(0)), scale=-v),%0A        pi.IntegralTerm(pi.Product(x1, psi1), limits=spat_domain.bounds, scale=c1),%0A        pi.IntegralTerm(pi.Product(x2, psi1), limits=spat_domain.bounds, scale=-c1),%0A    %5D,%0A    name=%22fluid temperature%22%0A)%0Aweak_form2 = pi.WeakFormulation(%0A    %5B%0A        pi.IntegralTerm(pi.Product(x2.derive(temp_order=1), psi2),limits=spat_domain.bounds),%0A        pi.IntegralTerm(pi.Product(x1, psi2), limits=spat_domain.bounds, scale=-c2),%0A        pi.IntegralTerm(pi.Product(x2, psi2), limits=spat_domain.bounds, scale=c2 + c3),%0A    %5D,%0A    name=%22wall temperature%22%0A)%0A%0Aics = %7Bweak_form1.name: %5Bpi.Function(lambda z: 0)%5D,%0A       weak_form2.name: %5Bpi.Function(lambda z: 0)%5D%7D%0Aspat_domains = %7Bweak_form1.name: spat_domain, weak_form2.name: spat_domain%7D%0Aevald1, evald2 = pi.simulate_systems(%5Bweak_form1, weak_form2%5D, ics, temp_domain,%0A                                     spat_domains)%0A%0Awin1 = pi.PgAnimatedPlot(%5Bevald1, evald2%5D, labels=dict(bottom='z'))%0Awin3 = pi.PgSurfacePlot(evald1, title=weak_form1.name)%0Awin4 = pi.PgSurfacePlot(evald2, title=weak_form2.name)%0Api.show()%0A

Add first example with 2 pde's