clear all;
options(1)=2;
% This option allows the display of convergence behavior
% X0 is the initial guess

x0 = 0.5;
%fsolve is called with the equations stored in a file cstreqn.m
%for p1=0:100:10000,
x = fsolve('p426nlfun',x0);
%end
%end

%Type X to display converged vaues

