function out = model
%
% heat2D.m
%
% Model exported on Feb 25 2011, 23:36 by COMSOL 4.1.0.154.

import com.comsol.model.*
import com.comsol.model.util.*

model = ModelUtil.create('Model');

model.modelPath('C:\COMSOL41');

model.modelNode.create('mod1');

model.geom.create('geom1', 2);

model.mesh.create('mesh1', 'geom1');

model.physics.create('ht', 'HeatTransfer', 'geom1', {'I1' 'I2' 'I3' 'I4' 'I5' 'I6' 'I7' 'I8' 'I9' 'I10' 'I11' 'I12'});

model.study.create('std1');
model.study('std1').feature.create('stat', 'Stationary');

model.geom('geom1').feature.create('sq1', 'Square');
model.geom('geom1').feature('sq1').set('base', 'corner');
model.geom('geom1').feature('sq1').set('pos', {'-0.6' '-5.5511151231258E-17'});
model.geom('geom1').feature('sq1').set('l', '0.4');
model.geom('geom1').run('sq1');
model.geom('geom1').run;

model.physics('ht').feature('solid1').set('k_mat', 1, 'userdef');
model.physics('ht').feature('solid1').set('k', {'1' '0' '0' '0' '1' '0' '0' '0' '1'});
model.physics('ht').feature('solid1').set('rho_mat', 1, 'userdef');
model.physics('ht').feature('solid1').set('rho', 1, '1');
model.physics('ht').feature('solid1').set('Cp_mat', 1, 'userdef');
model.physics('ht').feature('solid1').set('Cp', 1, '1');
model.physics('ht').feature.create('temp1', 'TemperatureBoundary', 1);
model.physics('ht').feature('temp1').selection.set([1]);
model.physics('ht').feature('temp1').set('T0', 1, '100');

model.mesh('mesh1').run;


model.sol.create('sol1');
model.sol('sol1').study('std1');
model.sol('sol1').feature.create('st1', 'StudyStep');
model.sol('sol1').feature('st1').set('study', 'std1');
model.sol('sol1').feature('st1').set('studystep', 'stat');
model.sol('sol1').feature.create('v1', 'Variables');
model.sol('sol1').feature.create('s1', 'Stationary');
model.sol('sol1').feature('s1').feature.create('fc1', 'FullyCoupled');
model.sol('sol1').feature('s1').feature('fc1').set('dtech', 'auto');
model.sol('sol1').feature('s1').feature('fc1').set('initstep', 0.01);
model.sol('sol1').feature('s1').feature('fc1').set('minstep', 1.0E-6);
model.sol('sol1').feature('s1').feature('fc1').set('maxiter', 50);
model.sol('sol1').feature('s1').feature.create('d1', 'Direct');
model.sol('sol1').feature('s1').feature('d1').set('linsolver', 'pardiso');
model.sol('sol1').feature('s1').feature('fc1').set('linsolver', 'd1');
model.sol('sol1').feature('s1').feature('fc1').set('dtech', 'auto');
model.sol('sol1').feature('s1').feature('fc1').set('initstep', 0.01);
model.sol('sol1').feature('s1').feature('fc1').set('minstep', 1.0E-6);
model.sol('sol1').feature('s1').feature('fc1').set('maxiter', 50);
model.sol('sol1').feature('s1').feature.remove('fcDef');
model.sol('sol1').attach('std1');

model.result.create('pg1', 2);
model.result('pg1').set('data', 'dset1');
model.result('pg1').feature.create('surf1', 'Surface');
model.result('pg1').feature('surf1').set('expr', 'T');
model.result('pg1').feature('surf1').set('descr', 'Temperature');

model.sol('sol1').runAll;

model.result('pg1').set('windowtitle', 'Graphics');
model.result('pg1').run;



mphplot(model,'pg1');

out = model;