function []=klclose_anal(cinit,winit,CTin);

%% klclose_anal(cinit,winit,CT);
%% Does time and nullcline analysis

[t,x]=ode45('klclosefun',[0:2:1500],[cinit;winit],[],CTin);

figure(1);clf;plot(t,x);
xlabel('time');ylabel('[Ca2+] and w');
title(['cinit = ' num2str(cinit) ' winit = ' num2str(winit) ' CT = ' num2str(CTin)]);
mess=' time analysis, hit return to continue', pause;
figure(2);clf;
[Cs,Ws,lams,CRAY,nullC,nullw]=klclosejac(CTin);
hold on;
plot(x(:,1),x(:,2),'mo-');
mess='phase plane plot, hit return to continue',pause;