init = asseminit(fem);
fem.sol=femstatic(fem, ...
'init',init, ...
'solcomp',{'phi0'}, ...
'outcomp',{'lm3','phi','phi0','pc','nc','lm4'});
fem0=fem;
fem.sol=femstatic(fem, ...
'init',fem0.sol, ...
'solcomp',{'lm3','phi','pc','nc','lm4'}, ...
'outcomp',{'lm3','phi','phi0','pc','lm4','nc'}, ...
'pname','Vd', ...
'plist',[0:0.02:0.2 0.25:0.05:1]);
fem0=fem;
fem.sol=femstatic(fem, ...
'init',init, ...
'solcomp',{'phi0'}, ...
'outcomp',{'lm3','phi','phi0','pc','nc','lm4'});
fem0=fem;
fem.sol=femstatic(fem, ...
'init',fem0.sol, ...
'solcomp',{'lm3','phi','pc','nc','lm4'}, ...
'outcomp',{'lm3','phi','phi0','pc','lm4','nc'}, ...
'pname','Vd', ...
'plist',[0:0.02:0.2 0.25:0.05:1]);
fem0=fem;
