clear all
cguess=100:10:9000;
tguess=390:-(390-298)/891:298;
csol=[];
tsol=[];
for i=1:(9000-100)/10+1
[res,fval]=fsolve('cstrpt',[cguess(i),tguess(i)],[]);
csol(i)=res(1);
tsol(i)=res(2);
cval(i)=fval(1);
tval(i)=fval(2);
end
plot(csol,tsol,'x')
figure
plot(csol,cval,'o')
figure
plot(tsol,tval,'+')