geometric compatibility factor or M prime by schmid factor

[mugikumar]

Hi everyone,
I was trying to construct the map of geometric compatibility factor or m prime in the misorientation space with help of Schmid factor but it was showing error in Matlab, but when i was using the demo set it was working fine. i think slip system for bcc have some issue for array operation. please somebody help me to resolve it and find the attachment
if i have a workspace of m prime and misorientation of two grains. is this possible to construct the map of misorientation vs m prime. for example image is attached below

[Tijmenvermeij]

I think you're going to need to post some more code to see what is going on before anyone can help you.