calcpsi

clc;
close all;
clear all;
psibetax=psi(0.3055,0)