sundialsTB, a Matlab Interface to SUNDIALS - Site Index Page ...

231 box on ;
232 g r i d on ;
233 x l a b e l ( ’ t ’ ) ;
234 y l a b e l ( ’ | Step s i z e | ’ ) ;
235
236 pl = pl +1;
237 subplot ( npl , 1 , pl )
238 p l o t ( t , q , ’− ’ ) ;
239 i f strcmp ( xaxis , ’ l o g ’ )
240 s e t ( gca , ’ XScale ’ , ’ l o g ’ ) ;
241 end
242 hold on ;
243 box on ;
244 g r i d on ;
245 x l a b e l ( ’ t ’ ) ;
246 y l a b e l ( ’ Order ’ ) ;
247 end
248
249 % Counters
250 i f cntr
251 pl = pl +1;
252 subplot ( npl , 1 , pl )
253 semilogy ( t , nst , ’ k− ’ ) ;
254 hold on ;
255 semilogy ( t , nfe , ’b− ’ ) ;
256 semilogy ( t , nni , ’ r− ’ ) ;
257 semilogy ( t , netf , ’ g− ’ ) ;
258 semilogy ( t , ncfn , ’ c− ’ ) ;
259 i f strcmp ( xaxis , ’ l o g ’ )
260 s e t ( gca , ’ XScale ’ , ’ l o g ’ ) ;
261 end
262 box on ;
263 g r i d on ;
264 x l a b e l ( ’ t ’ ) ;
265 y l a b e l ( ’ Counters ’ ) ;
266 end
267
268 % S o l u t i o n components
269 i f s o l
270 f i g u r e ( hf2 ) ;
271 map = colormap ;
272 n c o l s = s i z e (map , 1 ) ;
273 hold on ;
274 f o r i = 1 :N
275 hp = p l o t ( t , y ( i , : ) , ’− ’ ) ;
276 i c = 1+( i −1)∗ f l o o r ( n c o l s /N) ;
277 s e t ( hp , ’ Color ’ ,map( ic , : ) ) ;
278 end
279 i f strcmp ( xaxis , ’ l o g ’ )
280 s e t ( gca , ’ XScale ’ , ’ l o g ’ ) ;
281 end
282 box on ;
283 g r i d on ;
284 x l a b e l ( ’ t ’ ) ;
21